I have 2 bam files which I produced by mapping the reads with bowtie using hg19. This is a CHIP seq data for one of the histone marks, where one file is input(control) and other file is the chip sample. I want to show for certain genes which I am interested in, read plots around the Transcription Start Sites(TSS). Let's say I have 6 genes of interest, so in one TSS plot(for that histone mark) I want to see how histone marks are represented for different genes.
How can I do it?