1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
|
Author: Steffen Moeller <moeller@debian.org>
Description: Avoid division by 0 error, along idea behind PR#96
Index: ngmlr/src/AlignmentBuffer.cpp
===================================================================
--- ngmlr.orig/src/AlignmentBuffer.cpp
+++ ngmlr/src/AlignmentBuffer.cpp
@@ -1881,6 +1881,7 @@ int AlignmentBuffer::computeMappingQuali
verbose(0, true, "Computing mapping quality:");
readCoordsTree->findOverlapping(alignment.QStart, readLength - alignment.QEnd, results);
+ if (0 == results.size()) return(0);
int mqSum = 0;
int mqCount = 0;
for (int j = 0; j < results.size(); ++j) {
|
