def test_compute_quals_numerical_stability(self, probs, call, expected_gq): max_qual = round( genomics_math.ptrue_to_bounded_phred(1.0), postprocess_variants._QUAL_PRECISION) gq, qual = postprocess_variants.compute_quals(probs, call) self.assertEquals(expected_gq, gq) self.assertEquals(max_qual, qual)
def test_compute_quals(self, probs, call, expected_gq, expected_qual): gq, qual = postprocess_variants.compute_quals(probs, call) self.assertEquals(gq, expected_gq) self.assertAlmostEquals(qual, expected_qual, places=6)