[Bioperl-l] Slightly OT - how to create a blast archive file?

[Smithies, Russell]

Slightly off-topic perhaps, but does anyone know how to create a blast archive file from existing blast output?
I have 25GB of compressed blastp output that I need to trim to the top 50 descriptions and alignments from each of 2,500,000 queries as the application that's going to post-process the data (MEGAN) can't handle the volume.
I thought the blast_formatter would be ideal but as far as I can see, it only accepts blast archive files :(

Plan B is SearchIO to read then write each file but I suspect it might take a while.

Any better suggestions?

Russell Smithies



=======================================================================
Attention: The information contained in this message and/or attachments
from AgResearch Limited is intended only for the persons or entities
to which it is addressed and may contain confidential and/or privileged
material. Any review, retransmission, dissemination or other use of, or
taking of any action in reliance upon, this information by persons or
entities other than the intended recipients is prohibited by AgResearch
Limited. If you have received this message in error, please notify the
sender immediately.
=======================================================================