script to plot read counts around tss in bam file