context size correct

Schaudge · Schaudge · commit 65ac218eaad3 · 2022-07-10T22:55:39.000+08:00
diff --git a/src/lib/position/pileup_position.rs b/src/lib/position/pileup_position.rs
@@ -120,9 +120,9 @@ impl PileupPosition {
                 }
                 bam::pileup::Indel::Del(len) => {
                     let del_pos = alignment.qpos().unwrap();
-                    if len > 5 && del_pos >= context_size && del_pos + context_size <= record.seq_len() {
+                    if len > 5 && del_pos >= context_size && del_pos + context_size < record.seq_len() {
                         let mut serial_seq: String = String::from("");
-                        for ip in del_pos - context_size - 1..del_pos + context_size + 1 {
+                        for ip in del_pos - context_size + 1..del_pos + context_size + 1 {
                             serial_seq.push(record.seq()[ip] as char);
                         }
                         *self.del_seq_map.entry(serial_seq).or_insert(0) += 1;

Original file line number	Diff line number	Diff line change
`@@ -120,9 +120,9 @@ impl PileupPosition {`
`120`	`120`	`}`
`121`	`121`	`bam::pileup::Indel::Del(len) => {`
`122`	`122`	`let del_pos = alignment.qpos().unwrap();`
`123`		`- if len > 5 && del_pos >= context_size && del_pos + context_size <= record.seq_len() {`
	`123`	`+ if len > 5 && del_pos >= context_size && del_pos + context_size < record.seq_len() {`
`124`	`124`	`let mut serial_seq: String = String::from("");`
`125`		`- for ip in del_pos - context_size - 1..del_pos + context_size + 1 {`
	`125`	`+ for ip in del_pos - context_size + 1..del_pos + context_size + 1 {`
`126`	`126`	`serial_seq.push(record.seq()[ip] as char);`
`127`	`127`	`}`
`128`	`128`	`*self.del_seq_map.entry(serial_seq).or_insert(0) += 1;`