 Moments of the mutation count on a neutral coalescent tree

The independent mutation counts in coalescent epochs have the geometric mutation count in a coalescent epoch laws. Summing their expectations and variances gives the displayed formulas, where $H_m={\sum }_{i=1}^{m}1/i$ and $H_m^{(2)}={\sum }_{i=1}^{m}1/i^2$. The extra quadratic term records randomness of the ancestral tree, beyond the conditional Poisson distribution variance.