Question: Psi-Blast Commandline Version Using Biopython
gravatar for Molly
9.4 years ago by
Molly30 wrote:

Hi Niek,

I have been trying to get my program for PSI-Blast run, i have been using the wrapper from Bio.Blast.Applications.NcbipsiblastCommanline: I saw the code you posted at

i tried to use the same but it doesnt seem to work for me. can you suggest me why?

code that i used:

    psi_cline = NcbipsiblastCommandline('psiblast', db = 'refseq_protein',\
                                        query = queryID+".fasta", evalue =  10 , \
                                        out = queryID+"_psi.xml", outfmt = 7, \
                                        out_pssm = queryID+"_pssm")
    #p = subprocess.Popen(str(psi_cline),stdout=subprocess.PIPE,\
                        # stderr=subprocess.PIPE,shell=(sys.platform!="win32"))

here is the error that i get : Traceback (most recent call last): File "", line 116, in <module> psicline() File "/usr/local/lib/python2.6/dist-packages/biopython-1.57-py2.6-linux-x8664.egg/Bio/Application/", line 432, in call stdoutstr, stderrstr) Bio.Application.ApplicationError: Command 'psiblast -out NP012649psi.xml -outfmt 7 -query NP012649.fasta -db refseqprotein -evalue 10 -outpssm NP012649_pssm' returned non-zero exit status 127, '/bin/sh: psiblast: not found'

i do not get any XML as my output as the command is not found apparently. in this case queryID is a protein name and the queryID.fasta consists of the fasta file.

looking forward for your reply. Molly

biopython • 4.3k views
ADD COMMENTlink modified 9.4 years ago • written 9.4 years ago by Molly30

Hi Molly! Please don't cross-post on the BioPython Mailing List and here. Reference:

ADD REPLYlink modified 14 months ago by RamRS30k • written 9.4 years ago by Michael Schubert7.0k
gravatar for Peter
9.4 years ago by
Scotland, UK
Peter5.8k wrote:

As Michael pointed out, you've already been discussing this on the Biopython mailing list where people have tried to explained the error message probably means you have not installed the BLAST+ command line tools (in this case the binary psiblast), or the binaries are not on the Unix PATH (in which case update the PATH environment variable, or use the explicit full path to the binary).

ADD COMMENTlink modified 14 months ago by RamRS30k • written 9.4 years ago by Peter5.8k
gravatar for Molly
9.3 years ago by
Molly30 wrote:

Sorry people for the cross question. but i posted it here when nobody was replying me :( anyways, my problem is solved now.

Thanks for the help.

ADD COMMENTlink written 9.3 years ago by Molly30
Please log in to add an answer.


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