Commit 8266b1bb authored by Mathieu Giraud's avatar Mathieu Giraud
Browse files

core/segment.cpp: '-x', marks as UNSEG_TOO_FEW_ZERO if there is no second maximal k-mer

parent b52bafed
......@@ -241,6 +241,12 @@ KmerSegmenter::KmerSegmenter(Sequence seq, Germline *germline)
CountKmerAffectAnalyser ckaa(*(germline->index), sequence);
pair <KmerAffect, KmerAffect> max12 = ckaa.sortLeftRight(ckaa.max12(forbidden));
if (max12.first.isUnknown() || max12.second.isUnknown())
{
because = UNSEG_TOO_FEW_ZERO ;
return ;
}
strand = nb_strand[0] > nb_strand[1] ? -1 : 1 ;
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