Extract subject sequences from Blast result
I ran TBLASTN (query=protein, subject=nucleotide) and got the following result:
qseqid sseqid pident sstart send evalue MP_11765 contig_17 95.322 5460 5221 0
I noticed that sstart > ssend meaning this query is mapping on the reverse complement of the subject, correct? In this case how can I pull out the DNA sequence from my subject nucleotide file that corresponds to this hit?
Simply extracting the following range (5221-5460) from the subject nucleotide file does not seem to be the right way. Any ideas on a program that can do this?
• 64 views