def METHOD_NAME():
    tables = list(
        VcfReader(
            "tests/data/polyphasegenetic.test.parents.vcf",
            only_snvs=False,
            genotype_likelihoods=False,
            ploidy=4,
            mav=True,
        )
    )
    table = tables[0]
    param = PolyphaseGeneticParameter(4, 20, 0.06, 2, 0, True, True, False, "")
    vi = compute_phasable_variants(table, "Parent_A", "Parent_B", param)
    non_phasable = [33, 34, 36, 38, 96, 99, 106]
    phasable = vi.get_phasable()
    print(phasable)
    assert len(vi) == 135
    assert all([x not in phasable for x in non_phasable])
    assert [x for x in range(135) if x not in non_phasable] == phasable
    vi = compute_phasable_variants(table, "Parent_B", "Parent_A", param)
    phasable = vi.get_phasable()
    assert phasable == [10]