def expected_count(af, effective_mutation_rate):
    """Calculate the expected number of somatic variants
       greater than a given allele frequency given an effective mutation
       rate, according to the model of Williams et al. Nature 
       Genetics 2016"""
    return effective_mutation_rate * (1.0 / af - 1.0)