|
|
@@ -722,13 +722,13 @@ fn matched_range (cigar: &str, flag: &i32, matched_seq: &mut str) -> (Range<usiz
|
|
|
let mut current_add = num_acc.parse::<usize>().unwrap();
|
|
|
num_acc = "".to_string();
|
|
|
|
|
|
- // if n_head_to_rm > 0 && curr_op == 1 {
|
|
|
- // current_add = current_add - n_head_to_rm;
|
|
|
- // }
|
|
|
- // if n_trail_to_rm > 0 && curr_op == n_op {
|
|
|
- // current_add = current_add - n_trail_to_rm;
|
|
|
- // }
|
|
|
- // curr_op += 1;
|
|
|
+ if n_head_to_rm > 0 && curr_op == 1 {
|
|
|
+ current_add = current_add - n_head_to_rm;
|
|
|
+ }
|
|
|
+ if n_trail_to_rm > 0 && curr_op == n_op {
|
|
|
+ current_add = current_add - n_trail_to_rm;
|
|
|
+ }
|
|
|
+ curr_op += 1;
|
|
|
|
|
|
if f == "M" {
|
|
|
if first_m {
|
