#!/usr/bin/perl -w use strict; use Bio::DB::Query::GenBank; use Bio::DB::GenBank; use Bio::SeqIO; my $query_string = 'Oryza[Organism] AND EST[Keyword]'; my $query = Bio::DB::Query::GenBank->new(-db=>'nucleotide', -query=>$query_string, -mindate => '2002', -maxdate => '2002'); my $count = $query->count; warn "There were $count sequences from 2002 that met criteria: $query_string\n"; my $output = Bio::SeqIO->new(-format => 'genbank', -file => '>sequences.gbk'); my $output2 = Bio::SeqIO->new(-format => 'fasta', -file => '>sequences.fa'); my $gb = Bio::DB::GenBank->new; my $stream = $gb->get_Stream_by_query($query); my $i = 0; while (my $seq = $stream->next_seq) { $output->write_seq($seq); $output2->write_seq($seq); if( $i++ > 10 ) { last; } }