[Bioperl-l] HTMLResultWriter
lpaulet at ual.es
lpaulet at ual.es
Mon Jan 25 20:31:08 UTC 2010
Hi all,
I'm trying to generate a subroutine that performs a BLAST search and
returns the corresponding reports in txt, xml and html format. I´m
experiencing problems with the latter, as the program returns the
following error message:
"Can't call method "next_result" without a package or object reference at..."
sub blasting {
my ($query, $E_value) = @_;
my ($outputfilenameB, $outputfilenameX, $outputfilenameH);
$outputfilenameB=$query.".BLAST.txt";
$outputfilenameX=$query.".BLAST.xml";
$outputfilenameH=$query.".BLAST.html";
#legacy_blast.pl blastall -i query -d nr -o blast.out --path /opt/blast/bin
print qx(du -s /tmp);
my $blast_report =qx/$blast -p blastp -d $database -i $query -e
$E_value -b 20000 -o $outputfilenameB/;
my $XMLblast_report =qx/$blast -p blastp -d $database -i $query -e
$E_value -m 7 -b 20000 -o $outputfilenameX/;
my $writerhtml = new Bio::SearchIO::Writer::HTMLResultWriter();
my $outhtml = new Bio::SearchIO(-writer => $writerhtml,
-file => ">$outputfilenameH");
while( my $result = $blast_report->next_result ) { # get a result from
Bio::SearchIO parsing or build it up in memory
$outhtml->write_result($result);
}
}
Can anyone see where the problem is?
Cheers!
Lorenzo
More information about the Bioperl-l
mailing list