Commit 0703c4ec authored by Mathieu Giraud's avatar Mathieu Giraud

core/segment.cpp: MAX12, do not segment if the two most frequent kmers are not so frequent

parent 364e819d
......@@ -254,6 +254,9 @@ KmerSegmenter::KmerSegmenter(Sequence seq, Germline *germline)
strand = nb_strand[0] > nb_strand[1] ? -1 : 1 ;
computeSegmentation(strand, max12.first, max12.second);
if (!detected)
because = UNSEG_TOO_FEW_ZERO ;
// The pseudo-germline should never take precedence over the regular germlines
score = 1 ;
}
......
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