The mitochondrial free radical theory of aging (mFRTA) implicates Reactive Oxygen Species (ROS)-induced mutations of mitochondrial DNA (mtDNA) as a major cause of aging. However, fifty years after its inception, several of its premises are intensely debated. Much of this uncertainty is due to the large range of values in the reported experimental data, for example on oxidative damage and mutational burden in mtDNA. This is in part due to limitations with available measurement technologies. Here we show that sample preparations in some assays necessitating high dilution of DNA (single molecule level) may introduce significant statistical variability. Adding to this complexity is the intrinsically stochastic nature of cellular processes, which manifests in cells from the same tissue harboring varying mutation load. In conjunction, these random elements make the determination of the underlying mutation dynamics extremely challenging. Our in silico stochastic study reveals the effect of coupling the experimental variability and the intrinsic stochasticity of aging process in some of the reported experimental data. We also show that the stochastic nature of a de novo point mutation generated during embryonic development is a major contributor of different mutation burdens in the individuals of mouse population. Analysis of simulation results leads to several new insights on the relevance of mutation stochasticity in the context of dividing tissues and the plausibility of ROS vicious cycle hypothesis.