SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Parsing a VCF file for indels bnfoguy Bioinformatics 1 07-12-2012 10:19 AM
Dindel crash whilst parsing bam-sorted file michalkovac Bioinformatics 8 02-16-2012 06:56 PM
Blast > parsing result in Exel Giorgio C Bioinformatics 6 11-16-2011 12:16 AM
Parsing a Chromosome sequence with start and end coordinates empyrean Bioinformatics 13 09-10-2011 09:23 PM
Bias toward G in first nucleotide in sequence? sem Sample Prep / Library Generation 0 01-16-2009 11:54 AM

Reply
 
Thread Tools
Old 09-11-2012, 02:26 PM   #1
sunfuhui
Junior Member
 
Location: athens, ga

Join Date: Jul 2012
Posts: 8
Default Inquiry of Parsing tblastn result to get nucleotide sequence from the db file

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
sunfuhui is offline   Reply With Quote
Old 09-12-2012, 02:03 AM   #2
TiborNagy
Senior Member
 
Location: Budapest

Join Date: Mar 2010
Posts: 329
Default

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.
TiborNagy is offline   Reply With Quote
Old 09-12-2012, 05:39 AM   #3
sunfuhui
Junior Member
 
Location: athens, ga

Join Date: Jul 2012
Posts: 8
Default

Thank you, TiborNagy!
sunfuhui is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 02:51 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO