|
|
@@ -391,7 +391,7 @@ impl Sam {
|
|
|
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 - 1;
|
|
|
+ ref_range.start = ref_pos as usize + len;
|
|
|
} else {
|
|
|
let len = ref_range.len();
|
|
|
ref_range.start = ref_pos as usize;
|
|
|
@@ -677,7 +677,7 @@ fn matched_range (cigar: &str, flag: &i32, matched_seq: &mut str) -> (Range<usiz
|
|
|
if is_reverse(*flag) {
|
|
|
let query_len = range_query.len();
|
|
|
range_query.start = query_pos - range_query.end;
|
|
|
- range_query.end = range_query.start + query_len;
|
|
|
+ range_query.end = range_query.start + query_len - 1;
|
|
|
}
|
|
|
|
|
|
assert_eq!(range_ref.len(), ref_cigar_string.len());
|
