Commit d5359b7d authored by Li's avatar Li Committed by Haowen Zhang
Browse files

Replace N by T if reverse complemented to make sure the results is consistent... 

Replace N by T if reverse complemented to make sure the results is consistent with manually rc externally.
parent 59e405c3

src/sequence_batch.h

+3 −0
Original line number Diff line number Diff line
@@ -165,6 +165,9 @@ class SequenceBatch { 
          seed = ((seed << 2) | current_base) & mask;  // forward k-mer

        } else {

          seed = (seed << 2) & mask;  // N->A

          if (effective_range_[2] == -1) {

            seed |= 3ull;  // N->T if the sequence are reverse complemented

          }

        }

      } else {

        seed = (seed << 2) & mask;  // Pad A

CRA Git | Maintained and supported by SUSTech CRA and CCSE