[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