m

src/lib.rs

@@ -390,8 +390,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) - 1;
+            ref_range.start = ref_pos as usize + len - 1;
         } else {
             let len = ref_range.len(); 
             ref_range.start = ref_pos as usize;