Consider incorporating a quartet sampling approach (substituting parsimony for likelihood): https://bsapubs.onlinelibrary.wiley.com/doi/10.1002/ajb2.1016