forked from mdelcorvo/TOSCA
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Snakefile
37 lines (29 loc) · 1.14 KB
/
Snakefile
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
# Copyright 2021 Marcello Del Corvo.
# Licensed under the MIT license (http://opensource.org/licenses/MIT)
# This file may not be copied, modified, or distributed
# except according to those terms.
include: "rules/functions.smk"
# configfile: "./config.yaml"
outputdir = getpath(config["output"])
import pandas as pd
from snakemake.utils import validate
from snakemake.utils import min_version
samples = pd.read_table(config["meta"], dtype=str).set_index(["sample"], drop=False)
if "meta_N" in config:
normals = pd.read_table(config["meta_N"], dtype=str).set_index(["id"], drop=False)
normaldir = getpath(config["normDB"])
rule all:
input:
outputdir + "results/Report.html",
outputdir + "results/Somatic_Prediction.txt",
expand(outputdir + "PureCN/{sample}.rds",sample=samples["sample"]) if "meta_N" in config else []
##### Modules #####
include: "rules/genome.smk"
include: "rules/database.smk"
include: "rules/mapping.smk"
include: "rules/calling.smk"
include: "rules/filtering.smk"
include: "rules/qc.smk"
if "meta_N" in config:
include: "rules/normal.smk"
include: "rules/purecn.smk"