filter.rs: fix bugs in average quality score calculation (f64.round())

sharkLoc · Dec 20, 2023 · 3295c35 · 3295c35
1 parent e65486b
commit 3295c35
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/filter.rs b/src/filter.rs
@@ -181,5 +181,5 @@ fn phred_mean(
         .map(|x| { 10.0f64.powf((x - phred) as f64 / -10.0) })
         .sum::<f64>() / qual.len() as f64;
 
-    (-10.0f64 * ave_error.log10()) as u8
+    (-10.0f64 * ave_error.log10()).round() as u8
 }