From 1c7937a47fdf740e78f921089c1fc998b1580f3d Mon Sep 17 00:00:00 2001 From: Shifu Chen Date: Mon, 31 Oct 2016 15:11:21 +0800 Subject: [PATCH] update correction condition from Q27 to Q30 --- preprocesser.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/preprocesser.py b/preprocesser.py index 208df53..6581b9a 100755 --- a/preprocesser.py +++ b/preprocesser.py @@ -542,7 +542,7 @@ def run(self): if b1 != b2: # print(TOTAL_READS, o, b1, b2, q1, q2) this_is_corrected = False - if util.qualNum(q1) >= 27 and util.qualNum(q2) <= 16: + if util.qualNum(q1) >= 30 and util.qualNum(q2) <= 14: if b1!='N' and b2!='N': err_mtx[util.complement(b1)][util.complement(b2)] += 1 if not self.options.no_correction: @@ -550,7 +550,7 @@ def run(self): r2[3] = util.changeString(r2[3], -o-1, q1) corrected += 1 this_is_corrected = True - elif util.qualNum(q2) >= 27 and util.qualNum(q1) <= 16: + elif util.qualNum(q2) >= 30 and util.qualNum(q1) <= 14: if b1!='N' and b2!='N': err_mtx[b2][b1] += 1 if not self.options.no_correction: