compare to negative strand for dna_start_t_withstop

unipept · Sep 22, 2021 · d5c451a · d5c451a
1 parent 692c10c
commit d5c451a
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/viterbi.rs b/src/viterbi.rs
@@ -616,7 +616,7 @@ fn build_genes(
             dna.push(seq[t]);
             dna_start_t_withstop = t + 1;
             dna_start_t = t + 1;
-            if vpath[t] == hmm::State::M1 || vpath[t] == hmm::State::M4r {
+            if vpath[t] == hmm::State::M1r || vpath[t] == hmm::State::M4r {
                 if t > 2 {
                     dna_start_t_withstop = t - 2;
                 }