Function rustybam::trim_overlap::trim_overlapping_pafs
source · [−]pub fn trim_overlapping_pafs(
left: &mut PafRecord,
right: &mut PafRecord,
match_score: i32,
diff_score: i32,
indel_score: i32
)
Expand description
Example
use rustybam::trim_overlap::*;
use rustybam::paf::*;
let mut left = PafRecord::new("Q 10 0 10 + T 20 0 10 3 9 60 cg:Z:7=1X2=").unwrap();
left.aligned_pairs();
let mut right = PafRecord::new("Q 10 5 10 - T 20 10 15 3 9 60 cg:Z:3=1X1=").unwrap();
right.aligned_pairs();
trim_overlapping_pafs(&mut left, &mut right, 1 ,1 ,1);
assert_eq!(left.cigar.to_string(), "7=");
assert_eq!(right.cigar.to_string(), "3=");