def METHOD_NAME():
    reads = """
    1001 01001
      1010010  01
       100 10 0010
        010  100
          001 100 01
    """
    am = AlleleMatrix(string_to_readset(reads))
    s1 = am.extractInterval(0, 4, True)
    s2 = am.extractInterval(0, 4, False)
    s3 = am.extractInterval(10, 16, True)
    s4 = am.extractInterval(10, 16, False)
    assert len(am) == len(s2) == len(s4)
    assert len(s1) == 3
    assert len(s3) == 4
    assert s1.getNumPositions() == s2.getNumPositions() == 4
    assert s3.getNumPositions() == s4.getNumPositions() == 6
    for i in range(len(s3)):
        for j in range(s3.getNumPositions()):
            assert am.getAllele(i + 1, j + 10) == s3.getAllele(i, j)
            assert am.getAlleleDepths(j + 10) == s3.getAlleleDepths(j)
            assert am.localToGlobal(j + 10) == s3.localToGlobal(j)
    assert s4.getRead(0) == []
    assert s2.getRead(3) == s2.getRead(4) == []