[Bioperl-l] Downloading a sequence in genbank format

Diogo Tschoeke diogoat at gmail.com
Tue May 15 22:44:59 UTC 2007


Dear All,

I need to download a lot of sequence of Leishmania major in genbank
format...
But i can't download on the page of NCBI, because the downloaded file are
corrupted... when i use a browser to download this sequences
And them i looking for some script to download that`s file and fink
something like that:


#########################################################
use strict;
use warnings;

use Bio::Seq;
use Bio::SeqIO;
use Bio::DB::GenBank;

my $query = Bio::DB::Query::GenBank->new
                                (-query   =>'Leishmania major [Organism]',
                                -db      => 'nucleotide');
my $gb = new Bio::DB::GenBank;
my $seqio = $gb->get_Stream_by_query($query);

my $out = Bio::SeqIO->new(-format => 'genbank',
                          -file => '>>teste6.gb');
$out->write_seq($seqio);
#########################################################

And the system return me this erros
[diogo1 at genome perl]$ perl teste6.pl

-------------------- WARNING ---------------------
MSG:  Bio::SeqIO::genbank=HASH(0x96c0f08) is not a SeqI compliant module.
Attempting to dump, but may fail!
---------------------------------------------------
Can't locate object method "seq" via package "Bio::SeqIO::genbank" at
/usr/lib/perl5/site_perl/5.8.8/Bio/SeqIO/genbank.pm line 692.

Any Ideia?

Thank`s

Diogo Tschoeke
Laboratory of Molecular Biology of Trypanosomatides
Fundação Osvaldo Cruz - Fiocruz RJ, Brazil
http:biowebdb.org <http://www.ncbs.res.in/>




More information about the Bioperl-l mailing list