[Bioperl-l] extracting 'species' from ncbi records

Jason Stajich jason@cgt.mc.duke.edu
Fri Jan 17 16:49:26 EST 2003


$seq->species() returns and object (Bio::Species)

You want
next unless $seq->species->species() =~ /.../;


On Fri, 17 Jan 2003, Charles Hauser wrote:

> 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
>
>
>
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l@bioperl.org
> http://bioperl.org/mailman/listinfo/bioperl-l
>

--
Jason Stajich
Duke University
jason at cgt.mc.duke.edu



More information about the Bioperl-l mailing list