From cc97887069d26dd79ab1a28e3040ee7beb0a6f6d Mon Sep 17 00:00:00 2001
From: Tom Clegg <tom@curii.com>
Date: Thu, 30 Dec 2021 11:24:44 -0500
Subject: [PATCH] Fix coordinates in hgvs annotations.

refs #18438

Arvados-DCO-1.1-Signed-off-by: Tom Clegg <tom@curii.com>
---
 slicenumpy.go | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/slicenumpy.go b/slicenumpy.go
index 8696aa53f3..82cfff75ce 100644
--- a/slicenumpy.go
+++ b/slicenumpy.go
@@ -454,8 +454,10 @@ func (cmd *sliceNumpy) RunCommand(prog string, args []string, stdin io.Reader, s
 						continue
 					}
 					diffs, _ := hgvs.Diff(reftilestr, strings.ToUpper(string(tv.Sequence)), 0)
+					for i := range diffs {
+						diffs[i].Position += rt.pos
+					}
 					for _, diff := range diffs {
-						diff.Position += rt.pos
 						fmt.Fprintf(annow, "%d,%d,%d,%s:g.%s,%s,%d,%s,%s,%s\n", tag, outcol, v, rt.seqname, diff.String(), rt.seqname, diff.Position, diff.Ref, diff.New, diff.Left)
 					}
 					if *hgvsChunked {
-- 
2.30.2