[Bioperl-l] updating a reciprocal blast file
jvaughn7 at utk.edu
Thu Jun 5 11:53:54 EDT 2008
I have a large reciprocal blast file that contains 3 proteomes. I'd like to
integrate another proteome for downstream clustering. I imagine a
command-line script that takes as input the new proteome in fasta format,
the directory of the the old proteomes in fasta format, and the pre-existing
reciprocal blast file, and then performs the proper blasts and updates the
pre-existing reciprocal blast file accordingly. I am using blast locally
and the downstream processing is done by OrthoMCL. I assume this has been
handled before, but I can't track down the code. If anyone is familiar with
a pre-exisiting script or has pertinent advice, I'd be much obliged.
View this message in context: http://www.nabble.com/updating-a-reciprocal-blast-file-tp17673277p17673277.html
Sent from the Perl - Bioperl-L mailing list archive at Nabble.com.
More information about the Bioperl-l