Test-drive the GATK tools and Best Practices pipelines on Terra
Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.
DepthOfCoverage missing to report some locations
Dear GATK team
We recently noticed that some DepthOfCoverage output (no suffix: per locus coverage) is missing to report some locations of the BED file input. This was a bit more odd when the same locations were actually reported for a different base-quality cut-off.
Specifically, there are locations of the BED file that have output entries for -mbq 30 but don't for -mbq 20.
Irrespective of -mbq, from what it's read from the documentation (no suffix: per locus coverage) I would assume that the output contain all locations of the BED reported even if the coverage is zero. Right?
Would you please clarify how the actual engine works?