This is a test version of Biostars. For the public version, visit https://www.biostars.org.
psiblast to get pssm. How to get a pssm from a multiple alignment ?

Hi,

I have problems using psiblast this being my command line : psiblast -subject test.fa -in_msa aligment_test.fa -out_ascii_pssm pssm_test1.txt -ignore_msa_master

I get this error message whenever I have more that 4 sequences : BLAST options error: Repeated Seq-IDs detected in multiple sequence alignment file, please ensure all Seq-IDs are unique before proceeding. I am 100% sure I have no repeats and try also adding false sequences to my alignment to check that it would always crash when I have more that 4 sequences.

Do you have any idea of how I could proceed to get PSSMs from multiple alignments ?

Thank you very much in advance for your help, have a nice day

pssm
• 2,497 views • link
updated by Biostar Biostar • written by lagartija lagartija

Please post an example of your input data.

• 0 views • link
written by Joe Joe

Just headers needed. So grep "^>" test.fa | head -5 would do.

• 0 views • link
written by GenoMax GenoMax

0 answers

No answers yet.

Log in to answer this question.

More posts like this