## 3 Results
### Phylogeny  
The number of sequences obtained varied by morphospecies and gene; at least two morphospecies were covered for each gene except for _tuba1_. Only one _tuba1_ sequence was obtained from \(G\). nr. _gilesi_ and \(G\). sp. 2. The final concatenated dataset, pruned of potential rogue taxa, comprised 39 individuals across the seven morphospecies and four outgroup taxa in the Labiduridae and Pygidicranidae families (Table 4). The concatenated ML gene tree resolved all Anisolabididae morphospecies to clades with 100% bootstrap support, and so the tree is presented with morphospecies-internal branches reduced for clarity, opposite the *BEAST species tree (Figure 4). The one exception is \(G\). sp. 3 which, although reliably split from \(G\). sp. 2, has low bootstrap support for the node splitting the two individuals. Branch lengths are not shown to emphasise topology. Outgroups have been pruned in the final figures. The ML and Bayesian estimates are mostly identical; \(A\). sp. 1 is sister to all other taxa, followed by \(G\). _nr. gilesi_ (Figure 4). In both trees, \(G\). _forcipata_, \(G\). sp. 2, and \(G\). sp. 3 form a clade with high support (Figure 4) (referred to henceforth as the SA clade). The *BEAST tree resolves \(G\). sp. 1 and \(A\). sp. 2 to a highly supported clade, sister to the SA clade (Figure 4). The ML tree places \(G\). sp. 1 and sister to the SA clade and \(A\). sp. 2, and \(A\). sp. 2 as sister to the SA clade, but with low bootstrap support (Figure 4). While not annotated on the tree, the branches preceding the nodes with the lowest bootstrap support are also some of the shortest in the tree, both at 0.028. \begin{table}
\begin{tabular}{c c c c c c} \hline \hline
**Individual** & **Morphospecies** & **Family** & _cox1_ & _cob_ & _28S_ & _tuba1_ \\ \hline  
K399681 \\ \hline
308 & & Anisolabididae & MK399564 & MK399456 & MK399647 & MK399677 \\
309 & _Gonolabis_ nr. _gilesi_ & Anisolabididae & MK399565 & MK399457 & MK399468 & - \\
484 & & Anisolabididae & MK399610 & MK399471 & MK399658 & - \\ \hline
146 & & Anisolabididae & MK399500 & MK399433 & MK399632 & - \\
152 & & Anisolabididae & MK399501 & MK39944 & MK399633 & - \\
384 & & Anisolabididae & MK399508 & MK399467 & MK399655 & MK399682 \\
385 & & Anisolabididae & MK399599 & MK399468 & MK399656 & MK3
Regardless, the generic assignments following Steinmann [6] are paraphyletic and polyphyletic for _Gonolabis_ and _Anisolabis_ respectively according to the molecular phylogeny. The log-likelihoods of the best-known ML trees unconstrained and constrained by the monophylies of _Gonolabis_ and _Anisolabis_ were \(-8205.895\) and \(-8216.307\), respectively. The Shimodairao-Hasegawa test did not indicate that the constrained tree was significantly worse (\(p>0.05\)), although plotting it did reveal that the constraint caused many of the branches to reduce to almost zero. This would raise the constrained tree's \begin{table}
\begin{tabular}{c c c c c c} \hline \hline
**Individual** & **Morphospecies** & **Family** & _cox1_ & _cob_ & **28S** & _tuba1_ \\ \hline
232 & & Anisolabididae & MK399531 & MK399443 & MK399639 & - \\
236 & _Gonolabis_ sp. 2 & Anisolabididae & MK3995532 & MK399444 & MK399640 & MK399669 \\
239 & & Anisolabididae & MK399536 & MK399445 & MK399641 & - \\
274 & & Anisolabididae & MK399546 & MK399450 & - & - \\ \hline
251 & _Gonolabis_ sp. 3 & Anisolabididae & MK399539 & MK399446 & MK399642 & MK399670 \\
276 & _Gonolabis_ sp. 3 & Anisolabididae & MK399547 & MK399451 & MK399645 & MK399673 \\ \hline C\_f\_Betch & _Challia fletcher_ & Pygidciranidae & NC\_018538 & - & - \\ L\_ripar & _Labidura rijvaria_ & Labiduridae & KX06089 & - & AB119553 & KX069030 \\ N\_tenui & _Nala tunicornis_ & Labiduridae & KX060909 & - & - & KX069055 \\ N\_livid & _Nala tividipes_ & Labiduridae & KX060609 & - & - & KX069048 \\ \hline \hline \end{tabular}
\end{table} 
Table 4: _Cont._  
Figure 4: Results of Bayesian and maximum-likelihood phylogenetic analysis estimating the species tree of the Anisolabididae morphospecies.