Check for parsimony before looking for equivalent shift edges.
Check that "topologically equivalent" edges allocation produce same distribution (function qr.solve_exact).