This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How I can get percent identity matrix from clustalw or clustalo?

I have a fasta file containing multiple sequences.
I can't get percentage identity matrix, despite of the option I added.

clustalo --infile={my_fasta} -t DNA --percent-id -o {out_file}

clustalw -infile={my_fasta} -type=DNA -pim -outfile={out_file}

How can I get percentage identity matrix from this software in my server?

clustalo clustalw msa multiplesequencealignment pairwise
• 3,051 views • link
written by Jiun • Jiun

I also tried option: -pim=TRUE or -pim=True

• 0 views • link
written by Jiun • Jiun

1 answer

The option --percent-id will tell clustalo to calculate % identity instead of distance. Normally, it calculates the distance as [1 - % identity]. You still need to specify where the matrix will be saved (the -o switch specifies where to save the alignment), and also have to force full alignment.

clustalo --infile={my_fasta} -t DNA --percent-id --full --distmat-out={matrix_file} -o {out_file}
• 0 views • link
written by Mensur Dlakic Mensur Dlakic

Finally I can get the matrix of percent identity. Thank you so much!! :)

• 0 views • link
written by Jiun • Jiun

Log in to answer this question.

More posts like this