[Bioperl-l] RemoteBlast's save_output not saving properly for blasttable
Ulas Karaoz
UKaraoz at lbl.gov
Fri Aug 15 00:03:51 UTC 2008
Hi,
I found out that the save_output routine in RemoteBlast.pm doesn't
save the output properly into a file when the Blast output is a
hittable(blasttable). Might the reason be the fact that the tabular
output has a line that starts with a # while the parser is looking for
a line starting with just BLASTN, as in the section pasted below:
if( $l =~ /^(?:[T]?BLAST[NPX])\s*.+$/i ||
$l =~/^RPS-BLAST\s*.+$/i ) {
$seentop=1;
}
More information about the Bioperl-l
mailing list