# ABE_count_FIREs.sh was run first to count FIRE elements overlapping with peaks for each haplotype.

# ABE_peak_stats.py was run next to calculate raw and adjusted p values and to generate volcano plot.