Evaluate and compare base quality score recalibration (BQSR) tables
This tool generates plots to assess the quality of a recalibration run as part of the Base Quality Score Recalibration (BQSR) procedure.
The goal of this procedure is to correct for systematic bias that affects the assignment of base quality scores by the sequencer. The first pass consists of calculating error empirically and finding patterns in how error varies with basecall features over all bases. The relevant observations are written to a recalibration table. The second pass consists of applying numerical corrections to each individual basecall based on the patterns identified in the first step (recorded in the recalibration table) and writing out the recalibrated data to a new BAM or CRAM file.
The tool can take up to three different sets of recalibration tables. The resulting plots will be overlaid on top of each other to make comparisons easy.
|Original||-before||BEFORE||Maroon1||First pass recalibration tables obtained from applying org.broadinstitute.hellbender.transformers.BQSRReadTransformer on the original alignment.|
|Recalibrated||-after||AFTER||Blue||Second pass recalibration tables results from the application of org.broadinstitute.hellbender.transformers.BQSRReadTransformer on the alignment recalibrated using the first pass tables|
|Input||-bqsr||BQSR||Black||Any recalibration table without a specific role|
You need to specify at least one set. Multiple sets need to have the same values for the following parameters:
covariate (order is not important), no_standard_covs, run_without_dbsnp, solid_recal_mode, solid_nocall_strategy, mismatches_context_size, mismatches_default_quality, deletions_default_quality, insertions_default_quality, maximum_cycle_value, low_quality_tail, default_platform, force_platform, quantizing_levels and binary_tag_name
Currently this tool generates two outputs:
You need to specify at least one of them.
gatk AnalyzeCovariates \ -bqsr recal1.table \ -plots AnalyzeCovariates.pdf
gatk AnalyzeCovariates \ -before recal1.table \ -after recal2.table \ -plots AnalyzeCovariates.pdf
gatk AnalyzeCovariates \ -bqsr recal1.table \ -before recal2.table \ -after recal3.table \ -plots AnalyzeCovariates.pdf
This table summarizes the command-line arguments that are specific to this tool. For more details on each argument, see the list further down below the table or click on an argument name to jump directly to that entry in the list.
|Argument name(s)||Default value||Summary|
|Optional Tool Arguments|
|null||file containing the BQSR second-pass report file|
||||read one or more arguments files and add them to the command line|
|null||file containing the BQSR first-pass report file|
|null||Input covariates table file for on-the-fly base quality score recalibration|
|20||If the GCS bucket channel errors out, how many times it will attempt to re-initiate the connection|
|false||display the help message|
||false||do not emit warning messages related to suspicious last modification time order of inputs|
|null||location of the csv intermediate file|
|null||location of the output report|
||false||display the version number for this tool|
|Optional Common Arguments|
||null||A configuration file to use with the GATK.|
||false||Whether to suppress job-summary info on System.err.|
|false||Whether to use the JdkDeflater (as opposed to IntelDeflater)|
|false||Whether to use the JdkInflater (as opposed to IntelInflater)|
||INFO||Control verbosity of logging.|
||false||display hidden arguments|
Arguments in this list are specific to this tool. Keep in mind that other arguments are available that are shared with other tools (e.g. command-line GATK arguments); see Inherited arguments above.
file containing the BQSR second-pass report file
File containing the recalibration tables from the second pass.
read one or more arguments files and add them to the command line
file containing the BQSR first-pass report file
File containing the recalibration tables from the first pass.
Input covariates table file for on-the-fly base quality score recalibration
Enables recalibration of base qualities, intended primarily for use with BaseRecalibrator and ApplyBQSR (see Best Practices workflow documentation). The covariates tables are produced by the BaseRecalibrator tool. Please be aware that you should only run recalibration with the covariates file created on the same input bam(s).
A configuration file to use with the GATK.
If the GCS bucket channel errors out, how many times it will attempt to re-initiate the connection
int 20 [ [ -∞ ∞ ] ]
display the help message
do not emit warning messages related to suspicious last modification time order of inputs
If true, it won't show a warning if the last-modification time of the before and after input files suggest that they have been reversed.
location of the csv intermediate file
Output csv file name.
location of the output report
Output report file name.
Whether to suppress job-summary info on System.err.
display hidden arguments
Whether to use the JdkDeflater (as opposed to IntelDeflater)
Whether to use the JdkInflater (as opposed to IntelInflater)
Control verbosity of logging.
The --verbosity argument is an enumerated type (LogLevel), which can have one of the following values:
display the version number for this tool
GATK version 188.8.131.52 built at 09-43-2018 09:43:10.