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 and it also resulted in an error. Each of the files has 30 to 40 million mappings, so I wouldn't expect them to be too big. I'll stick to bplapply. -------------------------------------- Dario Strbenac University of Sydney Camperdown NSW 2050 Australia _______________________________________________ Bioc-devel@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/bioc-devel