diff --git a/algo/core/kmerstore.h b/algo/core/kmerstore.h index e6d386ad8ff29a5b3c4f616797ca73a802186987..514f4ae135c70a741b52323a37488b5cc85c50f3 100644 --- a/algo/core/kmerstore.h +++ b/algo/core/kmerstore.h @@ -333,7 +333,7 @@ string IKmerStore::getLabel(T kmer) const { for (typename list< pair >::const_iterator it = labels.begin(); it != labels.end(); ++it) if (it->first == kmer) return it->second ; - return "" ; + return "?" ; } // .getResults() diff --git a/algo/core/segment.cpp b/algo/core/segment.cpp index 69cc66ac3b020473f23948267fb08a6f09747dc3..7e6899fa427a3193a88255ff0f8add8b37a7ac9f 100644 --- a/algo/core/segment.cpp +++ b/algo/core/segment.cpp @@ -642,9 +642,10 @@ FineSegmenter::FineSegmenter(Sequence seq, Germline *germline, Cost segment_c) if (!germline->rep_5.size() || !germline->rep_3.size()) { - // We check whether this sequence is segmented with MAX12 (with default e-value parameters) + // We check whether this sequence is segmented with MAX12 or MAX1U (with default e-value parameters) KmerSegmenter *kseg = new KmerSegmenter(seq, germline, THRESHOLD_NB_EXPECTED, 1); - if (kseg->isSegmented() && (!strcmp(germline->code.c_str(), PSEUDO_GERMLINE_MAX12))) + if (kseg->isSegmented() && ((!strcmp(germline->code.c_str(), PSEUDO_GERMLINE_MAX12)) + || !strcmp(germline->code.c_str(), PSEUDO_GERMLINE_MAX1U))) { reversed = kseg->isReverse();