Skip to content

Commit

Permalink
fix: omit deeptools_qc also for scRNA-seq in case of single samples
Browse files Browse the repository at this point in the history
  • Loading branch information
steffenheyne committed Jul 25, 2017
1 parent 6eb5cb0 commit 5ee4675
Showing 1 changed file with 14 additions and 10 deletions.
24 changes: 14 additions & 10 deletions workflows/scRNAseq/Snakefile
Original file line number Diff line number Diff line change
Expand Up @@ -46,7 +46,19 @@ def run_Trimming(trim):
def run_bamCoverage():
return( expand("bamCoverage/{sample}.coverage.bw", sample = samples) )


def run_deeptools_qc():
file_list = [
expand("bamCoverage/{sample}.RPKM.bw", sample = samples),
expand("bamCoverage/{sample}.coverage.bw", sample = samples),
"deepTools_qc/plotEnrichment/plotEnrichment.png",
"deepTools_qc/plotEnrichment/plotEnrichment.tsv"]
if len(samples)>1:
file_list.append( ["deepTools_qc/multiBigwigSummary/coverage.bed.npz",
"deepTools_qc/plotCorrelation/correlation.pearson.bed_coverage.heatmap.png",
"deepTools_qc/plotCorrelation/correlation.spearman.bed_coverage.heatmap.png",
"deepTools_qc/plotPCA/PCA.bed_coverage.png"] )
return(file_list)

### execute before workflow starts #############################################
################################################################################
onstart:
Expand Down Expand Up @@ -93,15 +105,7 @@ rule all:
expand("Counts/{sample}.featureCounts_summary.txt",sample = samples),
"Results/all_samples.gencode_genomic.coutt_merged.csv",
"QC_report/QC_report.all_samples.libstats_reads.tsv",
expand("bamCoverage/{sample}.RPKM.bw", sample = samples),
expand("bamCoverage/{sample}.coverage.bw", sample = samples),
"deepTools_qc/plotEnrichment/plotEnrichment.png",
"deepTools_qc/plotEnrichment/plotEnrichment.tsv",
"deepTools_qc/multiBigwigSummary/coverage.bed.npz",
"deepTools_qc/plotCorrelation/correlation.pearson.bed_coverage.heatmap.png",
"deepTools_qc/plotCorrelation/correlation.spearman.bed_coverage.heatmap.png",
"deepTools_qc/plotPCA/PCA.bed_coverage.png",

run_deeptools_qc(),
rule annot:
input:
"Annotation/genes.filtered.bed",
Expand Down

0 comments on commit 5ee4675

Please sign in to comment.