For example, if I have a file with 3 sequences:
>1 gtccgatgat >2 gatggatcggacgttagcaa >3 acatttgaga
I want to generate sequences with 10% of the nucleotides randomly mutated.
>1 gtcTgatgat >2 gaAggatcggacgtGagcaa >3 acatCtgaga
Does anybody know a program or code to do that?
And if you know a way to do this in python would be great!