Zubair Khalid

Virologist/Molecular Biologist | Veterinarian | Bioinformatician

Conventional & Molecular Virology • Vaccine Development • Computational Biology

Dr. Zubair Khalid is a veterinarian and virologist specializing in conventional and molecular virology, vaccine development, and computational biology. Dedicated to advancing animal health through innovative research and multi-omics approaches.

Dr. Zubair Khalid - Veterinarian, Virologist, and Vaccine Development Researcher specializing in Computational Biology, Multi-omics, Animal Health, and Infectious Disease Research

Section: Molecular Diagnostics

How to Interpret DNA Sequencing Chromatograms: Peaks, Quality, and Heterozygotes

The Science Laboratory at the Aspatria Agricultural college
Image by Unknown author Unknown author, Wikimedia Commons, licensed under Public domain.

Sanger sequencing chromatogram interpretation is the systematic visual and computational analysis of fluorescence-based electropherogram data to determine DNA sequence identity, assess base-calling confidence, and identify heterozygous positions or mixed bases. This method is essential for validating CRISPR/Cas9-induced mutations, confirming plasmid constructs, genotyping single-nucleotide polymorphisms (SNPs), and troubleshooting failed sequencing reactions. Chromatogram interpretation bridges raw instrument output and biologically meaningful sequence information, enabling researchers to distinguish high-quality data from artifacts, call heterozygous variants accurately, and make informed decisions about follow-up experiments.

At a Glance

Aspect Key Information
Purpose Determine DNA sequence from Sanger trace data; identify heterozygous bases and sequence quality
Input .ab1 or .scf files from capillary electrophoresis instruments
Key features Four colored peaks (A=green, T=red, C=blue, G=black); quality scores (Phred QV); peak spacing
Quality threshold QV ≥ 20 (99% base call accuracy) for reliable single-nucleotide calls; QV ≥ 30 for variant confirmation
Heterozygote detection Overlapping dual peaks at single position with ≥30% secondary peak height relative to primary
Common software SnapGene, Benchling, Chromas, 4Peaks, FinchTV, Geneious, ApE
Time required 5–15 minutes per chromatogram for trained analyst
BSL level BSL-1 for routine sequencing of non-pathogenic samples

Scientific Principle of Sanger Sequencing Chromatograms

Sanger sequencing relies on chain-termination dideoxynucleotides (ddNTPs) labeled with distinct fluorophores. During capillary electrophoresis, terminated fragments are separated by size, and a laser excites the fluorophores as fragments pass a detector. The instrument records fluorescence intensity at each wavelength across time, producing an electropherogram—a series of colored peaks where each color corresponds to a specific nucleotide base [5].

The chromatogram displays four key properties that inform interpretation:

Peak height reflects the relative fluorescence intensity at each position. Uniform peak heights across a trace indicate even incorporation and good sequencing quality. Peak height typically decreases after 600–800 bases from the primer due to accumulated termination events and polymerase fall-off.

Peak spacing should be regular, with approximately 1–2 base pairs between adjacent peaks. Irregular spacing suggests compression artifacts (common in GC-rich regions) or electrophoresis issues.

Peak shape should be sharp and symmetrical. Broad, split, or shouldered peaks indicate poor resolution, often from incomplete extension products or salt contamination.

Background fluorescence appears as low-level signal between peaks. Elevated background suggests primer-dimer artifacts, excess template, or incomplete dye-terminator removal.

The Phred quality score (QV) quantifies base-calling confidence. A QV of 20 corresponds to 99% accuracy (1 error per 100 bases), while QV 30 represents 99.9% accuracy. Most downstream applications require QV ≥ 20 for reliable analysis [5].

Instrumentation and Software Choices

Capillary Electrophoresis Platforms

Applied Biosystems (Thermo Fisher) instruments dominate Sanger sequencing, including the 3130, 3500, and 3730 series. These systems use polymer-filled capillaries, laser-induced fluorescence detection, and proprietary base-calling algorithms. Alternative platforms include Beckman Coulter CEQ series and LI-COR DNA sequencers, though these are less common in contemporary laboratories.

The choice of instrument affects chromatogram quality through:

  • Capillary length: Longer capillaries (50–80 cm) improve resolution for longer reads but increase run time
  • Polymer type: POP-7 polymer provides higher resolution than POP-4 for fragments >500 bases
  • Detection sensitivity: Newer instruments (3500 series) offer improved signal-to-noise ratios compared to older 3130 models

