DNA was then processed as previously described (41 (link),42 (link)). DNA libraries were validated and processed with Illumina cBot for cluster generation on the flowcell, following the manufacturer's instructions. The libraries were sequenced on single-end 50 bp mode on HiSeq2500 (Illumina) at a depth of ∼30–50 million sequences per sample. CASAVA 1.8.2 version (Illumina pipeline) was used to processed raw data for both format conversion and de-multiplexing.