Takes as input a file with assembled reads from a single diploid individual and returns maximum likelihood estimates of the population mutation rate, , the sequencing error , the zygosity correlation, and the population recombination rate.