Stochastic simulation of mitochondrial heteroplasmy dynamics