def METHOD_NAME(tmp_path):
    """Ensure that the block cuts stay consistent with pre-phasing and that
    lowest sensitivity allows prephasings to bridge unconnected (by reads) blocks"""
    results = []
    for s in range(6):
        outvcf = tmp_path / "output{}.vcf".format(s)
        run_polyphase(
            phase_input_files=["tests/data/polyploid.cuts.bam"],
            variant_file="tests/data/polyploid.cuts.vcf",
            ploidy=4,
            ignore_read_groups=True,
            block_cut_sensitivity=s,
            output=outvcf,
            use_prephasing=True,
        )
        assert os.path.isfile(outvcf)
        tables = list(VcfReader(outvcf, phases=True))
        assert len(tables) == 1
        block_starts = set([i.block_id for i in tables[0].phases_of("Test_Cuts") if i is not None])
        results.append(block_starts)
        print(block_starts)
    for s in range(5):
        assert all(cut in results[s + 1] for cut in results[s])
    assert results[0] == {3}
    assert results[1] == {3, 18}
    assert results[2] in [{3, 9, 18, 27}, {3, 9, 18}]
    assert results[3] in [{3, 9, 18, 27}, {3, 9, 18}]
    assert results[4] == {3, 9, 18, 27}
    assert results[5] == {3, 9, 18, 27}