#To run this code

1. Download the code/utils script.
2. Edit consensus_corces.R file (line 254) and provide the path for your summit files which you got from MACS2 peak calling.
3. Run this then with the following code:

module load R/4.1.0  ### On valinor
module load R/4.2.1  ### On h4h
sbatch -p himem -J consensus --export=ALL -N 1 -n 1 --mem 30G -t 2-0 --wrap "Rscript consensus_corces.R"

#This code uses annotation of genome (Promoter, Distal, Exon, Intron) from this file:

https://ftp.ebi.ac.uk/pub/databases/gencode/Gencode_human/release_39/gencode.v39.annotation.gtf.gz

#You can visualize the file on IGV by unziping it and sort:

igvtools sort gencode.v39.annotation.gtf