m

src/lib.rs

@@ -385,12 +385,13 @@ impl Sam {
         
         if is_reverse(flag) {
             let len = ref_range.len();
-            ref_range.end = ref_pos as usize;
+            // +1 because the range exclude the last [ref_range.start, ref_range.end[
+            ref_range.end = (ref_pos as usize) + 1;
             ref_range.start = ref_pos as usize + len ;
         } else {
             let len = ref_range.len(); 
             ref_range.start = ref_pos as usize;
-            ref_range.end = ref_pos as usize + len;
+            ref_range.end = (ref_pos as usize + len) + 1;
         }
 
         let ref_sequence = fa.get_sequence(&ref_name, ref_range.start.clone(), ref_range.end.clone());