fix the ultra stupid bug for - strand in cdna

hasindu2008 · Feb 28, 2024 · e741017 · e741017
1 parent 4eaf31d
commit e741017
Showing 1 changed file with 10 additions and 0 deletions.
diff --git a/src/genread.c b/src/genread.c
@@ -298,6 +298,16 @@ static char *gen_read_rna(core_t *core, char **ref_id, int32_t *ref_len, int32_t
 
     }
 
+    if(cdna){
+        if(*c == '-'){
+            char *r = reverse_complement(seq);
+            r[*rlen] = '\0';
+            free(seq);
+            seq = r;
+        }
+
+    }
+
     return seq;
 
 }