[Bioperl-l] Help on retrieving NT contigs with Bio::DB::GenBank.
Kun Zhang
kzhang@sph.uth.tmc.edu
Fri, 30 Nov 2001 16:28:58 -0600
Hello!
I got a error message (attached below) when trying to retrieving some NT
contigs from GenBank with the Bio::DB::Genbank module. It looks like the
problem occurs only on NT sequence because the getGenBank.pl came with the
bioperl-0.9.0 distribution works fine. And my perl script works when I
replacing the "NT_001035" with "AF303112". Can anyone help me out? Thanks!
Kun Zhang
Human Genetics Center
University of Texas-Houston
------------------------My codes-----------------------------
my $gb = new Bio::DB::GenBank;
$gb->request_format('fasta') ;
$contigSeq = $gb->get_Seq_by_acc('NT_001035');
==============================ERROR MESSAGE==================================
-------------------- EXCEPTION --------------------
MSG: Attempting to set the sequence to [<HTML] which does not look healthy
STACK Bio::PrimarySeq::seq
/usr/local/lib/perl5/site_perl/5.6.1/Bio/PrimarySeq.pm:251
STACK Bio::PrimarySeq::new
/usr/local/lib/perl5/site_perl/5.6.1/Bio/PrimarySeq.pm:226
STACK Bio::Seq::new /usr/local/lib/perl5/site_perl/5.6.1/Bio/Seq.pm:132
STACK Bio::SeqIO::fasta::next_primary_seq
/usr/local/lib/perl5/site_perl/5.6.1/Bio/SeqIO/fasta.pm:130
STACK Bio::SeqIO::fasta::next_seq
/usr/local/lib/perl5/site_perl/5.6.1/Bio/SeqIO/fasta.pm:85
STACK Bio::DB::WebDBSeqI::get_Seq_by_acc
/usr/local/lib/perl5/site_perl/5.6.1/Bio/DB/WebDBSeqI.pm:159
STACK toplevel ./splitSeq.pl:26