m

src/lib.rs

@@ -407,8 +407,8 @@ impl Sam {
         // in bam the positions are inclusives [start, end]
         if is_reverse(flag) {
             let len = ref_range.len();
-            ref_range.end = (ref_pos as usize) + 1;
-            ref_range.start = ref_pos as usize + len;
+            ref_range.end = (ref_pos as usize) + 2;
+            ref_range.start = ref_pos as usize + len - 2;
             query_range.end = query_range.end;
         } else {
             let len = ref_range.len();