Update to version 1.1.1

PacificBiosciences · Jul 29, 2024 · 7bff22d · 7bff22d
1 parent 4b7576d
commit 7bff22d
Show file tree

Hide file tree

Showing 4 changed files with 6 additions and 4 deletions.
diff --git a/Cargo.lock b/Cargo.lock
diff --git a/Cargo.toml b/Cargo.toml
@@ -1,6 +1,6 @@
 [package]
 name = "trgt"
-version = "1.1.0"
+version = "1.1.1"
 edition = "2021"
 build = "build.rs"
 

diff --git a/README.md b/README.md
@@ -123,6 +123,8 @@ tandem repeats at genome scale. 2024](https://www.nature.com/articles/s41587-023
   - Added subsampling of regions with ultra-high coverage (`>MAX_DEPTH * 3`, by default 750); implemented via reservoir sampling.
   - Fixed a cluster genotyper bug that occurred when only a single read covered a locus.
   - Added new logic for filtering non-HiFi reads: remove up to 3% of lower quality reads that do not match the expected repeat sequence.
+- 1.1.1
+  - Hotfix: Read filtering logic no longer removes reads without RQ tags.
 
 ### DISCLAIMER
 

diff --git a/src/trgt/workflows/tr.rs b/src/trgt/workflows/tr.rs
@@ -297,7 +297,7 @@ fn extract_reads(
                     continue;
                 }
 
-                if get_rq_tag(&record).unwrap_or(0.0) < min_read_qual {
+                if get_rq_tag(&record).unwrap_or(1.0) < min_read_qual {
                     n_filt += 1;
                     continue;
                 }
@@ -322,7 +322,7 @@ fn extract_reads(
                         continue;
                     }
 
-                    if get_rq_tag(&record).unwrap_or(0.0) < min_read_qual {
+                    if get_rq_tag(&record).unwrap_or(1.0) < min_read_qual {
                         n_filt += 1;
                         continue;
                     }