About Blast Scores In Matrix Format.
3
0
Entering edit mode
13.1 years ago
Azam ▴ 60

Hi All, How Can I get local Blast score in matrix format??( I have score but I need in Matrix format) How can I run blast of 200 sequences in order to alignment sequences by writing a python programme.( comparing one by one all sequences one another by blast). regards, Azam

matrix blast • 3.9k views
ADD COMMENT
1
Entering edit mode
13.1 years ago

edit, to expand on Chris' idea:

You can perform a BLAST search on a multi-FASTA file on itself.

blastp -query my.fasta -subject my.fasta -outfmt 7 # or 6

which yields all pairwise scores that you just need to parse from, e.g. tabular file format or XML using BioPython.

-- old answer, which might still be the cleaner way to do it --

I suggest you look at the following resources:

And, most importantly, show us what you have tried already and where you got stuck. We won't write the program for you, you will have to do this yourself.

As for the "score in Matrix format": a score is a scalar (one value) and a matrix consists of n*m values. Are you asking for a matrix of all pairwise scores?

ADD COMMENT
1
Entering edit mode
13.1 years ago

You can also put all sequences in a single file and blast every individual sequence against that file. (Although you will do every comparison twice in that way).

ADD COMMENT
1
Entering edit mode

I just tested, you can actually perform an all-against-all with just one BLAST query (see my edit).

ADD REPLY
0
Entering edit mode

Have you done it before? It has same results by all align by all??

ADD REPLY
0
Entering edit mode
13.1 years ago

If your data looks like this:

column 1 : seq 1
column 2 : seq 2
column 3 : Score

then, you'll find many examples on the web by searching for "pivot table". For example see: http://www.unix.com/shell-programming-scripting/155448-create-pivot-table.html

To call blast from python: http://stackoverflow.com/search?q=python++call+program+external

ADD COMMENT

Login before adding your answer.

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