Sim, E. M.; Fong, W.; Suster, C.; Agius, J. E.; Chandra, S.; Suliman, B.; Wang, Q.; Ngo, C.; Finemore, C.; Chen, S. C.-A.; Basile, K.; Sintchenko, V.

Shiga Toxin (Stx) producing Escherichia coli (STEC) is a subset of pathogenic E. coli that can produce two types of Stx, Stx1 and Stx2, which can be further subtyped into four and 15 subtypes respectively. Not all subtypes, however, are equal in virulence potential, and the risk of severe disease including haemolytic uraemic syndrome has been linked to certain Stx2 subtypes e.g. Stx2a, Stx2d, highlighting the importance to survey stx subtypes. Previously, we developed a STEC virulence barcode to capture pertinent information on virulence genes to infer pathogenic potential. However, the process required multiple manual curation steps to determine the barcode. Here we introduce STECode, a bioinformatic tool to automate the STEC virulence barcode generation from sequencing reads or genomic assemblies. The development, and validation of STECode is described using a set of publicly available completed STEC genomes, along with their corresponding short reads. STECode was applied to interrogate the virulence landscape and molecular epidemiology of human STEC isolated during the period of the international border closures related to COVID-19 in the state of New South Wales, Australia.