[Bioperl-l] problem with Bio::DB::Genbank
shalabh sharma
shalabh.sharma7 at gmail.com
Tue Aug 13 19:57:11 UTC 2013
Hi All,
I am trying to download some sequences. I am using Bio::DB::Genbank.
But i am not getting the sequence , i am only getting the header part.
This the the script:
#!/usr/bin/perl
use Bio::DB::GenBank;
use Bio::SeqIO;
my $gb = new Bio::DB::GenBank;
my $seqout = new Bio::SeqIO(-fh => \*STDOUT, -format => 'fasta');
my $seq = $gb->get_Seq_by_id('116513137');
$seqout->write_seq($seq);
Output:
>YP_812044 30S ribosomal protein S18 [Lactococcus lactis subsp. cremoris
SK11].
I would really appreciate if anyone can help me out.
Thanks
Shalabh
--
More information about the Bioperl-l
mailing list