fix quality score calculation

wdecoster · Apr 4, 2023 · dfb31cf · dfb31cf
1 parent 4e7c162
commit dfb31cf
Show file tree

Hide file tree

Showing 3 changed files with 4 additions and 3 deletions.
diff --git a/Cargo.lock b/Cargo.lock
diff --git a/Cargo.toml b/Cargo.toml
@@ -1,6 +1,6 @@
 [package]
 name = "chopper"
-version = "0.4.0"
+version = "0.5.0"
 authors = ["wdecoster <decosterwouter@gmail.com>"]
 edition = "2018"
 

diff --git a/src/main.rs b/src/main.rs
@@ -79,7 +79,8 @@ where
                     let read_len = record.seq().len();
                     // If a read is shorter than what is to be cropped the read is dropped entirely (filtered out)
                     if args.headcrop + args.tailcrop < read_len {
-                        let average_quality = ave_qual(record.qual());
+                        let average_quality =
+                            ave_qual(&record.qual().iter().map(|i| i - 33).collect::<Vec<u8>>());
                         if average_quality >= args.minqual
                             && average_quality <= args.maxqual
                             && read_len >= args.minlength