[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