Last updated: 2022-03-21

Checks: 7 0

Knit directory: chipseq-cross-species/

This reproducible R Markdown analysis was created with workflowr (version 1.7.0). The Checks tab describes the reproducibility checks that were applied when the results were created. The Past versions tab lists the development history.


Great! Since the R Markdown file has been committed to the Git repository, you know the exact version of the code that produced these results.

Great job! The global environment was empty. Objects defined in the global environment can affect the analysis in your R Markdown file in unknown ways. For reproduciblity it’s best to always run the code in an empty environment.

The command set.seed(20220209) was run prior to running the code in the R Markdown file. Setting a seed ensures that any results that rely on randomness, e.g. subsampling or permutations, are reproducible.

Great job! Recording the operating system, R version, and package versions is critical for reproducibility.

Nice! There were no cached chunks for this analysis, so you can be confident that you successfully produced the results during this run.

Great job! Using relative paths to the files within your workflowr project makes it easier to run your code on other machines.

Great! You are using Git for version control. Tracking code development and connecting the code version to the results is critical for reproducibility.

The results in this page were generated with repository version 01af838. See the Past versions tab to see a history of the changes made to the R Markdown and HTML files.

Note that you need to be careful to ensure that all relevant files for the analysis have been committed to Git prior to generating the results (you can use wflow_publish or wflow_git_commit). workflowr only checks the R Markdown file, but you know if there are other scripts or data files that it depends on. Below is the status of the Git repository when the results were generated:


Ignored files:
    Ignored:    .RData
    Ignored:    .Rhistory
    Ignored:    analysis/.RData
    Ignored:    analysis/.Rhistory
    Ignored:    data/genomic_data/
    Ignored:    data/geo_submission/
    Ignored:    output/annotations/
    Ignored:    output/bam_files/
    Ignored:    output/filtered_peaks/
    Ignored:    output/logs/
    Ignored:    output/peaks/
    Ignored:    output/plots/
    Ignored:    output/qc/A-1_input.SeqDepthNorm.bw
    Ignored:    output/qc/A-1_input.SeqDepthNorm_dunnart_downSampled.bw
    Ignored:    output/qc/A-1_input.ccurve.txt
    Ignored:    output/qc/A-1_input.extrap.txt
    Ignored:    output/qc/A-1_input_R1_trimmed.fastq.gz
    Ignored:    output/qc/A-1_input_est_lib_complex_metrics.txt
    Ignored:    output/qc/A-2_H3K4me3.SeqDepthNorm.bw
    Ignored:    output/qc/A-2_H3K4me3.SeqDepthNorm_dunnart_downSampled.bw
    Ignored:    output/qc/A-2_H3K4me3.ccurve.txt
    Ignored:    output/qc/A-2_H3K4me3.extrap.txt
    Ignored:    output/qc/A-2_H3K4me3_R1_trimmed.fastq.gz
    Ignored:    output/qc/A-2_H3K4me3_est_lib_complex_metrics.txt
    Ignored:    output/qc/A-2_H3K4me3_vs_A-1_input.frip_default.txt
    Ignored:    output/qc/A-2_H3K4me3_vs_A-1_input.frip_default_dunnart_downSampled.txt
    Ignored:    output/qc/A-2_H3K4me3_vs_A-1_input.frip_p0.01_dunnart_downSampled.txt
    Ignored:    output/qc/A-2_H3K4me3_vs_B-1_input.frip_default.txt
    Ignored:    output/qc/A-2_H3K4me3_vs_B-1_input.frip_default_dunnart_downSampled.txt
    Ignored:    output/qc/A-2_H3K4me3_vs_B-1_input.frip_p0.01_dunnart_downSampled.txt
    Ignored:    output/qc/A-3_H3K27ac.SeqDepthNorm.bw
    Ignored:    output/qc/A-3_H3K27ac.SeqDepthNorm_dunnart_downSampled.bw
    Ignored:    output/qc/A-3_H3K27ac.ccurve.txt
    Ignored:    output/qc/A-3_H3K27ac.extrap.txt
    Ignored:    output/qc/A-3_H3K27ac_R1_trimmed.fastq.gz
    Ignored:    output/qc/A-3_H3K27ac_est_lib_complex_metrics.txt
    Ignored:    output/qc/A-3_H3K27ac_vs_A-1_input.frip_default.txt
    Ignored:    output/qc/A-3_H3K27ac_vs_A-1_input.frip_default_dunnart_downSampled.txt
    Ignored:    output/qc/A-3_H3K27ac_vs_A-1_input.frip_p0.01_dunnart_downSampled.txt
    Ignored:    output/qc/A-3_H3K27ac_vs_B-1_input.frip_default.txt
    Ignored:    output/qc/A-3_H3K27ac_vs_B-1_input.frip_default_dunnart_downSampled.txt
    Ignored:    output/qc/A-3_H3K27ac_vs_B-1_input.frip_p0.01_dunnart_downSampled.txt
    Ignored:    output/qc/B-1_input.SeqDepthNorm.bw
    Ignored:    output/qc/B-1_input.SeqDepthNorm_dunnart_downSampled.bw
    Ignored:    output/qc/B-1_input.ccurve.txt
    Ignored:    output/qc/B-1_input.extrap.txt
    Ignored:    output/qc/B-1_input_R1_trimmed.fastq.gz
    Ignored:    output/qc/B-1_input_est_lib_complex_metrics.txt
    Ignored:    output/qc/B-2_H3K4me3.SeqDepthNorm.bw
    Ignored:    output/qc/B-2_H3K4me3.SeqDepthNorm_dunnart_downSampled.bw
    Ignored:    output/qc/B-2_H3K4me3.ccurve.txt
    Ignored:    output/qc/B-2_H3K4me3.extrap.txt
    Ignored:    output/qc/B-2_H3K4me3_R1_trimmed.fastq.gz
    Ignored:    output/qc/B-2_H3K4me3_est_lib_complex_metrics.txt
    Ignored:    output/qc/B-2_H3K4me3_vs_A-1_input.frip_default.txt
    Ignored:    output/qc/B-2_H3K4me3_vs_A-1_input.frip_default_dunnart_downSampled.txt
    Ignored:    output/qc/B-2_H3K4me3_vs_A-1_input.frip_p0.01_dunnart_downSampled.txt
    Ignored:    output/qc/B-2_H3K4me3_vs_B-1_input.frip_default.txt
    Ignored:    output/qc/B-2_H3K4me3_vs_B-1_input.frip_default_dunnart_downSampled.txt
    Ignored:    output/qc/B-2_H3K4me3_vs_B-1_input.frip_p0.01_dunnart_downSampled.txt
    Ignored:    output/qc/B-3_H3K27ac.SeqDepthNorm.bw
    Ignored:    output/qc/B-3_H3K27ac.SeqDepthNorm_dunnart_downSampled.bw
    Ignored:    output/qc/B-3_H3K27ac.ccurve.txt
    Ignored:    output/qc/B-3_H3K27ac.extrap.txt
    Ignored:    output/qc/B-3_H3K27ac_R1_trimmed.fastq.gz
    Ignored:    output/qc/B-3_H3K27ac_est_lib_complex_metrics.txt
    Ignored:    output/qc/B-3_H3K27ac_vs_A-1_input.frip_default.txt
    Ignored:    output/qc/B-3_H3K27ac_vs_A-1_input.frip_default_dunnart_downSampled.txt
    Ignored:    output/qc/B-3_H3K27ac_vs_A-1_input.frip_p0.01_dunnart_downSampled.txt
    Ignored:    output/qc/B-3_H3K27ac_vs_B-1_input.frip_default.txt
    Ignored:    output/qc/B-3_H3K27ac_vs_B-1_input.frip_default_dunnart_downSampled.txt
    Ignored:    output/qc/B-3_H3K27ac_vs_B-1_input.frip_p0.01_dunnart_downSampled.txt
    Ignored:    output/qc/ENCFF011NFM_E12.5_H3K27ac.SeqDepthNorm.bw
    Ignored:    output/qc/ENCFF011NFM_vs_ENCFF058AUT_E12.5_H3K27ac_downSampled.frip.txt
    Ignored:    output/qc/ENCFF045IPK_E10.5_H3K4me3.SeqDepthNorm_mouse.bw
    Ignored:    output/qc/ENCFF045IPK_vs_ENCFF825AVI_E10.5_H3K4me3_downSampled.frip.txt
    Ignored:    output/qc/ENCFF124TAB_E13.5_H3K4me3.SeqDepthNorm_mouse.bw
    Ignored:    output/qc/ENCFF124TAB_vs_ENCFF248PGK_E13.5_H3K4me3_downSampled.frip.txt
    Ignored:    output/qc/ENCFF124UYX_E10.5_H3K4me3.SeqDepthNorm_mouse.bw
    Ignored:    output/qc/ENCFF124UYX_vs_ENCFF157KEH_E10.5_H3K4me3_downSampled.frip.txt
    Ignored:    output/qc/ENCFF182ZPF_E12.5_H3K4me3.SeqDepthNorm_mouse.bw
    Ignored:    output/qc/ENCFF182ZPF_vs_ENCFF203JQV_E12.5_H3K4me3_downSampled.frip.txt
    Ignored:    output/qc/ENCFF194ORC_E13.5_H3K27ac.SeqDepthNorm.bw
    Ignored:    output/qc/ENCFF194ORC_vs_ENCFF117QRC_E13.5_H3K27ac_downSampled.frip.txt
    Ignored:    output/qc/ENCFF213EBC_E10.5_H3K27ac.SeqDepthNorm.bw
    Ignored:    output/qc/ENCFF213EBC_vs_ENCFF157KEH_E10.5_H3K27ac_downSampled.frip.txt
    Ignored:    output/qc/ENCFF258KCR_E15.5_H3K4me3.SeqDepthNorm_mouse.bw
    Ignored:    output/qc/ENCFF258KCR_vs_ENCFF727QTS_E15.5_H3K4me3_downSampled.frip.txt
    Ignored:    output/qc/ENCFF290ZNF_E13.5_H3K27ac.SeqDepthNorm.bw
    Ignored:    output/qc/ENCFF290ZNF_vs_ENCFF248PGK_E13.5_H3K27ac_downSampled.frip.txt
    Ignored:    output/qc/ENCFF327VAO_E14.5_H3K27ac.SeqDepthNorm.bw
    Ignored:    output/qc/ENCFF327VAO_vs_ENCFF784ORI_E14.5_H3K27ac_downSampled.frip.txt
    Ignored:    output/qc/ENCFF394TZN_E12.5_H3K27ac.SeqDepthNorm.bw
    Ignored:    output/qc/ENCFF394TZN_vs_ENCFF203JQV_E12.5_H3K27ac_downSampled.frip.txt
    Ignored:    output/qc/ENCFF401BKM_E15.5_H3K4me3.SeqDepthNorm_mouse.bw
    Ignored:    output/qc/ENCFF401BKM_vs_ENCFF182XFG_E15.5_H3K4me3_downSampled.frip.txt
    Ignored:    output/qc/ENCFF485UDC_E13.5_H3K4me3.SeqDepthNorm_mouse.bw
    Ignored:    output/qc/ENCFF485UDC_vs_ENCFF117QRC_E13.5_H3K4me3_downSampled.frip.txt
    Ignored:    output/qc/ENCFF512SFE_E11.5_H3K27ac.SeqDepthNorm.bw
    Ignored:    output/qc/ENCFF512SFE_vs_ENCFF184CUE_E11.5_H3K27ac_downSampled.frip.txt
    Ignored:    output/qc/ENCFF515PKL_E11.5_H3K27ac.SeqDepthNorm.bw
    Ignored:    output/qc/ENCFF515PKL_vs_ENCFF376FGM_E11.5_H3K27ac_downSampled.frip.txt
    Ignored:    output/qc/ENCFF548BRR_E10.5_H3K27ac.SeqDepthNorm.bw
    Ignored:    output/qc/ENCFF548BRR_vs_ENCFF825AVI_E10.5_H3K27ac_downSampled.frip.txt
    Ignored:    output/qc/ENCFF584JFB_E15.5_H3K27ac.SeqDepthNorm.bw
    Ignored:    output/qc/ENCFF584JFB_vs_ENCFF727QTS_E15.5_H3K27ac_downSampled.frip.txt
    Ignored:    output/qc/ENCFF665QBJ_E14.5_H3K4me3.SeqDepthNorm_mouse.bw
    Ignored:    output/qc/ENCFF665QBJ_vs_ENCFF002HZV_E14.5_H3K4me3_downSampled.frip.txt
    Ignored:    output/qc/ENCFF707WKL_E15.5_H3K27ac.SeqDepthNorm.bw
    Ignored:    output/qc/ENCFF707WKL_vs_ENCFF182XFG_E15.5_H3K27ac_downSampled.frip.txt
    Ignored:    output/qc/ENCFF717QDV_E11.5_H3K4me3.SeqDepthNorm_mouse.bw
    Ignored:    output/qc/ENCFF717QDV_vs_ENCFF376FGM_E11.5_H3K4me3_downSampled.frip.txt
    Ignored:    output/qc/ENCFF724DMU_E14.5_H3K4me3.SeqDepthNorm_mouse.bw
    Ignored:    output/qc/ENCFF724DMU_vs_ENCFF784ORI_E14.5_H3K4me3_downSampled.frip.txt
    Ignored:    output/qc/ENCFF760QYZ_E11.5_H3K4me3.SeqDepthNorm_mouse.bw
    Ignored:    output/qc/ENCFF760QYZ_vs_ENCFF184CUE_E11.5_H3K4me3_downSampled.frip.txt
    Ignored:    output/qc/ENCFF902HAR_E14.5_H3K27ac.SeqDepthNorm.bw
    Ignored:    output/qc/ENCFF902HAR_vs_ENCFF002HZV_E14.5_H3K27ac_downSampled.frip.txt
    Ignored:    output/qc/ENCFF941QJZ_E12.5_H3K4me3.SeqDepthNorm_mouse.bw
    Ignored:    output/qc/ENCFF941QJZ_vs_ENCFF058AUT_E12.5_H3K4me3_downSampled.frip.txt
    Ignored:    output/qc/H3K27ac_multiBAM_fingerprint_metrics.txt
    Ignored:    output/qc/H3K27ac_multiBAM_fingerprint_metrics_mouse.txt
    Ignored:    output/qc/H3K27ac_multiBAM_fingerprint_rawcounts.txt
    Ignored:    output/qc/H3K27ac_multiBAM_fingerprint_rawcounts_mouse.txt
    Ignored:    output/qc/H3K27ac_multibamsum.npz
    Ignored:    output/qc/H3K27ac_multibamsum.tab
    Ignored:    output/qc/H3K27ac_multibamsum_mouse.npz
    Ignored:    output/qc/H3K27ac_multibamsum_mouse.tab
    Ignored:    output/qc/H3K27ac_pearsoncor_multibamsum_matrix.txt
    Ignored:    output/qc/H3K27ac_pearsoncor_multibamsum_matrix_mouse.txt
    Ignored:    output/qc/H3K27ac_plot_coverage_rawcounts.tab
    Ignored:    output/qc/H3K27ac_plot_coverage_rawcounts_mouse.tab
    Ignored:    output/qc/H3K4me3_multiBAM_fingerprint_metrics_mouse.txt
    Ignored:    output/qc/H3K4me3_multiBAM_fingerprint_rawcounts_mouse.txt
    Ignored:    output/qc/H3K4me3_pearsoncor_multibamsum_matrix_mouse.txt
    Ignored:    output/qc/bamPEFragmentSize_rawcounts.tab
    Ignored:    output/qc/bamPEFragmentSize_rawcounts_dunnart_downSampled.tab
    Ignored:    output/qc/multiBAM_fingerprint_metrics.txt
    Ignored:    output/qc/multiBAM_fingerprint_metrics_dunnart_downSampled.txt
    Ignored:    output/qc/multiBAM_fingerprint_rawcounts.txt
    Ignored:    output/qc/multiBAM_fingerprint_rawcounts_dunnart_downSampled.txt
    Ignored:    output/qc/multibamsum.npz
    Ignored:    output/qc/multibamsum.tab
    Ignored:    output/qc/multibamsum_dunnart_downSampled.npz
    Ignored:    output/qc/multibamsum_dunnart_downSampled.tab
    Ignored:    output/qc/multibamsum_mouse.npz
    Ignored:    output/qc/multibamsum_mouse.tab
    Ignored:    output/qc/pearsoncor_multibamsum_matrix.txt
    Ignored:    output/qc/pearsoncor_multibamsum_matrix_dunnart_downSampled.txt
    Ignored:    output/qc/pearsoncor_multibamsum_matrix_mouse.txt
    Ignored:    output/wga/

