WGS was performed as detailed elsewhere (Perez-Lago et al., 2014 (link)) using a HiSeq 2000 device and a Miseq device (Illumina), which generated 101-51–bp paired-end reads. We mapped the reads for each strain using the Burrows-Wheeler Aligner and the ancestral MTB genome as detailed elsewhere (Comas et al., 2013 (link)). SNP calls were made with SAMtools and VarScan (coverage of at least 10×, mean SNP mapping quality of 20). The genome was compared between strains using an in-house script written in R.
Free full text: Click here