Alignment & files

Pairwise Sequence Alignment

Align two DNA or protein sequences with Needleman–Wunsch (global) or Smith–Waterman (local) and affine gap penalties, with identity, similarity and gap statistics. Runs in your browser.

Pairwise Sequence AlignmentRuns in your browser; nothing is uploaded

Global (Needleman–Wunsch) alignment of protein sequences, BLOSUM62, gap open 11, extend 1.

Score286
Identity43.6%
Similarity60.4%
Gaps9
Length149
sp|P69905| 1 MV-LSPADKTNVKAAWGKVGAHAGEYGAEALERMFLSFPTTKTYFPHF-DLS-----HGS 53 || |:|.:|:.|.|.|||| :..|.|.|||.|:.:.:|.|:.:|..| ||| .|: sp|P68871| 1 MVHLTPEEKSAVTALWGKV--NVDEVGGEALGRLLVVYPWTQRFFESFGDLSTPDAVMGN 58 sp|P69905| 54 AQVKGHGKKVADALTNAVAHVDDMPNALSALSDLHAHKLRVDPVNFKLLSHCLLVTLAAH 113 .:||.|||||..|.::.:||:|::....:.||:||..||.|||.||:||.:.|:..||.| sp|P68871| 59 PKVKAHGKKVLGAFSDGLAHLDNLKGTFATLSELHCDKLHVDPENFRLLGNVLVCVLAHH 118 sp|P69905| 114 LPAEFTPAVHASLDKFLASVSTVLTSKYR 142 ...||||.|.|:..|.:|.|:..|..||. sp|P68871| 119 FGKEFTPPVQAAYQKVVAGVANALAHKYH 147

Global and local alignment

A global alignment (Needleman–Wunsch) aligns two sequences from end to end and suits sequences of similar length that are related along their whole length. A local alignment (Smith–Waterman) finds the best-scoring region shared by the two sequences and ignores the rest, which suits a domain inside a longer protein or a read against a longer reference.

Both are computed exactly by dynamic programming. This tool uses affine gap penalties, in which opening a gap costs more than extending it, so a single long indel is preferred to many short ones, as in real sequence evolution.

Scoring

DNA is scored with a match reward of +2 and a mismatch penalty of −3, with a gap-open penalty of 5 and a gap-extension penalty of 2. Proteins are scored with the BLOSUM62 substitution matrix, a gap-open penalty of 11 and an extension penalty of 1, similar to the BLASTP defaults. A gap of length L costs open + extend × (L − 1).

Reading the result

In the middle line, | marks an identical position, : a positive BLOSUM62 score (a conservative substitution in proteins), . a mismatch and a space a gap. Identity is the fraction of aligned columns with identical residues; the alignment score depends on the scoring scheme and is only comparable between alignments scored the same way.

Learn the method properly. Bioinformatics Foundations: Databases, file formats, alignment and the command line: the working toolkit of every bioinformatician. View the course.

FAQ

Frequently asked questions

Not answered here? Send us your question and we will reply by e-mail.

Up to a few thousand residues each: the alignment matrix grows with the product of the two lengths, and the browser limit is about 25 million cells. For database searches use BLAST.

No. BLAST is a heuristic that searches a database quickly; this tool computes the exact optimal alignment of two sequences you supply.

No. The calculation runs in your browser with JavaScript. Nothing you paste is sent to a server, which makes the tool safe for unpublished or confidential sequences.