Analysis Software

Software Platform Key Features Cost
SnapGene Windows/Mac Integrated chromatogram viewer, automatic base calling, heterozygote detection Commercial
Benchling Web-based Cloud storage, collaboration, CRISPR editing tools Freemium
Chromas Windows Simple viewer, quality scoring, export options Free
4Peaks Mac Lightweight viewer, peak height measurement Free
FinchTV Windows/Mac/Linux Multi-trace alignment, quality visualization Free
Geneious Windows/Mac/Linux Advanced assembly, variant calling, batch processing Commercial
ApE Windows/Mac/Linux Open-source, restriction site mapping Free

For CRISPR/Cas9 validation, ICE (Inference of CRISPR Edits) analysis software provides automated chromatogram decomposition to identify indel mutations and estimate editing efficiency [1][2]. This tool compares Sanger traces from edited samples against wild-type controls to infer mutation spectra.

Controls for Reliable Interpretation

Positive Controls

  • Wild-type reference sequence: A known, high-quality chromatogram from unedited or wild-type DNA provides baseline peak heights, spacing, and quality metrics. This control is essential for identifying heterozygous positions in CRISPR-edited samples [1].
  • Confirmed heterozygous sample: A sample with a known heterozygous SNP (e.g., from a validated cell line or individual) validates the software's heterozygote detection parameters.

Negative Controls

  • No-template control (NTC): Sequencing reaction without DNA template identifies primer-dimer artifacts and contamination. NTC chromatograms should show only primer peaks or flat baseline.
  • Water control: PCR amplification with water instead of template, followed by sequencing, confirms absence of environmental DNA contamination.

Quality Controls

  • Primer-only sequencing: Sequencing reaction with primer but no template reveals primer peaks and dye-blob artifacts.
  • Duplicate reactions: Independent sequencing of the same template from separate PCR reactions identifies stochastic errors versus reproducible variants.

Conceptual Workflow for Chromatogram Interpretation

Step 1: Visual Inspection of Raw Trace

Open the chromatogram file in your chosen software. Examine the overall trace quality:

  1. Check signal intensity: The initial 20–50 bases often show irregular peak heights due to primer interference. Reliable sequence begins after this region.
  2. Assess peak uniformity: High-quality traces show consistent peak heights (±20% variation) across the readable region.
  3. Identify dye blobs: Large, broad peaks at the beginning of the trace (typically around 50–100 bases) represent unincorporated dye terminators. These should be excluded from analysis.
  4. Evaluate baseline noise: The baseline between peaks should be flat. Elevated baseline suggests template contamination or incomplete purification.

Step 2: Quality Score Assessment

Most software displays quality values as colored bars or numerical scores beneath the trace. Apply these thresholds:

  • QV ≥ 30 (green): High confidence; suitable for variant confirmation and publication-quality data
  • QV 20–29 (yellow): Moderate confidence; acceptable for routine genotyping but requires confirmation for novel variants
  • QV < 20 (red): Low confidence; bases should be called manually or excluded from analysis

For CRISPR/Cas9 validation, focus on the region spanning the predicted cut site (typically 50–200 bases from the primer). Quality scores in this region must be ≥20 for reliable indel detection [1][2].

Step 3: Base Calling Verification

Compare the software's automated base calls with visual inspection of peak colors:

  1. Confirm each base call: For positions with QV ≥ 20, automated calls are usually correct. For lower-quality regions, manually assign bases by identifying the highest peak at each position.
  2. Resolve ambiguities: When two peaks of similar height appear at one position, the software may call "N" (any base) or the higher peak. Manual inspection determines whether this represents a true heterozygote or artifact.
  3. Check homopolymer runs: Tracts of identical bases (e.g., AAAAA) often show decreasing peak heights and increased spacing errors. Count peaks carefully; the software may under- or over-call repeat length.

Step 4: Heterozygote Detection

Heterozygous positions appear as two overlapping peaks at a single nucleotide position. Key criteria for calling a true heterozygote:

  • Secondary peak height ≥30% of primary peak height: Lower secondary peaks may represent background noise or dye artifacts
  • Consistent appearance in both forward and reverse reads: True heterozygotes should be detectable from both directions
  • Peak shape symmetry: Both peaks should show similar width and shape
  • Position within high-quality region: Heterozygote calls in regions with QV < 20 require confirmation

