[Bioperl-l] retreive all seq of a species

Paulo Almeida paulo.david at netvisao.pt
Wed Jun 9 05:45:45 EDT 2004


Something like this should work:

use Bio::DB::Query::GenBank;
use Bio::DB::GenBank;

open(OUT,">test");
my $query_string = '"Antheraea mylitta"[Organism]';
my $query = Bio::DB::Query::GenBank->new(
                                                                
-db=>'nucleotide',
                                                                 
-query=>$query_string,
                                                                 );
my $gb = new Bio::DB::GenBank;
my $stream = $gb->get_Stream_by_query($query);
while (my $seq = $stream->next_seq) {
     print OUT ">" , $seq->display_id , "\n" , $seq->seq , "\n\n";
}

 > Using bio perl modules how i can retreive all the
 > nucleotide sequences submitted in genbank for a
 > particular species (eg. Antheraea mylitta) .



More information about the Bioperl-l mailing list