We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent e65486b commit 3295c35Copy full SHA for 3295c35
1 file changed
src/filter.rs
@@ -181,5 +181,5 @@ fn phred_mean(
181
.map(|x| { 10.0f64.powf((x - phred) as f64 / -10.0) })
182
.sum::<f64>() / qual.len() as f64;
183
184
- (-10.0f64 * ave_error.log10()) as u8
+ (-10.0f64 * ave_error.log10()).round() as u8
185
}
0 commit comments