add plotEnrichment
from BAMs by coordinate index run https://deeptools.readthedocs.io/en/develop/content/tools/plotEnrichment.html
plotEnrichment --bamfiles results/merged/NeuNposLmx1bneg_sort.bam \
results/merged/NeuNposLmx1bpos_sort.bam \
results/merged/NeuNpos_sort.bam \
--BED results/genrich/NeuNposLmx1bneg.narrowPeak \
results/genrich/NeuNposLmx1bpos.narrowPeak \
results/genrich/NeuNpos.narrowPeak -p 4 \
-- samFlagInclude 64 \
--minMappingQuality 30 --centerReads -o enrichment_center_q30.png