In order to thoroughly examine the population and temporal structure of Bhutanese S. sonnei and the population from which they emerged, all sequences from Bhutan, including identical sequences removed in the previous phylogeny, were aligned with the six most closely related isolates from other countries and trimmed down to include only the 544 SNPs present in the alignment of 77 taxa. These data were then utilized for maximum-likelihood reconstruction using PhyML under the TVM+Γ4 substitution model (Guindon et al., 2010). In order to investigate whether temporal signal was embedded in the resulting maximum-likelihood phylogeny, we applied Path-O-Gen (version 1.4) to assess the linear relationship between root-to-tip divergence and date of isolation (in day, month, year).