[Bioperl-guts-l] [14786] bioperl-live/trunk/Bio/SearchIO/blasttable.pm: Previous changes killed RemoteBlast again; make a little more bulletproof
Christopher John Fields
cjfields at dev.open-bio.org
Tue Aug 5 18:00:01 EDT 2008
Revision: 14786
Author: cjfields
Date: 2008-08-05 18:00:00 -0400 (Tue, 05 Aug 2008)
Log Message:
-----------
Previous changes killed RemoteBlast again; make a little more bulletproof
Modified Paths:
--------------
bioperl-live/trunk/Bio/SearchIO/blasttable.pm
Modified: bioperl-live/trunk/Bio/SearchIO/blasttable.pm
===================================================================
--- bioperl-live/trunk/Bio/SearchIO/blasttable.pm 2008-08-05 01:17:28 UTC (rev 14785)
+++ bioperl-live/trunk/Bio/SearchIO/blasttable.pm 2008-08-05 22:00:00 UTC (rev 14786)
@@ -162,8 +162,6 @@
local $_;
my ($alg, $ver);
while( defined ($_ = $self->_readline) ) {
- # skip any HTML cruft (e.g. from RemoteBlast)
- next if (m{^<\/?PRE>});
# WU-BLAST -mformat 3 only
if(m{^#\s((?:\S+?)?BLAST[NPX])\s(\d+\.\d+.+\d{4}\])}) {
($alg, $ver) = ($1, $2);
@@ -183,6 +181,7 @@
next if /^#/ || /^\s*$/;
my @fields = split;
+ next if @fields == 1;
my ($qname,$hname, $percent_id, $hsp_len, $mismatches,$gapsm,
$qstart,$qend,$hstart,$hend,$evalue,$bits);
# WU-BLAST-specific
@@ -207,9 +206,6 @@
# we need total gaps in the alignment
$gapsm=$qgaps+$sgaps;
}
- else {
- $self->throw("Unknown BLAST tabular format");
- }
# Remember Jim's code is 0 based
if( defined $lastquery &&
More information about the Bioperl-guts-l
mailing list