Sequence Alignment: Score and Identity
Global Needleman–Wunsch alignment maximizes a sequence score; identity counts matching symbols in the aligned columns.
Free to browse · flashcards or a multiple-choice quiz · import from Essential
Formula
Alignment score in arbitrary score points · g and s in the same score points · identity in %
Variables & units
| Symbol | Meaning | Unit |
|---|---|---|
| Best score for prefixes ending at positions i and j | score points | |
| Substitution score for the symbol pair | score points | |
| Positive linear gap penalty per gap column | score points | |
| Number of identical symbol pairs in the final alignment | count | |
| Number of alignment columns, including gap columns | count | |
| Percent identity for the stated alignment | % |
LaTeX
F_{i,j}=\max\{F_{i-1,j-1}+s(x_i,y_j),\ F_{i-1,j}-g,\ F_{i,j-1}-g\},\quad I=\frac{N_{match}}{L_{align}}\cdot100\%
Original schematic by Quanta
Free with a Quanta account. Your download stays ready after registration.
Derivation & background
The recurrence implements global alignment with a linear gap score. A diagonal step scores a symbol pair; horizontal or vertical steps insert a gap. Initialize F00=0, Fi0=−ig and F0j=−jg. Identity is a separate descriptive statistic whose denominator must be declared; this page counts gap columns in alignment length. Local alignment, affine gap costs and substitution matrices require different settings. Similarity or identity alone does not establish shared function or ancestry.
Exam blueprint
Validity range
Global alignment of two sequences with a linear gap score; initialization and scoring must be stated. Identity here includes gap columns. Other alignment goals require different methods.
Derivation steps
For two prefixes, choose the best score from a symbol pair, a gap in sequence 1, or a gap in sequence 2.
- 1Initialize F00=0 and the borders with −i·g and −j·g.
- 2At each cell, maximize the diagonal step plus substitution score and the two gap steps.
- 3Trace back the best path; then count identical columns and divide by alignment length.
Rearrangements
Match count from identity and alignment length
Use I as a percent and keep the same gap convention.
Task variant
Align ACT and AGT globally: match +2, mismatch −1, gap −2. Find score and identity.
Without gaps: score=2−1+2=3.
Two of three columns match, so I=2/3×100%=66.7%.
Common mistakes
Equating alignment score with percent identity or computing identity without an alignment.
The score depends on substitution and gap values; identity counts identical columns under a stated alignment convention.
Exam context
- Compare global DNA or protein sequences, state scoring rules, and interpret metrics separately.
Formula cluster
Bioinformatics and sequences
Alignment complements composition and assembly metrics but does not replace a homology assessment.
Worked example
Align ACT against AGT globally with match +2, mismatch −1 and gap −2. The gap-free alignment scores 2−1+2=3 points. Two of three columns match, so identity=2/3×100%=66.7%. Check: adding a gap changes both the score and the denominator and must follow the same stated convention.
Applications
Compare global similarity of DNA or protein sequences computationally and document the scoring convention.
How do you calculate with Sequence Alignment: Score and Identity?
Here is how to work through a typical Sequence Alignment: Score and Identity (Fᵢⱼ = max(Diagonale, Lücke, Lücke); Identität = Treffer/Ausrichtung) task step by step:
- 1
Task
Align ACT and AGT globally: match +2, mismatch −1, gap −2. Find score and identity.
Solution path
Without gaps: score=2−1+2=3.
Two of three columns match, so I=2/3×100%=66.7%.
Frequently asked questions about Sequence Alignment: Score and Identity
What does an alignment score mean?+
An alignment score adds rewards for matching or substituted symbols and subtracts gap costs. For ACT versus AGT with match +2, mismatch −1 and gap −2, the gap-free alignment scores 2−1+2=3. This value is neither a probability nor a percentage; it depends on the scoring matrix and gap penalties. Compare scores only when method, parameters, sequence type and alignment goal are comparable. A higher score means a better computational alignment under that same convention, but it does not automatically establish biological function.
How does the Needleman–Wunsch recurrence work?+
Dynamic programming fills a table for every pair of prefixes. Each cell checks three paths: align two symbols diagonally and add their substitution score, or place a gap in either sequence and subtract the gap penalty. Border values account for aligning an entire prefix to gaps. After filling the table, trace back from the final cell to the origin to recover the optimal global alignment. This reuses subproblems instead of enumerating every possible alignment. Affine gap costs require an expanded state model.
Is sequence identity the same as sequence similarity?+
No. Identity counts exactly matching letters or amino acids in aligned columns. Similarity may also reward conservative amino-acid substitutions through a substitution matrix. A protein sequence can therefore obtain a good score without every position being identical. Identity also depends on the alignment and its denominator; this page includes gap columns in alignment length. Report percent identity together with alignment length, gap fraction and method. Identity or similarity alone proves neither homology nor shared function, because composition and database searching can mislead interpretation.
How do gaps affect identity?+
A gap is introduced when the method explains a sequence difference as an insertion or deletion. In the score, each linear gap costs the chosen value g; the remaining symbols may then align better. Identity requires a declared denominator convention. Here, gap columns remain in the denominator: ten identical positions across twelve columns yield 83.3 percent identity. Some software reports a different denominator. Avoid comparing percentages until you have checked gap conventions, terminal gaps and alignment type. State this convention in reports.
Does a high alignment prove a shared biological function?+
No. An alignment is first a statistically scored comparison of two strings. A high score can depend on sequence length, composition, database choice, scoring matrix or search space. A defensible homology claim needs suitable null models, significance measures, domain context and sometimes phylogenetic evidence. Even homology does not mean both proteins have the same present-day function. Use the score to prioritize further analysis and assess conserved motifs, structure and experiments separately. Document parameters, database version and selection criteria. A suitable null model tests whether a hit exceeds chance expectations.
Quanta exam set
Curated exam set for "Sequence Alignment: Score and Identity":
Question (front)
Which formula describes Sequence Alignment: Score and Identity?
Answer in your set
Question (front)
How do you rearrange Sequence Alignment: Score and Identity for Match count from identity and alignment length?
Answer in your set
Question (front)
Which common mistake happens with Sequence Alignment: Score and Identity?
Answer in your set
+ 8 more cards: units, variables, derivation, example, exam task
The formula set contains 11 cards. After importing, review them and start practicing; your ratings determine later review intervals.
Subject references and further reading
This list names subject references and further reading. A page or section number appears when it is part of the reference; these are not verbatim quotations.
Related formulas
More formulas: Biochemistry and molecular biology
Retain Sequence Alignment: Score and Identity for exams
Create a curated FSRS exam set for Fᵢⱼ = max(Diagonale, Lücke, Lücke); Identität = Treffer/Ausrichtung: formula recall, variables, derivation, rearrangement, worked example, common mistakes and exam context.
Curated formula set · LaTeX · FSRS spaced repetition · import from Essential