[Bioperl-l] Bio::Tools::Run::StandAloneBlastPlus - Illegal seek error running blastn

[Mark A. Jensen]

Mike - please try updating your bioperl-live (the core) to the latest code 
(revision 16761 or so).
CommandExts is a work in progress; from the stack errors it looks like you've 
got an older version.
Try it then ping us back, if you would--
Thanks
Mark
----- Original Message ----- 
From: "mike stubbington (BI)" <mike.stubbington at bbsrc.ac.uk>
To: <bioperl-l at lists.open-bio.org>
Sent: Thursday, January 28, 2010 10:41 AM
Subject: [Bioperl-l] Bio::Tools::Run::StandAloneBlastPlus - Illegal seek error 
running blastn


Dear all,

I am attempting to blast some primers against the mouse genome. I have created a 
local mouse genome blast database and I can search against it using 'blastn' at 
the command line.

I have perl code that creates an array of bioperl sequence objects called 
@primers

I then create a StandAloneBlastPlus factory using the following code…

my $blastFactory = Bio::Tools::Run::StandAloneBlastPlus->new(
-db_dir => '/Users/stubbing/localBlast/',
-db_name => 'MouseGenome'
);

and then attempt to blast my primers using this…

my @shortPrimers;
my $count=1;
foreach (@primers) {
my $currentSeq = $_;
print "Checking primer $count/$primerNumber ";
if ($_->length < 40) {
push(@shortPrimers,$_);
print "Too short!\n";
}
else {
print "BLASTing...";
my $blastResult = $blastFactory->blastn(-query => $currentSeq);
}
$count++;
}

This fails with the following error…

------------- EXCEPTION -------------
MSG: /usr/local/ncbi/blast/bin/blastn call crashed: There was a problem running 
/usr/local/ncbi/blast/bin/blastn : Illegal seek at 
/Library/Perl/5.10.0/Bio/Tools/Run/WrapperBase/CommandExts.pm line 989, <DATA> 
line 532.

STACK Bio::Tools::Run::WrapperBase::_run 
/Library/Perl/5.10.0/Bio/Tools/Run/WrapperBase/CommandExts.pm:994
STACK Bio::Tools::Run::StandAloneBlastPlus::AUTOLOAD 
/Library/Perl/5.10.0/Bio/Tools/Run/StandAloneBlastPlus.pm:1236
STACK Bio::Tools::Run::StandAloneBlastPlus::run 
/Library/Perl/5.10.0/Bio/Tools/Run/StandAloneBlastPlus/BlastMethods.pm:267
STACK Bio::Tools::Run::StandAloneBlastPlus::AUTOLOAD 
/Library/Perl/5.10.0/Bio/Tools/Run/StandAloneBlastPlus.pm:1233
STACK toplevel ./5CTest.pl:63
-------------------------------------

Line 63 in my code is (as you might expect) the one that calls blastn on my 
factory object.

I'd appreciate any help you might be able to provide to shed light on this.

Thanks in advance,

Mike



_______________________________________________
Bioperl-l mailing list
Bioperl-l at lists.open-bio.org
http://lists.open-bio.org/mailman/listinfo/bioperl-l