diff --git a/dunnart/Snakefile b/dunnart/Snakefile index 816e9a9163345ccc13bcddf47055bc9be6ccb81f..cd994e1b5abdec9bf6e1aadd76a48611c6f9cbe8 100644 --- a/dunnart/Snakefile +++ b/dunnart/Snakefile @@ -56,6 +56,10 @@ rule all: expand("results/bowtie2/{sample}_PPq30.sorted.dupmark.bam", sample=all_samples), expand("results/bowtie2/{sample}_PPq30.sorted.dedup.bam", sample=all_samples), expand("results/bowtie2/{sample}_PPq30.sorted.dedup.bai", sample=all_samples), + "results/bowtie2/H3K4me3_pooled_PPq30.sorted.dedup.bam", + "results/bowtie2/H3K27ac_pooled_PPq30.sorted.dedup.bam", + "logs/H3K4me3.mergeBAM", + "logs/H3K27ac.mergeBAM", expand("results/qc/{sample}.unfiltered.flagstat.qc", sample=all_samples), expand("results/qc/{sample}.dedup.flagstat.qc", sample=all_samples), expand("results/qc/{sample}.dupmark.flagstat.qc", sample=all_samples), @@ -191,8 +195,8 @@ rule indexBam: rule mergeBAMreplicates: input: - H3K4me3 = ["results/bowtie2/A-2_H3K4me3_PPq30.sorted.dedup.bam, results/bowtie2/B-2_H3K4me3_PPq30.sorted.dedup.bam"], - H3K27ac = ["results/bowtie2/A-3_H3K27ac_PPq30.sorted.dedup.bam, results/bowtie2/B-3_H3K27ac_PPq30.sorted.dedup.bam"] + H3K4me3 = ["results/bowtie2/A-2_H3K4me3_PPq30.sorted.dedup.bam", "results/bowtie2/B-2_H3K4me3_PPq30.sorted.dedup.bam"], + H3K27ac = ["results/bowtie2/A-3_H3K27ac_PPq30.sorted.dedup.bam", "results/bowtie2/B-3_H3K27ac_PPq30.sorted.dedup.bam"] output: H3K4me3 = "results/bowtie2/H3K4me3_pooled_PPq30.sorted.dedup.bam", H3K27ac = "results/bowtie2/H3K27ac_pooled_PPq30.sorted.dedup.bam"