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

Jason Stajich jason.stajich at gmail.com
Sat Jul 24 12:41:30 EDT 2010


Hi. You should ask questions to the list. There are many who can  
answer your question there.

In the older versions you can put the retrieve statement in an eval  
block. I think newer behavior may be only to issue a warning and keep  
going.

You may be able to also set verbose to -1 (-verbose => -1) in the  
initialization of the db obj and see whether it skips or fails on  
missing ID.
  -Jason.
Sent from my iPod

On Jul 23, 2010, at 6:48, Hongseok Tae <htae at vbi.vt.edu> wrote:

> Hi Jason Stajich,
>
> I am using Bioperl and it is very useful. I have a question about  
> Bio::DB::GenBank.
> This is a simple code.
>
> --------------------------
>
> $db = Bio::DB::GenBank->new();
> @ids = ("LOC441435", "X78121");
> foreach $id (@ids){
> $seqobj = $db->get_Seq_by_id($id);
> $seqstr   = $seqobj->seq();
> print "$>id\n$seqstr\n";
> }
>
> --------------------------
>
> If there is no entry for an ID, I would like to skip it and to get  
> next one. But this code stops running when an ID does not exist in  
> genbank.
> How can I solve this problem?
>
> Thanks.
> Hongseok Tae
>



More information about the Bioperl-l mailing list