Nonmetric multidimensional scaling (NMDS) biplot of a Robinson–Foulds symmetric absolute distance matrix. Ordination of topological differences observed between the topologies obtained under different types of searches and inferences used in this study, after rotation procedures using principal component analysis. The relationships among topologies were forced to be presented into only two dimensions, so that the procedure preserved the distances among topologies after iteratively trying to position the objects in the space in such a way to minimise a stress function (stress = 0.0314). The more different topologies found under different analyses, the higher the distance among the points on NMDS biplot. Abbreviations: A1, analysis 1; A2, analysis 2; A3, analysis 3; A4, analysis 4; A5, analysis 5; A6, analysis 6; SPR0, analysis with a combination of continuous characters producing a Subtree Pruning Regrafting distance equal to zero; SPR2, analysis with a combination of continuous characters producing the Subtree Pruning Regrafting distance equal to two; SPR6, analysis with a combination of continuous characters producing the Subtree Pruning Regrafting distance equal to six; MLDisc, maximum likelihood inference performed on the matrix of only discrete characters; BAYESDisc, Bayesian inference performed on the matrix with only discrete characters; BAYESDiscCont, Bayesian inference performed on the matrix with all datasets combined.