fix(call): vcf output with --respect-ref

davidlougheed · Jan 30, 2024 · 8f39a3d · 8f39a3d
1 parent ebe9e28
commit 8f39a3d
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/strkit/VERSION b/strkit/VERSION
@@ -1 +1 @@
-0.13.0
+0.13.1
diff --git a/strkit/call/output/vcf.py b/strkit/call/output/vcf.py
@@ -126,7 +126,7 @@ def _write_contig_vrs():
 
             # seq_alleles = (ref_start_anchor + ref_seq, *(seq_alts or (".",))) if call is not None else (".",)
 
-            start = result["start_adj"] - len(ref_start_anchor)
+            start = result.get("start_adj", result["start"]) - len(ref_start_anchor)
             vr: pysam.VariantRecord = vf.new_record(
                 contig=contig,
                 start=start,