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

James Wasmuth james.wasmuth at ed.ac.uk
Thu Jul 3 12:36:06 EDT 2003


Hi all,

I using the below script:


use Bio::DB::GenBank;
my @query_acc = qw/AF043104 CE24123 AF172351 AF435952/;
$embl=new Bio::DB::GenBank;
foreach $acc (@query_acc)       {
    print ":$acc:\n";
    $seq = $embl->get_Seq_by_acc($acc);
    ...
    ...
}

However I get the following error...
------------- EXCEPTION  -------------
MSG: acc does not exist
STACK Bio::DB::WebDBSeqI::get_Seq_by_acc 
/usr/local/lib/perl5/site_perl/5.8.0/Bio/DB/WebDBSeqI.pm:176


But AF043104 does exist.
If I add the version number so AF043104.1 then get_Seq_by_acc works!?!  

Am I doing something wrong or is this a bug?

thanks
James

-- 

Nematode Bioinformatics
Blaxter Nematode Genomics Group
Institute of Cell, Animal and Population Biology
Ashworth Labs
University of Edinburgh
King's Buildings
Edinburgh
EH9 3JT
UK

(+44)(0)131 650 7403




More information about the Bioperl-l mailing list