For CRISPR/Cas9-edited samples, heterozygous indels produce overlapping traces downstream of the cut site. This "mixed trace" pattern—where peaks become unreadable after the mutation—indicates successful editing [1][2]. ICE analysis quantifies this by comparing trace decomposition against wild-type controls.

Step 5: Sequence Alignment and Variant Calling

Export the consensus sequence (including ambiguous bases) and align to a reference sequence using alignment software:

  1. Identify single-nucleotide variants (SNVs): Positions where the sample sequence differs from reference
  2. Detect insertions/deletions (indels): Gaps in alignment indicate indels; confirm by examining chromatogram for shifted peaks
  3. Validate with reverse read: All variants should be confirmed in the complementary strand sequence

Quality Checks and Troubleshooting

Common Chromatogram Artifacts

Observation Likely Cause Discriminating Check
Broad, rounded peaks throughout Excess template DNA (100–500 ng recommended) Repeat with 1:5 and 1:10 template dilutions
Multiple peaks at single positions (not heterozygote) Primer-dimer or secondary PCR products Run PCR product on agarose gel; purify single band
Sudden drop in peak height after clean sequence Secondary structure (GC-rich region) Add DMSO (5%) or betaine (1 M) to sequencing reaction
Gradual peak height decrease from beginning Insufficient template or degraded DNA Quantify template; check DNA integrity on gel
Elevated baseline with small peaks Salt contamination from PCR Purify PCR product with column or ethanol precipitation
Peaks with shoulders or split peaks Incomplete extension products Increase extension time in cycle sequencing
No peaks except primer peak Failed sequencing reaction Verify primer Tm; check primer-template complementarity
All four colors present at every position Mixed template (contamination) Sequence individual colonies or re-purify template
Peaks spaced irregularly Electrophoresis problem or capillary issue Re-run sample; check instrument maintenance logs

Troubleshooting Table

Problem Observation Likely Cause Recommended Action
Weak signal Peak heights <200 RFU Insufficient template Increase template to 100–200 ng; verify DNA concentration
Strong signal but poor quality Peak heights >5000 RFU with high background Excess template Dilute template 1:5; reduce cycle number
Dye blob interference Large broad peak at 50–100 bases Incomplete dye-terminator removal Purify sequencing reaction with column or ethanol precipitation
Heterozygote not reproducible Secondary peak appears in only one direction PCR artifact or sequencing error Repeat from independent PCR; sequence both strands
Indel detection failure Clean trace but no editing confirmed Editing below detection threshold Increase sequencing read length; use ICE analysis with wild-type control
GC-rich region compression Peaks too close together Secondary structure Add 5% DMSO; sequence with GC-rich buffer; use higher denaturation temperature

Result Interpretation and Documentation

Reporting Heterozygous Variants

When documenting heterozygous positions, include:

  1. Position relative to reference: e.g., "chr7:g.55249018G>A"
  2. Peak height ratio: e.g., "G peak 65%, A peak 35%"
  3. Quality scores: QV for both primary and secondary peaks
  4. Confirmation method: e.g., "Confirmed in both forward and reverse reads"
  5. Biological context: e.g., "Heterozygous missense variant in exon 4 of TP53"

CRISPR/Cas9 Editing Validation

For CRISPR-edited samples, chromatogram interpretation follows specific criteria [1][2]:

  1. Wild-type control: Clean, single-peak trace throughout the target region
  2. Edited sample: Mixed trace pattern beginning at or near the predicted cut site (3–4 bases upstream of PAM)
  3. ICE analysis output: Indel percentage, mutation spectrum, and R² value (goodness of fit)
  4. Confirmation: Sanger sequencing of cloned PCR products to identify individual alleles

Document the following for each edited sample:

  • Target gene and exon
  • Guide RNA sequence and PAM location
  • Predicted cut site coordinates
  • Chromatogram file name and date
  • Software version and analysis parameters
  • Indel percentage and mutation types identified

Limitations of Chromatogram Interpretation

  1. Detection threshold: Heterozygous variants with allele frequency <15–20% may not be reliably detected by visual inspection
  2. Homopolymer regions: Tracts of >6 identical bases cannot be accurately sized by Sanger sequencing
  3. GC-rich templates: Regions with >65% GC content often produce poor-quality traces
  4. Large indels: Insertions or deletions >50 bases may not be resolved by capillary electrophoresis
  5. Mosaicism: Low-frequency mosaic mutations in F0 CRISPR animals may be missed [2]

