m

src/lib.rs

@@ -407,14 +407,13 @@ 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);
+            ref_range.end = (ref_pos as usize) + 1;
             ref_range.start = ref_pos as usize + len - 1;
             query_range.end = query_range.end;
         } else {
             let len = ref_range.len(); 
             ref_range.start = ref_pos as usize;
             ref_range.end = (ref_pos as usize + len) - 1;
-            query_range.end = query_range.end ;
         }
 
         let ref_sequence = fa.get_sequence(&ref_name, ref_range.start.clone(), ref_range.end.clone());