I have a large amount of protein sequences in the .fasta format and I would like to extract only the amino acid mutations from these sequences, so that, in the end, I want to have a list that looks something like this: I456L, W675T, etc Is there a program or way to do this? Thankful
$\begingroup$
$\endgroup$
-
$\begingroup$ Are all of your sequences of the same protein, are they all the same length, do you know which one is the WT / baseline? $\endgroup$ – Nicolai Mar 24 '20 at 11:40
-
$\begingroup$ Yes, the proteins are of the same, in the same length and I have a reference genome to compare $\endgroup$ – DaudAssassin Mar 24 '20 at 11:49
-
2$\begingroup$ I would suggest posting this question on bioinformatics.stackexchange.com/questions. $\endgroup$ – Maximilian Press Mar 25 '20 at 19:20
$\begingroup$
$\endgroup$
I'd blastp, and parse the alignment you get back.
