length distribution was off by 1 due to pysam.aend used as 0-base while 1-based