#!/bin/bash
set -euo pipefail

INDIR="/BRC/yan/heart/adt_heart_narrowPeak"
OUTDIR="/BRC/yan/heart/adt_heart_narrowPeak/peak_union"
mkdir -p ${OUTDIR}

OUT=${OUTDIR}/all_peak_union.bed

echo "[$(date)] Processing all peaks..."

cat ${INDIR}/*_Adule_*_peaks.narrowPeak \
    | awk 'BEGIN{OFS="\t"}{print $1,$2,$3}' \
    | grep -v "^chrM" \
    | sort -k1,1 -k2,2n \
    | bedtools merge -i stdin \
    > ${OUT}

echo "  → union: $(wc -l < ${OUT}) peaks  →  ${OUT}"
echo "[$(date)] Done."