Untracked files:
    Untracked:  .snakemake/
    Untracked:  Rplots.pdf
    Untracked:  code/.swp
    Untracked:  data/raw_reads/
    Untracked:  output/qc/A-1_input.downSampled.flagstat.qc
    Untracked:  output/qc/A-2_H3K4me3.downSampled.flagstat.qc
    Untracked:  output/qc/A-3_H3K27ac.downSampled.flagstat.qc
    Untracked:  output/qc/B-1_input.downSampled.flagstat.qc
    Untracked:  output/qc/B-2_H3K4me3.downSampled.flagstat.qc
    Untracked:  output/qc/B-3_H3K27ac.downSampled.flagstat.qc
    Untracked:  output/qc/E10.5_H3K27ac.pooled.macs2_downSampled
    Untracked:  output/qc/E10.5_H3K4me3.pooled.macs2_downSampled
    Untracked:  output/qc/E10.5_H3K4me3.pooled.macs2_downSampled.log
    Untracked:  output/qc/E11.5_H3K27ac.pooled.macs2_downSampled
    Untracked:  output/qc/E11.5_H3K4me3.pooled.macs2_downSampled
    Untracked:  output/qc/E11.5_H3K4me3.pooled.macs2_downSampled.log
    Untracked:  output/qc/E12.5_H3K27ac.pooled.macs2_downSampled
    Untracked:  output/qc/E12.5_H3K4me3.pooled.macs2_downSampled
    Untracked:  output/qc/E12.5_H3K4me3.pooled.macs2_downSampled.log
    Untracked:  output/qc/E13.5_H3K27ac.pooled.macs2_downSampled
    Untracked:  output/qc/E13.5_H3K4me3.pooled.macs2_downSampled
    Untracked:  output/qc/E13.5_H3K4me3.pooled.macs2_downSampled.log
    Untracked:  output/qc/E14.5_H3K27ac.pooled.macs2_downSampled
    Untracked:  output/qc/E14.5_H3K4me3.pooled.macs2_downSampled
    Untracked:  output/qc/E14.5_H3K4me3.pooled.macs2_downSampled.log
    Untracked:  output/qc/E15.5_H3K4me3.pooled.macs2_downSampled.log
    Untracked:  output/qc/ENCFF002HZV_E14.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF002HZV_E14.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF011NFM_E12.5_H3K27ac.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF011NFM_E12.5_H3K27ac.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF011NFM_E12.5_H3K27ac.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF011NFM_E12.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF011NFM_E12.5_H3K27ac.pbc.qc
    Untracked:  output/qc/ENCFF011NFM_E12.5_H3K27ac.q30.flagstat.qc
    Untracked:  output/qc/ENCFF011NFM_E12.5_H3K27ac.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF011NFM_E12.5_H3K27ac_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF011NFM_E12.5_H3K27ac_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF045IPK_E10.5_H3K4me3.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF045IPK_E10.5_H3K4me3.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF045IPK_E10.5_H3K4me3.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF045IPK_E10.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF045IPK_E10.5_H3K4me3.pbc.qc
    Untracked:  output/qc/ENCFF045IPK_E10.5_H3K4me3.q30.flagstat.qc
    Untracked:  output/qc/ENCFF045IPK_E10.5_H3K4me3.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF045IPK_E10.5_H3K4me3_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF045IPK_E10.5_H3K4me3_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF058AUT_E12.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF058AUT_E12.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF117QRC_E13.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF117QRC_E13.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF124TAB_E13.5_H3K4me3.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF124TAB_E13.5_H3K4me3.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF124TAB_E13.5_H3K4me3.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF124TAB_E13.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF124TAB_E13.5_H3K4me3.pbc.qc
    Untracked:  output/qc/ENCFF124TAB_E13.5_H3K4me3.q30.flagstat.qc
    Untracked:  output/qc/ENCFF124TAB_E13.5_H3K4me3.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF124TAB_E13.5_H3K4me3_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF124TAB_E13.5_H3K4me3_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF124UYX_E10.5_H3K4me3.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF124UYX_E10.5_H3K4me3.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF124UYX_E10.5_H3K4me3.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF124UYX_E10.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF124UYX_E10.5_H3K4me3.pbc.qc
    Untracked:  output/qc/ENCFF124UYX_E10.5_H3K4me3.q30.flagstat.qc
    Untracked:  output/qc/ENCFF124UYX_E10.5_H3K4me3.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF124UYX_E10.5_H3K4me3_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF124UYX_E10.5_H3K4me3_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF157KEH_E10.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF157KEH_E10.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF182XFG_E15.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF182XFG_E15.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF182ZPF_E12.5_H3K4me3.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF182ZPF_E12.5_H3K4me3.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF182ZPF_E12.5_H3K4me3.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF182ZPF_E12.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF182ZPF_E12.5_H3K4me3.pbc.qc
    Untracked:  output/qc/ENCFF182ZPF_E12.5_H3K4me3.q30.flagstat.qc
    Untracked:  output/qc/ENCFF182ZPF_E12.5_H3K4me3.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF182ZPF_E12.5_H3K4me3_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF182ZPF_E12.5_H3K4me3_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF184CUE_E11.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF184CUE_E11.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF194ORC_E13.5_H3K27ac.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF194ORC_E13.5_H3K27ac.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF194ORC_E13.5_H3K27ac.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF194ORC_E13.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF194ORC_E13.5_H3K27ac.pbc.qc
    Untracked:  output/qc/ENCFF194ORC_E13.5_H3K27ac.q30.flagstat.qc
    Untracked:  output/qc/ENCFF194ORC_E13.5_H3K27ac.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF194ORC_E13.5_H3K27ac_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF194ORC_E13.5_H3K27ac_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF203JQV_E12.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF203JQV_E12.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF213EBC_E10.5_H3K27ac.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF213EBC_E10.5_H3K27ac.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF213EBC_E10.5_H3K27ac.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF213EBC_E10.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF213EBC_E10.5_H3K27ac.pbc.qc
    Untracked:  output/qc/ENCFF213EBC_E10.5_H3K27ac.q30.flagstat.qc
    Untracked:  output/qc/ENCFF213EBC_E10.5_H3K27ac.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF213EBC_E10.5_H3K27ac_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF213EBC_E10.5_H3K27ac_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF248PGK_E13.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF248PGK_E13.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF258KCR_E15.5_H3K4me3.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF258KCR_E15.5_H3K4me3.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF258KCR_E15.5_H3K4me3.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF258KCR_E15.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF258KCR_E15.5_H3K4me3.pbc.qc
    Untracked:  output/qc/ENCFF258KCR_E15.5_H3K4me3.q30.flagstat.qc
    Untracked:  output/qc/ENCFF258KCR_E15.5_H3K4me3.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF258KCR_E15.5_H3K4me3_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF258KCR_E15.5_H3K4me3_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF290ZNF_E13.5_H3K27ac.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF290ZNF_E13.5_H3K27ac.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF290ZNF_E13.5_H3K27ac.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF290ZNF_E13.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF290ZNF_E13.5_H3K27ac.pbc.qc
    Untracked:  output/qc/ENCFF290ZNF_E13.5_H3K27ac.q30.flagstat.qc
    Untracked:  output/qc/ENCFF290ZNF_E13.5_H3K27ac.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF290ZNF_E13.5_H3K27ac_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF290ZNF_E13.5_H3K27ac_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF327VAO_E14.5_H3K27ac.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF327VAO_E14.5_H3K27ac.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF327VAO_E14.5_H3K27ac.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF327VAO_E14.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF327VAO_E14.5_H3K27ac.pbc.qc
    Untracked:  output/qc/ENCFF327VAO_E14.5_H3K27ac.q30.flagstat.qc
    Untracked:  output/qc/ENCFF327VAO_E14.5_H3K27ac.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF327VAO_E14.5_H3K27ac_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF327VAO_E14.5_H3K27ac_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF376FGM_E11.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF376FGM_E11.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF394TZN_E12.5_H3K27ac.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF394TZN_E12.5_H3K27ac.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF394TZN_E12.5_H3K27ac.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF394TZN_E12.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF394TZN_E12.5_H3K27ac.pbc.qc
    Untracked:  output/qc/ENCFF394TZN_E12.5_H3K27ac.q30.flagstat.qc
    Untracked:  output/qc/ENCFF394TZN_E12.5_H3K27ac.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF394TZN_E12.5_H3K27ac_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF394TZN_E12.5_H3K27ac_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF401BKM_E15.5_H3K4me3.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF401BKM_E15.5_H3K4me3.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF401BKM_E15.5_H3K4me3.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF401BKM_E15.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF401BKM_E15.5_H3K4me3.pbc.qc
    Untracked:  output/qc/ENCFF401BKM_E15.5_H3K4me3.q30.flagstat.qc
    Untracked:  output/qc/ENCFF401BKM_E15.5_H3K4me3.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF401BKM_E15.5_H3K4me3_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF401BKM_E15.5_H3K4me3_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF485UDC_E13.5_H3K4me3.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF485UDC_E13.5_H3K4me3.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF485UDC_E13.5_H3K4me3.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF485UDC_E13.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF485UDC_E13.5_H3K4me3.pbc.qc
    Untracked:  output/qc/ENCFF485UDC_E13.5_H3K4me3.q30.flagstat.qc
    Untracked:  output/qc/ENCFF485UDC_E13.5_H3K4me3.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF485UDC_E13.5_H3K4me3_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF485UDC_E13.5_H3K4me3_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF512SFE_E11.5_H3K27ac.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF512SFE_E11.5_H3K27ac.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF512SFE_E11.5_H3K27ac.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF512SFE_E11.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF512SFE_E11.5_H3K27ac.pbc.qc
    Untracked:  output/qc/ENCFF512SFE_E11.5_H3K27ac.q30.flagstat.qc
    Untracked:  output/qc/ENCFF512SFE_E11.5_H3K27ac.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF512SFE_E11.5_H3K27ac_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF512SFE_E11.5_H3K27ac_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF515PKL_E11.5_H3K27ac.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF515PKL_E11.5_H3K27ac.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF515PKL_E11.5_H3K27ac.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF515PKL_E11.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF515PKL_E11.5_H3K27ac.pbc.qc
    Untracked:  output/qc/ENCFF515PKL_E11.5_H3K27ac.q30.flagstat.qc
    Untracked:  output/qc/ENCFF515PKL_E11.5_H3K27ac.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF515PKL_E11.5_H3K27ac_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF515PKL_E11.5_H3K27ac_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF548BRR_E10.5_H3K27ac.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF548BRR_E10.5_H3K27ac.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF548BRR_E10.5_H3K27ac.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF548BRR_E10.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF548BRR_E10.5_H3K27ac.pbc.qc
    Untracked:  output/qc/ENCFF548BRR_E10.5_H3K27ac.q30.flagstat.qc
    Untracked:  output/qc/ENCFF548BRR_E10.5_H3K27ac.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF548BRR_E10.5_H3K27ac_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF548BRR_E10.5_H3K27ac_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF584JFB_E15.5_H3K27ac.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF584JFB_E15.5_H3K27ac.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF584JFB_E15.5_H3K27ac.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF584JFB_E15.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF584JFB_E15.5_H3K27ac.pbc.qc
    Untracked:  output/qc/ENCFF584JFB_E15.5_H3K27ac.q30.flagstat.qc
    Untracked:  output/qc/ENCFF584JFB_E15.5_H3K27ac.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF584JFB_E15.5_H3K27ac_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF584JFB_E15.5_H3K27ac_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF665QBJ_E14.5_H3K4me3.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF665QBJ_E14.5_H3K4me3.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF665QBJ_E14.5_H3K4me3.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF665QBJ_E14.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF665QBJ_E14.5_H3K4me3.pbc.qc
    Untracked:  output/qc/ENCFF665QBJ_E14.5_H3K4me3.q30.flagstat.qc
    Untracked:  output/qc/ENCFF665QBJ_E14.5_H3K4me3.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF665QBJ_E14.5_H3K4me3_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF665QBJ_E14.5_H3K4me3_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF707WKL_E15.5_H3K27ac.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF707WKL_E15.5_H3K27ac.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF707WKL_E15.5_H3K27ac.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF707WKL_E15.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF707WKL_E15.5_H3K27ac.pbc.qc
    Untracked:  output/qc/ENCFF707WKL_E15.5_H3K27ac.q30.flagstat.qc
    Untracked:  output/qc/ENCFF707WKL_E15.5_H3K27ac.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF707WKL_E15.5_H3K27ac_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF707WKL_E15.5_H3K27ac_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF717QDV_E11.5_H3K4me3.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF717QDV_E11.5_H3K4me3.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF717QDV_E11.5_H3K4me3.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF717QDV_E11.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF717QDV_E11.5_H3K4me3.pbc.qc
    Untracked:  output/qc/ENCFF717QDV_E11.5_H3K4me3.q30.flagstat.qc
    Untracked:  output/qc/ENCFF717QDV_E11.5_H3K4me3.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF717QDV_E11.5_H3K4me3_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF717QDV_E11.5_H3K4me3_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF724DMU_E14.5_H3K4me3.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF724DMU_E14.5_H3K4me3.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF724DMU_E14.5_H3K4me3.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF724DMU_E14.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF724DMU_E14.5_H3K4me3.pbc.qc
    Untracked:  output/qc/ENCFF724DMU_E14.5_H3K4me3.q30.flagstat.qc
    Untracked:  output/qc/ENCFF724DMU_E14.5_H3K4me3.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF724DMU_E14.5_H3K4me3_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF724DMU_E14.5_H3K4me3_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF727QTS_E15.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF727QTS_E15.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF760QYZ_E11.5_H3K4me3.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF760QYZ_E11.5_H3K4me3.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF760QYZ_E11.5_H3K4me3.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF760QYZ_E11.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF760QYZ_E11.5_H3K4me3.pbc.qc
    Untracked:  output/qc/ENCFF760QYZ_E11.5_H3K4me3.q30.flagstat.qc
    Untracked:  output/qc/ENCFF760QYZ_E11.5_H3K4me3.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF760QYZ_E11.5_H3K4me3_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF760QYZ_E11.5_H3K4me3_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF784ORI_E14.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF784ORI_E14.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF825AVI_E10.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF825AVI_E10.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF902HAR_E14.5_H3K27ac.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF902HAR_E14.5_H3K27ac.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF902HAR_E14.5_H3K27ac.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF902HAR_E14.5_H3K27ac.dupmark.qc
    Untracked:  output/qc/ENCFF902HAR_E14.5_H3K27ac.pbc.qc
    Untracked:  output/qc/ENCFF902HAR_E14.5_H3K27ac.q30.flagstat.qc
    Untracked:  output/qc/ENCFF902HAR_E14.5_H3K27ac.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF902HAR_E14.5_H3K27ac_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF902HAR_E14.5_H3K27ac_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/ENCFF941QJZ_E12.5_H3K4me3.dedup.flagstat.qc
    Untracked:  output/qc/ENCFF941QJZ_E12.5_H3K4me3.downSampled.flagstat.qc
    Untracked:  output/qc/ENCFF941QJZ_E12.5_H3K4me3.dupmark.flagstat.qc
    Untracked:  output/qc/ENCFF941QJZ_E12.5_H3K4me3.dupmark.qc
    Untracked:  output/qc/ENCFF941QJZ_E12.5_H3K4me3.pbc.qc
    Untracked:  output/qc/ENCFF941QJZ_E12.5_H3K4me3.q30.flagstat.qc
    Untracked:  output/qc/ENCFF941QJZ_E12.5_H3K4me3.unfiltered.flagstat.qc
    Untracked:  output/qc/ENCFF941QJZ_E12.5_H3K4me3_filt_15Mreads.SE.cc.plot.pdf
    Untracked:  output/qc/ENCFF941QJZ_E12.5_H3K4me3_filt_15Mreads.SE.cc.qc
    Untracked:  output/qc/H3K27ac_E10.5_overlap.frip
    Untracked:  output/qc/H3K27ac_E11.5_overlap.frip
    Untracked:  output/qc/H3K27ac_E12.5_overlap.frip
    Untracked:  output/qc/H3K27ac_E13.5_overlap.frip
    Untracked:  output/qc/H3K27ac_E14.5_overlap.frip
    Untracked:  output/qc/H3K27ac_E15.5_overlap.frip
    Untracked:  output/qc/H3K27ac_multiBAM_fingerprint.pdf
    Untracked:  output/qc/H3K27ac_pearsoncor_multibamsum.pdf
    Untracked:  output/qc/H3K27ac_plot_coverage.pdf
    Untracked:  output/qc/H3K4me3_overlap_default_dunnart_downSampled.frip
    Untracked:  output/qc/H3K4me3_overlap_p0.01_dunnart_downSampled.frip
    Untracked:  output/qc/ucsc_alignment/
    Untracked:  output/rnaseq/
    Untracked:  slurm-34229611.out
    Untracked:  slurm-34229694.out
    Untracked:  slurm-34229768.out
    Untracked:  slurm-34230768.out
    Untracked:  slurm-34230769.out
    Untracked:  slurm-34230770.out
    Untracked:  slurm-34230771.out
    Untracked:  slurm-34230772.out
    Untracked:  slurm-34230852.out
    Untracked:  slurm-34230853.out
    Untracked:  slurm-34230854.out
    Untracked:  slurm-34230855.out
    Untracked:  slurm-34231065.out
    Untracked:  slurm-34231066.out
    Untracked:  slurm-34231071.out
    Untracked:  slurm-34231170.out
    Untracked:  slurm-34231171.out
    Untracked:  slurm-34231172.out
    Untracked:  slurm-34231173.out
    Untracked:  slurm-34231174.out
    Untracked:  slurm-34231175.out
    Untracked:  slurm-34231176.out
    Untracked:  slurm-34231177.out
    Untracked:  slurm-34231178.out
    Untracked:  slurm-34231244.out
    Untracked:  slurm-34231245.out
    Untracked:  slurm-34231246.out
    Untracked:  slurm-34231247.out
    Untracked:  slurm-34231248.out
    Untracked:  slurm-34231249.out
    Untracked:  slurm-34231250.out
    Untracked:  slurm-34231253.out
    Untracked:  slurm-34231355.out
    Untracked:  slurm-34231356.out
    Untracked:  slurm-34231362.out
    Untracked:  slurm-34233668.out
    Untracked:  slurm-34233684.out
    Untracked:  slurm-34233685.out
    Untracked:  slurm-34233686.out
    Untracked:  slurm-34233687.out
    Untracked:  slurm-34233688.out
    Untracked:  slurm-34233778.out
    Untracked:  slurm-34233811.out
    Untracked:  slurm-34233812.out
    Untracked:  slurm-34233813.out
    Untracked:  slurm-34233910.out
    Untracked:  slurm-34233911.out
    Untracked:  slurm-34233921.out
    Untracked:  slurm-34235405.out
    Untracked:  slurm-34235417.out
    Untracked:  slurm-34235418.out
    Untracked:  slurm-34235419.out
    Untracked:  slurm-34235420.out
    Untracked:  slurm-34235421.out
    Untracked:  slurm-34235431.out
    Untracked:  slurm-34235432.out
    Untracked:  slurm-34235433.out
    Untracked:  slurm-34235434.out
    Untracked:  slurm-34235435.out
    Untracked:  slurm-34235469.out
    Untracked:  slurm-34235470.out
    Untracked:  slurm-34235471.out
    Untracked:  slurm-34235472.out
    Untracked:  slurm-34235473.out
    Untracked:  slurm-34235474.out
    Untracked:  slurm-34235475.out
    Untracked:  slurm-34235476.out
    Untracked:  slurm-34235477.out
    Untracked:  slurm-34235478.out
    Untracked:  slurm-34235509.out
    Untracked:  slurm-34235510.out
    Untracked:  slurm-34235511.out
    Untracked:  slurm-34235512.out
    Untracked:  slurm-34235513.out
    Untracked:  slurm-34235514.out
    Untracked:  slurm-34235515.out
    Untracked:  slurm-34235516.out
    Untracked:  slurm-34235544.out
    Untracked:  slurm-34235545.out
    Untracked:  slurm-34235546.out
    Untracked:  slurm-34235547.out
    Untracked:  slurm-34235559.out
    Untracked:  slurm-34235560.out
    Untracked:  slurm-34235561.out
    Untracked:  slurm-34235562.out
    Untracked:  slurm-34235563.out
    Untracked:  slurm-34235564.out
    Untracked:  slurm-34235565.out
    Untracked:  slurm-34235566.out
    Untracked:  slurm-34235567.out
    Untracked:  slurm-34235568.out
    Untracked:  slurm-34235569.out
    Untracked:  slurm-34235571.out
    Untracked:  slurm-34235592.out
    Untracked:  slurm-34235593.out
    Untracked:  slurm-34235594.out
    Untracked:  slurm-34235595.out
    Untracked:  slurm-34235598.out
    Untracked:  slurm-34235610.out
    Untracked:  slurm-34235611.out

