[Bioperl-l] Downloading a sequence in genbank format - related problem

Diogo Tschoeke diogoat at gmail.com
Wed May 16 20:02:44 UTC 2007


Dear all,

The script wich i wrote with your helps is working very good ( I paste the
script in the end of e-mail).
But I have another problem now, all the times wich I use the script im every
all the file have a diferent size...
Any ideia? what is the problem..? My conection? Problem on Ncbi? The script
maybe?

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

#############################################################
use strict;
use warnings;
use Bio::Seq;
use Bio::SeqIO;
use Bio::DB::GenBank;
my $query = Bio::DB::Query::GenBank->new
                                (-query   =>'Trypanosoma cruzi [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 => '>>Trypanosoma_cruzi1.gb');
while (my $seq = $seqio->next_seq){
         $out->write_seq($seq);
                        }
#########################################################




More information about the Bioperl-l mailing list