blast in python
Hey guys, i want make a blast with two fasta sequences in python i am using this code:
import os
blastcommand = 'blastp -query a.fasta -subject b.fasta'
os.system(blastcommand)
the output is always "1"
how can i fix that?
• 1,403 views
•
link
updated
by
Wayne
•
written
by
brunofelicianodeomena •
•
written
by
brunofelicianodeomena •
1 answer
my guess, use
blastcommand = '/full/path/to/blastp -query /path/to/a.fasta -subject /path/to/b.fasta'
• 0 views
•
link
written
by
Pierre Lindenbaum
Log in to answer this question.
More posts like this
-
seaborn dispersion graph
written by brunofelicianodeomena •
-
make a heat map with python
written by brunofelicianodeomena •
-
keep getting error with psiblast
written by Mo
-
What's the most efficient way to load k-mers into dict in python?
written by Ray •
-
Variant Calling pipeline to be run parallel on multiple cores
written by AR •
-
blast using python
written by santhoshhegde278 •
-
Different results with BLAST depending on if subject is formatted database or FASTA-file
written by driliwyr •
-
NcbiblastnCommandline not detecting hits
written by horea.christ •
-
Protein Sequence Alignment Using Blast
written by p.s •
-
How To Force Blastp To Display All The Possible Combination Of Alignments?
written by Rein •
That's just the system exit code indicating an error occurred. To capture the output you need to write it to a file and parse it thereafter. See https://www.w3schools.com/python/ref_os_system.asp
I cover sending BLAST output to parseable files & parsing them with Python & Pandas in a series of Jupyter Notebooks available here. I even bring up alternatives like awk for parsing the output. I have a script that I use there that converts the output to a Pandas dataframe that is useful for sorting through data in a Python context. In the context of Jupyter, you can even skip the file output intermediates for the results. You can work through this without installing anything all by going there and clicking 'launch binder' badge near the top or click here to launch a session directly. Everything is installed and loaded in the temporary Jupyter session served via the MyBInder service to work through the examples in the available notebooks that come up listed when the session starts. I suggest working through the available notebooks in the order they are listed there. You'll probably primarily be interested in the earlier ones in the series for now.