diff --git a/python/obitools3/commands/uniq.pyx b/python/obitools3/commands/uniq.pyx index 36ff184..6ec3e7e 100644 --- a/python/obitools3/commands/uniq.pyx +++ b/python/obitools3/commands/uniq.pyx @@ -23,7 +23,9 @@ from cpython.exc cimport PyErr_CheckSignals __title__="Group sequence records together" - + + +REVERSE_QUALITY_COLUMN_NAME = b"REVERSE_QUALITY" # TODO from ngsfilter, move to C def addOptions(parser): @@ -491,9 +493,11 @@ cdef uniq_sequences(View_NUC_SEQS view, View_NUC_SEQS o_view, ProgressBar pb, li o_idx += 1 - # Deletes quality column if there is one because the matching between sequence and quality will be broken (quality set to NA when sequence not) + # Deletes quality columns if there is one because the matching between sequence and quality will be broken (quality set to NA when sequence not) if QUALITY_COLUMN in view: o_view.delete_column(QUALITY_COLUMN) + if REVERSE_QUALITY_COLUMN_NAME in view: + o_view.delete_column(REVERSE_QUALITY_COLUMN_NAME) if taxonomy is not None: print("") # TODO because in the middle of progress bar. Better solution?