Compare phred scores to phred threshold, not MAPQ threshold! Also removed a comment. The answer to the question asked there is that we retrieve a set of alignments that overlap a region, but if any of them overhang it then the resulting positions that we iterate over may not always be in the region.