WebHopper
Posts: 7 Joined: Feb. 2017

Wesley, if I understand you well the expected correct bases after mutation is
E = C + (u * (L  C) / K)  (u * C * (K  1) / K)
where C = number of correct bases K = alphabet length L = number of base pairs u = mutation rate
I guess with "expected" E you refer to the mean. So if C=L we have
E = L (K + u  K u)/K
This is a linear law with respect to L if the other parameters are held constant. However, this does not match a Monte Carlo simulation, which shows that E is exponential with respect to L. So there must be some error in your calculation or I don't understand you well.
DiEb: Thanks for the link, I will look into this. Seems to be complicated however...
