[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