Reputation: 21
So I'm following the NCBI instructions available here: https://www.ncbi.nlm.nih.gov/books/NBK52640/
and I can't for the life of me understand what's wrong here.
And here's my Python code:
from Bio.Blast.Applications import NcbipsiblastCommandline
import subprocess
psi_cline = NcbipsiblastCommandline('psiblast', db = 'refseq_protein.00',\
query = "results.fasta", evalue = 10 , \
out = "out_psi.xml", outfmt = 7, \
out_pssm ="pssm-results_pssm")
print(psi_cline)
Upvotes: 0
Views: 5888
Reputation: 21
I had this very problem today and it turns out that my blast-database was corrupted in some manner. I recreated my database with makeblastdb and this error went away.
Upvotes: 2