Question: Blastall or blast+ (conversion tool?)
0
gravatar for jomaco
4.9 years ago by
jomaco190
European Union
jomaco190 wrote:

Following on from AHRD functional annotation program (error), the functional annotation program AHRD (https://github.com/groupschoof/AHRD) appears to require the results from blastp searches of proteins to various databases (sprot, trembl etc.) to be in blastall pairwise format instead of blast+ pairwise format. 

Does it make a difference whether I use blastall or blast+? It seems the use of blast+ is recommended by NCBI - would the results be better using blast+?

Is there a script already available which can convert blast+ pairwise format files to blastall pairwise format files? (I have lots of blast+ generated files which I would like to use if possible).

Thanks

ahrd blast+ blastall pairwise • 1.8k views
ADD COMMENTlink modified 4.9 years ago by pld4.8k • written 4.9 years ago by jomaco190

Blast+ outputs 11 different formats + unlimited number of custom formats. Which particular format are you referring to? BLAST archive format (ASN.1)?

ADD REPLYlink modified 4.9 years ago • written 4.9 years ago by 5heikki8.5k

Hi, I am referring to the "pairwise" output format (output format zero) which appears to be slightly different between blast+ blastp (2.2.29+) and blastall blastp (2.2.18)

ADD REPLYlink written 4.9 years ago by jomaco190
0
gravatar for pld
4.9 years ago by
pld4.8k
United States
pld4.8k wrote:

As far as I know they didn't change the pairwise output format, but I could be wrong. At least if they did change it, the BLAST team doesn't seem to have documented it.

From the module in BioJava that is throwing the exception, it seems that an outdated version of BioJava will throw this exception when it encounters a version of BLAST it isn't sure about.

http://www.biojava.org/pipermail/biojava-l/2003-January/003427.html

You may want to see if you can update your version of BioJava, or if you're feeling brave you can set a lazy mode that won't throw an exception if BioJava can't tell which BLAST version you're using.

I'm sorta confused as to why AHRD would want to use the text pairwise output as your desired input format. XML, tabular and so on are much easier to read and much less sensitive to small changes that could show up between BLAST versions.

ADD COMMENTlink written 4.9 years ago by pld4.8k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1028 users visited in the last hour