Commit a9b4ecc4 authored by Mathieu Giraud's avatar Mathieu Giraud
Browse files

core/segment.cpp: if a sequence is AMBIGUOUS, do no try the remaining germlines

Getting UNSEG_AMBIGUOUS is not a normal situation. When it happens, we want
to exit and to report this sequence to the user, and not to segment by chance
the sequence on another germline.
parent a8d2c873
......@@ -233,6 +233,11 @@ KmerSegmenter::KmerSegmenter(Sequence seq, MultiGermline *multigermline)
system = germline->code;
return ;
// If a read is AMBIGUOUS, do not test other germlines
if (because == UNSEG_AMBIGUOUS)
return ;
} // end for (Germlines)
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