From 4f4994099157799cabc863fc8d056d4d3fedc4c6 Mon Sep 17 00:00:00 2001 From: Yan Gao Date: Thu, 31 Jan 2019 22:52:43 -0500 Subject: [PATCH] fix bug with N base --- abPOA | 2 +- src/seq.c | 6 +++--- 2 files changed, 4 insertions(+), 4 deletions(-) diff --git a/abPOA b/abPOA index 5c1b560..3cac111 160000 --- a/abPOA +++ b/abPOA @@ -1 +1 @@ -Subproject commit 5c1b5600f6b0e55eab2a1a1446e75804f8545345 +Subproject commit 3cac11133fc8686ee4671d697caee840faedb7b7 diff --git a/src/seq.c b/src/seq.c index 79ed1e5..6708dd5 100644 --- 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; }