If you happen to see a question you know the answer to, please do chime in and help your fellow community members. We encourage our fourm members to be more involved, jump in and help out your fellow researchers with their questions. GATK forum is a community forum and helping each other with using GATK tools and research is the cornerstone of our success as a genomics research community.We appreciate your help!
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.
GATK4 Mutect2 bundle resources
Is it possible to get a description of how the gnomAD files (for Mutect2 and CalculateContamination) were generated? Additionally is possible to get some guideline on how they should be generated for non-human species (since gnomAD is human only)?
I saw the wdl script on GitHub. It's very helpful, but incomplete. For example, there is a "minimum_allele_frequency" variable that is expected as an input, so its value is unknown. Also, there is some regex that expects a specifically formatted VCF. Since gnomAD provides only chromosome VCFs, it's possible that the whole genome VCF would differ depending on how exactly those are merged. The VCF can be acquired from other sources (such as Ensembl) and that will also have slightly different fields.
More generally, af-only-gnomad.hg38.vcf.gz contains ~260M records, but small_exac_common_3.hg38.vcf.gz has only ~62K. That is very heavily filtered. That seems like extremely heavy filtering.