forked from epigen/genome_tracks
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathsnakefile
18 lines (15 loc) · 958 Bytes
/
snakefile
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
import os
import pandas as pd
import pathlib
from snakemake.utils import Paramspace
samples = Paramspace(pd.read_csv("sample_metadata.csv"))
sample_names = set(samples.sample_name)
configfile: "config.yml"
# TODO: define expand parameters from config file
rule all:
input:
bigwig = expand("results/UCSCGenomeBrowser/{species}/{genome}/bw/{sample}_{species}_{direction}.bw", species = config["genomes"].keys(), genome = config["genomes"]['cmv']['genbank'], sample = sample_names, direction = ['for', 'rev']),
genomes_file = expand("results/UCSCGenomeBrowser/{species}/genomes.txt", species = config["genomes"].keys()),
hub_file = expand("results/UCSCGenomeBrowser/{species}/hub.txt", species = config["genomes"].keys()),
trackdb_file = expand("results/UCSCGenomeBrowser/{species}/{genome}/trackDb.txt", species = config["genomes"].keys(), genome = config["genomes"]['cmv']['genbank'])
include: "rules/visualization.smk"