Sequence analysis
Pairwise sequence alignment
Compare two sequences directly, preserve the aligned strings and scoring settings, and separate correspondence from similarity scores.
Inputs
1 required
Methods
2 connected
- 01MAFFT · Pairwise Alignment
- 02StringZilla v5 · Global Score
MAFFT returns aligned sequences while StringZilla independently returns a Needleman–Wunsch global score; the workflow keeps those outputs distinct.
Use this templateWhat is pairwise sequence alignment?
Pairwise sequence alignment is the process of arranging two biological sequences so corresponding residues or nucleotides appear in shared positions. An explicit substitution and gap model determines the alignment path and score. Those outputs are related but not interchangeable: a score ranks similarity under one parameterization, while the aligned strings expose matches, substitutions, insertions, deletions, and terminal gaps.
Use pairwise alignment when the scientific question concerns one sequence versus one other sequence: two alleles, an engineered construct and its parent, two domains, or a query and a known homolog. Choose global or local behavior according to whether the full lengths or only the best-matching region should correspond.
Check molecular type, orientation, sequence boundaries, gap penalties, and substitution model before interpreting percent identity. The same pair can receive different alignments and scores under different settings, especially around repeats, low-complexity regions, long insertions, and weakly conserved termini.
When to use pairwise sequence alignment
- Best fit. Two homologs, alleles, constructs, or domains
- Required input. Exactly two compatible FASTA sequences and a global-versus-local decision
Benefits of pairwise sequence alignment
- Clear correspondence. Makes residue correspondence explicit
- Connected evidence. Supports direct construct comparison
- Reusable output. Produces compact reviewable outputs
Primary limitations
- Method dependence. Results depend on scoring choices
- Input dependence. Repeats can create alternative alignments
- Interpretive limit. Similarity does not prove shared function
Pairwise sequence alignment methods
Dynamic-programming methods optimize a defined objective across a score matrix. End-to-end algorithms penalize terminal differences, while local algorithms allow the comparison to begin and end inside the sequences.
Progressive multiple-alignment programs can align a two-record FASTA, but their output should not be described as an exact Needleman–Wunsch path unless that is the algorithm actually run. Preserve method names and settings with every result.
Pairwise sequence alignment applications
Pairwise sequence alignment is best suited to two homologs, alleles, constructs, or domains. The result can support comparative review, sequence curation, annotation, profile construction, phylogenetic preparation, structural interpretation, or experimental planning when those downstream uses match the alignment scope.
Keep the alignment as evidence rather than a conclusion. Downstream claims should remain tied to sequence provenance, coverage, method agreement, relevant biological context, and any independent structural, evolutionary, or experimental support.
How to run pairwise sequence alignment online
Use the connected workflow to keep input records, method settings, native outputs, warnings, and exports together. Review every stage before using the result for annotation, phylogeny, variant interpretation, or experimental decisions.
- Define scope. Decide whether the full sequences or only a shared region should be compared.
- Prepare the pair. Verify molecule type, orientation, headers, boundaries, and ambiguous symbols.
- Set scoring. Record the substitution and gap settings rather than relying on an unlabeled score.
- Run comparison. Generate aligned strings and, when useful, an independent pairwise score.
- Interpret and export. Inspect identity, coverage, gaps, terminal effects, and biologically important positions.
How to interpret pairwise sequence alignment results
Report identity together with aligned length, coverage, gap count, and the treatment of ambiguous residues. Identity calculated over the aligned region answers a different question from identity divided by either full sequence length.
Inspect substitutions at known motifs or functional sites directly. A high overall identity can hide a critical local change, while a modest global identity can coexist with a strongly conserved domain.
How pairwise sequence alignment works
MAFFT returns aligned sequences while StringZilla independently returns a Needleman–Wunsch global score; the workflow keeps those outputs distinct.
- Define scope. Decide whether the full sequences or only a shared region should be compared.
- Prepare the pair. Verify molecule type, orientation, headers, boundaries, and ambiguous symbols.
- Set scoring. Record the substitution and gap settings rather than relying on an unlabeled score.
- Run comparison. Generate aligned strings and, when useful, an independent pairwise score.
- Interpret and export. Inspect identity, coverage, gaps, terminal effects, and biologically important positions.
Inputs and outputs
Check formats before running, then inspect and download the result from every workflow step.
Inputs
- Alignment input.
FASTAPDBmmCIFTwo protein, DNA, or RNA sequences in FASTA format.
Outputs
- Alignment outputs.
FASTACSVTSVPDBJSONAligned FASTA plus method-native score tables and downloadable settings.
Tools for pairwise sequence alignment
Use these methods to prepare inputs, run the core analysis, inspect outputs, and validate the evidence described in this workflow.
MAFFT
Create configurable protein, DNA, or RNA alignments
StringZilla v5
Calculate pairwise global, local, or edit-distance score matrices
Clustal Omega
Create scalable protein or nucleotide multiple-sequence alignments
MUSCLE5
Generate conventional or ensemble multiple-sequence alignments
HMMER
Find homologs with profile hidden Markov models
MMseqs2
Search and cluster large protein or nucleotide sequence sets
USAlign
Align macromolecular structures and derive residue correspondence
FoldSeek
Find and compare structurally similar proteins
PDB to FASTA converter
Extract sequences from structures before comparison
FastTree
Estimate trees from large sequence alignments
IQ-TREE
Infer maximum-likelihood phylogenies from alignments
RAxML-NG
Run maximum-likelihood phylogenetic analysis
Other sequence analysis workflows
Compare related approaches based on the molecular system, available evidence, required inputs, and decision you need to support.
Multiple sequence alignment
Aligns three or more homologous sequences to identify shared positions, insertions, deletions, and conserved regions.
Global sequence alignment
Compares sequences end to end, including terminal differences and gaps across their full lengths.
Protein sequence alignment
Aligns amino-acid sequences using substitution-aware methods suited to protein evolution and function.
DNA sequence alignment
Aligns nucleotide sequences to compare homologous genes, amplicons, loci, transcripts, or constructs.
Local sequence alignment
Finds or scores the best-matching subsequences without forcing unrelated flanks into the comparison.
Whole genome alignment
Maps large homologous regions between genome assemblies and reports coordinates, rearrangements, and sequence differences.
Structure-based sequence alignment
Uses three-dimensional correspondence to align residues whose sequence similarity alone may be weak.
Frequently asked questions
Two protein, DNA, or RNA sequences in FASTA format.
Aligned FASTA plus method-native score tables and downloadable settings.
Start from the scientific scope: global or local, pairwise or multiple, sequence or structure, and conventional or genome scale. Then record the method, substitution model, gap settings, sequence type, and any filtering rather than relying on defaults without provenance.
No. Scores and identities quantify similarity under a defined model. Homology is an evolutionary interpretation, and shared function requires additional evidence such as domain context, conserved residues, structure, phylogeny, experiments, or curated annotation.
Preserve both raw sequences, algorithm, substitution and gap settings, aligned strings, and score definition.
A complete pairwise sequence alignment project is usually quote-based because providers scope sequence curation, method selection, alignment review, interpretation, and downstream analysis together. Harvard’s FY26 bioinformatics core first defines deliverables and a time estimate, then charges $180–$265 per hour; MSU lists $84–$110 per hour and expects at least eight consultant hours for custom analysis.
The total depends on sequence count and length, input cleanup, molecular type, the number of methods compared, manual review, genome scale, figures, phylogenetic or structural follow-up, and whether the deliverable includes interpretation or only alignment files.
ProteinIQ self-service starts at $29 per month for academic Plus and $99 per month for commercial Pro, with the configured pairwise sequence alignment run estimated in credits before submission. Done-for-you analysis is scoped separately and can include data preparation, method comparison, interpretation, and a reproducible handoff.
Start with a workflow you can inspect and edit
Add your inputs, review the settings, and keep every structure, score, table, and file connected to the step that produced it.