Inquiry of Parsing tblastn result to get nucleotide sequence from the db file

[sunfuhui]

Hi all,

I met a problem when I tried to parse tblastn result:

Using either bioperl Bio::SearchIO ($hsp->hit_string) or -outfmt parameter of tblastn can only return the amino acid sequence of the hit.

But I want to get the original nucleotide sequence before translation (not the whole one, but only the part which matches the $hsp->query_string).

How could I do that? Any idea?

Thanks a lot in advance!!!!!!!!!

All the best,
Ying

[TiborNagy]

Hi,

tblastn won't give you the nucleotide sequences, however you get the subject start and end positions. You can use an Emboss program (seqret) or Bioperl to get the sequence from the original file.

[sunfuhui]

Thank you, TiborNagy!