fix: deal with zero length blocks getting added to segment

neherlab · Mar 18, 2022 · f6c39c3 · f6c39c3
1 parent 36f72d5
commit f6c39c3
Showing 1 changed file with 3 additions and 4 deletions.
diff --git a/src/block.jl b/src/block.jl
@@ -146,20 +146,19 @@ function partition(alignment; minblock=500)
 
     # ----------------------------
     # internal operators
-
     function finalize_block!()
         length(segment) == 0 && @goto advance
 
         push!(block, (
             range   = (
-                qry = Pos(qry.start,qry.stop-1), 
-                ref = Pos(ref.start,ref.stop-1)
+                qry = (qry.stop-1 ≥ qry.start) ? Pos(qry.start,qry.stop-1) : nothing,
+                ref = (ref.stop-1 ≥ ref.start) ? Pos(ref.start,ref.stop-1) : nothing,
             ),
             segment = segment
         ))
 
         segment = PosPair[]
-        
+
         @label advance
         advance!(qry)
         advance!(ref)