“…Sample sequences were identified using the Basic Local Alignment Search Tool (BLAST) (Mount, 2007) for the local similarity of 18S, 28S rDNA and cox1 from the nucleotide databases of The National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov). For other flatworms and outgroup species, nucleotide sequences were searched from these databases (accession numbers from Cribb et al., 2001; Lockyer et al., 2003; Van Steenkiste et al., 2013; Fraija-Fernandez et al., 2015; Vdacny et al., 2011; Laumer and Giribet, 2014; Ngamniyon et al., 2014; Tkach et al., 2003; Devi et al., 2010; Vdacny et al., 2012; Hoyal-Cuthill et al., 2016; Leung et al., 2009; Keeney et al., 2009; Ng et al., 2018). Sequences were aligned with the Clustal Omega multiple sequence alignment programme (http://www.genome.jp/tools/clustalw/), and phylogenetic trees inferred using the maximum-likelihood approach in MEGA 5.1 (Tamura et al., 2011).…”