Ensure that query range used is limited to contig start/stop. (#31)

* Ensure that query range used is limited to contig start/stop.
fulcrumgenomics · Oct 18, 2023 · 70519fd · 70519fd
1 parent 204a5ee
commit 70519fd
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 4 deletions.
diff --git a/docs/tools/index.md b/docs/tools/index.md
@@ -4,7 +4,7 @@ title: fgsv tools
 
 # fgsv tools
 
-The following tools are available in fgsv version 0.1.0-fd73568.
+The following tools are available in fgsv version 0.1.0-dcaa891.
 ## All tools
 
 All tools.

diff --git a/src/main/scala/com/fulcrumgenomics/sv/tools/AggregateSvPileup.scala b/src/main/scala/com/fulcrumgenomics/sv/tools/AggregateSvPileup.scala
@@ -351,8 +351,8 @@ case class AggregatedBreakpointPileup(id: String,
           source       = src,
           contig       = contig,
           breakends    = pileups,
-          queryMin     = minPos - flank,
-          queryMax     = maxPos + flank,
+          queryMin     = Math.max(1, minPos - flank),
+          queryMax     = Math.min(maxPos + flank, src.dict(contig).length),
           minTotal     = minTotal,
           minFrequency = minFrequency,
         )
@@ -461,7 +461,6 @@ case class AggregatedBreakpointPileup(id: String,
   private def overlapsTarget(contig: String, start: Int, end: Int, targets: OverlapDetector[BEDFeature]): Boolean = {
     targets.overlapsAny(new Interval(contig, start, end))
   }
-
 }