[Bioperl-l] producing an organism report from BLAST results?
Jason Stajich
jason.stajich at duke.edu
Mon Nov 21 16:26:23 EST 2005
Well isn't encoded in the hit name?
my ($gi) = ($hit->name =~ /gi\|(\d+)/ );
If you run WU-BLAST you need to provide the -gi option on the command
line.
If you run NCBI-BLAST you need to provide the -I option
From the blastall help:
-I Show GI's in deflines [T/F]
default = F
-jason
On Nov 21, 2005, at 4:22 PM, Tim Reid wrote:
>
> You need to look up the taxonomid ID from the gi number for a record.
>
> thanks for your reply, but i am unable to find the gi number of a
> blast hit.
>
> i have looked at the each_accession_number method in
> Bio::Search::HitI, but it does not appear to yield a gi number.
>
> your help is appreciated.
--
Jason Stajich
Duke University
http://www.duke.edu/~jes12
More information about the Bioperl-l
mailing list