psiblast to get pssm. How to get a pssm from a multiple alignment ?
0
0
Entering edit mode
2.5 years ago
lagartija ▴ 130

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 • 1.3k views
0
Entering edit mode

0
Entering edit mode

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