Nextera XT libraries were prepared for each of the samples. The barcoded libraries were pooled into equimolar concentrations following manufacturer’s guidelines (Illumina, San Diego, CA) using the Mid-Output Kit for paired-end sequencing (2×150 bp) on an Illumina NextSeq500 sequencing platform. Raw data (mean virtual coverage 361x) was demultiplexed and subsequently clipped of adapters using Trimmomatic v0.38 with default parameters46 (link). Quality control passing read-pairs were aligned against reference genome/plasmids (Accession numbers: NC_016810.1, NC_017718.1, NC_017719.1, NC_017720.1) with bwa v0.7.1747 (link). Genomic variant were called using Pilon v1.23 48 (link). with the following parameters: (i) minimum coverage 10x; (ii) minimum quality score = 20; (iii) minimum read mapping quality = 10. SnpEff v4.3 was used to annotate variants according to NCBI and predict their effect on genes49 (link).