HMMtBroadPeak()
Call broad peak by Hidden Markov Models with t emissions
check_identical_header()
helper function to check the bam file header
plotPeaks()
Plot tool for called peaks