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)
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...