Ah yes, I had kind of assumed that MACS2 had multithread support but according to this feature request on their GitHub it isn't supported yet. Guess I'll just request one core and a ton of memory. Thanks!
How much memory / how many cores should I request for peak calling?
I've never run ChIPseq before, so I've been writing a pipeline from scratch to perform all the necessary read trimming, alignment, etc. I've successfully gotten to the peak calling step, which I will be performing using MACS2. My question is, how much memory and how many cores should I request from our HPC cluster for each peak call? I have no idea how computationally intensive the process is, and I'd like it to be completed fairly quickly, while not requesting so many resources that the jobs I submit take forever to get picked up. We use slurm as a job scheduler, if that information makes any difference (doubt it does).
1 answer
It will depend on the peak caller you are planning to use. If the peak caller is not thread-aware/multi-threaded you can just ask for a single core. Asking for more cores is not going to speed the process up or allow the program to use multiple cores, if it is not designed/written to do so.
You shouldn't need more than 8 GB of memory for MACS.
• 0 views
•
link
written
by
jared.andrews07
Log in to answer this question.
More posts like this
-
Unsupported attributes (CPU/Memory) for pbcromwell workflow manager
written by Zeng Hao •
-
Cromwell slurm singularity call caching
written by asmariyaz23 •
-
How to submit multithreaded jobs on SLURM and determine number of threads?
written by steve
-
HISAT index incomplete even with HUGE RAM
written by emcc •
-
Blog: "No-coding" Bioinformatics - learn by analyzing public domain data
written by elia.brodsky
-
Tutorial: InsideDNA: Running approximate Bayesian computation inferences about population history f…
written by shilparaopradeep •
