“…However, as demonstrated here and elsewhere, variant calling pipelines and the various mapping/alignment, SNP calling, and SNP filtering practices that they employ (e.g., removal of recombination and clustered SNPs) can influence the identification of SNPs in WGS data and, thus, the topology of a resulting phylogeny (Pightling et al, 2014, 2015; Croucher et al, 2015; Hwang et al, 2015; Katz et al, 2017; Sandmann et al, 2017). This can be particularly problematic for outbreak and cluster detection in bacterial pathogen surveillance: pairwise SNP thresholds are currently widely used to make initial decisions regarding the inclusion or exclusion of isolates in a given outbreak (Taylor et al, 2015; Gymoese et al, 2017; Mair-Jenkins et al, 2017; McCloskey and Poon, 2017; Walker et al, 2018). In such scenarios, just a few SNPs can be the deciding factor in whether a bacterial pathogen is included or excluded as part of an outbreak or cluster (Katz et al, 2017), rendering the choice of variant calling method as non-trivial.…”