def METHOD_NAME():
    reads = """
    1001 01001
      1010010  01
       100 10 0010
        010  100
          001 100 01
    """
    am = AlleleMatrix(string_to_readset(reads))
    pos1 = list(range(16))
    pos2 = [0, 1, 9, 10, 13, 14]
    reads = [1, 2, 3]
    s1 = am.extractSubMatrix(pos1, reads, True)
    s2 = am.extractSubMatrix(pos2, reads, True)
    assert len(s1) == 3
    assert len(s2) == 2
    assert s1.getNumPositions() == 16
    assert s2.getNumPositions() == 6
    for i in range(len(s1)):
        for j in range(s1.getNumPositions()):
            assert am.getAllele(i + 1, j) == s1.getAllele(i, j)
    for i in range(len(s2)):
        for j in range(s2.getNumPositions()):
            assert am.getAllele(i + 2, pos2[j]) == s2.getAllele(i, j)