Tutorials › Statistical Genetics › Variant Records and Allele Normalization

Genotyping Data and Quality Control · Tutorial 80 of 989

Variant Records and Allele Normalization

Learn to read small-variant records and recognize when different-looking insertion or deletion records describe the same alternate sequence.

Beginner 10 min read

What You'll Learn

  • Define a genomic locus and reference and alternate alleles
  • Distinguish biallelic records from multiallelic records
  • Interpret small variants as changes to a reference sequence
  • Trim redundant bases to obtain minimal allele representations
  • Left-align insertions and deletions in repetitive sequence
  • Explain why consistent normalization matters when comparing variant records

Why Variant Representation Matters

A small genomic change can be described in more than one way. This is especially common for insertions and deletions in repeated sequence: deleting one of several consecutive identical bases can produce the same resulting DNA sequence, even when the stated position differs. Before comparing variants across samples or data sets, it is therefore important to distinguish a biological sequence change from the particular record used to describe it.

This tutorial defines the basic parts of a small-variant record and introduces normalization: a set of representation rules used to make equivalent records consistent. The focus is on substitutions and short insertions and deletions, not on assigning genotypes to individuals.

Definition: A locus is a genomic location or region at which a sequence difference is described. A reference allele is the sequence present in the specified reference genome at that record. An alternate allele is a sequence used to describe a possible difference from the reference allele. A variant record identifies the reference context and one or more alternate alleles; it does not, by itself, say which allele an individual carries.

Reference and Alternate Alleles

A variant description depends on a reference genome: a specified genome assembly and its sequence. For the examples here, positions are numbered from 1 within a chromosome, and a record’s position is the coordinate of the first base in its reference allele. The reference allele, written REF, is the sequence beginning at that position; an alternate allele, written ALT, replaces that sequence. The unchanged sequence before and after the record is retained.

For example, suppose position 104 in a hypothetical reference has the base C. A record with position 104, REF C, and ALT T describes replacing C with T. This is a single-nucleotide substitution. The record says that T is an alternative sequence at that locus; it does not say how common T is, whether it is associated with a trait, or which individuals carry it.

For a deletion, REF is longer than ALT; for an insertion, ALT is longer than REF. In common anchored representations, both allele strings contain at least one base, even when the biological change is an insertion or deletion. That retained base provides sequence context for locating the change. The exact representation convention matters when comparing records, so the reference genome and coordinate convention must be held fixed.

Worked Example: Reading a Substitution Record

In a hypothetical reference genome, position 104 is C. A record gives position 104, REF C, ALT T. What sequence change does it describe?

The reference sequence at the locus is C, and the alternate sequence specified by the record is T. Replacing the one-base REF string with the one-base ALT string gives a C-to-T substitution at position 104. Both alleles have length 1, so neither an insertion nor a deletion is described.

This interpretation is limited to the record itself: it identifies a possible alternate sequence, not a carrier, a genotype, or a biological effect.

Biallelic and Multiallelic Records

The number of alternate alleles listed in a record determines whether the record is biallelic or multiallelic. “Biallelic” counts the reference allele as one of the possible alleles: a biallelic record has one REF allele and one ALT allele. A multiallelic record has one REF allele and at least two distinct ALT alleles. These terms describe the alleles listed at a locus, not the number of individuals observed or the number of copies carried by any one individual.

Definition: A biallelic record lists one reference allele and one alternate allele. A multiallelic record lists one reference allele and two or more alternate alleles. The distinct allele strings are interpreted relative to the same reference context.

Worked Example: Interpreting a Multiallelic Record

At a hypothetical position, the reference base is A. A record lists REF A and ALT C,G. Identify the alleles and classify the record.

The listed alleles are A (the reference allele), C (one alternate allele), and G (a second alternate allele). There are two ALT alleles, so this is a multiallelic record. It describes two possible substitutions relative to A: A-to-C and A-to-G.

The record does not assert that both alternate alleles occur in the same individual, or that either is present in any particular sample. Those questions require genotype data, which are separate from the variant description.

Normalization: A Consistent Representation

Normalization puts a variant record into a standard form so that equivalent descriptions can be compared. For small variants, two useful operations are trimming redundant shared bases and, for indels in repeated sequence, left alignment. A minimal representation removes shared sequence from the allele strings when doing so leaves valid nonempty strings. Left alignment chooses the leftmost equivalent placement on the same reference sequence.

A basic check for equivalence is to compare the alternate sequence produced in the same reference context. If the unchanged sequence to the left is \(L\), the REF string is \(R\), the ALT string is \(A\), and the unchanged sequence to the right is \(S\), the alternate sequence is \(L A S\). Two records are equivalent in that context if applying either one yields the same alternate sequence. A difference in written position alone does not establish that two indel records describe different changes.

Key procedure: First remove redundant common bases while keeping the allele strings nonempty. For an insertion or deletion in a repeat, shift the minimal representation left as far as possible while preserving the same alternate sequence relative to the reference. Record the reference assembly and representation convention when reporting the result.

Worked Example: Trimming a Redundant Base

At position 20 of a hypothetical reference, REF is AC and ALT is AT. Find a shorter equivalent representation.

