def getNumberOfAlignments(bamfile): ''' return number of alignments in bamfile. ''' samfile = pysam.AlignmentFile(bamfile) return samfile.mapped