I'm able to generate a PSSM and a sequence from combinatorial experimental data. I want to BLAST the results to find proteins displaying these motifs. However, I have information about each position's sensitivity to substitutions, in addition to a consensus sequence. I know how to create a PSSM from this information, but I don't know how to incorporate this into a BLAST search. Does anyone know a cmdline argument to blast or a way of hacking a PSSM into a PSI-BLAST search to add position information to a sequence?
Question: Applying A Pssm To Psi-Blast
1
white.d.andrew • 10 wrote:
ADD COMMENT
• link
•
modified 8.6 years ago
by
Leonor Palmeira ♦ 3.8k
•
written
8.6 years ago by
white.d.andrew • 10
0
Leonor Palmeira ♦ 3.8k wrote:
You might find some information in this post : http://www.biostars.org/post/show/1699/how-to-scale-a-pssm-for-ncbi-blast/
As well as here : http://etutorials.org/Misc/blast/Part+V+BLAST+Reference/Chapter+13.+NCBI-BLAST+Reference/13.8+blastpgp+Parameters+PSI-BLAST+and+PHI-BLAST/ and here http://gogarten.uconn.edu/MCB372/Laboratories/assign3.html
What you probably need to do:
- generate a .chk checkpoint file similar to the one generated by 'blastpgp' containing your PSSM
feed it to blastpgp with something similar to:
blastpgp -i test1.fa -d genome.faa -R test1.chk -o genome.test1.br -I T -a 2
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: 990 users visited in the last hour