Extract variant annotation using SnpSift

Author: ahmig

workflow/rules/vaf.smk

@@ -137,6 +137,26 @@ rule tsv_to_vcf:
         "../scripts/tsv_to_vcf.py"
 
 
+rule extract_vcf_fields:
+    threads: 1
+    conda: "../envs/snpeff.yaml"
+    params:
+        extract_columns = [
+            "CHROM", "REF", "POS", "ALT", "DP",
+            '"GEN[*].ALT_DP"', '"GEN[*].ALT_RV"', '"GEN[*].ALT_FREQ"',
+            '"GEN[*].ALT_QUAL"', '"ANN[*].GENE"', '"ANN[*].HGVS_P"'
+        ],
+        sep = ","
+    input:
+        vcf = OUTDIR/f"{OUTPUT_NAME}.vcf"
+    output:
+        tsv = OUTDIR/f"{OUTPUT_NAME}.vcf_fields.tsv"
+    log:
+        LOGDIR / "tsv_to_vcf" / "log.txt"
+    shell:
+        'SnpSift extractFields -s {params.sep:q} {input.vcf:q} {params.extract_columns} >{output.tsv:q} 2>{log:q}'
+
+
 rule variants_effect:
     threads: 1
     shadow: "minimal"