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