Bioperl: Two newbie questions
philipp.pagel@cmp.yale.edu
philipp.pagel@cmp.yale.edu
Sun, 16 Apr 2000 22:31:26 -0400 (EDT)
Hi!
Disclaimer: I am not at all experienced with perl in general and just start
with bioperl...
And I am trying to find my way around in the examples and documentation,
because I want to write a litte script to do something usefull for me :-)
I have two questions that I could not figure out with the documentation:
1) I read from a file:
$in = Bio::SeqIO->new( '-format' => $informat , -file => $file);
while ($seq = $in->next_seq();){
...
}
This works fine for the fasta files included in the bioperl distribution.
But my Sequence happens to come in gcg format - so i just set $informat to
"gcg". Ok - my files include only one single sequence, but anyway... I
expected it to work nevertheless. But this turns out to be an endless loop
delivering the same sequence over and over again. Why is that?
2) I played with run_blast_remote.pl and copuld not get it to work. I just
tried the first example the script returns when evoked with -eg. Here is what
I get:
------------------------snip--------------------------
Use of uninitialized value at ./run_blast_remote.pl line 114.
-------------------- WARNING --------------------
MSG: /usr/lib/perl5/site_perl/5.005/Bio/Tools/Blast/Run/Webblast.pm:1205
Seq::type - deprecated method. You should use $obj->moltype in preference
(notice that
moltype returns lowercase strings)
-------------------------------------------------
Use of uninitialized value at
/usr/lib/perl5/site_perl/5.005/Bio/PrimarySeqI.pm
line 847.
Use of uninitialized value at
/usr/lib/perl5/site_perl/5.005/Bio/PrimarySeqI.pm
line 848.
Use of uninitialized value at
/usr/lib/perl5/site_perl/5.005/Bio/PrimarySeqI.pm
line 849.
-------------------- WARNING --------------------
MSG: /usr/lib/perl5/site_perl/5.005/Bio/Tools/Blast/Run/Webblast.pm:1206
Seq::seq_len - deprecated method. You should use $obj->length in preference
-------------------------------------------------
Use of uninitialized value at
/usr/lib/perl5/site_perl/5.005/Bio/PrimarySeqI.pm
line 717.
Use of uninitialized value at
/usr/lib/perl5/site_perl/5.005/Bio/Tools/Blast/Run/Webblast.pm line 1228.
The sequence length of (0) too short for blasting: min length = 15
*** Trouble creating Blast object:
-------------------- EXCEPTION --------------------
MSG: Blast failed: No file created or a error occured upon submission.
CONTEXT: Error in uNKNOWN CONTEXT
SCRIPT: ./run_blast_remote.pl
STACK:
Bio::Tools::Blast::Run::Webblast::_blast(1147)
Bio::Tools::Blast::Run::Webblast::blast_remote(560)
Bio::Tools::Blast::_run_remote(1244)
Bio::Tools::Blast::run(1174)
Bio::Tools::SeqAnal::_initialize(279)
Bio::Root::Object::new(462)
main::create_blast(555)
main::./run_blast_remote.pl(131)
---------------------------------------------------
main::create_blast called at ./run_blast_remote.pl line 131
eval {...} called at ./run_blast_remote.pl line 130
----------------------snip---------------------------
What am I doing wrong?
cu and thanks for any hints
Philipp
--
Dr. Philipp Pagel
Cellular and Molecular Physiology phone: (203) 785-6835
SHM BE-30
Yale University
333 Cedar ST
New Haven, CT 06520
=========== Bioperl Project Mailing List Message Footer =======
Project URL: http://bio.perl.org/
For info about how to (un)subscribe, where messages are archived, etc:
http://www.techfak.uni-bielefeld.de/bcd/Perl/Bio/vsns-bcd-perl.html
====================================================================