[Bioperl-l] problems with Bio::Tools::BPbl2seq
Arye Shemesh
ARYES at wicc.weizmann.ac.il
Mon Dec 22 06:28:44 EST 2003
Hi,
I have a problem with Bio::Tools::BPbl2seq. When I try to run bl2seq on 2
Bio::Seq objects, i get these error messages:
Use of uninitialized value in print at
/usr/local/lib/perl5/site_perl/5.6.1/Bio/Root/IO.pm line 305, <GEN0> line 7
(#1)
The weird thing is that i do get the ruslts i need from my script!
Can you please help me find out what's going on here?
This is the context in which I use the bl2seq:
sub BlastSeqs
{
my $factory = Bio::Tools::Run::StandAloneBlast->new('program' =>
'blastp');
$factory->F("F");
my $StructSeqObj = Bio::Seq->new(-seq=>$SeqFromStruct);
foreach my $seqName (keys %{$DBSeqHashRef})
#iterating over the data hash
{
my $alignedSeqObj = Bio::Seq->new(
-id=>$seqName,
-seq=>$DBSeqHashRef->{$seqName});
my $blast_report;
eval
{
#THIS IS THE LINE THAT CAUSES THE ERROR:
$blast_report = $factory->bl2seq ($StructSeqObj, alignedSeqObj);
}; #of eval
print STDERR $@ if ($@);
}
} # of sub BlastSeqs
Thanks,
Arye
More information about the Bioperl-l
mailing list