[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