[Bioperl-l] Fix for Bio::Tools::Blat

Brian Osborne brian_osborne at cognia.com
Wed Nov 10 08:09:49 EST 2004


I've added your patch, please check it to see that I've put it in the proper

Pardon my ignorance but why do you use Tools::Blat rather than

Brian O.

-----Original Message-----
From: bioperl-l-bounces at portal.open-bio.org
[mailto:bioperl-l-bounces at portal.open-bio.org]On Behalf Of Neil Saunders
Sent: Wednesday, November 10, 2004 1:16 AM
To: bioperl-l at portal.open-bio.org
Subject: [Bioperl-l] Fix for Bio::Tools::Blat

dear all,

I mentioned in a previous post that Bio::Tools::Blat has a problem with
"tblastx-like" blat output (i.e. using the -t=dnax -q=dnax options).
This generates PSL files with 2 strands in the 'strand' column - one
for query, one for target (e.g. ++).  If you use the next_result()
method of Bio::Tools::Blat.pm on such files, you get errors like:

Argument "++" isn't numeric in numeric ne (!=) at
/usr/local/share/perl/5.8.4/Bio/Location/Atomic.pm line 170, <GEN0>
line 228.

I have added the following lines to my Blat.pm:

if($strand =~/([+-])[+-]/) {
$strand = $1;

which seems to fix it up.  Perhaps we could make this an 'official'

 School of Biotechnology and Biomolecular Sciences,
 The University of New South Wales,
 Sydney 2052,

Bioperl-l mailing list
Bioperl-l at portal.open-bio.org

More information about the Bioperl-l mailing list