Extract subject sequences from Blast result

1

Hi all,

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?

Thank you.


sequence


blast


homology

• 64 views



Source link