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