scripts for bam analysis