The two strings share the first base A. Removing this common prefix advances the start coordinate by one and leaves REF C and ALT T at position 21. The original replacement, AC with AT, and the shortened replacement, C with T, produce the same alternate sequence: the A is unchanged, and the following C is replaced by T.

Thus the minimal representation is position 21, REF C, ALT T. Both allele strings remain nonempty. Removing the common A without advancing the coordinate would describe a different location, so coordinate adjustment is part of trimming a shared prefix.

Equivalent Insertions in a Repeat

When several neighboring reference bases are identical, the same inserted base can be placed at more than one point in the run without changing the resulting sequence. A left-aligned representation resolves this ambiguity by choosing the leftmost placement. In an anchored record, the base immediately before the insertion is retained in both allele strings.

Worked Example: Left-Aligning an Insertion

Consider the hypothetical reference sequence G C A A A T, with positions 1 through 6 in that order. One A is inserted into the run of A bases. Show equivalent anchored descriptions and determine the left-aligned record.

An insertion after position 2 can be written at position 2 as REF C, ALT CA: the unchanged C is retained and an A is added after it. The same resulting sequence can be described by inserting after position 3, 4, or 5, using respectively position 3, REF A, ALT AA; position 4, REF A, ALT AA; or position 5, REF A, ALT AA.

Each record produces the sequence G C A A A A T. For example, the position-2 record replaces the C in G C A A A T with CA, giving G C A A A A T. In the other records, inserting A elsewhere within the existing run gives that same sequence. The leftmost insertion point is after position 2, so the left-aligned representation is position 2, REF C, ALT CA.

The several records are not evidence for several distinct insertion events. They are equivalent descriptions of one resulting alternate sequence in this reference context.

Equivalent Deletions in a Repeat

The same representational ambiguity occurs for deletions. In an anchored representation, the base before the deleted sequence is retained in both alleles. If a run contains identical bases, deleting any one of them may yield the same sequence. Left alignment again selects the leftmost equivalent representation.

Worked Example: Left-Aligning a Deletion

Consider the hypothetical reference sequence T G C C C A, with positions 1 through 6. One C is deleted from the run of three Cs. Write equivalent anchored descriptions and find the left-aligned record.

Deleting the C at position 3 can be represented at position 2 as REF GC, ALT G: the G at position 2 is retained, while the following C is removed. Deleting the C at position 4 can be written at position 3 as REF CC, ALT C. Deleting the C at position 5 can be written at position 4 as REF CC, ALT C.

Applying each record yields T G C C A. For the first record, replacing GC with G in T G C C C A leaves T G C C A. The other records remove a different written C but yield the same resulting sequence because the deleted bases are identical and adjacent. The leftmost representation is therefore position 2, REF GC, ALT G.

Why Normalization Matters

Without normalization, equivalent indels can appear at different coordinates or have different allele strings. A direct comparison might then treat them as distinct records, even though applying each record to the same reference produces the same alternate sequence. Normalization makes record matching and downstream summaries more reliable; it does not determine whether a variant is real, harmful, or carried by a particular person.

The reference assembly is essential. Coordinates and reference sequences can differ between assemblies, so a position-and-allele description without its reference context may be ambiguous. Normalization also has a limit: it chooses a consistent representation under a stated convention, but it cannot make records comparable if they use incompatible reference sequences or inconsistent coordinate systems.

1
Confirm the reference context.
Check the reference sequence and coordinate convention used by the record.
2
Check the allele strings.
Interpret REF as the sequence being replaced and ALT as the sequence that replaces it.
3
Trim shared sequence.
Remove redundant common bases where possible, keeping the allele strings nonempty and adjusting the position when a shared prefix is removed.
4
Left-align a repetitive indel.
Shift the representation left only when the resulting alternate sequence remains unchanged.

Common Mistake

A common error is to conclude that records at different positions must describe different indels. That reasoning ignores repeated sequence. For example, deleting one C from a run of three Cs can be written at different positions, yet all such descriptions can yield the same alternate sequence. The correction is to compare the resulting sequences in the same reference context, then use minimal, left-aligned representation to make the description consistent.

A second error is to interpret “multiallelic” as meaning that one individual carries multiple alternate alleles. The term describes the alternatives listed at a locus. Whether an individual carries a particular allele is a separate question about genotype data.

Key takeaway: A variant record describes a change relative to a reference sequence. Biallelic and multiallelic refer to the number of listed alleles; normalization makes equivalent small-variant descriptions consistent by trimming redundant sequence and left-aligning repetitive indels.

Check Your Understanding

Use the definitions and sequence comparisons above to answer these questions.

  1. A record lists REF G and ALT A,T. Is it biallelic or multiallelic, and how many alternate alleles are listed?
  2. Why must a coordinate be adjusted when a shared prefix is trimmed from both REF and ALT?
  3. In the reference sequence C A A A G, where can an A be inserted to produce the same final sequence, and which anchored position is leftmost?
  4. Two deletion records have different positions in a run of identical bases. What sequence-based check can determine whether they are equivalent?
  5. Does a multiallelic record establish that a particular individual carries more than one alternate allele? Explain.
Clinical Biostats