Illumina MiSeq paired end sequencing
The concentration of the amplicons was estimated using a Gel Doc� XR+ System coupled with Image Lab� Software (BioRad, Hercules, CA USA) with 3 �l of O'GeneRulerTM 100 bp Plus DNA Ladder (Thermo Fisher Scientific, Inc., Waltham, MA, USA) as the internal standard for band intensity measurement. The samples of individual gels were pooled into approximately equimolar subpools as indicated by band intensity and measured with the Qubit dsDNA br Assay Kit (Life Technologies GmbH, Darmstadt, Germany). Subpools were mixed in an equimolar fashion and stored at -20�C until sequencing