I am looking for a way to take a pairwise alignment and output the amino acid modification. For example, the sequences could be:
>seq1 MYPROT >seq2 MYYROT
The desired output would be: P3Y
I have hundreds of protein sequences so I'd like something high-throughput because doing this by eye is tedious and error-prone.
Any help is greatly appreciated!