The current GATK version is 3.7-0
Examples: Monday, today, last week, Mar 26, 3/26/04

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Powered by Vanilla. Made with Bootstrap.
GATK 3.7 is here! Be sure to read the Version Highlights and optionally the full Release Notes.
Register now for the upcoming GATK Best Practices workshop, Feb 20-22 in Leuven, Belgium. Open to all comers! More info and signup at

Recommendation on performance using scatter/gather

LaurentLaurent Member, Collaborator Posts: 43 ✭✭
edited April 2013 in Ask the GATK team

Dear all,

I am currently running an analysis using the HaplotypeCaller on 300 large BAM files on our cluster and decided to chunk the the genome in 3MB bins in order for them to be processed in a decent time. I'm however experiencing very long runtimes as more and more jobs get scheduled to run in parallel on the same files. Looking at the GATK options, I saw these 2 that I thought could be of help and was wondering what were the recommendation for using them:

More precisely, does the num_bam_file_handles increase processing time by a lot? and what is the default value for --read_buffer_size ?

Thanks a lot,


Best Answer


  • LaurentLaurent Member, Collaborator Posts: 43 ✭✭

    Hi Mark,

    Thanks for the explanation! I will digg more into the problem, but at the moment what I am reporting is only "observation" of my runtimes getting extremely high when running in parallel. I have observed similar problems when running other walkers using scatter/gather in the past on our cluster. I'll give a shot at extracting the regions using PrintReads to a local scratch beforehand and let you know if this helps! I'll also use the nightly build to benefit from the latest improvements.

Sign In or Register to comment.