Dear all,

I have a list of genes, for which I have retrieved the UniProt protein sequences. I want to retrieve the sequence of particular mutations to the protein sequences, for example W155L, and save it in fasta format.

I could create my own function but I wonder if this can be done using out of the box Bioconductor functions or packages, which would be more robust and accommodate different kinds of mutations.

Any hints would be greatly appreciated.