Unstaged changes:
    Modified:   analysis/mouse_dunnart_peak_features.Rmd
    Modified:   code/configs/cluster.json
    Deleted:    output/qc/H3K4me3_overlap_default.frip

Note that any generated files, e.g. HTML, png, CSS, etc., are not included in this status report because it is ok for generated content to have uncommitted changes.


These are the previous versions of the repository in which changes were made to the R Markdown (analysis/mouse_dunnart_data_processing_for_comparison.Rmd) and HTML (docs/mouse_dunnart_data_processing_for_comparison.html) files. If you’ve configured a remote Git repository (see ?wflow_git_remote), click on the hyperlinks in the table below to view the files as they were in that past version.

File Version Author Date Message
Rmd 01af838 lecook 2022-03-21 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html 8d8ca2b lecook 2022-03-21 Build site.
Rmd 8bd045b lecook 2022-03-21 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html cf8c198 lecook 2022-03-21 Build site.
Rmd 46f85ad lecook 2022-03-21 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html 512e82e lecook 2022-03-21 Build site.
html dae2727 lecook 2022-03-21 Build site.
Rmd fc98159 lecook 2022-03-21 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html 11f4efd lecook 2022-03-21 Build site.
html b7d3bab lecook 2022-03-21 Build site.
html 4a02b9b lecook 2022-03-21 Build site.
Rmd 82786d9 lecook 2022-03-21 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html 9ff818e lecook 2022-03-21 Build site.
Rmd 5eada99 lecook 2022-03-21 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html 31faf07 lecook 2022-03-18 Build site.
html 5ec5742 lecook 2022-03-18 Build site.
html 6281cd3 lecook 2022-03-18 Build site.
html 5762794 lecook 2022-03-18 Build site.
html 352deb8 lecook 2022-03-18 Build site.
Rmd 582184a lecook 2022-03-18 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html 6afb369 lecook 2022-03-18 Build site.
html fd09f75 lecook 2022-03-18 Build site.
html 0766c62 lecook 2022-03-18 Build site.
html b67e70b lecook 2022-03-17 Build site.
html eaf5612 lecook 2022-03-17 Build site.
Rmd af396d2 lecook 2022-03-17 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html a04c949 lecook 2022-03-17 Build site.
Rmd 6062643 lecook 2022-03-17 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html 7c3a344 lecook 2022-03-17 Build site.
Rmd eca7e88 lecook 2022-03-17 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html 80957c9 lecook 2022-03-17 Build site.
Rmd 6fc8343 lecook 2022-03-17 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html 56ecaa0 lecook 2022-03-16 Build site.
Rmd 7f26be0 lecook 2022-03-16 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html 76c6eee lecook 2022-03-16 Build site.
Rmd 9c12dfe lecook 2022-03-16 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html 1456cf5 lecook 2022-03-16 Build site.
Rmd 9bf4173 lecook 2022-03-16 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html e662415 lecook 2022-03-16 Build site.
Rmd 2720905 lecook 2022-03-16 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html 7767732 lecook 2022-03-16 Build site.
Rmd c447893 lecook 2022-03-16 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html a90d5a0 lecook 2022-03-16 Build site.
Rmd 87c4146 lecook 2022-03-16 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html 93879c2 lecook 2022-03-16 Build site.
html a552f17 lecook 2022-03-16 Build site.
Rmd 5a17b44 lecook 2022-03-16 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html 11fbc81 lecook 2022-03-16 Build site.
Rmd dffebdd lecook 2022-03-16 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html 4b71068 lecook 2022-03-16 Build site.
Rmd a287d1b lecook 2022-03-16 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
Rmd eaeea84 lecook 2022-03-16 updated
Rmd 81d5677 lecook 2022-03-16 fix downsampled
html c5d4b95 lecook 2022-03-16 Build site.
html f7966c8 lecook 2022-03-16 Build site.
Rmd 040a6a7 lecook 2022-03-16 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html f593bed lecook 2022-03-16 Build site.
Rmd 039df2a lecook 2022-03-16 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html 9d125bb lecook 2022-03-16 Build site.
html 9bcce4a lecook 2022-03-16 Build site.
Rmd c0a1798 lecook 2022-03-16 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
html 841e0f4 lecook 2022-03-16 Build site.
Rmd 275510f lecook 2022-03-16 wflow_publish(“analysis/mouse_dunnart_data_processing_for_comparison.Rmd”)
Rmd faabd95 lecook 2022-03-03 updated
Rmd 2103ae7 lecook 2022-03-03 updates
Rmd 323d8b6 lecook 2022-03-03 copied scripts over
Rmd ae523af lecook 2022-03-02 fix code formatting for bam coverage rule
Rmd a97d3c5 lecook 2022-03-01 first commit
html 275d229 lecook 2022-03-01 Build site.
Rmd 784c4ec lecook 2022-03-01 update
Rmd bf8fac9 lecook 2022-03-01 updated
Rmd 5fb3661 lecook 2022-03-01 clean up
Rmd c2eeed7 lecook 2022-02-28 debugged

