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$ 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

I'd blastp, and parse the alignment you get back.


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.