Sanity check sequence length after primer stripping

neufeld · Nov 6, 2014 · f1c79dd · f1c79dd
1 parent 447a4d3
commit f1c79dd
Showing 1 changed file with 5 additions and 0 deletions.
diff --git a/assembler.c b/assembler.c
@@ -70,6 +70,11 @@ static bool align(
 	double rquality = 0;
 	ptrdiff_t len;
 
+	if (assembler->minoverlap + result->forward_offset >= result->forward_length || assembler->minoverlap + result->reverse_offset >= result->reverse_length) {
+		LOG(PANDA_DEBUG_BUILD, PANDA_CODE_NEGATIVE_SEQUENCE_LENGTH);
+		return false;
+	}
+
 	if (assembler->maxoverlap == 0) {
 		maxoverlap = result->forward_length < result->reverse_length ? result->forward_length : result->reverse_length;
 	} else if (maxoverlap > assembler->maxoverlap) {