[Bioc-devel] readGAlignmentPairs Fails if Used Inside mclapply Loop

[Dario Strbenac]

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