[Bioperl-l] timing out a blast in StandAloneBlast.pm

Peter Kos kos@rite.or.jp" <kos@rite.or.jp
Wed, 3 Jul 2002 02:26:54 +0900


Hi Bonnie,

Wait a minute. This is a completely different issue.
I usually suggest first playing with the parameters, since people 
tend to forget about them and encounter problems, which can be easily 
solved this way.
But if the problem is only the memory limitation of the list of HSPs, 
why don't you just chop up the genome to appropriately sized and 
perhaps appropriately overlapping slices for running blast on them 
separately? Then by simply adjusting the coordinates you may generate 
exactly the same fragmented list that you wanted, and you can 
collect/store it in the separate file for a final sorting and later 
reprocessing.
(The same thing, just in another dimension/direction/axis/whatever in 
the time-sequence-memory space.)

It may not be as elegant as adding Megablast to StandAloneBlast.pm, 
but at least you may keep all the possible HSPs that you originally 
wanted. And, moreover, it seems to be more simple, so I would be able 
to do this, while I am not good enough to mess around in the modules. 
You can be much better, of course.

Best regards
Peter

> Thank you so much for your response.  I completely agree with you.
>  Since
> BLAST is not retrieving HSPs in any particular order timing out a
> BLAST
> does not make scientific sense.  Our problem is that our linux 
boxes
> seem
> to crash on sequences that have a large number of HSPs (yes, we do
> need
> better hardware...).  My thought was to time these out and capture
> them in
> a separate file for more "specialized" processing later.  But, I
> think
> perhaps a better solution is to play with the BLAST parameters to
> allow
> less HSPs through and use this set of parameters for the whole set,
> as you
> suggested.  Unfortunately, BLAST doesn't have options like 
Megablast
> does
> for limiting searches on "-p percent identity" and "-s score", so I
> am
> working on adding Megablast to the existing StandAloneBlast.pm
> program,
> which will hopefully help me to limit the HSPs a little better.
>  After
> looking at the code for StandAloneBlast.pm it is just a wrapper for
> BLAST,
> so there is no magic going on there that isn't in the regular BLAST
> program.  Thank you for your all of your help!
>
> Kind Regards,
> Bonnie