def METHOD_NAME():
    reads = """
    1 101 01
      01   01
       00 1001
        1 00 1
    """
    am = AlleleMatrix(string_to_readset(reads))
    ad = [[0, 1], [1, 1], [2, 1], [1, 2], [2, 1], [3, 1], [1, 1], [0, 2]]
    for i in range(am.getNumPositions()):
        assert am.getAlleleDepths(i) == ad[i]