File tree 6 files changed +18
-18
lines changed
6 files changed +18
-18
lines changed Original file line number Diff line number Diff line change 1
1
# our samples are pre-configured
2
- SAMPLES = ["A", "B"]
2
+ samples = ["A", "B"]
3
3
4
4
rule all:
5
5
input:
@@ -38,8 +38,8 @@ rule samtools_index:
38
38
rule bcftools_call:
39
39
input:
40
40
fa="data/genome.fa",
41
- bam=expand("sorted_reads/{sample}.bam", sample=SAMPLES ),
42
- bai=expand("sorted_reads/{sample}.bam.bai", sample=SAMPLES )
41
+ bam=expand("sorted_reads/{sample}.bam", sample=samples ),
42
+ bai=expand("sorted_reads/{sample}.bam.bai", sample=samples )
43
43
output:
44
44
"calls/all.vcf"
45
45
shell:
Original file line number Diff line number Diff line change 1
1
# our samples are pre-configured
2
- SAMPLES = ["A", "B"]
2
+ samples = ["A", "B"]
3
3
4
4
5
5
rule all:
@@ -40,8 +40,8 @@ rule samtools_index:
40
40
rule bcftools_call:
41
41
input:
42
42
fa="data/genome.fa",
43
- bam=expand("sorted_reads/{sample}.bam", sample=SAMPLES ),
44
- bai=expand("sorted_reads/{sample}.bam.bai", sample=SAMPLES )
43
+ bam=expand("sorted_reads/{sample}.bam", sample=samples ),
44
+ bai=expand("sorted_reads/{sample}.bam.bai", sample=samples )
45
45
output:
46
46
"calls/all.vcf"
47
47
shell:
Original file line number Diff line number Diff line change 1
1
# our samples are pre-configured
2
- SAMPLES = ["A", "B"]
2
+ samples = ["A", "B"]
3
3
4
4
configfile: "config/config.yaml"
5
5
@@ -41,8 +41,8 @@ rule samtools_index:
41
41
rule bcftools_call:
42
42
input:
43
43
fa="data/genome.fa",
44
- bam=expand("sorted_reads/{sample}.bam", sample=SAMPLES ),
45
- bai=expand("sorted_reads/{sample}.bam.bai", sample=SAMPLES )
44
+ bam=expand("sorted_reads/{sample}.bam", sample=samples ),
45
+ bai=expand("sorted_reads/{sample}.bam.bai", sample=samples )
46
46
output:
47
47
"calls/all.vcf"
48
48
shell:
Original file line number Diff line number Diff line change 1
1
# our samples are pre-configured
2
- SAMPLES = ["A", "B"]
2
+ samples = ["A", "B"]
3
3
4
4
configfile: "config/config.yaml"
5
5
@@ -41,8 +41,8 @@ rule samtools_index:
41
41
rule bcftools_call:
42
42
input:
43
43
fa="data/genome.fa",
44
- bam=expand("sorted_reads/{sample}.bam", sample=SAMPLES ),
45
- bai=expand("sorted_reads/{sample}.bam.bai", sample=SAMPLES )
44
+ bam=expand("sorted_reads/{sample}.bam", sample=samples ),
45
+ bai=expand("sorted_reads/{sample}.bam.bai", sample=samples )
46
46
output:
47
47
"calls/all.vcf"
48
48
shell:
Original file line number Diff line number Diff line change 1
1
# our samples are pre-configured
2
- SAMPLES = ["A", "B"]
2
+ samples = ["A", "B"]
3
3
4
4
#TASK: Figure out, which your final target is.
5
5
# This is the final target of our workflow, so far.
@@ -44,8 +44,8 @@ rule samtools_index:
44
44
rule bcftools_call:
45
45
input:
46
46
fa="data/genome.fa",
47
- bam=expand("sorted_reads/{sample}.bam", sample=SAMPLES ),
48
- bai=expand("sorted_reads/{sample}.bam.bai", sample=SAMPLES )
47
+ bam=expand("sorted_reads/{sample}.bam", sample=samples ),
48
+ bai=expand("sorted_reads/{sample}.bam.bai", sample=samples )
49
49
output:
50
50
"calls/all.vcf"
51
51
shell:
Original file line number Diff line number Diff line change 1
1
# our samples are pre-configured
2
- SAMPLES = ["A", "B"]
2
+ samples = ["A", "B"]
3
3
4
4
#TASKS: This is an almost working Snakefile!
5
5
# But it is scrambled and contains a bug.
@@ -36,8 +36,8 @@ rule samtools_sort:
36
36
rule bcftools_call:
37
37
input:
38
38
fa="data/genome.fa",
39
- bam=expand("sorted_reads/{sample}.bam", sample=SAMPLES ),
40
- bai=expand("sorted_reads/{sample}.bam.bai", sample=SAMPLES )
39
+ bam=expand("sorted_reads/{sample}.bam", sample=samples ),
40
+ bai=expand("sorted_reads/{sample}.bam.bai", sample=samples )
41
41
output:
42
42
"calls/all.vcf"
43
43
shell:
You can’t perform that action at this time.
0 commit comments