Hello,
I fixed the suggested test case, since the command didn't specify the
connection and produced an error. It appears to work without a problem.
> bamFile <- mappedReadsGenome[2]
> length(serialize(readGAlignmentPairs(bamFile, strandMode=2), NULL))
[1] 1329295005
I also tried mc.cores = 2
Good day,
I found that readGAlignmentPairs fails when used inside an mclapply loop but
not an sapply loop. I haven't had such problems with other functions when using
mclapply.
> class(mappedToGenomeFiles)
[1] "character"
> length(mappedToGenomeFiles)
[1] 13
> mappedReadsGenome <-