From 3acd6facb8508acb15f588c25beb07a00d55a9ac Mon Sep 17 00:00:00 2001 From: BenLangmead Date: Sat, 15 Apr 2017 23:12:41 -0400 Subject: [PATCH] bugfix to avoid calling rdid on invalid read --- aligner.h | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/aligner.h b/aligner.h index 534d4be..2e5ac30 100644 --- a/aligner.h +++ b/aligner.h @@ -279,7 +279,7 @@ class MixedMultiAligner { // Get a new read pair ret = ps->nextReadPair(); saw_last_read = ret.second; - if(ps->rdid() < qUpto_ && ret.first) { + if(ret.first && ps->rdid() < qUpto_) { if(ps->paired()) { // Read currently in buffer is paired-end (*alignersPE_)[0]->setQuery(ps);