add 1 to read count to avoid MathError (zeroDivison)

nf-core · Oct 24, 2018 · 3c398e5 · 3c398e5
1 parent 81b1508
commit 3c398e5
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/bin/normalizedReadCount b/bin/normalizedReadCount
@@ -244,6 +244,8 @@ if __name__ == "__main__":
     writeBam(inbam=BAM1, outbam=outbam1, commonReads=commonReads)
     writeBam(inbam=BAM2, outbam=outbam2, commonReads=commonReads)
 
+    nb1 += 1
+    nb2 += 1
     NormalizedReadRatio = log2((nb1 / gs1) / (nb2 / gs2))
 
     # Template output file structure