feat: added config file to tutorial folder

setup/solutions/08_Snakefile

@@ -1,25 +1,34 @@
-SAMPLES = ['A', 'B']
+configfile: "config/config.yaml"
 
 rule all:
     input:
         "plots/quals.svg",
         "calls/all.vcf"
 
+
+def get_bwa_map_input_fastqs(wildcards):
+    return config["samples"][wildcards.sample]
+
 rule bwa_map:
     input:
         "data/genome.fa",
-        "data/samples/{sample}.fastq"
+        get_bwa_map_input_fastqs
     output:
-        "mapped_reads/{sample}.bam"
+        temp("mapped_reads/{sample}.bam")
+    params:
+        rg=r"@RG\tID:{sample}\tSM:{sample}"
+    log: "logs/bwa_mem/{sample}.log"
+    threads: 8
     shell:
-        "bwa mem {input} | "
-        "samtools view -Sb - > {output}"
+        "(bwa mem -R '{params.rg}' -t {threads} {input} | "
+        "samtools view -Sb - > {output}) 2> {log}"
+
 
 rule samtools_sort:
     input:
         "mapped_reads/{sample}.bam"
     output:
-        "sorted_reads/{sample}.bam"
+        temp("sorted_reads/{sample}.bam")
     shell:
         "samtools sort -T sorted_reads/{wildcards.sample} "
         "-O bam {input} > {output}"
@@ -33,6 +42,7 @@ rule samtools_index:
     shell:
         "samtools index {input}"
 
+
 rule bcftools_call:
     input:
         fa="data/genome.fa",

setup/tutorial/config/config.yaml

@@ -0,0 +1,6 @@
+samples:
+  A: A.fastq
+  B: B.fastq
+
+bwa_map:
+  rg=r"@RG\tID:{sample}\tSM:{sample}"