gatk_print_reads

gatk_print_reads recalibrates base quality scores in sequencing reads to improve the accuracy of downstream variant calling.

Key Features:

• GATK component: Implemented as a component of the Genome Analysis Toolkit (GATK).
• Two-pass framework: Operates during the second pass of a two-pass processing framework for base quality score recalibration (BQSR).
• Per-base traversal: Traverses each read by base to evaluate and adjust quality information.
• Covariate calculation: Calculates user-specified covariates including read group, reported quality score, cycle position, and dinucleotide context.
• Hashmap-based empirical scoring: Uses covariates as keys in a large hashmap to derive empirical base quality scores.
• Quality score replacement: Overwrites existing base quality scores in reads with the derived empirical quality scores.
• Integration with GATK workflows: Fits into GATK analytic steps including read mapping, local realignment around indels, BQSR, SNP discovery and genotyping, and machine learning approaches to distinguish true variants from sequencing artifacts.
• Sequencing dataset applicability: Applicable to high-throughput sequencing datasets such as deep whole-genome sequencing, whole-exome captures, and multi-sample low-pass sequencing.

Scientific Applications:

• Variant calling improvement: Enhances fidelity of SNP and indel calls by correcting systematic errors in base quality scores.
• Comparative genomics and population studies: Supports analyses relevant to human disease, ancestry, and evolutionary studies by improving underlying read quality metrics.
• High-throughput sequencing projects: Used in pipelines for deep WGS, WES, and multi-sample low-pass sequencing to produce higher-quality variant datasets.

Methodology:

In the second pass of a two-pass framework, the tool traverses each read by base, calculates covariates (read group, reported quality score, cycle position, dinucleotide context), uses those covariates as keys in a large hashmap to derive empirical base quality scores, and overwrites the reads' quality scores with the empirical values.