[Bioperl-l] Can't find Carp/Heavy.pm using BPbl2seq module
Gisel, Andreas
Andreas.Gisel@fmi.ch
Wed, 19 Jun 2002 14:28:17 +0200
Dear Bioperl users
I am new in the Perl/Bioperl world and run into the following problem I don't know how to solve:
I am using Bio::Tools::Run::StandAloneBlast and Bio::Tools::BPbl2seq for alignments between two DNA sequences (see below part of the script). The modules run in a loop and work fin until the 60th loop and the I get an error 'can't find Carp/Heavy.pm'. I checked and Carp/Heavy.pm it is at the place where it should be!!
What is the problem and how can I solve it??
Thousand Thanks
Andreas Gisel
my $seqio = Bio::SeqIO->new (-file => $queryin);
my $seq2 = $seqio->next_seq();
do
{
$primer = substr($genome, $start, $long);
my $seq1 = Bio::Seq->new ( -seq => $primer,
-id => 'my seq');
my $factory = Bio::Tools::Run::StandAloneBlast->new( 'prog' => 'blastn',
'outfile' => 'bl2seq.out');
my $bl2seq_report = $factory->bl2seq($seq1, $seq2);
#$bl2seq_report->next_feature;
#my $report = Bio::Tools::BPbl2seq->new(-file => 'bl2seq.out');
#$report->sbjctName;
#$report->sbjctLength;
while(my $hsp = $bl2seq_report->next_feature)
{
print "start: $start ",
"P value: ", $hsp->P,
" length: ", $hsp->length,
" (1st seq length : ", $seq1->length, " 2nd seq length: ",$seq2->length,")",
"\n";
}
$start++;
}
while ($start+$long < $que_leng);
_____________________________________________________________
Novartis Forschungsstiftung
Zweigniederlassung FMI
Maulbeerstrasse 66
4058 Basel
Switzerland
Tel.: ++41 61 697 8099
Fax: ++41 61 697 3976