A Python utility for calculating statistics against genome positions based on sequence alignments from a SAM or BAM file.