Search both the forward and revcomp of the ref and alt (#1)

fulcrumgenomics · Oct 13, 2021 · 0514930 · 0514930
1 parent 8bc0046
commit 0514930
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/main.rs b/src/main.rs
@@ -366,8 +366,12 @@ fn main() -> Result<()> {
                 .into_par_iter()
                 .map(|(r1, r2)| {
                     let ref_match = r1.seq.find(ref_forward_seq).is_some()
+                        || r1.seq.find(&ref_revcomp_seq).is_some()
+                        || r2.seq.find(ref_forward_seq).is_some()
                         || r2.seq.find(&ref_revcomp_seq).is_some();
                     let alt_match = r1.seq.find(alt_forward_seq).is_some()
+                        || r1.seq.find(&alt_revcomp_seq).is_some()
+                        || r2.seq.find(alt_forward_seq).is_some()
                         || r2.seq.find(&alt_revcomp_seq).is_some();
                     let m = match (ref_match, alt_match) {
                         (true, true) => Matches::Both,