intervene_figure.sh

# Make overlap of enhancer data


cp ../H3K4me1_filtered_peaks.bed H3K4me1_enhancers.bed
cp ../H3K27ac_filtered_peaks.bed H3K27ac_enhancers.bed
cp  ../SUNP_filtered_peaks.bed .
cp  ../SUNP_filtered_peaks.bed .
cut -f1,2,3  SUNP_filtered_peaks.bed > t
cut -f1,2,3 SSP_filtered_peaks_merged.bed >> t
sort-bed t > GRO-seq_filtered_peaks.bed

intervene upset -i *enhancers.bed --filename --order degree  --sxlabel "No. of Unique Enhancers" --dpi 1200  --figsize 10 8 --mbcolor "#108DCC"  --sbcolor "#E62128"


# All called example:
chr10:51,534,904-51,592,022
chr12:131,332,043-131,362,365

Maybe:
chr13:28,405,807-28,514,065