fix for snps

brentp · May 7, 2019 · 913d1c1 · 913d1c1
1 parent 89327be
commit 913d1c1
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/duphold.nim b/src/duphold.nim
@@ -522,7 +522,7 @@ proc annotate*(snps:snpset, variant:Variant, sample_i:int) =
   var ci = innerCI(variant)
   var i = lowerBound(snps.starts, ci.left.int32, system.cmp)
   var n = 0
-  while snps.starts[i] < ci.right:
+  while i < snps.starts.len and snps.starts[i] < ci.right:
     dhgt[5 * sample_i + snps.nalts[i]] += 1
     n += 1
     i += 1