DISOPRED3 Blastpgp Failed 256

Question

[euduca]

Hi, all =)

I'm trying to use de DISOPRED3 program, but this error appears:

$  ../run_disopred.pl ../examples/example.fasta 
Running PSI-BLAST search ...

[../run_disopred.pl] ERROR: /home/eac16/local/ncbi/blast/legacy/2.2.26/bin/blastpgp -i /home/eac16/local/DISOPRED/3.16/examples/example.fasta -d /scratch/eac16/dbs/blast/legacy/uniprot/uniref90 -a 2 -b 0 -j 3 -h 0.001 -C /home/eac16/local/DISOPRED/3.16/examples/example_16870_8d0a0101.chk -o /home/eac16/local/DISOPRED/3.16/examples/example_16870_8d0a0101.blast 
 **failed: 256**

The DISOPRED3 uses the legacy blast, so I used the command "formatdb" (2.2.26 version ) to create the database.

$ formatdb -i uniref90.fasta -p T -o T

The paths are correct and the variables at the run_disopred.pl's script too:

my $NCBI_DIR = "/home/eac16/local/ncbi/blast/legacy/2.2.26/bin/";
my $SEQ_DB = "/scratch/eac16/dbs/blast/legacy/uniprot/uniref90";
my $PSIBLAST_PAR = "-a 1 -b 0 -j 3 -h 0.001";

I installed and used the DISOPRED3 in a HPC once and use to predict (proteins) an entire genome and it was done. But now I moved to another lab, and I don't know what is wrong with the software/instalation.

Can someone help me?

Thanks. Eduardo.

[fishgolden]

Try blastpgp with the same command. Then blast pgp will give you more information.

/home/eac16/local/ncbi/blast/legacy/2.2.26/bin/blastpgp -i /home/eac16/local/DISOPRED/3.16/examples/example.fasta -d /scratch/eac16/dbs/blast/legacy/uniprot/uniref90 -a 2 -b 0 -j 3 -h 0.001 -C /home/eac16/local/DISOPRED/3.16/examples/example_16870_8d0a0101.chk -o /home/eac16/local/DISOPRED/3.16/examples/example_16870_8d0a0101.blast

However, I'm afraid that

my $SEQ_DB = "/scratch/eac16/dbs/blast/legacy/uniprot/uniref90";

should be

my $SEQ_DB = "/scratch/eac16/dbs/blast/legacy/uniprot/uniref90.fasta";

but not checked.

[euduca]

Hi Yamule. You're right about the variable $SEQ_DB. I fixed this.

But the error still persisted, so I figured out that variable $BLASTDB was pointed, by default, to the NR database in the HPC. So, before to run the script I set the variable $BLASTDB to:

export BLASTDB=/scratch/eac16/dbs/blast/legacy/uniprot/uniref90

After this the blastpgp worked properly.

I wrongly presumed that variable $SEQ_DB would be enough to set the database.

Thanks for your help.