Chapter 3: Step 1 — Quality Control (QC) & Preprocessing

Chapter 3: Step 1 — Quality Control (QC) & Preprocessing

Every clinical diagnostic assay starts with raw binary reads emitted by sequencing platforms. In this chapter, we master the mathematics of Phred quality scores, implement automated quality reporting using FastQC, and deploy high-performance trimming algorithms using FastP and Trimmomatic to eliminate technical artifacts.

The 4-Line Anatomy of a FASTQ Record

Each sequencing read occupies exactly four lines in a standard FASTQ file:

@A00123:456:H5TY2DSX7:1:1101:1024:1000 1:N:0:GATCGT  <- Line 1: Header (Instrument, Run, Flowcell, Tile, X:Y coordinates, Index)
ATGCGATCGATCGTAGCTAGCTAGGCTAGCTAACTAGCTAGCTA         <- Line 2: Nucleotide Sequence Call (A, C, G, T, N)
+                                                    <- Line 3: Separator (optional repeat of header)
FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF         <- Line 4: ASCII Phred Quality Characters (ASCII 33 offset)

Phred Quality Score Mathematics

The Phred score ($Q$) is a logarithmic transformation of the estimated base-calling error probability ($P_e$):

Production Preprocessing Pipelines: FastP vs Trimmomatic

In modern high-throughput labs, FastP has largely replaced legacy Java-based tools due to its multithreaded C++ architecture, automatic adapter detection, and built-in poly-G tail removal (crucial for two-color Illumina chemistry where no-signal corresponds to G):

# Modern High-Performance Trimming with FastP
fastp     --in1 sample_R1.fastq.gz     --in2 sample_R2.fastq.gz     --out1 clean_R1.fastq.gz     --out2 clean_R2.fastq.gz     --detect_adapter_for_pe     --trim_poly_g     --cut_front --cut_tail     --cut_window_size 4     --cut_mean_quality 20     --length_required 50     --thread 16     --html qc_report.html     --json qc_metrics.json

# Legacy Trimmomatic Pipeline Execution
trimmomatic PE -threads 16     sample_R1.fastq.gz sample_R2.fastq.gz     clean_paired_R1.fq.gz clean_unpaired_R1.fq.gz     clean_paired_R2.fq.gz clean_unpaired_R2.fq.gz     ILLUMINACLIP:adapters/TruSeq3-PE.fa:2:30:10:2:keepBothReads     LEADING:3 TRAILING:3     SLIDINGWINDOW:4:20     MINLEN:50