[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