A script to calculate reads per kilobase from bam file