Blast in Biopython
0
0
Entering edit mode
7.9 years ago
midox ▴ 290

hello, I have two sequences and i need to use Blast under Biopython to align these two sequences and retrieve the results in list for example. But i can't do this. can you explain me using commande line please? thanks

Blast alignment biopython • 1.8k views
ADD COMMENT
0
Entering edit mode

Can you show us what you tried and what didn't work? People will be more eager to help you when you show some effort from your side.

ADD REPLY
0
Entering edit mode

here is a simple description of my program:

from Bio import SeqIO
import os
import sys
from Bio.Blast.Applications import NcbiblastnCommandline

com = NcbiblastnCommandline(query = "GGCTAGCTTT", subject = "GGCTTTCTTT", out = 'blast2')
myStdout, myStderr = com()

but here the problem is query and subject must a files. can i use directly a sequences? Thanks

ADD REPLY
0
Entering edit mode

I think this page is quite clear to what the input should be:

ADD REPLY
0
Entering edit mode

You posted earlier with an indication that you had already done this: calculate the e-value of alignment Since your question about e-values has not bee answered are you retrying to get an answer with this new question?

ADD REPLY
0
Entering edit mode

no it's not the same problem. the e-value problem is : i have created the SW matrix and i would determinate the e-value according to the values used in BLAST but the problem is my sequences is very long so the creation of SW matrix and recuperate the results is too long. Here, i try another solution using BLAST directly in Biopython, i would do an alignment between two sequences using BLAST under Biopython(standalone or in internet). if it's rapid so used this method but I do not know how to use BLAST under python. I try to use it.

ADD REPLY

Login before adding your answer.

Traffic: 1905 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