fix bug with N base

Xinglab · Feb 1, 2019 · 4f49940 · 4f49940
1 parent c12c7ae
commit 4f49940
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 4 deletions.
diff --git a/abPOA b/abPOA
diff --git a/src/seq.c b/src/seq.c
@@ -78,9 +78,9 @@ uint8_t *get_bseq(char *seq, int seq_len) {
     int i;
     uint8_t *bseq = (uint8_t*)_err_malloc(seq_len * sizeof(uint8_t));
     for (i = 0; i < seq_len; ++i) {
-        // TODO N(ambiguous base)
-        bseq[i] = hash_nt4_table[nst_nt4_table[(int)seq[i]]];
-        // bseq[i] = nst_nt4_table[(int)seq[i]];
+        // N(ambiguous base)
+        // bseq[i] = hash_nt4_table[nst_nt4_table[(int)seq[i]]];
+        bseq[i] = nst_nt4_table[(int)seq[i]];
     }
     return bseq;
 }