cd /BRC/yan/heart/adt_heart_narrowPeak/peak_union

FEMALE=$(wc -l < female_peak_union.bed)
MALE=$(wc -l < male_peak_union.bed)

OVERLAP=$(bedtools intersect -a female_peak_union.bed -b male_peak_union.bed -u | wc -l)

bedtools intersect -a female_peak_union.bed -b male_peak_union.bed -v > female_only.bed
bedtools intersect -a male_peak_union.bed -b female_peak_union.bed -v > male_only.bed

FEMALE_ONLY=$(wc -l < female_only.bed)
MALE_ONLY=$(wc -l < male_only.bed)

echo "=============================="
echo "Female total : ${FEMALE}"
echo "Male total   : ${MALE}"
echo "------------------------------"
echo "Overlap      : ${OVERLAP}"
printf "  %% of Female : %.1f%%\n" $(echo "${OVERLAP} ${FEMALE}" | awk '{print $1/$2*100}')
printf "  %% of Male   : %.1f%%\n" $(echo "${OVERLAP} ${MALE}"   | awk '{print $1/$2*100}')
echo "------------------------------"
echo "Female-only  : ${FEMALE_ONLY}"
printf "  %% of Female : %.1f%%\n" $(echo "${FEMALE_ONLY} ${FEMALE}" | awk '{print $1/$2*100}')
echo "Male-only    : ${MALE_ONLY}"
printf "  %% of Male   : %.1f%%\n" $(echo "${MALE_ONLY} ${MALE}"     | awk '{print $1/$2*100}')
echo "=============================="