[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