If you happen to see a question you know the answer to, please do chime in and help your fellow community members. 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.
Release notes for GATK version 2.4
GATK 2.4 was released on February 26, 2013. Highlights are listed below. Read the detailed version history overview here: http://www.broadinstitute.org/gatk/guide/version-history
Important note 1 for this release: with this release comes an updated licensing structure for the GATK. Different files in our public repository are protected with different licenses, so please see the text at the top of any given file for details as to its particular license.
Important note 2 for this release: the GATK team spent a tremendous amount of time and engineering effort to add extensive tests for many of our core tools (a process that will continue into future releases). Unsurprisingly, as part of this process many small (and some not so small) bugs were uncovered during testing that we subsequently fixed. While we usually attempt to enumerate in our release notes all of the bugs fixed during a given release, that would entail quite a Herculean effort for release 2.4; so please just be aware that there were many smaller fixes that may be omitted from these notes.
Base Quality Score Recalibration
- The underlying calculation of the recalibration has been improved and generalized so that the empirical quality is now calculated through a Bayesian estimate. This radically improves the accuracy in particular for bins with small numbers of observations.
- Added many run time improvements so that this tool now runs much faster.
- Print Reads writes a header when used with the -BQSR argument.
- Added a check to make sure that BQSR is not being run on a reduced bam (which would be bad).
- The --maximum_cycle_value argument can now be specified during the Print Reads step to prevent problems when running on bams with extremely long reads.
- Fixed bug where reads with an existing BQ tag and soft-clipped bases could cause the tool to error out.
- Fixed the QUAL calculation for monomorphic (homozygous reference) sites (the math for previous versions was not correct).
- Biased downsampling (i.e. contamination removal) values can now be specified as per-sample fractions.
- Fixed bug where biased downsampling (i.e. contamination removal) was not being performed correctly in the presence of reduced reads.
- The indel likelihoods calculation had several bugs (e.g. sometimes the log likelihoods were positive!) that manifested themselves in certain situations and these have all been fixed.
- Small run time improvements were added.
- Extensive performance improvements were added to the Haplotype Caller. This includes run time enhancements (it is now much faster than previous versions) plus improvements in accuracy for both SNPs and indels. Internal assessment now shows the Haplotype Caller calling variants more accurately than the Unified Genotyper. The changes for this tool are so extensive that they cannot easily be enumerated in these notes.
- The QD annotation is now divided by the average length of the alternate allele (weighted by the allele count); this does not affect SNPs but makes the calculation for indels much more accurate.
- Fixed Fisher Strand annotation where p-values sometimes summed to slightly greater than 1.0.
- Fixed Fisher Strand annotation for indels where reduced reads were not being handled correctly.
- The Haplotype Score annotation no longer applies to indels.
- Added the Variant Type annotation (not enabled by default) to annotate the VCF record with the variant type.
- The DepthOfCoverage annotation has been renamed to Coverage.
- Several small run time improvements were added to make this tool slightly faster.
- By default this tool now uses a downsampling value of 40x per start position.
- Fixed bug where some reads with soft clipped bases were not be realigned.
- Run time performance improvements added where one uses the PRIORITIZE or REQUIRE_UNIQUE options.
- The --regenotype functionality has been removed from SelectVariants and transferred into its own tool: RegenotypeVariants.
- Removed the GenotypeConcordance evaluation module (which had many bugs) and converted it into its own tested, standalone tool (called GenotypeConcordance).
- The VariantContext and related classes have been moved out of the GATK codebase and into Picard's public repository. The GATK now uses the variant.jar as an external library.
- Added a new Read Filter to reassign just a particular mapping quality to another one (see the ReassignOneMappingQualityFilter).
- Added the Regenotype Variants tool that allows one to regenotype a VCF file (which must contain likelihoods in the PL field) after samples have been added/removed.
- Added the Genotype Concordance tool that calculates the concordance of one VCF file against another.
- Bug fix for VariantsToVCF for records where old dbSNP files had '-' as the reference base.
- The GATK now automatically converts IUPAC bases in the reference to Ns and errors out on other non-standard characters.
- Fixed bug for the DepthOfCoverage tool which was not counting deletions correctly.
- Added Cat Variants, a standalone tool to quickly combine multiple VCF files whose records are non-overlapping (e.g. as produced during scatter-gather).
- The Somatic Indel Detector has been removed from our codebase and moved to the Broad Cancer group's private repository.
- Fixed Validate Variants rsID checking which wasn't working if there were multiple IDs.
- Picard jar updated to version 1.84.1337.
- Tribble jar updated to version 1.84.1337.
- Variant jar updated to version 1.85.1357.