I just wanted to ask what parameters did you use for the clustal alignment? As I can see from the input sequences, it aligned according to the input sequences. But how to optimize it for aligning as asked in the question, i.e. "align multiple phosphorylated peptides around their phosphosites "