Pipeline for mouse and dunnart peak calling

Download mouse unfiltered alignments from ENCODE

See mouse_data_ENCODE for details on these samples.

xargs -L 1 curl -O -J -L < ENCODE_files.txt
xargs -L 1 curl -O -J -L < control_files.txt

Filter mouse alignments

Uses samtools and picards MarkDuplicates to remove low quality reads. This is based on and in accordance with the ENCODE Guidelines and pipeline (https://github.com/ENCODE-DCC/chip-seq-pipeline2)

See full snakemake script for details: code/mouse_peak_calling/mouse_H3K4me3 and code/mouse_peak_calling/mouse_H3K27ac

Table. Read alignment quality metrics before and after filtering steps
file processing num_reads
A-1_input duplicates removed 47645260
A-1_input subsampled 10 million reads 9994998
A-1_input remove orphan reads and MAPQ < 30 87165286
A-1_input mapped reads 130587134
A-2_H3K4me3 duplicates removed 35617708
A-2_H3K4me3 subsampled 10 million reads 9995100
A-2_H3K4me3 remove orphan reads and MAPQ < 30 62678360
A-2_H3K4me3 mapped reads 103324742
A-3_H3K27ac duplicates removed 45700530
A-3_H3K27ac subsampled 10 million reads 10000116
A-3_H3K27ac remove orphan reads and MAPQ < 30 86503078
A-3_H3K27ac mapped reads 131071676
B-1_input duplicates removed 44433864
B-1_input subsampled 10 million reads 9995472
B-1_input remove orphan reads and MAPQ < 30 72303316
B-1_input mapped reads 111574640
B-2_H3K4me3 duplicates removed 42357180
B-2_H3K4me3 subsampled 10 million reads 10003240
B-2_H3K4me3 remove orphan reads and MAPQ < 30 67048994
B-2_H3K4me3 mapped reads 114146802
B-3_H3K27ac duplicates removed 43939460
B-3_H3K27ac subsampled 10 million reads 10003238
B-3_H3K27ac remove orphan reads and MAPQ < 30 66615704
B-3_H3K27ac mapped reads 104714846
ENCFF011NFM_E12.5_H3K27ac duplicates removed 28442356
ENCFF011NFM_E12.5_H3K27ac subsampled 10 million reads 10000583
ENCFF011NFM_E12.5_H3K27ac remove orphan reads and MAPQ < 30 32754035
ENCFF011NFM_E12.5_H3K27ac mapped reads 42905369
ENCFF045IPK_E10.5_H3K4me3 duplicates removed 27327155
ENCFF045IPK_E10.5_H3K4me3 subsampled 10 million reads 10004541
ENCFF045IPK_E10.5_H3K4me3 remove orphan reads and MAPQ < 30 34283317
ENCFF045IPK_E10.5_H3K4me3 mapped reads 40151871
ENCFF124TAB_E13.5_H3K4me3 duplicates removed 33020017
ENCFF124TAB_E13.5_H3K4me3 subsampled 10 million reads 9998371
ENCFF124TAB_E13.5_H3K4me3 remove orphan reads and MAPQ < 30 38139564
ENCFF124TAB_E13.5_H3K4me3 mapped reads 50188754
ENCFF124UYX_E10.5_H3K4me3 duplicates removed 47429840
ENCFF124UYX_E10.5_H3K4me3 subsampled 10 million reads 10002064
ENCFF124UYX_E10.5_H3K4me3 remove orphan reads and MAPQ < 30 65480729
ENCFF124UYX_E10.5_H3K4me3 mapped reads 80931343
ENCFF182ZPF_E12.5_H3K4me3 duplicates removed 26489691
ENCFF182ZPF_E12.5_H3K4me3 subsampled 10 million reads 9996266
ENCFF182ZPF_E12.5_H3K4me3 remove orphan reads and MAPQ < 30 30140587
ENCFF182ZPF_E12.5_H3K4me3 mapped reads 38675670
ENCFF194ORC_E13.5_H3K27ac duplicates removed 23235076
ENCFF194ORC_E13.5_H3K27ac subsampled 10 million reads 9998367
ENCFF194ORC_E13.5_H3K27ac remove orphan reads and MAPQ < 30 39115139
ENCFF194ORC_E13.5_H3K27ac mapped reads 51326572
ENCFF213EBC_E10.5_H3K27ac duplicates removed 61061524
ENCFF213EBC_E10.5_H3K27ac subsampled 10 million reads 9999609
ENCFF213EBC_E10.5_H3K27ac remove orphan reads and MAPQ < 30 74296920
ENCFF213EBC_E10.5_H3K27ac mapped reads 88370370
ENCFF258KCR_E15.5_H3K4me3 duplicates removed 38278946
ENCFF258KCR_E15.5_H3K4me3 subsampled 10 million reads 10002892
ENCFF258KCR_E15.5_H3K4me3 remove orphan reads and MAPQ < 30 44138244
ENCFF258KCR_E15.5_H3K4me3 mapped reads 57759464
ENCFF290ZNF_E13.5_H3K27ac duplicates removed 26624266
ENCFF290ZNF_E13.5_H3K27ac subsampled 10 million reads 9997991
ENCFF290ZNF_E13.5_H3K27ac remove orphan reads and MAPQ < 30 37416425
ENCFF290ZNF_E13.5_H3K27ac mapped reads 50570025
ENCFF327VAO_E14.5_H3K27ac duplicates removed 19228162
ENCFF327VAO_E14.5_H3K27ac subsampled 10 million reads 10002150
ENCFF327VAO_E14.5_H3K27ac remove orphan reads and MAPQ < 30 30104464
ENCFF327VAO_E14.5_H3K27ac mapped reads 38601448
ENCFF394TZN_E12.5_H3K27ac duplicates removed 29223815
ENCFF394TZN_E12.5_H3K27ac subsampled 10 million reads 10001895
ENCFF394TZN_E12.5_H3K27ac remove orphan reads and MAPQ < 30 35669429
ENCFF394TZN_E12.5_H3K27ac mapped reads 47079762
ENCFF401BKM_E15.5_H3K4me3 duplicates removed 29687861
ENCFF401BKM_E15.5_H3K4me3 subsampled 10 million reads 10001562
ENCFF401BKM_E15.5_H3K4me3 remove orphan reads and MAPQ < 30 34005723
ENCFF401BKM_E15.5_H3K4me3 mapped reads 43924122
ENCFF485UDC_E13.5_H3K4me3 duplicates removed 33094082
ENCFF485UDC_E13.5_H3K4me3 subsampled 10 million reads 9998225
ENCFF485UDC_E13.5_H3K4me3 remove orphan reads and MAPQ < 30 38287309
ENCFF485UDC_E13.5_H3K4me3 mapped reads 50049213
ENCFF512SFE_E11.5_H3K27ac duplicates removed 11879316
ENCFF512SFE_E11.5_H3K27ac subsampled 10 million reads 9999949
ENCFF512SFE_E11.5_H3K27ac remove orphan reads and MAPQ < 30 12827968
ENCFF512SFE_E11.5_H3K27ac mapped reads 18655610
ENCFF515PKL_E11.5_H3K27ac duplicates removed 11961286
ENCFF515PKL_E11.5_H3K27ac subsampled 10 million reads 10000378
ENCFF515PKL_E11.5_H3K27ac remove orphan reads and MAPQ < 30 13753883
ENCFF515PKL_E11.5_H3K27ac mapped reads 20764363
ENCFF548BRR_E10.5_H3K27ac duplicates removed 38069875
ENCFF548BRR_E10.5_H3K27ac subsampled 10 million reads 9999586
ENCFF548BRR_E10.5_H3K27ac remove orphan reads and MAPQ < 30 43572898
ENCFF548BRR_E10.5_H3K27ac mapped reads 50325403
ENCFF584JFB_E15.5_H3K27ac duplicates removed 35551034
ENCFF584JFB_E15.5_H3K27ac subsampled 10 million reads 9997290
ENCFF584JFB_E15.5_H3K27ac remove orphan reads and MAPQ < 30 40125566
ENCFF584JFB_E15.5_H3K27ac mapped reads 54217191
ENCFF665QBJ_E14.5_H3K4me3 duplicates removed 19804325
ENCFF665QBJ_E14.5_H3K4me3 subsampled 10 million reads 9999703
ENCFF665QBJ_E14.5_H3K4me3 remove orphan reads and MAPQ < 30 22484958
ENCFF665QBJ_E14.5_H3K4me3 mapped reads 29707720
ENCFF707WKL_E15.5_H3K27ac duplicates removed 27522173
ENCFF707WKL_E15.5_H3K27ac subsampled 10 million reads 10000248
ENCFF707WKL_E15.5_H3K27ac remove orphan reads and MAPQ < 30 30334974
ENCFF707WKL_E15.5_H3K27ac mapped reads 40139064
ENCFF717QDV_E11.5_H3K4me3 duplicates removed 13504654
ENCFF717QDV_E11.5_H3K4me3 subsampled 10 million reads 9998305
ENCFF717QDV_E11.5_H3K4me3 remove orphan reads and MAPQ < 30 14136366
ENCFF717QDV_E11.5_H3K4me3 mapped reads 21310481
ENCFF724DMU_E14.5_H3K4me3 duplicates removed 17278836
ENCFF724DMU_E14.5_H3K4me3 subsampled 10 million reads 10000466
ENCFF724DMU_E14.5_H3K4me3 remove orphan reads and MAPQ < 30 18964886
ENCFF724DMU_E14.5_H3K4me3 mapped reads 26143918
ENCFF760QYZ_E11.5_H3K4me3 duplicates removed 16515484
ENCFF760QYZ_E11.5_H3K4me3 subsampled 10 million reads 9999147
ENCFF760QYZ_E11.5_H3K4me3 remove orphan reads and MAPQ < 30 17547752
ENCFF760QYZ_E11.5_H3K4me3 mapped reads 26619353
ENCFF902HAR_E14.5_H3K27ac duplicates removed 15203030
ENCFF902HAR_E14.5_H3K27ac subsampled 10 million reads 10000944
ENCFF902HAR_E14.5_H3K27ac remove orphan reads and MAPQ < 30 22156787
ENCFF902HAR_E14.5_H3K27ac mapped reads 28628975
ENCFF941QJZ_E12.5_H3K4me3 duplicates removed 26459778
ENCFF941QJZ_E12.5_H3K4me3 subsampled 10 million reads 9997132
ENCFF941QJZ_E12.5_H3K4me3 remove orphan reads and MAPQ < 30 29985012
ENCFF941QJZ_E12.5_H3K4me3 mapped reads 39588024

Library complexity

Library Complexity ChIP-seq Standards:

Table. Library complexity metrics
PBC1 PBC2 Bottlenecking level NRF Complexity Flag colors
< 0.5 < 1 Severe < 0.5 Concerning Orange
0.5 ≤ PBC1 < 0.8 1 ≤ PBC2 < 3 Moderate 0.5 ≤ NRF < 0.8 Acceptable Yellow
0.8 ≤ PBC1 < 0.9 3 ≤ PBC2 < 10 Mild 0.8 ≤ NRF < 0.9 Compliant None
≥ 0.9 ≥ 10 None > 0.9 Ideal None
file TotalReadPairs DistinctReadPairs OneReadPair TwoReadPairs NRF=Distinct/Total PBC1=OnePair/Distinct PBC2=OnePair/TwoPair
A-1_input 15253416 11800078 9139380 2050439 0.773602 0.774519 4.457280
A-2_H3K4me3 10971609 8664376 6841920 1438893 0.789709 0.789661 4.754989
A-3_H3K27ac 15138750 11532208 8783937 2093115 0.761768 0.761687 4.196586
B-1_input 12653570 10379893 8511240 1534850 0.820313 0.819974 5.545324
B-2_H3K4me3 11731352 9749128 8099076 1373138 0.831032 0.830749 5.898224
B-3_H3K27ac 11653321 9873338 8362931 1280857 0.847255 0.847022 6.529168
ENCFF011NFM_E12.5_H3K27ac 32749172 28539951 24838523 3292599 0.871471 0.870307 7.543744
ENCFF045IPK_E10.5_H3K4me3 34281533 27360053 22206158 3901376 0.798099 0.811627 5.691878
ENCFF124TAB_E13.5_H3K4me3 38137879 33884965 30508506 2743682 0.888486 0.900355 11.119549
ENCFF124UYX_E10.5_H3K4me3 65477344 47549452 35813376 7957026 0.726197 0.753182 4.500849
ENCFF182ZPF_E12.5_H3K4me3 30138204 26530683 23539279 2536145 0.880301 0.887247 9.281519
ENCFF194ORC_E13.5_H3K27ac 39109371 24688776 15163850 6093811 0.631275 0.614200 2.488402
ENCFF213EBC_E10.5_H3K27ac 74283125 61117294 50208502 9085217 0.822761 0.821511 5.526395
ENCFF258KCR_E15.5_H3K4me3 44135436 39270102 35533421 2950620 0.889764 0.904847 12.042696
ENCFF290ZNF_E13.5_H3K27ac 37412029 27631957 20118736 5726925 0.738585 0.728097 3.513008
ENCFF327VAO_E14.5_H3K27ac 30101456 19803731 12562311 4968458 0.657899 0.634341 2.528412
ENCFF394TZN_E12.5_H3K27ac 35664003 29335171 24004004 4520699 0.822543 0.818267 5.309799
ENCFF401BKM_E15.5_H3K4me3 34003038 30398350 27571247 2274666 0.893989 0.906998 12.121009
ENCFF485UDC_E13.5_H3K4me3 38285412 33973468 30522248 2825345 0.887374 0.898414 10.803016
ENCFF512SFE_E11.5_H3K27ac 12825706 11914376 11059801 803865 0.928945 0.928274 13.758282
ENCFF515PKL_E11.5_H3K27ac 13750546 12050668 10525979 1367478 0.876377 0.873477 7.697366
ENCFF548BRR_E10.5_H3K27ac 43570876 38093302 33332547 4153213 0.874284 0.875024 8.025725
ENCFF584JFB_E15.5_H3K27ac 40113267 35979441 32286974 3361753 0.896946 0.897373 9.604208
ENCFF665QBJ_E14.5_H3K4me3 22484295 20039171 18001895 1705759 0.891252 0.898335 10.553598
ENCFF707WKL_E15.5_H3K27ac 30326867 27756562 25422521 2162883 0.915247 0.915910 11.753997
ENCFF717QDV_E11.5_H3K4me3 14135672 13565058 13040722 488272 0.959633 0.961347 26.707905
ENCFF724DMU_E14.5_H3K4me3 18964071 17429599 16101305 1158889 0.919085 0.923791 13.893742
ENCFF760QYZ_E11.5_H3K4me3 17546903 16596513 15740663 779061 0.945837 0.948432 20.204660
ENCFF902HAR_E14.5_H3K27ac 22153805 15213820 10058815 3762364 0.686736 0.661163 2.673536
ENCFF941QJZ_E12.5_H3K4me3 29982687 26500387 23589168 2501027 0.883856 0.890144 9.431793

Normalise dunnart and mouse BAM files to 10M reads

This is so we can compare to mouse peaks.

Based on https://davemcg.github.io/post/easy-bam-downsampling/ script for doing this.

Run script in the directory with the files you want to subsample.

cd output/bam_files/
bash subsample.sh

Read QC with deepTools

Plot correlation between BAM files

plotCorrelation computes the overall similarity between two or more files based on read coverage (or other scores) within genomic regions.

This helps to determine whether the different sample types can be separated, i.e., samples of different conditions are expected to be more dissimilar to each other than replicates within the same condition.

Fingerprint plots

This tool samples indexed BAM files and plots a profile of cumulative read coverages for each. All reads overlapping a window (bin) of the specified length are counted; these counts are sorted and the cumulative sum is finally plotted.

It determines how well the signal in the ChIP-seq sample can be differentiated from the background distribution of reads in the control sample.

An ideal input with perfect uniform distribution of reads along the genome (i.e. without enrichments in open chromatin etc.) and infinite sequencing coverage should generate a straight diagonal line. A very specific and strong ChIP enrichment will be indicated by a prominent and steep rise of the cumulative sum towards the highest rank. This means that a big chunk of reads from the ChIP sample is located in few bins which corresponds to high, narrow enrichments typically seen for transcription factors.

Cross-correlation analysis

Information from: https://docs.google.com/document/d/1lG_Rd7fnYgRpSIqrIfuVlAz2dW1VaSQThzk836Db99c/edit

This set of programs operate on mapped Illumina single-end read datasets in tagAlign or BAM format. Because my data is paired-end I need to only use the forward read.

A high-quality ChIP-seq experiment will produce significant clustering of enriched DNA sequence tags/reads at locations bound by the protein of interest; the expectation is that we can observe a bimodal enrichment of reads (sequence tags) on both the forward and the reverse strands.

Cross-correlation analysis is done on a filtered (but not-deduped) and subsampled BAM. There is a special fastq trimming for cross-correlation analysis. Read1 fastq is trimmed to 50bp first using trimfastq.py (last modified 2017/11/08, https://github.com/ENCODE-DCC/chip-seq-pipeline2/blob/master/src/trimfastq.py). And then it is separately mapped as SE. Reads are filtered but duplicates are not removed. Then 15 million reads are randomly sampled and used for cross-correlation analysis.

ENCODE standards:

col. abbreviation description
1 Filename tagAlign/BAM filename
2 numReads effective sequencing depth i.e. total number of mapped reads in input file
3 estFragLen comma separated strand cross-correlation peak(s) in decreasing order of correlation.
4 corr_estFragLen comma separated strand cross-correlation value(s) in decreasing order (COL2 follows the same order)
5 phantomPeak Read length/phantom peak strand shift
6 corr_phantomPeak Correlation value at phantom peak
7 argmin_corr strand shift at which cross-correlation is lowest
8 min_corr minimum value of cross-correlation
9 NSC Normalized strand cross-correlation coefficient (NSC) = COL4 / COL8
10 RSC Relative strand cross-correlation coefficient (RSC) = (COL4 - COL8) / (COL6 - COL8)
11 QualityTag Quality tag based on thresholded RSC (codes= -2:veryLow, -1:Low, 0:Medium, 1:High, 2:veryHigh)

NSC; NSC>1.1 (higher values indicate more enrichment; 1 = no enrichment)

RSC; RSC>0.8 (0 = no signal; <1 low quality ChIP; >1 high enrichment

Quality tag based on thresholded RSC (codes: -2:veryLow,-1:Low,0:Medium,1:High; 2:veryHigh)

Table.

Table. Peak calling quality metrics for mouse and dunnart ChIP-seq peaks after peak calling with MACS2
file numReads estFragLen corr_estFragLen PhantomPeak corr_phantomPeak argmin_corr min_corr phantomPeakCoef relPhantomPeakCoef QualityTag
A-1_input 15000000 95,135,150 0.123657819405274,0.123466870849122,0.123426393106245 50 0.1253390 1500 0.1220062 1.013537 0.4955580 -1
A-2_H3K4me3 15000000 205,225,240 0.264180702500576,0.264046178090335,0.263856306475591 50 0.2533387 1500 0.1496324 1.765532 1.1045450 1
A-3_H3K27ac 15000000 210,220,235 0.22968659727008,0.229482796607544,0.229305284391192 55 0.2184643 1500 0.1570855 1.462175 1.1828360 1
B-1_input 15000000 115,130,150 0.128496809011662,0.128442767656697,0.12834049936991 50 0.1304347 1500 0.1268980 1.012599 0.4520538 -1
B-2_H3K4me3 15000000 210,225,245 0.316307968017767,0.316149524309277,0.315180730244617 50 0.3014805 1500 0.1675044 1.888356 1.1106720 1
B-3_H3K27ac 15000000 220,230,255 0.229222829886832,0.228995326028105,0.228262435658256 55 0.2203839 1500 0.1693505 1.353541 1.1731980 1
ENCFF011NFM_E12.5_H3K27ac 15000000 180,190,410 0.194079083135316,0.194062831994441,0.186974057921796 50 0.1893703 1500 0.1796170 1.080516 1.4827870 1
ENCFF045IPK_E10.5_H3K4me3 15000000 175 0.449054683973344 55 0.4162880 1500 0.1825270 2.460209 1.1401720 1
ENCFF124TAB_E13.5_H3K4me3 15000000 195 0.412343721762235 55 0.3704569 1500 0.1887255 2.184886 1.2304870 1
ENCFF124UYX_E10.5_H3K4me3 15000000 155,165 0.437234654565518,0.437141396190917 55 0.4121142 1500 0.1925528 2.270726 1.1144120 1
ENCFF182ZPF_E12.5_H3K4me3 15000000 185,195 0.346452943050889,0.34639588058122 55 0.3186158 1500 0.1799956 1.924786 1.2008160 1
ENCFF194ORC_E13.5_H3K27ac 15000000 190,395,425 0.175800485522065,0.159059423243546,0.158420421820687 50 0.1610650 1500 0.1480548 1.187401 2.1326190 2
ENCFF213EBC_E10.5_H3K27ac 15000000 150,160,445 0.204867659440794,0.204776302717579,0.197792501901694 55 0.2025141 1500 0.1858585 1.102277 1.1413040 1
ENCFF258KCR_E15.5_H3K4me3 15000000 190 0.403909589945092 55 0.3545045 1500 0.1796852 2.247874 1.2826070 1
ENCFF290ZNF_E13.5_H3K27ac 15000000 185,395,440 0.180342953287366,0.168514796709703,0.167835635046527 50 0.1702597 1500 0.1606540 1.122555 2.0497020 2
ENCFF327VAO_E14.5_H3K27ac 15000000 190,415,425 0.172586720457292,0.159861157605353,0.159494736213282 40 0.1618822 1500 0.1452528 1.188182 1.6437050 2
ENCFF394TZN_E12.5_H3K27ac 15000000 175,375,395 0.18623772611897,0.179380277089497,0.179074307889335 50 0.1814191 1500 0.1731371 1.075666 1.5818280 2
ENCFF401BKM_E15.5_H3K4me3 15000000 195 0.438425258985157 55 0.3745745 1500 0.1781736 2.460664 1.3251040 1
ENCFF485UDC_E13.5_H3K4me3 15000000 195 0.40296925574699 55 0.3623074 1500 0.1868044 2.157172 1.2316870 1
ENCFF512SFE_E11.5_H3K27ac 12827968 190,355,375 0.166257322984659,0.160998712648616,0.160715915726815 35 0.1617791 1500 0.1543163 1.077380 1.6000620 2
ENCFF515PKL_E11.5_H3K27ac 13753883 200,405,425 0.161626774030533,0.154448376422035,0.154114551276721 35 0.1555308 1500 0.1492136 1.083190 1.9649750 2
ENCFF548BRR_E10.5_H3K27ac 15000000 185,495,505 0.223064425357532,0.208818555781315,0.208444603188354 55 0.2166170 1500 0.1875370 1.189442 1.2217140 1
ENCFF584JFB_E15.5_H3K27ac 15000000 170,340,375 0.199001232367767,0.192417898412095,0.191830737999882 50 0.1942060 1500 0.1843538 1.079453 1.4867190 1
ENCFF665QBJ_E14.5_H3K4me3 15000000 195 0.489610077518068 40 0.4293305 1500 0.1933835 2.531809 1.2554800 1
ENCFF707WKL_E15.5_H3K27ac 15000000 175,350,400 0.200317801866649,0.192867946044151,0.191945224097497 55 0.1944893 1500 0.1836933 1.090502 1.5398700 2
ENCFF717QDV_E11.5_H3K4me3 14136366 195 0.325676561742523 40 0.2982427 1500 0.1802414 1.806891 1.2324880 1
ENCFF724DMU_E14.5_H3K4me3 15000000 210,225 0.450702967917296,0.449650385690451 40 0.3975410 1500 0.1929002 2.336457 1.2597820 1
ENCFF760QYZ_E11.5_H3K4me3 15000000 195,205 0.34817213254642,0.34812492238547 40 0.3161701 1500 0.1872628 1.859270 1.2482560 1
ENCFF902HAR_E14.5_H3K27ac 15000000 185,315,375 0.162598991064155,0.155567973364012,0.152921683350782 40 0.1533997 1500 0.1404980 1.157305 1.7130270 2
ENCFF941QJZ_E12.5_H3K4me3 15000000 185 0.311302296801822 55 0.2904834 1500 0.1864295 1.669812 1.2000780 1

Peak QC

Number of peaks

Dunnart

echo "Number of peaks A3 H3K27ac dunnart default"
less output/peaks/A-3_H3K27ac_vs_A-1_input_macs2_default_dunnart_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks A2 H3K4me3 dunnart default"
less output/peaks/A-2_H3K4me3_vs_A-1_input_macs2_default_dunnart_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks B3 H3K27ac dunnart default"
less output/peaks/B-3_H3K27ac_vs_B-1_input_macs2_default_dunnart_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks B2 H3K4me3 dunnart default"
less output/peaks/B-2_H3K4me3_vs_B-1_input_macs2_default_dunnart_downSampled_peaks.narrowPeak | wc -l

echo "Number of peaks A3 H3K27ac dunnart P < 0.01"
less output/peaks/A-3_H3K27ac_vs_A-1_input_macs2_p0.01_dunnart_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks A2 H3K4me3 dunnart P < 0.01"
less output/peaks/A-2_H3K4me3_vs_A-1_input_macs2_p0.01_dunnart_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks B3 H3K27ac dunnart P < 0.01"
less output/peaks/B-3_H3K27ac_vs_B-1_input_macs2_p0.01_dunnart_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks B2 H3K4me3 dunnart P < 0.01"
less output/peaks/B-2_H3K4me3_vs_B-1_input_macs2_p0.01_dunnart_downSampled_peaks.narrowPeak | wc -l
Number of peaks A3 H3K27ac dunnart default
74012
Number of peaks A2 H3K4me3 dunnart default
51122
Number of peaks B3 H3K27ac dunnart default
78381
Number of peaks B2 H3K4me3 dunnart default
52680
Number of peaks A3 H3K27ac dunnart P < 0.01
123095
Number of peaks A2 H3K4me3 dunnart P < 0.01
98023
Number of peaks B3 H3K27ac dunnart P < 0.01
127978
Number of peaks B2 H3K4me3 dunnart P < 0.01
96055

Mouse

H3K4me3

echo "Number of peaks rep 1 E10.5 H3K4me3"
less output/peaks/ENCFF124UYX_vs_ENCFF157KEH_E10.5_H3K4me3_macs2_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks rep 2 E10.5 H3K4me3"
less output/peaks/ENCFF045IPK_vs_ENCFF825AVI_E10.5_H3K4me3_macs2_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks rep 1 E11.5 H3K4me3"
less output/peaks/ENCFF760QYZ_vs_ENCFF184CUE_E11.5_H3K4me3_macs2_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks rep 2 E11.5 H3K4me3"
less output/peaks/ENCFF717QDV_vs_ENCFF376FGM_E11.5_H3K4me3_macs2_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks rep 1 E12.5 H3K4me3"
less output/peaks/ENCFF941QJZ_vs_ENCFF058AUT_E12.5_H3K4me3_macs2_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks rep 2 E12.5 H3K4me3"
less output/peaks/ENCFF182ZPF_vs_ENCFF203JQV_E12.5_H3K4me3_macs2_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks rep 1 E13.5 H3K4me3"
less output/peaks/ENCFF485UDC_vs_ENCFF117QRC_E13.5_H3K4me3_macs2_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks rep 2 E13.5 H3K4me3"
less output/peaks/ENCFF124TAB_vs_ENCFF248PGK_E13.5_H3K4me3_macs2_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks rep 1 E14.5 H3K4me3"
less output/peaks/ENCFF665QBJ_vs_ENCFF002HZV_E14.5_H3K4me3_macs2_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks rep 2 E14.5 H3K4me3"
less output/peaks/ENCFF724DMU_vs_ENCFF784ORI_E14.5_H3K4me3_macs2_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks rep 1 E15.5 H3K4me3"
less output/peaks/ENCFF401BKM_vs_ENCFF182XFG_E15.5_H3K4me3_macs2_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks rep 2 E15.5 H3K4me3"
less output/peaks/ENCFF258KCR_vs_ENCFF727QTS_E15.5_H3K4me3_macs2_downSampled_peaks.narrowPeak | wc -l
Number of peaks rep 1 E10.5 H3K4me3
24597
Number of peaks rep 2 E10.5 H3K4me3
22489
Number of peaks rep 1 E11.5 H3K4me3
21033
Number of peaks rep 2 E11.5 H3K4me3
20656
Number of peaks rep 1 E12.5 H3K4me3
22192
Number of peaks rep 2 E12.5 H3K4me3
21495
Number of peaks rep 1 E13.5 H3K4me3
22183
Number of peaks rep 2 E13.5 H3K4me3
22626
Number of peaks rep 1 E14.5 H3K4me3
22150
Number of peaks rep 2 E14.5 H3K4me3
21532
Number of peaks rep 1 E15.5 H3K4me3
20829
Number of peaks rep 2 E15.5 H3K4me3
20529

H3K27ac

echo "Number of peaks rep 1 E10.5 H3K27ac"
less output/peaks/ENCFF213EBC_vs_ENCFF157KEH_E10.5_H3K27ac_macs2_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks rep 2 E10.5 H3K27ac"
less output/peaks/ENCFF548BRR_vs_ENCFF825AVI_E10.5_H3K27ac_macs2_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks rep 1 E11.5 H3K27ac"
less output/peaks/ENCFF512SFE_vs_ENCFF184CUE_E11.5_H3K27ac_macs2_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks rep 2 E11.5 H3K27ac"
less output/peaks/ENCFF515PKL_vs_ENCFF376FGM_E11.5_H3K27ac_macs2_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks rep 1 E12.5 H3K27ac"
less output/peaks/ENCFF011NFM_vs_ENCFF058AUT_E12.5_H3K27ac_macs2_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks rep 2 E12.5 H3K27ac"
less output/peaks/ENCFF394TZN_vs_ENCFF203JQV_E12.5_H3K27ac_macs2_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks rep 1 E13.5 H3K27ac"
less output/peaks/ENCFF194ORC_vs_ENCFF117QRC_E13.5_H3K27ac_macs2_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks rep 2 E13.5 H3K27ac"
less output/peaks/ENCFF290ZNF_vs_ENCFF248PGK_E13.5_H3K27ac_macs2_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks rep 1 E14.5 H3K27ac"
less output/peaks/ENCFF902HAR_vs_ENCFF002HZV_E14.5_H3K27ac_macs2_downSampled_peaks.narrowPeak | wc -l 
echo "Number of peaks rep 2 E14.5 H3K27ac"
less output/peaks/ENCFF327VAO_vs_ENCFF784ORI_E14.5_H3K27ac_macs2_downSampled_peaks.narrowPeak | wc -l
echo "Number of peaks rep 1 E15.5 H3K27ac"
less output/peaks/ENCFF707WKL_vs_ENCFF182XFG_E15.5_H3K27ac_macs2_downSampled_peaks.narrowPeak | wc -l 
echo "Number of peaks rep 2 E15.5 H3K27ac"
less output/peaks/ENCFF584JFB_vs_ENCFF727QTS_E15.5_H3K27ac_macs2_downSampled_peaks.narrowPeak | wc -l
Number of peaks rep 1 E10.5 H3K27ac
22728
Number of peaks rep 2 E10.5 H3K27ac
31810
Number of peaks rep 1 E11.5 H3K27ac
10821
Number of peaks rep 2 E11.5 H3K27ac
9967
Number of peaks rep 1 E12.5 H3K27ac
19496
Number of peaks rep 2 E12.5 H3K27ac
10618
Number of peaks rep 1 E13.5 H3K27ac
33386
Number of peaks rep 2 E13.5 H3K27ac
25417
Number of peaks rep 1 E14.5 H3K27ac
24114
Number of peaks rep 2 E14.5 H3K27ac
24708
Number of peaks rep 1 E15.5 H3K27ac
17820
Number of peaks rep 2 E15.5 H3K27ac
16197

Consensus peaksets for biological replicates

Number of consensus peaks

Dunnart
echo "Number of consensus H3K4me3 dunnart peaks"
less output/peaks/H3K4me3_overlap_default_dunnart_downSampled_peaks.narrowPeak | wc -l
echo "Number of consensus H3K27ac dunnart peaks"
less output/peaks/H3K27ac_overlap_default_dunnart_downSampled_peaks.narrowPeak | wc -l
Number of consensus H3K4me3 dunnart peaks
27902
Number of consensus H3K27ac dunnart peaks
42072
Mouse
echo "Number of consensus H3K4me3 mouse E10.5 peaks"
less output/peaks/H3K4me3_E10.5_overlap.narrowPeak | wc -l
echo "Number of consensus H3K4me3 mouse E11.5 peaks"
less output/peaks/H3K4me3_E11.5_overlap.narrowPeak | wc -l
echo "Number of consensus H3K4me3 mouse E12.5 peaks"
less output/peaks/H3K4me3_E12.5_overlap.narrowPeak | wc -l
echo "Number of consensus H3K4me3 mouse E13.5 peaks"
less output/peaks/H3K4me3_E13.5_overlap.narrowPeak | wc -l
echo "Number of consensus H3K4me3 mouse E14.5 peaks"
less output/peaks/H3K4me3_E14.5_overlap.narrowPeak | wc -l
echo "Number of consensus H3K4me3 mouse E15.5 peaks"
less output/peaks/H3K4me3_E15.5_overlap.narrowPeak | wc -l
echo "Number of consensus H3K27ac mouse E10.5 peaks"
less output/peaks/H3K27ac_E10.5_overlap.narrowPeak | wc -l
echo "Number of consensus H3K27ac mouse E11.5 peaks"
less output/peaks/H3K27ac_E11.5_overlap.narrowPeak | wc -l
echo "Number of consensus H3K27ac mouse E12.5 peaks"
less output/peaks/H3K27ac_E12.5_overlap.narrowPeak | wc -l
echo "Number of consensus H3K27ac mouse E13.5 peaks"
less output/peaks/H3K27ac_E13.5_overlap.narrowPeak | wc -l
echo "Number of consensus H3K27ac mouse E14.5 peaks"
less output/peaks/H3K27ac_E14.5_overlap.narrowPeak | wc -l
echo "Number of consensus H3K27ac mouse E15.5 peaks"
less output/peaks/H3K27ac_E15.5_overlap.narrowPeak | wc -l
Number of consensus H3K4me3 mouse E10.5 peaks
output/peaks/H3K4me3_E10.5_overlap.narrowPeak: No such file or directory
0
Number of consensus H3K4me3 mouse E11.5 peaks
output/peaks/H3K4me3_E11.5_overlap.narrowPeak: No such file or directory
0
Number of consensus H3K4me3 mouse E12.5 peaks
output/peaks/H3K4me3_E12.5_overlap.narrowPeak: No such file or directory
0
Number of consensus H3K4me3 mouse E13.5 peaks
output/peaks/H3K4me3_E13.5_overlap.narrowPeak: No such file or directory
0
Number of consensus H3K4me3 mouse E14.5 peaks
output/peaks/H3K4me3_E14.5_overlap.narrowPeak: No such file or directory
0
Number of consensus H3K4me3 mouse E15.5 peaks
output/peaks/H3K4me3_E15.5_overlap.narrowPeak: No such file or directory
0
Number of consensus H3K27ac mouse E10.5 peaks
output/peaks/H3K27ac_E10.5_overlap.narrowPeak: No such file or directory
0
Number of consensus H3K27ac mouse E11.5 peaks
output/peaks/H3K27ac_E11.5_overlap.narrowPeak: No such file or directory
0
Number of consensus H3K27ac mouse E12.5 peaks
output/peaks/H3K27ac_E12.5_overlap.narrowPeak: No such file or directory
0
Number of consensus H3K27ac mouse E13.5 peaks
output/peaks/H3K27ac_E13.5_overlap.narrowPeak: No such file or directory
0
Number of consensus H3K27ac mouse E14.5 peaks
output/peaks/H3K27ac_E14.5_overlap.narrowPeak: No such file or directory
0
Number of consensus H3K27ac mouse E15.5 peaks
output/peaks/H3K27ac_E15.5_overlap.narrowPeak: No such file or directory
0

Grouping peaks into putative enhancers and promoters

Find uniquely H3K4me3 sites (i.e. peaks that don’t overlap with H3K27ac):

conda activate chip
bedtools intersect -v -a output/peaks/H3K4me3_E10.5_overlap.narrowPeak -b output/peaks/H3K27ac_E10.5_overlap.narrowPeak > output/peaks/H3K4me3_E10.5_only.narrowPeak
bedtools intersect -v -a output/peaks/H3K4me3_E11.5_overlap.narrowPeak -b output/peaks/H3K27ac_E11.5_overlap.narrowPeak > output/peaks/H3K4me3_E11.5_only.narrowPeak
bedtools intersect -v -a output/peaks/H3K4me3_E12.5_overlap.narrowPeak -b output/peaks/H3K27ac_E12.5_overlap.narrowPeak > output/peaks/H3K4me3_E12.5_only.narrowPeak
bedtools intersect -v -a output/peaks/H3K4me3_E13.5_overlap.narrowPeak -b output/peaks/H3K27ac_E13.5_overlap.narrowPeak > output/peaks/H3K4me3_E13.5_only.narrowPeak
bedtools intersect -v -a output/peaks/H3K4me3_E14.5_overlap.narrowPeak -b output/peaks/H3K27ac_E14.5_overlap.narrowPeak > output/peaks/H3K4me3_E14.5_only.narrowPeak
bedtools intersect -v -a output/peaks/H3K4me3_E15.5_overlap.narrowPeak -b output/peaks/H3K27ac_E15.5_overlap.narrowPeak > output/peaks/H3K4me3_E15.5_only.narrowPeak


bedtools intersect -v -a output/peaks/H3K4me3_overlap_default_dunnart_downSampled_peaks.narrowPeak -b output/peaks/H3K27ac_overlap_default_dunnart_downSampled_peaks.narrowPeak > output/peaks/H3K4me3_only_dunnart_downSampled.narrowPeak

Find unique H3K27ac sites (i.e. peaks that don’t overlap with H3K4me3):

bedtools intersect -v -a output/peaks/H3K27ac_E10.5_overlap.narrowPeak -b output/peaks/H3K4me3_E10.5_overlap.narrowPeak > output/peaks/H3K27ac_E10.5_only.narrowPeak
bedtools intersect -v -a output/peaks/H3K27ac_E11.5_overlap.narrowPeak -b output/peaks/H3K4me3_E11.5_overlap.narrowPeak > output/peaks/H3K27ac_E11.5_only.narrowPeak
bedtools intersect -v -a output/peaks/H3K27ac_E12.5_overlap.narrowPeak -b output/peaks/H3K4me3_E12.5_overlap.narrowPeak > output/peaks/H3K27ac_E12.5_only.narrowPeak
bedtools intersect -v -a output/peaks/H3K27ac_E13.5_overlap.narrowPeak -b output/peaks/H3K4me3_E13.5_overlap.narrowPeak > output/peaks/H3K27ac_E13.5_only.narrowPeak
bedtools intersect -v -a output/peaks/H3K27ac_E14.5_overlap.narrowPeak -b output/peaks/H3K4me3_E14.5_overlap.narrowPeak > output/peaks/H3K27ac_E14.5_only.narrowPeak
bedtools intersect -v -a output/peaks/H3K27ac_E15.5_overlap.narrowPeak -b output/peaks/H3K4me3_E15.5_overlap.narrowPeak > output/peaks/H3K27ac_E15.5_only.narrowPeak


bedtools intersect -v -a output/peaks/H3K27ac_overlap_default_dunnart_downSampled_peaks.narrowPeak -b output/peaks/H3K4me3_overlap_default_dunnart_downSampled_peaks.narrowPeak > output/peaks/H3K27ac_only_dunnart_downSampled.narrowPeak

Find peaks common between H3K27ac & H3K4me3 with a reciprocal overlap of at least 50%.

bedtools intersect -f 0.5 -r -a output/peaks/H3K4me3_E10.5_overlap.narrowPeak -b output/peaks/H3K27ac_E10.5_overlap.narrowPeak > output/peaks/H3K4me3_and_H3K27ac_E10.5.narrowPeak
bedtools intersect -f 0.5 -r -a output/peaks/H3K4me3_E11.5_overlap.narrowPeak -b output/peaks/H3K27ac_E11.5_overlap.narrowPeak > output/peaks/H3K4me3_and_H3K27ac_E11.5.narrowPeak
bedtools intersect -f 0.5 -r -a output/peaks/H3K4me3_E12.5_overlap.narrowPeak -b output/peaks/H3K27ac_E12.5_overlap.narrowPeak > output/peaks/H3K4me3_and_H3K27ac_E12.5.narrowPeak
bedtools intersect -f 0.5 -r -a output/peaks/H3K4me3_E13.5_overlap.narrowPeak -b output/peaks/H3K27ac_E13.5_overlap.narrowPeak > output/peaks/H3K4me3_and_H3K27ac_E13.5.narrowPeak
bedtools intersect -f 0.5 -r -a output/peaks/H3K4me3_E14.5_overlap.narrowPeak -b output/peaks/H3K27ac_E14.5_overlap.narrowPeak > output/peaks/H3K4me3_and_H3K27ac_E14.5.narrowPeak
bedtools intersect -f 0.5 -r -a output/peaks/H3K4me3_E15.5_overlap.narrowPeak -b output/peaks/H3K27ac_E15.5_overlap.narrowPeak > output/peaks/H3K4me3_and_H3K27ac_E15.5.narrowPeak

bedtools intersect -f 0.5 -r -a output/peaks/H3K4me3_overlap_default_dunnart_downSampled_peaks.narrowPeak -b output/peaks/H3K27ac_overlap_default_dunnart_downSampled_peaks.narrowPeak > output/peaks/H3K4me3_H3K27ac_dunnart_downSampled.narrowPeak

Combine peak files

cat output/peaks/H3K4me3_E10.5_only.narrowPeak output/peaks/H3K4me3_and_H3K27ac_E10.5.narrowPeak | sort -k 4,4 > output/peaks/E10.5_promoter_peaks.narrowPeak
cat output/peaks/H3K4me3_E11.5_only.narrowPeak output/peaks/H3K4me3_and_H3K27ac_E11.5.narrowPeak | sort -k 4,4 > output/peaks/E11.5_promoter_peaks.narrowPeak
cat output/peaks/H3K4me3_E12.5_only.narrowPeak output/peaks/H3K4me3_and_H3K27ac_E12.5.narrowPeak | sort -k 4,4 > output/peaks/E12.5_promoter_peaks.narrowPeak
cat output/peaks/H3K4me3_E13.5_only.narrowPeak output/peaks/H3K4me3_and_H3K27ac_E13.5.narrowPeak | sort -k 4,4 > output/peaks/E13.5_promoter_peaks.narrowPeak
cat output/peaks/H3K4me3_E14.5_only.narrowPeak output/peaks/H3K4me3_and_H3K27ac_E14.5.narrowPeak | sort -k 4,4 > output/peaks/E14.5_promoter_peaks.narrowPeak
cat output/peaks/H3K4me3_E15.5_only.narrowPeak output/peaks/H3K4me3_and_H3K27ac_E15.5.narrowPeak | sort -k 4,4 > output/peaks/E15.5_promoter_peaks.narrowPeak

cat output/peaks/H3K4me3_only_dunnart_downSampled.narrowPeak output/peaks/H3K4me3_H3K27ac_dunnart_downSampled.narrowPeak| sort -k 4,4 > output/peaks/dunnart_downSampled_promoter_peaks.narrowPeak


cat output/peaks/H3K27ac_E10.5_only.narrowPeak | sort -k 4,4 > output/peaks/E10.5_enhancer_peaks.narrowPeak
cat output/peaks/H3K27ac_E11.5_only.narrowPeak | sort -k 4,4 > output/peaks/E11.5_enhancer_peaks.narrowPeak
cat output/peaks/H3K27ac_E12.5_only.narrowPeak | sort -k 4,4 > output/peaks/E12.5_enhancer_peaks.narrowPeak
cat output/peaks/H3K27ac_E13.5_only.narrowPeak | sort -k 4,4 > output/peaks/E13.5_enhancer_peaks.narrowPeak
cat output/peaks/H3K27ac_E14.5_only.narrowPeak | sort -k 4,4 > output/peaks/E14.5_enhancer_peaks.narrowPeak
cat output/peaks/H3K27ac_E15.5_only.narrowPeak | sort -k 4,4 > output/peaks/E15.5_enhancer_peaks.narrowPeak
cat output/peaks/H3K27ac_only_dunnart_downSampled.narrowPeak | sort -k 4,4 > output/peaks/dunnart_downSampled_enhancer_peaks.narrowPeak
cat: output/peaks/H3K4me3_E10.5_only.narrowPeak: No such file or directory
cat: output/peaks/H3K4me3_and_H3K27ac_E10.5.narrowPeak: No such file or directory
cat: output/peaks/H3K4me3_E11.5_only.narrowPeak: No such file or directory
cat: output/peaks/H3K4me3_and_H3K27ac_E11.5.narrowPeak: No such file or directory
cat: output/peaks/H3K4me3_E12.5_only.narrowPeak: No such file or directory
cat: output/peaks/H3K4me3_and_H3K27ac_E12.5.narrowPeak: No such file or directory
cat: output/peaks/H3K4me3_E13.5_only.narrowPeak: No such file or directory
cat: output/peaks/H3K4me3_and_H3K27ac_E13.5.narrowPeak: No such file or directory
cat: output/peaks/H3K4me3_E14.5_only.narrowPeak: No such file or directory
cat: output/peaks/H3K4me3_and_H3K27ac_E14.5.narrowPeak: No such file or directory
cat: output/peaks/H3K4me3_E15.5_only.narrowPeak: No such file or directory
cat: output/peaks/H3K4me3_and_H3K27ac_E15.5.narrowPeak: No such file or directory
cat: output/peaks/H3K4me3_only_dunnart_downSampled.narrowPeak: No such file or directory
cat: output/peaks/H3K4me3_H3K27ac_dunnart_downSampled.narrowPeak: No such file or directory
cat: output/peaks/H3K27ac_E10.5_only.narrowPeak: No such file or directory
cat: output/peaks/H3K27ac_E11.5_only.narrowPeak: No such file or directory
cat: output/peaks/H3K27ac_E12.5_only.narrowPeak: No such file or directory
cat: output/peaks/H3K27ac_E13.5_only.narrowPeak: No such file or directory
cat: output/peaks/H3K27ac_E14.5_only.narrowPeak: No such file or directory
cat: output/peaks/H3K27ac_E15.5_only.narrowPeak: No such file or directory
cat: output/peaks/H3K27ac_only_dunnart_downSampled.narrowPeak: No such file or directory

sessionInfo()
R version 4.1.0 (2021-05-18)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Red Hat Enterprise Linux

Matrix products: default
BLAS/LAPACK: /usr/local/easybuild-2019/easybuild/software/compiler/gcc/10.2.0/openblas/0.3.12/lib/libopenblas_haswellp-r0.3.12.so

locale:
 [1] LC_CTYPE=en_AU.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_AU.UTF-8        LC_COLLATE=en_AU.UTF-8    
 [5] LC_MONETARY=en_AU.UTF-8    LC_MESSAGES=en_AU.UTF-8   
 [7] LC_PAPER=en_AU.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_AU.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] knitr_1.33        forcats_0.5.1     readr_1.4.0       tidyr_1.1.3      
 [5] tibble_3.1.2      ggplot2_3.3.3     tidyverse_1.3.1   purrr_0.3.4      
 [9] stringr_1.4.0     dplyr_1.0.8       data.table_1.14.0 workflowr_1.7.0  

loaded via a namespace (and not attached):
 [1] Rcpp_1.0.8        lubridate_1.7.10  getPass_0.2-2     ps_1.6.0         
 [5] assertthat_0.2.1  rprojroot_2.0.2   digest_0.6.27     utf8_1.2.1       
 [9] R6_2.5.0          cellranger_1.1.0  backports_1.2.1   reprex_2.0.0     
[13] evaluate_0.14     highr_0.9         httr_1.4.2        pillar_1.6.1     
[17] rlang_1.0.2       readxl_1.3.1      rstudioapi_0.13   whisker_0.4      
[21] callr_3.7.0       jquerylib_0.1.4   rmarkdown_2.8     munsell_0.5.0    
[25] broom_0.7.6       compiler_4.1.0    httpuv_1.6.1      modelr_0.1.8     
[29] xfun_0.23         pkgconfig_2.0.3   htmltools_0.5.1.1 tidyselect_1.1.1 
[33] fansi_0.5.0       crayon_1.4.1      dbplyr_2.1.1      withr_2.4.2      
[37] later_1.2.0       grid_4.1.0        jsonlite_1.7.2    gtable_0.3.0     
[41] lifecycle_1.0.1   DBI_1.1.1         git2r_0.28.0      magrittr_2.0.1   
[45] scales_1.1.1      cli_2.5.0         stringi_1.6.2     fs_1.5.0         
[49] promises_1.2.0.1  xml2_1.3.2        bslib_0.2.5.1     ellipsis_0.3.2   
[53] generics_0.1.0    vctrs_0.3.8       tools_4.1.0       glue_1.4.2       
[57] hms_1.1.0         processx_3.5.2    yaml_2.2.1        colorspace_2.0-1 
[61] rvest_1.0.0       haven_2.4.1       sass_0.4.0