diff --git a/example/templates/old/stable/hg19/README b/example/templates/old/stable/hg19/README
deleted file mode 100644
index 7127c06b8941dbcc978579039164c66b8ae1fc29..0000000000000000000000000000000000000000
--- a/example/templates/old/stable/hg19/README
+++ /dev/null
@@ -1,5 +0,0 @@
-The config file is now in a new and more suitable format: yaml instead of the previous json.
-
-For parameter explanations, see the yaml file and the comments above each parameter.
-
-Feel free to add as many comments to the yaml file as you want
diff --git a/example/templates/old/stable/hg19/config.yaml b/example/templates/old/stable/hg19/config.yaml
deleted file mode 100644
index 6f7d44fe47ee122710d0d669dbf99e76fe2427b9..0000000000000000000000000000000000000000
--- a/example/templates/old/stable/hg19/config.yaml
+++ /dev/null
@@ -1,230 +0,0 @@
-################################################
-################################################
-# CONFIGURATION FILE FOR THE ATAC-SEQ PIPELINE #
-################################################
-################################################
-
-# This format allows comments and is therefore easier to work with than the json format we used before.
-# Quotation marks are optional for strings. Make sure to put ": " (that is, hyphen space) as separator
-
-
-##################
-# SECTION output #
-##################
-output:
-
-  # STRING. Absolute path to the output directory. Will be created if not yet present.
-  outdir: "/FULL/PATH/TO/OUTPUT/SAME/AS/IN/CLUSTER.json"
-
-  # BOOLEAN. “true” or “false”. Default “true”. Should the pipeline use MACS2 to call peaks? If yes, peaks will be called in 3 different flavors (ENCODE, stringent, non-stringent). See the section “par_peakCalling” in the configuration file for details.
-  doPeakCalling: true
-
-  # BOOLEAN. “true” or “false”. Default “false”. Should, in addition to treating all input files separately, the pipeline merge all replicate files and do all downstream analysis in addition? If set to true, for each individual as specified in the sample table, all samples belonging to this individual ID will be merged to one file after the post-processing. This makes only sense if you have more than 1 replicate per individual.
-  alsoMergeReplicates: false
-
-  # BOOLEAN. “true” or “false”. Default “false”. Should GC bias be assessed and corrected for the final BAM files in addition to the non-corrected ones? If set to true, the GC bias will be assessed and corrected using deepTools. Additionally, all downstream steps of the pipeline will be done for both GC-corrected and original files (including peak calling, PCA, coverage, etc). This greatly increases running time of th pipeline and essentially doubles the number of files. We recommend setting this to false unless the GC bias is important for downstream analyses
-  # !!! Important: Leave at false for now, I have to fix an issue in the Snakefile that currently causes an error !!
-  correctGCBias: false
-
-  # BOOLEAN. “true” or “false”. Default “false”. Should base qualities be recalibrated using GATK to detect and correct systematic errors in base quality scores? This step can be time-consuming and needs the following other parameters: GATK_jar, knownSNPs, knownIndels. I recommend turning this off for now.
-  doBaseRecalibration: false
-
-  # BOOLEAN. “true” or “false”. Default “true”. Should the pipeline produce coverage files and diagnostic plots? If set to true, coverage files for the final BAM files in bigwig and bedgraph.gz format will be produced as well as a coverage plot using deepTools.
-  generateCoverageFiles: false
-
-  # BOOLEAN. “true” or “false”. Default “true”. Not implemented yet.
-  doIDR: false
-
-  # BOOLEAN. "true" or "false". Default "false". 
-  encodePeaks: false
-
-  # BOOLEAN. "true" or "false". Default "false". Leave as "false" for now, this is not fully integrated yet
-  annotatePeaks: false
-
-###################
-# SECTION general #
-###################
-general:
-
-  # INTEGER. Default 12. Maximum number of cores per rule. For local computation, the minimum of this value and the --cores parameter will define the number of CPUs per rule, while in a cluster setting, the minimum of this value and the number of cores on the node the jobs runs is used.
-  maxCoresPerRule: 12
-
-###################
-# SECTION samples #
-###################
-samples:
-
-  # STRING. No default. Absolute path to the sample summary file. See section 2.4 for details.
-  summaryFile: "/FULL/PATH/TO/samples.csv"
-
-  # BOOLEAN. “true” or “false”. Default “true”. Paired-end data? Single-end ATAC-Seq data is not yet supported with this pipeline. If set to "false", the Snakemake pipeline will abort in the beginning.
-  pairedEnd: true
-
-###########################
-# SECTION additionalInput #
-###########################
-additionalInput:
-
-  # STRING. Absolute path to the scripts directory. Either "/g/scb2/zaugg/carnold/Projects/AtacSeq/src/Snakemake/R" or "/g/scb2/zaugg/carnold/Projects/AtacSeq/src/Snakemake/dev/R"
-  scriptsDir: "/g/scb2/zaugg/carnold/Projects/AtacSeq/src/Snakemake/R"
-
-  # STRING. Absolute path to the adapters file for Trimmomatic in fasta format. Default “/g/scb2/zaugg/zaugg_shared/Programs/Trimmomatic-0.33/adapters/NexteraPE-PE.fa”. There is usually no need to change this unless for your experiment, this adapter file is not suited.
-  trimmomatic_adapters: "/g/scb2/zaugg/zaugg_shared/Programs/Trimmomatic-0.33/adapters/NexteraPE-PE.fa"
-
-  # STRING. Absolute path to a BED file that contains the genomic regions that should be filtered from the peaks. The default depends on the genome assembly, see the templates for details. Only needed if doPeakCalling is set to true.
-  blacklistRegions: "/g/scb2/zaugg/zaugg_shared/annotations/hg19/blacklisted/hg19-blacklist.v2.bed"
-
-  # STRING. Absolute path to a database of known polymorphic sites (SNPs and Indels, respectively). This is only needed if (1) doBaseRecalibration is set to true and (2) the genome is either hg19 or hg38 and ignored otherwise. The default depends on the genome assembly, see the templates for details. Supported formats from GATK: BCF2, BEAGLE, BED, BEDTABLE, EXAMPLEBINARY, GELITEXT, RAWHAPMAP, REFSEQ, SAMPILEUP, SAMREAD, TABLE, VCF, VCF3.
-  knownSNPs: "/g/scb2/zaugg/zaugg_shared/annotations/hg19/GATK_bundle/dbsnp_138.hg19.vcf.gz"
-  knownIndels: "/g/scb2/zaugg/zaugg_shared/annotations/hg19/GATK_bundle/Mills_and_1000G_gold_standard.indels.hg19.sites.vcf.gz"
-
-  # STRING. The default depends on the genome assembly, see the templates for details. Absolute path to the reference genome in fasta and 2bit format, respectively, both of which have to correspond to the same genome assembly version as used for the alignment as well as the database of polymorphic sites (knownSNPs and knownIndels, if applicable).
-  refGenome_fasta: "/g/scb2/zaugg/zaugg_shared/annotations/hg19/GATK_bundle/ucsc.hg19.onlyRefChr.fasta"
-  refGenome_2bit: "/g/scb2/zaugg/zaugg_shared/annotations/hg19/GATK_bundle/ucsc.hg19.2bit"
-
-  # STRING. The default depends on the genome assembly, see the templates for details. Absolute path to an genome annotation file in GTF format.
-  annotationGTF: "/g/scb2/zaugg/zaugg_shared/annotations/hg19/Gencode_v19/gencode.v19.annotation.gtf"
-
-
-#######################
-# SECTION executables #
-#######################
-
-# If run with Singularity, this section can be ignored
-
-executables:
-
-  # STRING. Default “java”. Name of the executable for Java. Java version must be at least 1.8!
-  java: "java"
-
-  # STRING. Default "/g/scb2/zaugg/carnold/Projects/AtacSeq/src/Snakemake/tools/GenomeAnalysisTK.jar". Absolute path to a JAR file for the GATK suite.
-  GATK_jar: "/g/scb2/zaugg/carnold/Projects/AtacSeq/src/Snakemake/tools/GenomeAnalysisTK.jar"
-
-  # STRING. Default "/g/scb2/zaugg/zaugg_shared/Programs/Picardtools/picardOld.jar". Absolute path to a JAR file for the Picard suite.
-  PICARD_jar: "/g/scb2/zaugg/carnold/Projects/AtacSeq/src/Snakemake/tools/picardOld.jar"
-
-
-########################
-# SECTION par_trimming #
-########################
-par_trimming:
-
-  # STRING. Default "1:30:4:5:true". ILLUMINACLIP value. See trimmomatic manual
-  trimmomatic_ILLUMINACLIP: "1:30:4:1:true"
-
-  # INTEGER. Default 3. TRAILING value. See trimmomatic manual
-  trimmomatic_trailing: 3
-
-  # INTEGER. Default 20. MINLEN value. See trimmomatic manual
-  trimmomatic_minlen: 20
-
-  # STRING. Default "phred33". Phred type. See trimmomatic manual. The "-" is added automatically by the Snakemake pipeline.
-  trimmomatic_phredType: "phred33"
-
-#####################
-# SECTION par_align #
-#####################
-par_align:
-
-  # STRING. Default "--very-sensitive". Sensitivity. Leave empty for the default sensitivity. See bowtie2 manual.
-  bowtie2_sensitivity: "--very-sensitive"
-
-  # INTEGER. Default 2000. Value for parameter X. See bowtie2 manual.
-  bowtie2_maxFragmentLength: 2000
-
-  # STRING. Default "/g/scb/zaugg/zaugg_shared/annotations/hg19/referenceGenome/Bowtie2/hg19". Absolute path to the reference genome + the prefix of the file with the index. In the default case, we use those files within /g/scb/zaugg/zaugg_shared/annotations/hg19/referenceGenome/Bowtie2 that start with hg19. See bowtie2 manual for more details.
-  bowtie2_refGenome: "/g/scb2/zaugg/zaugg_shared/annotations/hg19/referenceGenome/Bowtie2/hg19"
-
-  # STRING. Default “hg19”. Reference genome assembly version. Must match the one used by the alignment program.
-  assemblyVersion: "hg19"
-
-
-#########################
-# SECTION par_postalign #
-#########################
-par_postalign:
-
-  # INTEGER. Default 10. Minimum MAPQ score. Reads with a lower MAPQ quality will be removed during the post-alignment processing.
-  minMAPQscore: 10
-
-  # STRING. Default “LENIENT”. Value of the VALIDATION_STRINGENCY from SortSam (Picard tools). See the manual for details.
-  ValidationStringencySortSam: "LENIENT"
-
-  # STRING. Default “SILENT”. Value of the VALIDATION_STRINGENCY from MarkDuplicates (Picard tools). See the manual for details.
-  ValidationStringencyMarkDuplicates: "SILENT"
-
-  # STRING. Defauled “ID”. Used for filtering reads. Relates to the one letter abbreviations for CIGAR strings such as I for insertion and D for deletion.  Specify all the one letter abbreviations in the CIGAR string of a read here that should be filtered. "ID" keeps a read only if the CIGAR string does not contain the letters "I" and "D" (e.g., only M for example)
-  CIGAR: "ID"
-
-  # INTEGER. Default 4 and -5, respectively. Adjustment of the read start positions on the forward and reverse strand and should be a positive or negative number, respectively. See the Buenrostro paper for details.
-  adjustRSS_forward: 4
-  adjustRSS_reverse: -5
-
-#######################
-# SECTION par_scripts #
-#######################
-
-# NOTE: There is usually no need to modify the parameters here.
-
-par_scripts:
-
-  # INTEGER. Default 4000. The region size in bp that specifies what is considered within a TSS. The STATS script does a TSS enrichment test to test whether or not ATAC-Seq reads are primarily located within annotated TSS as opposed to outside of TSS regions. A value of 4000 means the region from -2kb up to +2kb of annotated TSS.
-  STATS_script_withinThr: 4000
-
-  # INTEGER. Default 1000. The size of the region adjacent to the within TSS region that is considered outside of a TSS. A value of 1000 therefore denotes the 1kb region up- and downstream of the within TSS region (from -3 to -2kb upstream and from +2 to +3 kb downstream of annotated TSS.)
-  STATS_script_outsideThr: 1000
-
-  # STRING. Default "protein_coding". Gene type to keep / do the analayses for. Allowed are gene types as specified by GENCODE. The default is "protein_coding"
-  STATS_script_geneTypesToKeep: "protein_coding"
-
-  # INTEGER. Default 600. Fragment length cutoff. All reads with a fragment length less than this value will be filtered for the purpose of this script.
-  FL_distr_script_cutoff: 600
-
-
-
-###########################
-# SECTION par_peakCalling #
-###########################
-par_peakCalling:
-
-  # STRING. Default "--nolambda –nomodel". Peak calling model for non-stringent peak calling.
-  modelNonStringent: "--nolambda --nomodel"
-
-  # STRING. Default "—nomodel”. Peak calling model for stringent peak calling.
-  modelStringent: "--nomodel"
-
-  # NUMERIC [0, 1]. Default 0.01. Minimum q-value threshold for stringent peak calling.
-  modelStringent_minQValue: 0.01
-
-  # NUMERIC [0, 1]. Default 0.1. Minimum q-value threshold for non-stringent peak calling.
-  modelNonStringent_minQValue: 0.1
-
-  # INTEGER. Default 10000. Value for slocal parameter for non-stringent peak calling.
-  modelNonStringent_slocal: 10000
-
-  # NUMERIC [0, 1]. Default 0.1. p-value threshold for ENCODE peak calling.
-  Encode_pValThreshold: 0.1
-
-  # STRING. Default "--nomodel --shift -75 --extsize 150 --broad --keep-dup all". Model for Encode peak calling (broad and gapped mode)
-  Encode_modelBroadAndGapped: "--nomodel --shift -75 --extsize 150 --broad --keep-dup all"
-
-  # STRING. Default "--nomodel --shift -75 --extsize 150 -B --SPMR --keep-dup all –call-summits". Model for Encode peak calling (narrow mode)
-  Encode_modelNarrow: "--nomodel --shift -75 --extsize 150 -B --SPMR --keep-dup all --call-summits"
-
-
-#########################
-# SECTION par_deepTools #
-#########################
-par_deepTools:
-
-  # INTEGER. The default depends on the genome assembly, see the templates for details. Length of the “mappable” genome in bp as defined by deepTools (see https://deeptools.readthedocs.io/en/develop/content/feature/effectiveGenomeSize.html).
-  effectiveGenomeSize: 2736124973
-
-  # STRING. Default “normalizeTo1x”. Either "normalizeTo1x NUMBER" OR "normalizeUsingRPKM" (note the missing (!) leading "--"), where NUMBER denotes the number of base pairs that are mappable. Beware of the mapping dependence on the read length for this number: The reported numbers on the websites are for 30bp reads, and we now have much longer reads usually. See Koehler et al. (2011) for numbers1.
-  bamCoverage_normalizationCoverage: "RPGC"
-
-  # INTEGER. Default 10. Size of the bins, in bases
-  bamCoverage_binSize: 10
-
-  # STRING. Default "--extendReads  --centerReads". Additional options that are supported by bamCoverage. Note that the "--" or "-" has to be present here.
-  bamCoverage_otherOptions: "--extendReads --centerReads"
diff --git a/example/templates/old/stable/hg38/README b/example/templates/old/stable/hg38/README
deleted file mode 100644
index 7127c06b8941dbcc978579039164c66b8ae1fc29..0000000000000000000000000000000000000000
--- a/example/templates/old/stable/hg38/README
+++ /dev/null
@@ -1,5 +0,0 @@
-The config file is now in a new and more suitable format: yaml instead of the previous json.
-
-For parameter explanations, see the yaml file and the comments above each parameter.
-
-Feel free to add as many comments to the yaml file as you want
diff --git a/example/templates/old/stable/hg38/config.yaml b/example/templates/old/stable/hg38/config.yaml
deleted file mode 100644
index baa83d0c1bbc46388cf22be3b405a734afa50314..0000000000000000000000000000000000000000
--- a/example/templates/old/stable/hg38/config.yaml
+++ /dev/null
@@ -1,230 +0,0 @@
-################################################
-################################################
-# CONFIGURATION FILE FOR THE ATAC-SEQ PIPELINE #
-################################################
-################################################
-
-# This format allows comments and is therefore easier to work with than the json format we used before.
-# Quotation marks are optional for strings. Make sure to put ": " (that is, hyphen space) as separator
-
-
-##################
-# SECTION output #
-##################
-output:
-
-  # STRING. Absolute path to the output directory. Will be created if not yet present.
-  outdir: "/FULL/PATH/TO/OUTPUT/SAME/AS/IN/CLUSTER.json"
-
-  # BOOLEAN. “true” or “false”. Default “true”. Should the pipeline use MACS2 to call peaks? If yes, peaks will be called in 3 different flavors (ENCODE, stringent, non-stringent). See the section “par_peakCalling” in the configuration file for details.
-  doPeakCalling: true
-
-  # BOOLEAN. “true” or “false”. Default “false”. Should, in addition to treating all input files separately, the pipeline merge all replicate files and do all downstream analysis in addition? If set to true, for each individual as specified in the sample table, all samples belonging to this individual ID will be merged to one file after the post-processing. This makes only sense if you have more than 1 replicate per individual.
-  alsoMergeReplicates: false
-
-  # BOOLEAN. “true” or “false”. Default “false”. Should GC bias be assessed and corrected for the final BAM files in addition to the non-corrected ones? If set to true, the GC bias will be assessed and corrected using deepTools. Additionally, all downstream steps of the pipeline will be done for both GC-corrected and original files (including peak calling, PCA, coverage, etc). This greatly increases running time of th pipeline and essentially doubles the number of files. We recommend setting this to false unless the GC bias is important for downstream analyses
-  # !!! Important: Leave at false for now, I have to fix an issue in the Snakefile that currently causes an error !!
-  correctGCBias: false
-
-  # BOOLEAN. “true” or “false”. Default “false”. Should base qualities be recalibrated using GATK to detect and correct systematic errors in base quality scores? This step can be time-consuming and needs the following other parameters: GATK_jar, knownSNPs, knownIndels. I recommend turning this off for now.
-  doBaseRecalibration: false
-
-  # BOOLEAN. “true” or “false”. Default “true”. Should the pipeline produce coverage files and diagnostic plots? If set to true, coverage files for the final BAM files in bigwig and bedgraph.gz format will be produced as well as a coverage plot using deepTools.
-  generateCoverageFiles: false
-
-  # BOOLEAN. “true” or “false”. Default “true”. Not implemented yet.
-  doIDR: false
-
-  # BOOLEAN. "true" or "false". Default "false".
-  encodePeaks: false
-
-  # BOOLEAN. "true" or "false". Default "false". Leave as "false" for now, this is not fully integrated yet
-  annotatePeaks: false
-
-###################
-# SECTION general #
-###################
-general:
-
-  # INTEGER. Default 12. Maximum number of cores per rule. For local computation, the minimum of this value and the --cores parameter will define the number of CPUs per rule, while in a cluster setting, the minimum of this value and the number of cores on the node the jobs runs is used.
-  maxCoresPerRule: 12
-
-###################
-# SECTION samples #
-###################
-samples:
-
-  # STRING. No default. Absolute path to the sample summary file. See section 2.4 for details.
-  summaryFile: "/FULL/PATH/TO/samples.csv"
-
-  # BOOLEAN. “true” or “false”. Default “true”. Paired-end data? Single-end ATAC-Seq data is not yet supported with this pipeline. If set to "false", the Snakemake pipeline will abort in the beginning.
-  pairedEnd: true
-
-###########################
-# SECTION additionalInput #
-###########################
-additionalInput:
-
-  # STRING. Absolute path to the scripts directory. Either "/g/scb2/zaugg/carnold/Projects/AtacSeq/src/Snakemake/R" or "/g/scb2/zaugg/carnold/Projects/AtacSeq/src/Snakemake/dev/R"
-  scriptsDir: "/g/scb2/zaugg/carnold/Projects/AtacSeq/src/Snakemake/R"
-
-  # STRING. Absolute path to the adapters file for Trimmomatic in fasta format. Default “/g/scb2/zaugg/zaugg_shared/Programs/Trimmomatic-0.33/adapters/NexteraPE-PE.fa”. There is usually no need to change this unless for your experiment, this adapter file is not suited.
-  trimmomatic_adapters: "/g/scb2/zaugg/zaugg_shared/Programs/Trimmomatic-0.33/adapters/NexteraPE-PE.fa"
-
-  # STRING. Absolute path to a BED file that contains the genomic regions that should be filtered from the peaks. The default depends on the genome assembly, see the templates for details. Only needed if doPeakCalling is set to true.
-  blacklistRegions: "/g/scb2/zaugg/zaugg_shared/annotations/hg38/blacklisted/hg38-blacklist.v2.bed"
-
-  # STRING. Absolute path to a database of known polymorphic sites (SNPs and Indels, respectively). This is only needed if (1) doBaseRecalibration is set to true and (2) the genome is either hg19 or hg38 and ignored otherwise. The default depends on the genome assembly, see the templates for details. Supported formats from GATK: BCF2, BEAGLE, BED, BEDTABLE, EXAMPLEBINARY, GELITEXT, RAWHAPMAP, REFSEQ, SAMPILEUP, SAMREAD, TABLE, VCF, VCF3.
-  knownSNPs: "/g/scb2/zaugg/zaugg_shared/annotations/hg38/GATK/Homo_sapiens_assembly38.dbsnp138.vcf.gz"
-  knownIndels: "/g/scb2/zaugg/zaugg_shared/annotations/hg38/GATK/Mills_and_1000G_gold_standard.indels.hg38.vcf.gz"
-
-  # STRING. The default depends on the genome assembly, see the templates for details. Absolute path to the reference genome in fasta and 2bit format, respectively, both of which have to correspond to the same genome assembly version as used for the alignment as well as the database of polymorphic sites (knownSNPs and knownIndels, if applicable).
-  refGenome_fasta: "/g/scb2/zaugg/zaugg_shared/annotations/hg38/GATK/Homo_sapiens_assembly38.fasta"
-  refGenome_2bit: "/g/scb2/zaugg/zaugg_shared/annotations/hg38/GATK/Homo_sapiens_assembly38.2bit"
-
-  # STRING. The default depends on the genome assembly, see the templates for details. Absolute path to an genome annotation file in GTF format.
-  annotationGTF: "/g/scb2/zaugg/zaugg_shared/annotations/hg38/Gencode_v34/gencode.v34.annotation.gtf"
-
-
-#######################
-# SECTION executables #
-#######################
-
-# If run with Singularity, this section can be ignored
-
-executables:
-
-  # STRING. Default “java”. Name of the executable for Java. Java version must be at least 1.8!
-  java: "java"
-
-  # STRING. Default "/g/scb2/zaugg/carnold/Projects/AtacSeq/src/Snakemake/tools/GenomeAnalysisTK.jar". Absolute path to a JAR file for the GATK suite.
-  GATK_jar: "/g/scb2/zaugg/carnold/Projects/AtacSeq/src/Snakemake/tools/GenomeAnalysisTK.jar"
-
-  # STRING. Default "/g/scb2/zaugg/zaugg_shared/Programs/Picardtools/picardOld.jar". Absolute path to a JAR file for the Picard suite.
-  PICARD_jar: "/g/scb2/zaugg/carnold/Projects/AtacSeq/src/Snakemake/tools/picardOld.jar"
-
-
-########################
-# SECTION par_trimming #
-########################
-par_trimming:
-
-  # STRING. Default "1:30:4:5:true". ILLUMINACLIP value. See trimmomatic manual
-  trimmomatic_ILLUMINACLIP: "1:30:4:1:true"
-
-  # INTEGER. Default 3. TRAILING value. See trimmomatic manual
-  trimmomatic_trailing: 3
-
-  # INTEGER. Default 20. MINLEN value. See trimmomatic manual
-  trimmomatic_minlen: 20
-
-  # STRING. Default "phred33". Phred type. See trimmomatic manual. The "-" is added automatically by the Snakemake pipeline.
-  trimmomatic_phredType: "phred33"
-
-#####################
-# SECTION par_align #
-#####################
-par_align:
-
-  # STRING. Default "--very-sensitive". Sensitivity. Leave empty for the default sensitivity. See bowtie2 manual.
-  bowtie2_sensitivity: "--very-sensitive"
-
-  # INTEGER. Default 2000. Value for parameter X. See bowtie2 manual.
-  bowtie2_maxFragmentLength: 2000
-
-  # STRING. Default "/g/scb/zaugg/zaugg_shared/annotations/hg19/referenceGenome/Bowtie2/hg19". Absolute path to the reference genome + the prefix of the file with the index. In the default case, we use those files within /g/scb/zaugg/zaugg_shared/annotations/hg19/referenceGenome/Bowtie2 that start with hg19. See bowtie2 manual for more details.
-  bowtie2_refGenome: "/g/scb2/zaugg/zaugg_shared/annotations/hg38/referenceGenome/Bowtie2/hg38"
-
-  # STRING. Default “hg38”. Reference genome assembly version. Must match the one used by the alignment program.
-  assemblyVersion: "hg38"
-
-
-#########################
-# SECTION par_postalign #
-#########################
-par_postalign:
-
-  # INTEGER. Default 10. Minimum MAPQ score. Reads with a lower MAPQ quality will be removed during the post-alignment processing.
-  minMAPQscore: 10
-
-  # STRING. Default “LENIENT”. Value of the VALIDATION_STRINGENCY from SortSam (Picard tools). See the manual for details.
-  ValidationStringencySortSam: "LENIENT"
-
-  # STRING. Default “SILENT”. Value of the VALIDATION_STRINGENCY from MarkDuplicates (Picard tools). See the manual for details.
-  ValidationStringencyMarkDuplicates: "SILENT"
-
-  # STRING. Defauled “ID”. Used for filtering reads. Relates to the one letter abbreviations for CIGAR strings such as I for insertion and D for deletion.  Specify all the one letter abbreviations in the CIGAR string of a read here that should be filtered. "ID" keeps a read only if the CIGAR string does not contain the letters "I" and "D" (e.g., only M for example)
-  CIGAR: "ID"
-
-  # INTEGER. Default 4 and -5, respectively. Adjustment of the read start positions on the forward and reverse strand and should be a positive or negative number, respectively. See the Buenrostro paper for details.
-  adjustRSS_forward: 4
-  adjustRSS_reverse: -5
-
-#######################
-# SECTION par_scripts #
-#######################
-
-# NOTE: There is usually no need to modify the parameters here.
-
-par_scripts:
-
-  # INTEGER. Default 4000. The region size in bp that specifies what is considered within a TSS. The STATS script does a TSS enrichment test to test whether or not ATAC-Seq reads are primarily located within annotated TSS as opposed to outside of TSS regions. A value of 4000 means the region from -2kb up to +2kb of annotated TSS.
-  STATS_script_withinThr: 4000
-
-  # INTEGER. Default 1000. The size of the region adjacent to the within TSS region that is considered outside of a TSS. A value of 1000 therefore denotes the 1kb region up- and downstream of the within TSS region (from -3 to -2kb upstream and from +2 to +3 kb downstream of annotated TSS.)
-  STATS_script_outsideThr: 1000
-
-  # STRING. Default "protein_coding". Gene type to keep / do the analayses for. Allowed are gene types as specified by GENCODE. The default is "protein_coding"
-  STATS_script_geneTypesToKeep: "protein_coding"
-
-  # INTEGER. Default 600. Fragment length cutoff. All reads with a fragment length less than this value will be filtered for the purpose of this script.
-  FL_distr_script_cutoff: 600
-
-
-
-###########################
-# SECTION par_peakCalling #
-###########################
-par_peakCalling:
-
-  # STRING. Default "--nolambda –nomodel". Peak calling model for non-stringent peak calling.
-  modelNonStringent: "--nolambda --nomodel"
-
-  # STRING. Default "—nomodel”. Peak calling model for stringent peak calling.
-  modelStringent: "--nomodel"
-
-  # NUMERIC [0, 1]. Default 0.01. Minimum q-value threshold for stringent peak calling.
-  modelStringent_minQValue: 0.01
-
-  # NUMERIC [0, 1]. Default 0.1. Minimum q-value threshold for non-stringent peak calling.
-  modelNonStringent_minQValue: 0.1
-
-  # INTEGER. Default 10000. Value for slocal parameter for non-stringent peak calling.
-  modelNonStringent_slocal: 10000
-
-  # NUMERIC [0, 1]. Default 0.1. p-value threshold for ENCODE peak calling.
-  Encode_pValThreshold: 0.1
-
-  # STRING. Default "--nomodel --shift -75 --extsize 150 --broad --keep-dup all". Model for Encode peak calling (broad and gapped mode)
-  Encode_modelBroadAndGapped: "--nomodel --shift -75 --extsize 150 --broad --keep-dup all"
-
-  # STRING. Default "--nomodel --shift -75 --extsize 150 -B --SPMR --keep-dup all –call-summits". Model for Encode peak calling (narrow mode)
-  Encode_modelNarrow: "--nomodel --shift -75 --extsize 150 -B --SPMR --keep-dup all --call-summits"
-
-
-#########################
-# SECTION par_deepTools #
-#########################
-par_deepTools:
-
-  # INTEGER. The default depends on the genome assembly, see the templates for details. Length of the “mappable” genome in bp as defined by deepTools (see https://deeptools.readthedocs.io/en/develop/content/feature/effectiveGenomeSize.html).
-  effectiveGenomeSize: 2747877777
-
-  # STRING. Default “normalizeTo1x”. Either "normalizeTo1x NUMBER" OR "normalizeUsingRPKM" (note the missing (!) leading "--"), where NUMBER denotes the number of base pairs that are mappable. Beware of the mapping dependence on the read length for this number: The reported numbers on the websites are for 30bp reads, and we now have much longer reads usually. See Koehler et al. (2011) for numbers1.
-  bamCoverage_normalizationCoverage: "RPGC"
-
-  # INTEGER. Default 10. Size of the bins, in bases
-  bamCoverage_binSize: 10
-
-  # STRING. Default "--extendReads  --centerReads". Additional options that are supported by bamCoverage. Note that the "--" or "-" has to be present here.
-  bamCoverage_otherOptions: "--extendReads --centerReads"
diff --git a/example/templates/old/stable/mm10/README b/example/templates/old/stable/mm10/README
deleted file mode 100644
index 7127c06b8941dbcc978579039164c66b8ae1fc29..0000000000000000000000000000000000000000
--- a/example/templates/old/stable/mm10/README
+++ /dev/null
@@ -1,5 +0,0 @@
-The config file is now in a new and more suitable format: yaml instead of the previous json.
-
-For parameter explanations, see the yaml file and the comments above each parameter.
-
-Feel free to add as many comments to the yaml file as you want
diff --git a/example/templates/old/stable/mm10/config.yaml b/example/templates/old/stable/mm10/config.yaml
deleted file mode 100644
index 2a1934aa5d0b199dcb88a05e9e11d4fb33efa683..0000000000000000000000000000000000000000
--- a/example/templates/old/stable/mm10/config.yaml
+++ /dev/null
@@ -1,230 +0,0 @@
-################################################
-################################################
-# CONFIGURATION FILE FOR THE ATAC-SEQ PIPELINE #
-################################################
-################################################
-
-# This format allows comments and is therefore easier to work with than the json format we used before.
-# Quotation marks are optional for strings. Make sure to put ": " (that is, hyphen space) as separator
-
-
-##################
-# SECTION output #
-##################
-output:
-
-  # STRING. Absolute path to the output directory. Will be created if not yet present.
-  outdir: "/FULL/PATH/TO/OUTPUT/SAME/AS/IN/CLUSTER.json"
-
-  # BOOLEAN. “true” or “false”. Default “true”. Should the pipeline use MACS2 to call peaks? If yes, peaks will be called in 3 different flavors (ENCODE, stringent, non-stringent). See the section “par_peakCalling” in the configuration file for details.
-  doPeakCalling: true
-
-  # BOOLEAN. “true” or “false”. Default “false”. Should, in addition to treating all input files separately, the pipeline merge all replicate files and do all downstream analysis in addition? If set to true, for each individual as specified in the sample table, all samples belonging to this individual ID will be merged to one file after the post-processing. This makes only sense if you have more than 1 replicate per individual.
-  alsoMergeReplicates: false
-
-  # BOOLEAN. “true” or “false”. Default “false”. Should GC bias be assessed and corrected for the final BAM files in addition to the non-corrected ones? If set to true, the GC bias will be assessed and corrected using deepTools. Additionally, all downstream steps of the pipeline will be done for both GC-corrected and original files (including peak calling, PCA, coverage, etc). This greatly increases running time of th pipeline and essentially doubles the number of files. We recommend setting this to false unless the GC bias is important for downstream analyses
-  # !!! Important: Leave at false for now, I have to fix an issue in the Snakefile that currently causes an error !!
-  correctGCBias: false
-
-  # BOOLEAN. “true” or “false”. Default “false”. Should base qualities be recalibrated using GATK to detect and correct systematic errors in base quality scores? This step can be time-consuming and needs the following other parameters: GATK_jar, knownSNPs, knownIndels. I recommend turning this off for now.
-  doBaseRecalibration: false
-
-  # BOOLEAN. “true” or “false”. Default “true”. Should the pipeline produce coverage files and diagnostic plots? If set to true, coverage files for the final BAM files in bigwig and bedgraph.gz format will be produced as well as a coverage plot using deepTools.
-  generateCoverageFiles: false
-
-  # BOOLEAN. “true” or “false”. Default “true”. Not implemented yet.
-  doIDR: false
-
-  # BOOLEAN. "true" or "false". Default "false".
-  encodePeaks: false
-
-  # BOOLEAN. "true" or "false". Default "false". Leave as "false" for now, this is not fully integrated yet
-  annotatePeaks: false
-
-###################
-# SECTION general #
-###################
-general:
-
-  # INTEGER. Default 12. Maximum number of cores per rule. For local computation, the minimum of this value and the --cores parameter will define the number of CPUs per rule, while in a cluster setting, the minimum of this value and the number of cores on the node the jobs runs is used.
-  maxCoresPerRule: 12
-
-###################
-# SECTION samples #
-###################
-samples:
-
-  # STRING. No default. Absolute path to the sample summary file. See section 2.4 for details.
-  summaryFile: "/FULL/PATH/TO/samples.csv"
-
-  # BOOLEAN. “true” or “false”. Default “true”. Paired-end data? Single-end ATAC-Seq data is not yet supported with this pipeline. If set to "false", the Snakemake pipeline will abort in the beginning.
-  pairedEnd: true
-
-###########################
-# SECTION additionalInput #
-###########################
-additionalInput:
-
-  # STRING. Absolute path to the scripts directory. Either "/g/scb2/zaugg/carnold/Projects/AtacSeq/src/Snakemake/R" or "/g/scb2/zaugg/carnold/Projects/AtacSeq/src/Snakemake/dev/R"
-  scriptsDir: "/g/scb2/zaugg/carnold/Projects/AtacSeq/src/Snakemake/R"
-
-  # STRING. Absolute path to the adapters file for Trimmomatic in fasta format. Default “/g/scb2/zaugg/zaugg_shared/Programs/Trimmomatic-0.33/adapters/NexteraPE-PE.fa”. There is usually no need to change this unless for your experiment, this adapter file is not suited.
-  trimmomatic_adapters: "/g/scb2/zaugg/zaugg_shared/Programs/Trimmomatic-0.33/adapters/NexteraPE-PE.fa"
-
-  # STRING. Absolute path to a BED file that contains the genomic regions that should be filtered from the peaks. The default depends on the genome assembly, see the templates for details. Only needed if doPeakCalling is set to true.
-  blacklistRegions: "/g/scb2/zaugg/zaugg_shared/annotations/mm10/blacklisted/mm10-blacklist.v2.bed"
-
-  # STRING. Absolute path to a database of known polymorphic sites (SNPs and Indels, respectively). This is only needed if (1) doBaseRecalibration is set to true and (2) the genome is either hg19 or hg38 and ignored otherwise. Thus, not needed for mm10
-  knownSNPs: ""
-  knownIndels: ""
-
-  # STRING. The default depends on the genome assembly, see the templates for details. Absolute path to the reference genome in fasta and 2bit format, respectively, both of which have to correspond to the same genome assembly version as used for the alignment as well as the database of polymorphic sites (knownSNPs and knownIndels, if applicable).
-  refGenome_fasta: "/g/scb2/zaugg/zaugg_shared/annotations/mm10/reference/mm10.fa"
-  refGenome_2bit: "/g/scb2/zaugg/zaugg_shared/annotations/mm10/reference/mm10.2bit"
-
-  # STRING. The default depends on the genome assembly, see the templates for details. Absolute path to an genome annotation file in GTF format.
-  annotationGTF: "/g/scb2/zaugg/zaugg_shared/annotations/mm10/Gencode_M16/gencode.vM16.annotation.gtf"
-
-
-#######################
-# SECTION executables #
-#######################
-
-# If run with Singularity, this section can be ignored
-
-executables:
-
-  # STRING. Default “java”. Name of the executable for Java. Java version must be at least 1.8!
-  java: "java"
-
-  # STRING. Default "/g/scb2/zaugg/carnold/Projects/AtacSeq/src/Snakemake/tools/GenomeAnalysisTK.jar". Absolute path to a JAR file for the GATK suite.
-  GATK_jar: "/g/scb2/zaugg/carnold/Projects/AtacSeq/src/Snakemake/tools/GenomeAnalysisTK.jar"
-
-  # STRING. Default "/g/scb2/zaugg/zaugg_shared/Programs/Picardtools/picardOld.jar". Absolute path to a JAR file for the Picard suite.
-  PICARD_jar: "/g/scb2/zaugg/carnold/Projects/AtacSeq/src/Snakemake/tools/picardOld.jar"
-
-
-########################
-# SECTION par_trimming #
-########################
-par_trimming:
-
-  # STRING. Default "1:30:4:5:true". ILLUMINACLIP value. See trimmomatic manual
-  trimmomatic_ILLUMINACLIP: "1:30:4:1:true"
-
-  # INTEGER. Default 3. TRAILING value. See trimmomatic manual
-  trimmomatic_trailing: 3
-
-  # INTEGER. Default 20. MINLEN value. See trimmomatic manual
-  trimmomatic_minlen: 20
-
-  # STRING. Default "phred33". Phred type. See trimmomatic manual. The "-" is added automatically by the Snakemake pipeline.
-  trimmomatic_phredType: "phred33"
-
-#####################
-# SECTION par_align #
-#####################
-par_align:
-
-  # STRING. Default "--very-sensitive". Sensitivity. Leave empty for the default sensitivity. See bowtie2 manual.
-  bowtie2_sensitivity: "--very-sensitive"
-
-  # INTEGER. Default 2000. Value for parameter X. See bowtie2 manual.
-  bowtie2_maxFragmentLength: 2000
-
-  # STRING. Default "/g/scb/zaugg/zaugg_shared/annotations/hg19/referenceGenome/Bowtie2/hg19". Absolute path to the reference genome + the prefix of the file with the index. In the default case, we use those files within /g/scb/zaugg/zaugg_shared/annotations/hg19/referenceGenome/Bowtie2 that start with hg19. See bowtie2 manual for more details.
-  bowtie2_refGenome: "/g/scb2/zaugg/zaugg_shared/annotations/mm10/Bowtie2_new_useThis/mm10"
-
-  # STRING. Default “mm10”. Reference genome assembly version. Must match the one used by the alignment program.
-  assemblyVersion: "mm10"
-
-
-#########################
-# SECTION par_postalign #
-#########################
-par_postalign:
-
-  # INTEGER. Default 10. Minimum MAPQ score. Reads with a lower MAPQ quality will be removed during the post-alignment processing.
-  minMAPQscore: 10
-
-  # STRING. Default “LENIENT”. Value of the VALIDATION_STRINGENCY from SortSam (Picard tools). See the manual for details.
-  ValidationStringencySortSam: "LENIENT"
-
-  # STRING. Default “SILENT”. Value of the VALIDATION_STRINGENCY from MarkDuplicates (Picard tools). See the manual for details.
-  ValidationStringencyMarkDuplicates: "SILENT"
-
-  # STRING. Defauled “ID”. Used for filtering reads. Relates to the one letter abbreviations for CIGAR strings such as I for insertion and D for deletion.  Specify all the one letter abbreviations in the CIGAR string of a read here that should be filtered. "ID" keeps a read only if the CIGAR string does not contain the letters "I" and "D" (e.g., only M for example)
-  CIGAR: "ID"
-
-  # INTEGER. Default 4 and -5, respectively. Adjustment of the read start positions on the forward and reverse strand and should be a positive or negative number, respectively. See the Buenrostro paper for details.
-  adjustRSS_forward: 4
-  adjustRSS_reverse: -5
-
-#######################
-# SECTION par_scripts #
-#######################
-
-# NOTE: There is usually no need to modify the parameters here.
-
-par_scripts:
-
-  # INTEGER. Default 4000. The region size in bp that specifies what is considered within a TSS. The STATS script does a TSS enrichment test to test whether or not ATAC-Seq reads are primarily located within annotated TSS as opposed to outside of TSS regions. A value of 4000 means the region from -2kb up to +2kb of annotated TSS.
-  STATS_script_withinThr: 4000
-
-  # INTEGER. Default 1000. The size of the region adjacent to the within TSS region that is considered outside of a TSS. A value of 1000 therefore denotes the 1kb region up- and downstream of the within TSS region (from -3 to -2kb upstream and from +2 to +3 kb downstream of annotated TSS.)
-  STATS_script_outsideThr: 1000
-
-  # STRING. Default "protein_coding". Gene type to keep / do the analayses for. Allowed are gene types as specified by GENCODE. The default is "protein_coding"
-  STATS_script_geneTypesToKeep: "protein_coding"
-
-  # INTEGER. Default 600. Fragment length cutoff. All reads with a fragment length less than this value will be filtered for the purpose of this script.
-  FL_distr_script_cutoff: 600
-
-
-
-###########################
-# SECTION par_peakCalling #
-###########################
-par_peakCalling:
-
-  # STRING. Default "--nolambda –nomodel". Peak calling model for non-stringent peak calling.
-  modelNonStringent: "--nolambda --nomodel"
-
-  # STRING. Default "—nomodel”. Peak calling model for stringent peak calling.
-  modelStringent: "--nomodel"
-
-  # NUMERIC [0, 1]. Default 0.01. Minimum q-value threshold for stringent peak calling.
-  modelStringent_minQValue: 0.01
-
-  # NUMERIC [0, 1]. Default 0.1. Minimum q-value threshold for non-stringent peak calling.
-  modelNonStringent_minQValue: 0.1
-
-  # INTEGER. Default 10000. Value for slocal parameter for non-stringent peak calling.
-  modelNonStringent_slocal: 10000
-
-  # NUMERIC [0, 1]. Default 0.1. p-value threshold for ENCODE peak calling.
-  Encode_pValThreshold: 0.1
-
-  # STRING. Default "--nomodel --shift -75 --extsize 150 --broad --keep-dup all". Model for Encode peak calling (broad and gapped mode)
-  Encode_modelBroadAndGapped: "--nomodel --shift -75 --extsize 150 --broad --keep-dup all"
-
-  # STRING. Default "--nomodel --shift -75 --extsize 150 -B --SPMR --keep-dup all –call-summits". Model for Encode peak calling (narrow mode)
-  Encode_modelNarrow: "--nomodel --shift -75 --extsize 150 -B --SPMR --keep-dup all --call-summits"
-
-
-#########################
-# SECTION par_deepTools #
-#########################
-par_deepTools:
-
-  # INTEGER. The default depends on the genome assembly, see the templates for details. Length of the “mappable” genome in bp as defined by deepTools (see https://deeptools.readthedocs.io/en/develop/content/feature/effectiveGenomeSize.html).
-  effectiveGenomeSize: 2407883318
-
-  # STRING. Default “normalizeTo1x”. Either "normalizeTo1x NUMBER" OR "normalizeUsingRPKM" (note the missing (!) leading "--"), where NUMBER denotes the number of base pairs that are mappable. Beware of the mapping dependence on the read length for this number: The reported numbers on the websites are for 30bp reads, and we now have much longer reads usually. See Koehler et al. (2011) for numbers1.
-  bamCoverage_normalizationCoverage: "RPGC"
-
-  # INTEGER. Default 10. Size of the bins, in bases
-  bamCoverage_binSize: 10
-
-  # STRING. Default "--extendReads  --centerReads". Additional options that are supported by bamCoverage. Note that the "--" or "-" has to be present here.
-  bamCoverage_otherOptions: "--extendReads --centerReads"
diff --git a/src/R/._consensusPeaks.R b/src/R/._consensusPeaks.R
deleted file mode 100755
index d70a0fb74ebb41a558cd74fdff0ad6c439b47412..0000000000000000000000000000000000000000
Binary files a/src/R/._consensusPeaks.R and /dev/null differ
diff --git a/src/R/.dropbox.attr b/src/R/.dropbox.attr
deleted file mode 100755
index 9e26dfeeb6e641a33dae4961196235bdb965b21b..0000000000000000000000000000000000000000
--- a/src/R/.dropbox.attr
+++ /dev/null
@@ -1 +0,0 @@
-{}
\ No newline at end of file