Runs DA pipeline and makes it ready for enrichment analyses
annotatePeaks(cp, reference.genome, show.annotation.pie = FALSE, verbose)
bed formatted consensus peak matrix: CHR, START, STOP and raw peak counts (peaks by 3+samples)
genome of interested species. It should be 'hg38', 'hg19' or 'mm10'.
shows the annotation pie chart produced with ChipSeeker
prints messages through running the pipeline
DApeaks returns DA peaks