[Bioperl-l] Bio::DB::GenBank batch retrieval question

Chris Chalfant CHALFANT_CHRIS_M@Lilly.com
Thu, 07 Mar 2002 08:21:20 -0500


When I use Bio::DB::GenBank::get_Stream_by_acc to retrieve a set of
accession numbers and include an invalid accession number in with
several valid accession numbers, I get no sequences back from Entrez.
For example, the code below returns no output (though I get output if I
remove the "bogus" accession).

Is this the expected behavior or am I using the code incorrectly?  If
this is the correct behavior, how would you suggest requesting genbank
records for a list which may include invalid accession numbers?

CODE:

my @accessions = qw(AB000095 AB000220 bogus);
my $gb = new Bio::DB::GenBank;
my $seqio = $gb->get_Stream_by_acc(\@accessions);

while (my $record = $seqio->next_seq) {
  print $record->primary_id, "\n";
}

Thanks,

Chris

Chris Chalfant
Bioinformatics
Eli Lilly and Company