def log_prob_correct_from_qual(q): """Computes the probability of no error given a phred quality.""" return np.log1p(- 10**(-0.1 * q))