Following on from AHRD functional annotation program (error), the functional annotation program 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
1 answer
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.
Log in to answer this question.
-
How to handle TrEMBL proteins without gene annotation in plasma proteomics?
written by Luwell •
-
comparison e-value from pairwise alignment to blastp
written by Olympia •
-
mRNA / Protein Annotation questions - etiquette and consensus
written by Biogeek
-
Blastp e-value, bit score, identify importance to identify homolog sequences
written by Mehmet
-
Blastp - through R
written by maahpishanu •
-
Blastp to know the name of the protein
written by utkarsh.sood •
-
Result difference between blastall -p blastp (legacy blast) and blastp -task blastp (blast+)
written by ahsan.chattha15 •
-
Blast2GO (csv output)
written by jomaco
-
AHRD functional annotation program (error)
written by jomaco
-
Some Questions About Using Orthomcl To Find Orthologs Within Many Species
written by User 7478 •
Blast+ outputs 11 different formats + unlimited number of custom formats. Which particular format are you referring to? BLAST archive format (ASN.1)?
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)