Commit 888d0c69 authored by Mathieu Giraud's avatar Mathieu Giraud

core/segment.cpp, MAX12: use sortLeftRight() to sort the two most frequent KmerAffects

parent 10b0e66e
...@@ -243,7 +243,7 @@ KmerSegmenter::KmerSegmenter(Sequence seq, Germline *germline) ...@@ -243,7 +243,7 @@ KmerSegmenter::KmerSegmenter(Sequence seq, Germline *germline)
forbidden.insert(KmerAffect::getUnknown()); forbidden.insert(KmerAffect::getUnknown());
CountKmerAffectAnalyser ckaa(*(germline->index), sequence); CountKmerAffectAnalyser ckaa(*(germline->index), sequence);
pair <KmerAffect, KmerAffect> max12 = ckaa.max12(forbidden); pair <KmerAffect, KmerAffect> max12 = ckaa.sortLeftRight(ckaa.max12(forbidden));
strand = nb_strand[0] > nb_strand[1] ? -1 : 1 ; strand = nb_strand[0] > nb_strand[1] ? -1 : 1 ;
computeSegmentation(strand, max12.first, max12.second); computeSegmentation(strand, max12.first, max12.second);
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment