Analysis of EMSeq data for project Boulard - Stork
Until count of methylation on tiles, promoters and CpG Islands, analysis is performed using snakemake pipeline in snake-make/Snakefile.
Snakefile is run on the EMBL cluster using pipeline_wrapper.sh in src/sh
Snakefile uses conda environments in env/conda and singularity containers built with recipes in singularity/recipes
config file for Snakefile and for SLURM are in config/
important for reproducibility:
conda version 4.8.3
singularity version used to build containers is 3.5.3
snakemake version to run Snakefile is 5.9.1
Using the Rdata files produced in Snakefile, I then continue the analysis in the Rmd files, which can run both on the cluster and locally on a personal computer (having at least 16Gb of RAM).
Steps outside of the Snakemake pipeline and the Rmd files
Definition of promoters as HCG or LCG and creation of 'HCG_transcripts.txt' and 'LCG_transcripts.txt' files is done in R script 'src/R/TSSs_CpG.R'.
The bed file containing genome tiles of 200 and 100 consecutive CpGs which are used in the 'CpG_tiles_counts.R' script (run in a step of the Snakemake pipeline) are produced using SeqMonk Read Position Probe Generator, asking for a minimum of 1 read count per cytocine over all sample.
The bed file containing coverage outliers used in the 'CpG_tiles_counts.R' script were also produced using SeqMonk, by quantifying methylation over 25 kb windows and picking those regions with methylation level 10x above median.