The BioNumerics software was used to generate a minimum spanning tree (MST) or a UPGMA hierarchical clustering as described previously [3 (link)]. The categorical coefficient was used to calculate the MST and the MST was based on in-house E. coli and K. pneumoniae wgMLST schemes. The NGS data of the K. pneumoniae and E. coli isolates were used for classical MLST and wgMLST analyses using in-house wgMLST schemes made in the SeqSphere software version 6.0.2 (Ridom). The in-house K. pneumoniae wgMLST scheme comprised 4978 genes (3471 core-genome and 1507 accessory-genome targets) using K. pneumoniae MGH 78,578 (NC_009648.1) as a reference genome. The in-house E. coli wgMLST scheme comprised 4503 genes (3199 core-genome and 1304 accessory-genome targets) using E. coli 536 (CP000247.1) as a reference genome.