sarf13 has asked for the wisdom of the Perl Monks concerning the following question:
Hi.. I was trying the following script for extracting sequences in fasta format from a list of geneID from genbank. I got that error message when i tried large number of sequence download(some 1000-2000) at a time. the script does not work for large number of sequence to be extracted at a time.
use Bio::Perl; $database="genbank"; $format="fasta"; $pipe ="\\|"; $space = " "; open(INPUTFILE, "<test.txt"); while(<INPUTFILE>) { my($line) = $_; chomp($line); $line=~ s/$space/:/; $line=~ s/$pipe/$space/; $line=~ s/g/G/; $line=~ s/i/I/; $id= "$line"; #print "$id"; #print "\n"; $sequence = get_sequence($database, $id); $test = write_sequence( ">>sequences_1.txt", $format, $sequence); open (CHK , ">>checking.txt"); print CHK <<HERE; $test HERE close CHK; } exit;
Thanks
|
|---|
| Replies are listed 'Best First'. | |
|---|---|
|
Re: more then 2000 file download from database error.
by marto (Cardinal) on May 19, 2012 at 13:56 UTC | |
|
Re: more then 2000 file download from database error.
by sauoq (Abbot) on May 19, 2012 at 14:02 UTC |