Base-wise coverage depth calculations (Table S4a in Additional file 2) were performed using collapsed alignment (forward and reverse reads merged). Amplicon-wise coverage depth calculations (Table S4b in Additional file 2) were performed using reads aligning to each specific amplicon (based on start site coordinate). This way, amplicons that overlap are kept separate for the coverage calculation and do not artificially inflate the coverage of the overlapping bases.
