[Bioperl-l] extracting 'species' from ncbi records
Charles Hauser
chauser@duke.edu
Fri Jan 17 16:28:17 EST 2003
All,
I am trying to parse species names from GenBank records.
the following code fails:
my $seqio = Bio::SeqIO->new('-format' => 'GenBank', -fh => \*FH);
while( my $seq = $seqio->next_seq ) {
next unless $seq->species() =~/reinhardtii/;
foreach my $feat ( $seq->top_SeqFeatures ) {
............
Do I need to call Bio::Species directly?
my $species = Bio::Species->new();
Chuck
More information about the Bioperl-l
mailing list