Hi,

I have a question regarding readGAlignments in the GenomicAlignments package. When we try to use readGAlignments within bplapply with multiple cores, we noticed that multiple number of cpus used within each thread, which caused some problem for us. Does readGAlignments automatically use parallel setting? Is it possible for users to customize it by supplying some arguments, tried to look into the scanBamParam and can't seem to find anything?

Any help is appreciated.

Many thanks