Create independent Snakemake language

lib/linguist/languages.yml

@@ -5389,7 +5389,6 @@ Python:
   - ".pyt"
   - ".pyw"
   - ".rpy"
-  - ".smk"
   - ".spec"
   - ".tac"
   - ".wsgi"
@@ -5399,7 +5398,6 @@ Python:
   - DEPS
   - SConscript
   - SConstruct
-  - Snakefile
   - wscript
   interpreters:
   - python
@@ -6509,6 +6507,22 @@ Smithy:
   extensions:
   - ".smithy"
   language_id: 1027892786
+Snakemake:
+  type: programming
+  group: Python
+  tm_scope: source.python
+  ace_mode: python
+  codemirror_mode: python
+  codemirror_mime_type: text/x-python
+  color: "#419179"
+  extensions:
+  - ".smk"
+  - ".snakefile"
+  filenames:
+  - Snakefile
+  aliases:
+  - snakefile
+  language_id: 151241392
 Solidity:
   type: programming
   color: "#AA6746"

samples/Snakemake/filenames/Snakefile

@@ -0,0 +1,49 @@
+#example from https://github.com/snakemake/snakemake/edit/main/examples/hello-world/Snakefile
+configfile: "config.yaml"
+
+
+rule all:
+    input:
+        expand(
+            "plots/{country}.hist.pdf",
+            country=config["countries"]
+        )
+
+
+rule select_by_country:
+    input:
+        "data/worldcitiespop.csv"
+    output:
+        "by-country/{country}.csv"
+    conda:
+        "envs/xsv.yaml"
+    shell:
+        "xsv search -s Country '{wildcards.country}' "
+        "{input} > {output}"
+
+
+rule plot_histogram:
+    input:
+        "by-country/{country}.csv"
+    output:
+        "plots/{country}.hist.svg"
+    container:
+        "docker://faizanbashir/python-datascience:3.6"
+    script:
+        "scripts/plot-hist.py"
+
+
+rule convert_to_pdf:
+    input:
+        "{prefix}.svg"
+    output:
+        "{prefix}.pdf"
+    wrapper:
+        "0.47.0/utils/cairosvg"
+
+
+rule download_data:
+    output:
+        "data/worldcitiespop.csv"
+    shell:
+        "curl -L https://burntsushi.net/stuff/worldcitiespop.csv > {output}"

samples/Snakemake/template.snakefile

@@ -0,0 +1,70 @@
+configfile: "config.yaml"
+
+
+rule all:
+    input:
+        "plots/quals.svg"
+
+
+def get_bwa_map_input_fastqs(wildcards):
+    return config["samples"][wildcards.sample]
+
+
+rule bwa_map:
+    input:
+        "data/genome.fa",
+        get_bwa_map_input_fastqs
+    output:
+        temp("mapped_reads/{sample}.bam")
+    params:
+        rg=r"@RG\tID:{sample}\tSM:{sample}"
+    log:
+        "logs/bwa_mem/{sample}.log"
+    threads: 8
+    shell:
+        "(bwa mem -R '{params.rg}' -t {threads} {input} | "
+        "samtools view -Sb - > {output}) 2> {log}"
+
+
+rule samtools_sort:
+    input:
+        "mapped_reads/{sample}.bam"
+    output:
+        protected("sorted_reads/{sample}.bam")
+    shell:
+        "samtools sort -T sorted_reads/{wildcards.sample} "
+        "-O bam {input} > {output}"
+
+
+rule samtools_index:
+    input:
+        "sorted_reads/{sample}.bam"
+    output:
+        "sorted_reads/{sample}.bam.bai"
+    shell:
+        "samtools index {input}"
+
+
+rule bcftools_call:
+    input:
+        fa="data/genome.fa",
+        bam=expand("sorted_reads/{sample}.bam", sample=config["samples"]),
+        bai=expand("sorted_reads/{sample}.bam.bai", sample=config["samples"])
+    output:
+        "calls/all.vcf"
+    params:
+        rate=config["prior_mutation_rate"]
+    log:
+        "logs/bcftools_call/all.log"
+    shell:
+        "(bcftools mpileup -f {input.fa} {input.bam} | "
+        "bcftools call -mv -P {params.rate} - > {output}) 2> {log}"
+
+
+rule plot_quals:
+    input:
+        "calls/all.vcf"
+    output:
+        "plots/quals.svg"
+    script:
+        "scripts/plot-quals.py"

vendor/README.md

@@ -505,6 +505,7 @@ This is a list of grammars that Linguist selects to provide syntax highlighting
 - **Smalltalk:** [tomas-stefano/smalltalk-tmbundle](https://github.com/tomas-stefano/smalltalk-tmbundle)
 - **Smarty:** [textmate/php-smarty.tmbundle](https://github.com/textmate/php-smarty.tmbundle)
 - **Smithy:** [awslabs/smithy-vscode](https://github.com/awslabs/smithy-vscode)
+- **Snakemake:** [MagicStack/MagicPython](https://github.com/MagicStack/MagicPython)
 - **Solidity:** [davidhq/SublimeEthereum](https://github.com/davidhq/SublimeEthereum)
 - **Soong:** [flimberger/android-system-tools](https://github.com/flimberger/android-system-tools)
 - **SourcePawn:** [Dreae/sourcepawn-vscode](https://github.com/Dreae/sourcepawn-vscode)