Browse code

Tweak findOverlaps() in read.bismark()

Peter Hickey authored on 13/06/2018 00:25:08
Showing 1 changed files

... ...
@@ -215,7 +215,7 @@
215 215
 
216 216
     # Construct 'M' and 'Cov' matrices -----------------------------------------
217 217
 
218
-    ol <- findOverlaps(loci_from_this_sample, loci)
218
+    ol <- findOverlaps(loci_from_this_sample, loci, type = "equal")
219 219
     M <- matrix(rep(0L, length(loci)), ncol = 1)
220 220
     Cov <- matrix(rep(0L, length(loci)), ncol = 1)
221 221
     M[subjectHits(ol)] <- dt[queryHits(ol), ][["M"]]