|
|
@@ -875,4 +875,19 @@ mod tests {
|
|
|
let n_ins = reads.iter().filter(|(_, b)| *b == b'I').count();
|
|
|
assert_eq!(n_ins, 5);
|
|
|
}
|
|
|
+
|
|
|
+ #[test]
|
|
|
+ fn del() {
|
|
|
+ init();
|
|
|
+ let id = "SALICETTO";
|
|
|
+ let result_dir = "/data/longreads_basic_pipe";
|
|
|
+ let diag_bam_path = format!("{result_dir}/{id}/diag/{id}_diag_hs1.bam");
|
|
|
+ let ins_pos = ("chr7", 511_886);
|
|
|
+ let mut bam = IndexedReader::from_path(diag_bam_path).unwrap();
|
|
|
+
|
|
|
+ let reads = qnames_at_base(&mut bam, ins_pos.0, ins_pos.1, false).unwrap();
|
|
|
+
|
|
|
+ let n_del = reads.iter().filter(|(_, b)| *b == b'D').count();
|
|
|
+ assert_eq!(n_del, 1);
|
|
|
+ }
|
|
|
}
|
