Merge a56bb6d into 9789a2d

bihealth · May 16, 2024 · 7a27eb9 · 7a27eb9
2 parents 9789a2d + a56bb6d
commit 7a27eb9
Show file tree

Hide file tree

Showing 2 changed files with 15 additions and 5 deletions.
diff --git a/snappy_wrappers/wrappers/mehari/annotate_seqvars/wrapper.py b/snappy_wrappers/wrappers/mehari/annotate_seqvars/wrapper.py
@@ -58,21 +58,31 @@
 
 # Run actual tools --------------------------------------------------------------------------------
 
-# Extract around BED file, if given.
+# Extract around BED file, if given.  Otherwise, "just" normalize.
 if [[ -n "{export_config[path_exon_bed]}" ]] && [[ "{export_config[path_exon_bed]}" != "None" ]]; then
     set -e
     bcftools view \
         -R {export_config[path_exon_bed]} \
         {snakemake.input.vcf} \
+    | bcftools norm \
+        -m -any \
+        --force \
+        --fasta-ref {snakemake.config[static_data_config][reference][path]} \
     | bcftools sort -T $TMPDIR \
-    | bcftools norm -d all \
     | bgzip -c \
     > $TMPDIR/tmp.vcf.gz
     tabix -f $TMPDIR/tmp.vcf.gz
 else
     set -e
-    ln -sr {snakemake.input.vcf} $TMPDIR/tmp.vcf.gz
-    ln -sr {snakemake.input.vcf}.tbi $TMPDIR/tmp.vcf.gz.tbi
+    bcftools norm \
+        -m -any \
+        --force \
+        --fasta-ref {snakemake.config[static_data_config][reference][path]} \
+        {snakemake.input.vcf} \
+    | bcftools sort -T $TMPDIR \
+    | bgzip -c \
+    > $TMPDIR/tmp.vcf.gz
+    tabix -f $TMPDIR/tmp.vcf.gz
 fi
 
 # Perform Mehari sequence variant annotation.

diff --git a/snappy_wrappers/wrappers/mehari/environment.yaml b/snappy_wrappers/wrappers/mehari/environment.yaml
@@ -5,7 +5,7 @@ dependencies:
   - bcftools==1.19
   - htslib==1.19.1
   - samtools==1.19.2
-  - mehari==0.21.1
+  - mehari==0.25.5
   - jq==1.7.1
   - vcfpy==0.13.8
   - pysam==0.22.0