Result difference between blastall -p blastp (legacy blast) and blastp -task blastp (blast+)
0
5
Entering edit mode
6.8 years ago

I am new to bioinformatics

I am running pancoreplot in cmg biotools.

I am trying to replace BLASTALL with BLAST+.

Here is the command of BLASTALL:

blastall -F 0 -i a.fsa -p blastp -e 1e-5 -m 7 -d b.fsa | tigrcut | gawk '{print $1"\t"$2}' > blastresult_ab.txt

I have tried to replace this command with BLAST+ command as follows:

blastp -task blastp -query a.fsa -evalue 1e-5 -outfmt 5 -db b.fsa | tigrcut | gawk '{print $1"\t"$2}' >output_ab.txt

But the end results are wrong. I am unable to correctly calculate the core and pan genomes after replacing BLASTALL with BLAST+.

Any help regarding this would be very helpful.

Thanks in advance.

software-error blast • 1.8k views
ADD COMMENT

Login before adding your answer.

Traffic: 1670 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6