Commit 8062d7cb authored by Kevin Kunzmann's avatar Kevin Kunzmann

made snakefile more robust to filename changes

parent fbe14c11
...@@ -30,13 +30,14 @@ rule download_imputed_genotype_chromosome: ...@@ -30,13 +30,14 @@ rule download_imputed_genotype_chromosome:
shell: shell:
""" """
mkdir -p output/imputed-genotypes mkdir -p output/imputed-genotypes
gsutil cp gs://fimm-horizon-outgoing-data/CENTER_TBI_data_freeze_190829/Imputed_data/CENTER_TBI_imputed_3695_1K_MAC1_freeze_190829_chr{wildcards.i}.vcf.gz \ gsutil cp gs://fimm-horizon-outgoing-data/CENTER_TBI_data_freeze_191014/Imputed_data/CENTER.TBI_imputed_N3741_hg19_HRC_MAC1_freeze_191014_chr1.vcf.gz
output/imputed-genotypes/CENTER_TBI_imputed_3695_1K_MAC1_freeze_190829_chr{wildcards.i}.vcf.gz \
output/imputed-genotypes/CENTER_TBI_imputed_genotype_chr{wildcards.i}.vcf.gz
""" """
rule download_imputed_genotypes: rule download_imputed_genotypes:
input: input:
expand("output/imputed-genotypes/CENTER_TBI_imputed_3695_1K_MAC1_freeze_190829_chr{i}.vcf.gz", expand("output/imputed-genotypes/CENTER_TBI_imputed_genotype_chr{i}.vcf.gz",
i = list(map(str, range(1, 23))) + ['X'] i = list(map(str, range(1, 23))) + ['X']
) )
...@@ -53,7 +54,7 @@ rule download_container: ...@@ -53,7 +54,7 @@ rule download_container:
rule vcf_to_dosages: rule vcf_to_dosages:
input: input:
"container.sif", "container.sif",
vcf_gz_file = "output/imputed-genotypes/CENTER_TBI_imputed_3695_1K_MAC1_freeze_190829_chr{i}.vcf.gz" vcf_gz_file = "output/imputed-genotypes/CENTER_TBI_imputed_genotype_chr{i}.vcf.gz"
output: output:
"{output_dir}/dosages/chr{i}.dosage.txt.gz" "{output_dir}/dosages/chr{i}.dosage.txt.gz"
singularity: singularity:
...@@ -79,7 +80,7 @@ rule vcf_to_dosages: ...@@ -79,7 +80,7 @@ rule vcf_to_dosages:
rule generate_samples_file: rule generate_samples_file:
input: input:
"container.sif", "container.sif",
vcf_gz_file = "output/imputed-genotypes/CENTER_TBI_imputed_3695_1K_MAC1_freeze_190829_chr1.vcf.gz" vcf_gz_file = "output/imputed-genotypes/CENTER_TBI_imputed_genotype_chr1.vcf.gz"
output: output:
"{output_dir}/dosages/samples.txt" "{output_dir}/dosages/samples.txt"
singularity: singularity:
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment