In a first step, we initialised simulations with one index case. For each primary case, we generated secondary cases according to a negative-binomial offspring distribution with mean R0 and dispersion k [7,8]. The dispersion parameter k quantifies the variability in the number of secondary cases, and can be interpreted as a measure of the impact of superspreading events (the lower the value of k, the higher the impact of superspreading). The generation time interval D was assumed to be gamma-distributed with a shape parameter of 2, and a mean that varied between 7 and 14 days. We explored a wide range of parameter combinations (Table) and ran 1,000 stochastic simulations for each individual combination. This corresponds to a total of 3.52 million one-index-case simulations that were run on UBELIX (http://www.id.unibe.ch/hpc), the high performance computing cluster at the University of Bern, Switzerland.