Biosafety Considerations

Sanger sequencing of routine DNA samples from non-pathogenic organisms falls under BSL-1 containment as defined by the CDC/NIH Biosafety in Microbiological and Biomedical Laboratories (BMBL) guidelines [3]. Standard microbiological practices apply:

  • Sample handling: Use gloves when handling DNA samples and sequencing reagents
  • Waste disposal: Discard sequencing reaction plates and capillary arrays as biohazardous waste if samples originated from BSL-2 organisms
  • Chemical safety: Formamide used in sample preparation is a reproductive toxin; handle in chemical fume hood
  • Instrument decontamination: Follow manufacturer's protocols for capillary and polymer block cleaning

For research involving recombinant DNA (including CRISPR/Cas9-edited organisms), follow NIH Guidelines for Research Involving Recombinant or Synthetic Nucleic Acid Molecules [4]. Institutional Biosafety Committee (IBC) approval is required for experiments involving:

  • Pathogenic organisms
  • Toxin production
  • Gene drive constructs
  • Select agent sequences

Frequently Asked Questions

Q1: How do I distinguish a true heterozygous SNP from a sequencing artifact?

A true heterozygous SNP shows two overlapping peaks at a single position with the secondary peak height ≥30% of the primary peak. The pattern should be reproducible in both forward and reverse sequencing reads. Artifacts typically appear as broad, irregular peaks or as secondary peaks that change position between reads. Confirm heterozygotes by sequencing an independent PCR product and checking for the variant in both directions.

Q2: What should I do when my chromatogram shows mixed peaks after the CRISPR cut site?

Mixed peaks downstream of the predicted cut site indicate successful CRISPR/Cas9 editing with multiple alleles present in the sample. Use ICE analysis software to quantify editing efficiency and identify the most common mutations. For clonal analysis, PCR-amplify the target region, clone into a plasmid vector, and sequence individual colonies to determine exact mutation sequences [1][2].

Q3: Why does my sequence quality drop dramatically after 500 bases?

Quality degradation after 500–800 bases is normal in Sanger sequencing due to polymerase fall-off and accumulated termination events. To extend read length, use fresh sequencing reagents, increase template amount, or design internal primers for overlapping reads. For targets requiring >1000 bases of continuous sequence, consider primer walking or next-generation sequencing approaches.

Q4: Can I trust automated base calls in low-quality regions?

Automated base calls in regions with QV < 20 should be treated as provisional. Manually inspect the chromatogram and assign bases by identifying the highest peak at each position. For critical applications (variant confirmation, clinical reporting), require QV ≥ 30 for all base calls. When manual inspection cannot resolve ambiguous positions, repeat the sequencing reaction with optimized conditions.

References and Further Reading

  1. Varshney GK, Carrington B, Pei W, et al. A high-throughput functional genomics workflow based on CRISPR/Cas9-mediated targeted mutagenesis in zebrafish. Nature Protocols. 2016. https://pubmed.ncbi.nlm.nih.gov/27809318/

    • Provides detailed methods for Sanger sequencing validation of CRISPR-induced mutations, including chromatogram analysis for indel detection.
  2. Jordan DC, Rivers ML, Belmont L, et al. Gene editing without a genome: generation and validation of F0 CRISPR mutants in gastropod mollusc Crepidula fornicata. bioRxiv. 2025. https://doi.org/10.1101/2025.11.26.690768

    • Demonstrates chromatogram interpretation for CRISPR validation in non-model organisms using ICE analysis.
  3. CDC and NIH. Biosafety in Microbiological and Biomedical Laboratories (BMBL), 6th Edition. U.S. Department of Health and Human Services. 2020. https://www.cdc.gov/labs/bmbl/index.html

    • Authoritative reference for laboratory biosafety levels and containment practices.
  4. National Institutes of Health. NIH Guidelines for Research Involving Recombinant or Synthetic Nucleic Acid Molecules. https://osp.od.nih.gov/policies/biosafety-and-biosecurity-policy/nih-guidelines-for-research-involving-recombinant-or-synthetic-nucleic-acid-molecules/

    • Regulatory framework for recombinant DNA research, including CRISPR/Cas9 applications.
  5. National Center for Biotechnology Information. Molecular Biology and Laboratory Methods. NCBI Bookshelf. https://www.ncbi.nlm.nih.gov/books/

    • Comprehensive collection of molecular biology protocols and methods references.

Related Articles