Fix mismatchsize paramter being ignored

Oshlack · Aug 14, 2022 · 572430d · 572430d
1 parent fa0c453
commit 572430d
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/pseudoaligner.rs b/src/pseudoaligner.rs
@@ -517,7 +517,7 @@ pub fn match_read<K: Kmer + Sync + Send>(optional: Option<(Vec<u32>, usize, usiz
 
 	match optional {
         	Some((eq_class,coverage,mismatches,readlen)) => {
-                                    if coverage >= seq.len() && mismatches == 0 &&  eq_class.len() == 1 {
+                                    if coverage >= seq.len() && mismatches <= mismatchsize &&  eq_class.len() == 1 {
                                     //if coverage >= READ_COVERAGE_THRESHOLD && eq_class.len() == 1 {
                                         debug!("{:?}",&seq);
                                       if trim { // trim check