/usr/share/doc/bioperl/examples/searchio/htmlwriter.pl is in bioperl 1.6.901-2.
This file is owned by root:root, with mode 0o755.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 | #!/usr/bin/perl
# Demonstrates the use of a SearchIO Blast parser and a SearchWriterI object
# for producing HTML Blast output from a Blast report input stream.
#
# Usage:
# STDIN: none; supply filename of BLAST report on command-line
# STDOUT: none; generates an output file "searchio.html"
# containing HTML-formatted Blast Report
# STDERR: Any errors that occurred.
#
# For more documentation about the writer, including
# a complete list of columns, see the docs for
# Bio::SearchIO::Writer::HTMLResultWriter.
#
# For more documentation about working with Blast result objects,
# see docs for these modules:
# Bio::Search::Result::BlastResult
# Bio::Search::Iteration::IterationI
# Bio::Search::Hit::BlastHit
# Bio::Search::HSP::BlastHSP
#
# For more documentation about the Blast parser, see docs for
# Bio::SearchIO
#
# Author: Steve Chervitz <sac@bioperl.org>
use strict;
use lib '../../';
use Bio::SearchIO;
use Bio::SearchIO::Writer::HTMLResultWriter;
my $outfile = "searchio.html";
my $file = shift or die "Usage: $0 <BLAST-report-file>\n HTML output is saved to $outfile\n";
my $in = Bio::SearchIO->new( -format => 'blast',
-file => $file, #comment this out to read STDIN
#-fh => \*ARGV, #uncomment this to read from STDIN
-verbose => 0 );
my $writer = new Bio::SearchIO::Writer::HTMLResultWriter();
my $out = new Bio::SearchIO(-writer => $writer,
-file => ">$outfile");
while ( my $result = $in->next_result() ) {
eval {
# printf STDERR "Report %d: $result\n", $in->result_count;
$out->write_result($result, 1);
};
if($@) {
warn "Warning: Blast parsing or writing exception caught for $result:\n$@\n";
}
}
printf STDERR "\n%d Blast report(s) processed.\n", $in->result_count;
printf STDERR "Output sent to file: %s\n", $out->file if $out->file;
|