The current GATK version is 3.2-2

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.
Bug Bulletin: The GenomeLocPArser error in SplitNCigarReads has been fixed; if you encounter it, use the latest nightly build.

Haplotype Caller -dcov flag

Posts: 6Member

Is there any way to override/re-implement the -dcov flag (which has been removed from the haplotype caller). I have a region that is almost 700 base pairs deep and continually causes a memory overload error. How can I restrict the coverage depth effectively?

Thanks, -Paul

Tagged:

Answers

• Posts: 6,176Administrator, GATK Developer admin

Hi Paul,

What version are you using? The -dcov option is available in the current version, which is 2.6-4.

Geraldine Van der Auwera, PhD

• Posts: 6Member

I am using GATK v2.5-2, maybe it is just the online documentation that is missing the -dcov discription, but I was unable to discern whether or not -dcov refers to the number of reads starting at a single location (i.e. -dcov = 2 would only allow 2 reads to begin at the same place), or if it refers to the number of reads at any location (i.e. -dcov =2 would only allow 2 reads to appear at a specific location in the genome).

Thanks, -Paul

• Posts: 6Member

Perfect, thanks for the quick responses!

Sign In or Register to comment.