Merge pull request #91 from jts/fix_ref_fasta_index

Add reference indexing, support different primer name prefixes

Author: rdeborja

workflow/rules/common.smk

@@ -245,19 +245,31 @@ rule make_amplicon_bed:
     params:
         script="primers_to_amplicons.py",
         offset=get_primer_offset,
-        bed_type_opt=get_primer_bed_type_opt
+        bed_type_opt=get_primer_bed_type_opt,
+        primer_prefix=get_primer_prefix
     shell:
-        "{params.script} --primers {input.primers} --offset {params.offset} --bed_type {params.bed_type_opt} --output {output}"
+        "{params.script} --primers {input.primers} --offset {params.offset} --bed_type {params.bed_type_opt} --output {output} --primer_prefix {params.primer_prefix}"
 
 rule make_amplicon_full_bed:
     input:
         primers=get_primer_bed
     output:
         "bed/amplicon_full.bed"
     params:
-        script="primers_to_amplicons.py"
+        script="primers_to_amplicons.py",
+        primer_prefix=get_primer_prefix
+    shell:
+        "{params.script} --primers {input.primers} --bed_type full --output {output} --primer_prefix {params.primer_prefix}"
+
+rule index_reference_genome:
+    input:
+        get_reference_genome
+    output:
+        expand(config["reference_genome"] + ".fai")
+    params:
+        exec="samtools faidx"
     shell:
-        "{params.script} --primers {input.primers} --bed_type full --output {output}"
+        "{params.exec} {input}"
 
 # make a bed file for the entire reference genome as a single record
 # from: https://bioinformatics.stackexchange.com/questions/91/how-to-convert-fasta-to-bed

workflow/rules/defaults.smk

@@ -42,6 +42,10 @@ def get_primer_offset(wildcards):
 def get_primer_bed_type_opt(wildcards):
     return config.get("bed_type", "unique_amplicons")
 
+# get the primer name prefix from the config.yaml file
+def get_primer_prefix(wildcards):
+    return config.get("primer_prefix", "nCoV-2019")
+
 def get_snp_tree_flag(wildcards=None):
     return config.get("build_snp_tree", True)