[Bioperl-l] accessing genbank
S.Paul
s.paul at surrey.ac.uk
Thu Jun 3 17:50:34 EDT 2004
Hi Everybody:
I have been trying to access the Genbank for the following accession numbers and outfut it to a file and am getting the following error message:
****************************************************************************
MSG: acc does not exist
STACK Bio::DB::WebDBSeqI::get_Seq_by_acc C:/Perl/site/lib/Bio/DB/WebDBSeqI.pm:1
7
STACK Bio::DB::GenBank::get_Seq_by_acc C:/Perl/site/lib/Bio/DB/GenBank.pm:216
STACK toplevel Genbank_seq_dopamine.pl:11
******************************************************************************************************
Not sure what I am doing wrong. I can, however, retrieve individual accession number. I am enclosing the code.
Thanks for the help in advance
Sujoy Paul
********************************************************************************************
use Bio::DB::GenBank;
use Bio::SeqIO;
use Bio::DB::WebDBSeqI;
my $gb = new Bio::DB::GenBank;
my $seqio = $gb->get_Seq_by_acc(['AAH21195','AAH38978']);
my $seqout = new Bio::SeqIO(-fh =>'dopamin_human.gbk', -format => 'genbank');
***********************************************************************************************
Sujoy Paul, PRISE Centre, UniS, s.paul at surrey.ac.uk
More information about the Bioperl-l
mailing list