more work on unconditional sampling of substitution histories on a tree