This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Global alignment of protein sequences all x all to identify similarities between them

I have a fasta file with the protein sequence of 48 genes. I would like to compare them all against each other and see the similarity between them. For example, Gene1 against all 48, Gene2 against all 48. I would like to make a matrix with the similarity between the genes. I tried to use the EMBOSS Needle tool, with the following code:

needle -outfile result.needle -asequence all-proteins.fasta -bsequence all-proteins.fasta -gapopen 10 -gapextend 0.5

But what it does is take the sequence of the first protein and compare it against the 48. Is there a way to do all against all at once, without having to put: Gene1 x all-proteins, Gene2 x all-proteins?

And another issue I noticed was that in one of the cases I don't get the same similarity between the genes. Example: Gene1 x Gene2 (9.6% similarity), but Gene2 x Gene1 (9.7% similarity).

Is there any way to calculate the similarity of all against all so that I get a single result and can do it in one go?

Thanks

needle protein global-alignment emboss

0 answers

No answers yet.

Log in to answer this question.