“…Supplemental to traditional clade support metrics, we used classical multidimensional scaling (MDS) to plot treespace (Gower, 1966; Hillis, Heath and St. John, 2005; St. John, 2017; Wright and Lloyd, 2020), with the goal of identifying the distribution of trees resolving key clades formed with KUMIP 314087 ( Table 1 ). Our R script inputs the unrooted post-burnin posterior samples (resultant from BI) and MPTs (resultant from MP) using ape v.5.3 (Paradis and Schliep, 2019), and employs phangorn v.2.5.5 (Schliep, 2011) to calculate pairwise unweighted Robinson-Foulds distances (RF, the proportion of bipartitions defined by a branch in one tree that is lacking in another tree) (Robinson and Foulds, 1981) for the total set of trees resulting from all analyses.…”