EstimateMutRates function

Estimate mutation rates