Site icon Ampersand Tutorials

How to Interpret BLAST Results: E-Value, Bit Score, Identity, and Coverage

BLAST results table showing E-value, bit score, percent identity and query coverage columns against a DNA sequence

Quick answer: No single number in a BLAST output decides whether a hit is biologically meaningful. The E-value tells you how statistically surprising a match is, the bit score tells you how good the alignment is on a common scale, percent identity tells you how similar the residues are, and query coverage tells you how much of your sequence the hit explains. A truthful BLAST interpretation reads these four together, plus the biology of your question.

Complements the two earlier BLAST guides on this blog: What Is BLAST in Bioinformatics and Choosing the BLAST Program. Here we focus on the result table itself.

Where your attention should start

Run a BLAST search and the results table greets you with columns of numbers: Query Cover, E value, Percent Identity, Score, Accession, and description text. Beginners often scan straight for the top row of percent identity. That instinct is half right: percent identity is useful, but it is not the number BLAST ranks by, so the top row of the table is not automatically your best answer.

BLAST orders hits by best score first, using the bit score. A long weakly-similar alignment can therefore out-rank a shorter, nearly perfect one, and a low-E-value hit for your real gene can sit below a flashy repeat match. So the skill you need is knowing which number answers which question:

The E-value: how often you would expect this by chance

The E-value is the expected number of chance alignments with a score at least this good in a database of this size. Two properties follow directly from the definition, and they explain most confusion:

A common rule of thumb for beginners: below 1e-5 is often considered statistically meaningful, but this threshold exists because someone tested it in a specific context. It is not biology. Two proteins from different families can share short motifs that trap you at loose E-values, and two sequences from recently duplicated genes can drift apart at genuine E-values. The E-value is a statistics layer; biology enters next.

The bit score: one scale across searches

The bit score rescales the raw alignment score (match, mismatch, gap penalties) into log-odds units, roughly meaning twice as many bits = four times as unlikely by chance. Practically, three facts make it useful:

When you write up BLAST results, give the bit score and the E-value together: they describe the same alignment quality from two directions (a purely statistical one and a normalized strength one).

Percent identity: how similar the matched residues are

Percent identity is the fraction of aligned residue positions that match exactly, usually reported over the aligned region. It answers how closely related the aligned slice is — but only for that slice, which is why it pairs with coverage:

A mismatch column and indel counts also matter, because transitions in a conserved domain can change function even when identity stays high. This is why a confident annotation often goes past the table into looking at the alignment itself.

Query coverage: how much of your sequence the hit explains

Query coverage is the fraction of your input sequence that the alignment spans. It matters because a hit can be born statistical and still be biologically wrong-scale:

Coverage pairs with a second metric very similar to it: the alignment length relative to the query length. When a hit’s coverage is low but E-value is strong, ask whether your query contains a region that genuinely matches something common in the database.

Reading them together: a worked example

HitQuery CoverE valueIdentityWhat to think
A68%1e-4599%A confident partial match: most of the gene aligns extremely well; the uncovered third may be absent from the subject or a weakly conserved tail
B97%3e-8092%A strong candidate homolog across nearly the whole query — usually the annotation-grade hit
C12%1e-60100%Short perfect match: likely a repeat, primer region, or small shared domain; careful before calling orthology
D92%0.834%It spans the query but sits at noise-level E-value — treat as no confident homology without supporting evidence

None of the four rows above is automatically good or bad until you state the question. For primer design you may want row type three. For annotating a transcriptome, type two is the workhorse. That is what interpreting means: metric combined with intent.

Practical thresholds without pretending they are laws

Rules that work better than any single cutoff:

BLAST itself also warns that short queries need adjusted word sizes and statistics. If your sequence is short, generic defaults distort everything above, so this applies most to longer sequences analyzed with standard parameters.

Common mistakes that cost hours

# What to record when you summarize a BLAST hit (pseudo-notebook cell)
record = {
  'accession': 'NP_000508.2',   # so anyone can reproduce
  'identity': 93.0,             # percent over aligned region
  'coverage': 96.0,             # percent of MY query
  'evalue': 2e-79,              # database-dependent chance count
  'bit_score': 412,             # normalized alignment strength
}

Key takeaways

Challenge

Run any sequence through NCBI BLAST, deliberately pick the fourth-ranked hit, and state for it: bit score, E-value, identity, coverage, and one sentence about which biological explanation (ortholog, shared domain, repeat, or merely weak similarity) fits those four numbers. If you cannot, pick the metrics you did not consider and re-read the hit before deciding.

Want help interpreting results in a course or research project? Ampersand Academy offers one-to-one mentoring in bioinformatics and data analysis — part of Mahadhi, a development and digital marketing company.

—

FAQ

Coverage, identity, and significance answer different questions; no single metric decides. In practice the combination of a strong bit score, a very low E-value, and coverage high enough to match the region of interest is what most people mean by a good hit.

The often-quoted cutoff near 0.05 is a convention for treatments, not a biological threshold. Interpretation depends on database size, query length, and whether the aligned region has biological sense.

The E-value depends on database size and search parameters. A hit against a small trial set will show different values from the giant NCBI nt database even with the identical alignment.

Usually not, or only partially. A perfect but tiny alignment often marks a repeat, a primer region, or a shared domain; check the matched region before concluding.

Both describe the surprise level, but the bit score is normalized and comparable across database sizes better than the E-value, and it is what determines the ranking order.

Ranking is by score, not by biology. The top hit can be a fragmented database entry or a different isoform; read descriptions and relationships, not just the first row.

Identity alone hides where differences cluster. A conserved protein can hide one mismatch inside an active-site motif, so inspect the alignment view before making a functional claim.

Which metric matters most in a BLAST result?

Coverage, identity, and significance answer different questions; no single metric decides. In practice the combination of a strong bit score, a very low E-value, and coverage high enough to match the region of interest is what most people mean by a good hit.

Is an E-value below 0.05 always meaningful?

The often-quoted cutoff near 0.05 is a convention for treatments, not a biological threshold. Interpretation depends on database size, query length, and whether the aligned region has biological sense.

Why do E-values change between two searches of the same hit?

The E-value depends on database size and search parameters. A hit against a small trial set will show different values from the giant NCBI nt database even with the identical alignment.

A hit has 100 percent identity but low query coverage. Is it my gene?

Usually not, or only partially. A perfect but tiny alignment often marks a repeat, a primer region, or a shared domain; check the matched region before concluding.

What does the bit score add beyond the E-value?

Both describe the surprise level, but the bit score is normalized and comparable across database sizes better than the E-value, and it is what determines the ranking order.

Should I always pick the top-ranked hit?

Ranking is by score, not by biology. The top hit can be a fragmented database entry or a different isoform; read descriptions and relationships, not just the first row.

How do mismatches and gaps affect interpretation?

Identity alone hides where differences cluster. A conserved protein can hide one mismatch inside an active-site motif, so inspect the alignment view before making a functional claim.

Exit mobile version