Open main menu

Biolecture.org β

Changes

Alignment

102,868 bytes added, 03:04, 2 December 2018
no edit summary
<p><span style="font-size:24px"><u><strongh1>Sequence Alignment</strong></u></span>alignment</ph1>
<p>&nbsp;From Wikipedia, the free encyclopedia</p>
<p>&nbsp;<a href="https://en.wikipedia.org/wiki/Sequence_alignment#mw-head">Jump to navigation</a><a href="https://en.wikipedia.org/wiki/Sequence_alignment#p-search">Jump to search</a></p>
<table> <tbody> <tr> <td> <p><span a href="https://en.wikipedia.org/wiki/File:Question_book-new.svg"><img alt="" src="https://upload.wikimedia.org/wikipedia/en/thumb/9/99/Question_book-new.svg/50px-Question_book-new.svg.png" style="font-sizeheight:39px; width:50px" /></a></p> </td> <td> <p>This article&nbsp;<strong>needs additional citations for&nbsp;<a href="https://en.wikipedia.org/wiki/Wikipedia:Verifiability" title="Wikipedia:20pxVerifiability">verification</a></strong>What is sequence alignment.&nbsp;Please help&nbsp;<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit">improve this article</a>&nbsp;by&nbsp;<a href="https://en.wikipedia.org/wiki/Help:Introduction_to_referencing_with_Wiki_Markup/1" title="Help:Introduction to referencing with Wiki Markup/1">adding citations to reliable sources</a>. Unsourced material may be challenged and removed.&nbsp;<small><em>(March 2009)</strongem></small><small><em>&nbsp;(<a href="https://en.wikipedia.org/wiki/Help:Maintenance_template_removal" title="Help:Maintenance template removal">Learn how and when to remove this template message</a>)</em></spansmall></p> </td> </tr> </tbody></table>
<p>In&nbsp;<a href="https://en.wikipedia.org/wiki/Bioinformatics" title="Bioinformatics">bioinformatics</a>, a&nbsp;<strong>sequence alignment</strong>&nbsp;is a way of arranging the sequences of&nbsp;<a href="https://en.wikipedia.org/wiki/DNA" title="DNA">DNA</a>,&nbsp;<a href="https://en.wikipedia.org/wiki/RNA" title="RNA">RNA</a>, or protein to identify regions of similarity that may be a consequence of functional,&nbsp;<a href="https://en.wikipedia.org/wiki/Structural_biology" title="Structural biology">structural</a>, or&nbsp;<a href="https://en.wikipedia.org/wiki/Evolution" title="Evolution">evolutionary</a>&nbsp;relationships between the sequences.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-mount-1">[1]</a></sup>&nbsp;Aligned sequences of&nbsp;<a href="https://en.wikipedia.org/wiki/Nucleotide" title="Nucleotide">nucleotide</a>&nbsp;or&nbsp;<a href="https://en.wikipedia.org/wiki/Amino_acid" title="Amino acid">amino acid</a>&nbsp;residues are typically represented as rows within a&nbsp;<a href="https://en.wikipedia.org/wiki/Matrix_(mathematics)" title="Matrix (mathematics)">matrix</a>. Gaps are inserted between the&nbsp;<a href="https://en.wikipedia.org/wiki/Residue_(chemistry)" title="Residue (chemistry)">residues</a>&nbsp;so that identical or similar characters are aligned in successive columns. Sequence alignments are also used for non-biological sequences, such as calculating the&nbsp;<a href="https://en.wikipedia.org/wiki/Edit_distance" title="Edit distance">edit distance cost</a>&nbsp;between strings in a&nbsp;<a href="https://en.wikipedia.org/wiki/Natural_language" title="Natural language">natural language</a>&nbsp;or in financial data.</p>
<p><span a href="https://en.wikipedia.org/wiki/File:Histone_Alignment.png"><img alt="" src="https://upload.wikimedia.org/wikipedia/commons/thumb/b/b5/Histone_Alignment.png/595px-Histone_Alignment.png" style="font-sizeheight:205px; width:14px595px"/>A sequence alignment is a way of arranging the sequences of DNA, RNA, or protein to identify regions of similarity that may be a consequence of functional, structural, or evolutionary relationships between the sequences.</spana></p>
<p>A sequence alignment, produced by&nbsp;<a href="https://en.wikipedia.org/wiki/ClustalO" title="ClustalO">ClustalO</a>, of mammalian&nbsp;<a href="https://en.wikipedia.org/wiki/Histone" title="Histone">histone</a>&nbsp;proteins.&nbsp;<br />Sequences are the&nbsp;<a href="https://en.wikipedia.org/wiki/Amino_acid#Table_of_standard_amino_acid_abbreviations_and_side_chain_properties" title="Amino acid">amino acids</a>&nbsp;for residues 120-180 of the proteins. Residues that are conserved across all sequences are highlighted in grey. Below the protein sequences is a key denoting&nbsp;<a href="https://en.wikipedia.org/wiki/Conserved_sequence" title="Conserved sequence">conserved sequence</a>&nbsp;(*),&nbsp;<a href="https://en.wikipedia.org/wiki/Conservative_mutation" title="Conservative mutation">conservative mutations</a>&nbsp;(:), semi-conservative mutations (.), and&nbsp;<a href="https://en.wikipedia.org/wiki/Segregating_site" title="Segregating site">non-conservative mutations</a>&nbsp;( ).<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-2">[2]</a></sup></p>
<p><span styleinput type="font-size:20pxcheckbox"><strong>How sequence is aligned?</strong></span></p>
<ph2>&nbsp;Contents</ph2>
<pul> <span styleli><a href="font-sizehttps://en.wikipedia.org/wiki/Sequence_alignment#Interpretation">1Interpretation</a></li> <li><a href="https:14px//en.wikipedia.org/wiki/Sequence_alignment#Alignment_methods">Human knowledge is applied in constructing algorithms to produce high-quality sequence alignments, and occasionally in adjusting the final results to reflect patterns that are difficult to represent algorithmically (especially in the case of nucleotide sequences)2Alignment methods</a></li> <li><a href="https://en.wikipedia. Computational approaches to sequence alignment generally fall into two categoriesorg/wiki/Sequence_alignment#Representations">3Representations</a></li> <li><a href="https: global alignments //en.wikipedia.org/wiki/Sequence_alignment#Global_and_local_alignments">4Global and local alignments</a></li> <li><a href="https://en.wikipedia.org/wiki/Sequence_alignment#Pairwise_alignment">5Pairwise alignment</a> <ul> <li><a href="https://en.wikipedia.org/wiki/Sequence_alignment#Dot-matrix_methods">5.1Dot-matrix methods</a></li> <li><a href="https://en.wikipedia.org/wiki/Sequence_alignment#Dynamic_programming">5. Calculating 2Dynamic programming</a global alignment is ></li> <li><a form of global optimization that &quot;forces&quot; the alignment to span the entire length of all query sequenceshref="https://en.wikipedia.org/wiki/Sequence_alignment#Word_methods">5. By contrast, local alignments identify regions of similarity within long sequences that are often widely divergent overall3Word methods</a></li> </ul> </li> <li><a href="https://en. Local alignments are often preferable, but can be more difficult to calculate because of the additional challenge of identifying the regions of similaritywikipedia. A variety of computational algorithms have been applied to the org/wiki/Sequence_alignment#Multiple_sequence_alignment">6Multiple sequence alignment problem</a> <ul> <li><a href="https://en.wikipedia.org/wiki/Sequence_alignment#Dynamic_programming_2">6. These include slow but formally correct methods like dynamic 1Dynamic programming</a></li> <li><a href="https://en. These also include efficient, heuristic algorithms or probabilistic wikipedia.org/wiki/Sequence_alignment#Progressive_methods">6.2Progressive methods designed for large-scale database searching, that do not guarantee to find best matches</a></li> <li><a href="https://en.wikipedia.org/wiki/Sequence_alignment#Iterative_methods">6. Hybrid 3Iterative methods, known as semi-global or &quot;glocal&quot; (short for global-local) methods, attempt to find the best possible </a></li> <li><a href="https://en.wikipedia.org/wiki/Sequence_alignment#Motif_finding">6.4Motif finding</a></li> <li><a href="https://en.wikipedia.org/wiki/Sequence_alignment#Techniques_inspired_by_computer_science">6.5Techniques inspired by computer science</a></li> </ul> </li> <li><a href="https://en.wikipedia.org/wiki/Sequence_alignment#Structural_alignment">7Structural alignment that includes the start and end </a> <ul> <li><a href="https://en.wikipedia.org/wiki/Sequence_alignment#DALI">7.1DALI</a></li> <li><a href="https://en.wikipedia.org/wiki/Sequence_alignment#SSAP">7.2SSAP</a></li> <li><a href="https://en.wikipedia.org/wiki/Sequence_alignment#Combinatorial_extension">7.3Combinatorial extension</a></li> </ul> </li> <li><a href="https://en.wikipedia.org/wiki/Sequence_alignment#Phylogenetic_analysis">8Phylogenetic analysis</a> <ul> <li><a href="https://en.wikipedia.org/wiki/Sequence_alignment#Assessment_of_significance">8.1Assessment of one or the other sequencesignificance</a></li> <li><a href="https://en.wikipedia.org/wiki/Sequence_alignment#Assessment_of_credibility">8. This can be especially useful when the downstream part 2Assessment of one sequence overlaps with the upstream part of the other sequencecredibility</a></li> <li><a href="https://en.wikipedia.org/wiki/Sequence_alignment#Scoring_functions">8.3Scoring functions</a></li> </ul> </li> <li><a href="https://en.wikipedia.org/wiki/Sequence_alignment#Other_biological_uses">9Other biological uses</a></li> <li><a href="https://en.wikipedia.org/wiki/Sequence_alignment#Non-biological_uses">10Non-biological uses</a></li> <li><a href="https://en.wikipedia.org/wiki/Sequence_alignment#Software">11Software</a></li> <li><a href="https://en.wikipedia. In this case, neither global nor local alignment is entirely appropriateorg/wiki/Sequence_alignment#See_also">12See also</a></li> <li><a href="https: //en.wikipedia.org/wiki/Sequence_alignment#References">13References</a global alignment would attempt to force the alignment to extend beyond the region of overlap, while ></li> <li><a local alignment might not fully cover the region of overlaphref="https://en.wikipedia.org/wiki/Sequence_alignment#External_links">14External links</a></spanli></pul>
<ph2>Interpretation[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&nbspamp;action=edit&amp;section=1" title="Edit section: Interpretation">edit</a>]</ph2>
<p>If two sequences in an alignment share a common ancestor, mismatches can be interpreted as&nbsp;<a href="https://en.wikipedia.org/wiki/Point_mutation" title="Point mutation">point mutations</a>&nbsp;and gaps as&nbsp;<span stylea href="https://en.wikipedia.org/wiki/Indel" title="font-sizeIndel">indels</a>&nbsp;(that is, insertion or deletion mutations) introduced in one or both lineages in the time since they diverged from one another. In sequence alignments of proteins, the degree of similarity between&nbsp;<a href="https:20px//en.wikipedia.org/wiki/Amino_acid" title="Amino acid">amino acids<strong/a>How &nbsp;occupying a particular position in the sequence alignment can be used?interpreted as a rough measure of how&nbsp;<a href="https:/strong/en.wikipedia.org/wiki/Conservation_(genetics)" title="Conservation (genetics)">conserved</a>&nbsp;a particular region or&nbsp;<a href="https://en.wikipedia.org/wiki/Sequence_motif" title="Sequence motif">sequence motif</a>&nbsp;is among lineages. The absence of substitutions, or the presence of only very conservative substitutions (that is, the substitution of amino acids whose&nbsp;<a href="https://en.wikipedia.org/wiki/Side_chain" title="Side chain">side chains</a>&nbsp;have similar biochemical properties) in a particular region of the sequence, suggest&nbsp;<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-predict-3">[3]</a></sup>&nbsp;that this region has structural or functional importance. Although DNA and RNA&nbsp;<a href="https://en.wikipedia.org/wiki/Nucleotide" title="Nucleotide">nucleotide</spana>&nbsp;bases are more similar to each other than are amino acids, the conservation of base pairs can indicate a similar functional or structural role.</p>
<ph2>Alignment methods[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&nbspamp;action=edit&amp;section=2" title="Edit section: Alignment methods">edit</a>]</ph2>
<p>Very short or very similar sequences can be aligned by hand. However, most interesting problems require the alignment of lengthy, highly variable or extremely numerous sequences that cannot be aligned solely by human effort. Instead, human knowledge is applied in constructing algorithms to produce high-quality sequence alignments, and occasionally in adjusting the final results to reflect patterns that are difficult to represent algorithmically (especially in the case of nucleotide sequences). Computational approaches to sequence alignment generally fall into two categories:&nbsp;<em>global alignments</em>&nbsp;and&nbsp;<em>local alignments</em>. Calculating a global alignment is a form of&nbsp;<span stylea href="font-sizehttps:14px//en.wikipedia.org/wiki/Global_optimization" title="Global optimization">Phylogenetics and sequence global optimization</a>&nbsp;that &quot;forces&quot; the alignment are closely related fields due to span the shared necessity entire length of evaluating sequence relatednessall query sequences. By contrast, local alignments identify regions of similarity within long sequences that are often widely divergent overall. Sequence alignment Local alignments are often preferable, but can be used for construction and interpretation more difficult to calculate because of phylogenetic trees, which are used to classify the evolutionary relationships between homologous genes represented in additional challenge of identifying the genomes regions of divergent speciessimilarity.<sup><a href="https://en.wikipedia. Roughly speaking, high org/wiki/Sequence_alignment#cite_note-Polyanovsky2011-4">[4]</a></sup>&nbsp;A variety of computational algorithms have been applied to the sequence identity suggests that the sequences in question have alignment problem. These include slow but formally correct methods like&nbsp;<a href="https://en.wikipedia.org/wiki/Dynamic_programming" title="Dynamic programming">dynamic programming</a comparatively young most recent common ancestor>. These also include efficient, while low identity suggests that the divergence is more ancient&nbsp;<a href="https://en.wikipedia.org/wiki/Heuristic_algorithm" title="Heuristic algorithm">heuristic algorithms</a>&nbsp;or&nbsp;<a href="https://en.wikipedia.org/wiki/Probability" title="Probability">probabilistic</spana>&nbsp;methods designed for large-scale database search, that do not guarantee to find best matches.</p>
<ph2>Representations[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&nbspamp;action=edit&amp;section=3" title="Edit section: Representations">edit</a>]</ph2>
<p>Alignments are commonly represented both graphically and in text format. In almost all sequence alignment representations, sequences are written in rows arranged so that aligned residues appear in successive columns. In text formats, aligned columns containing identical or similar characters are indicated with a system of conservation symbols. As in the image above, an asterisk or pipe symbol is used to show identity between two columns; other less common symbols include a colon for conservative substitutions and a period for semiconservative substitutions. Many sequence visualization programs also use color to display information about the properties of the individual sequence elements; in DNA and RNA sequences, this equates to assigning each nucleotide its own color. In protein alignments, such as the one in the image above, color is often used to indicate amino acid properties to aid in judging the&nbsp;<span stylea href="https://en.wikipedia.org/wiki/Conservation_(genetics)" title="Conservation (genetics)">conservation</a>&nbsp;of a given amino acid substitution. For multiple sequences the last row in each column is often the&nbsp;<a href="https://en.wikipedia.org/wiki/Consensus_sequence" title="Consensus sequence">consensus sequence</a>&nbsp;determined by the alignment; the consensus sequence is also often represented in graphical format with a&nbsp;<a href="https://en.wikipedia.org/wiki/Sequence_logo" title="Sequence logo">sequence logo</a>&nbsp;in which the size of each nucleotide or amino acid letter corresponds to its degree of conservation.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-Schneider-5">[5]</a></sup></p> <p>Sequence alignments can be stored in a wide variety of text-based file formats, many of which were originally developed in conjunction with a specific alignment program or implementation. Most web-based tools allow a limited number of input and output formats, such as&nbsp;<a href="https://en.wikipedia.org/wiki/FASTA_format" title="FASTA format">FASTA format</a>&nbsp;and&nbsp;<a href="https://en.wikipedia.org/wiki/GenBank" title="GenBank">GenBank</a>&nbsp;format and the output is not easily editable. Several conversion programs that provide graphical and/or command line interfaces are available<sup>[<em><a href="https://en.wikipedia.org/wiki/Wikipedia:Link_rot" title="Wikipedia:Link rot">dead link</a></em>]</sup>, such as&nbsp;<a href="https://web.archive.org/web/20071024223546/http://bioweb.pasteur.fr/seqanal/interfaces/readseq.html" rel="nofollow">READSEQ</a>&nbsp;and&nbsp;<a href="https://en.wikipedia.org/wiki/EMBOSS" title="EMBOSS">EMBOSS</a>. There are also several programming packages which provide this conversion functionality, such as&nbsp;<a href="https://en.wikipedia.org/wiki/BioPython" title="BioPython">BioPython</a>,&nbsp;<a href="https://en.wikipedia.org/wiki/BioRuby" title="BioRuby">BioRuby</a>&nbsp;and&nbsp;<a href="https://en.wikipedia.org/wiki/BioPerl" title="BioPerl">BioPerl</a>. The&nbsp;<a href="https://en.wikipedia.org/wiki/SAM_(file_format)" title="SAM (file format)">SAM/BAM files</a>&nbsp;use the CIGAR (Compact Idiosyncratic Gapped Alignment Report) string format to represent an alignment of a sequence to a reference by encoding a sequence of events (e.g. match/mismatch, insertions, deletions).<sup><a href="fonthttps://en.wikipedia.org/wiki/Sequence_alignment#cite_note-6">[6]</a></sup></p> <h2>Global and local alignments[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=4" title="Edit section: Global and local alignments">edit</a>]</h2> <p>Global alignments, which attempt to align every residue in every sequence, are most useful when the sequences in the query set are similar and of roughly equal size. (This does not mean global alignments cannot start and/or end in gaps.) A general global alignment technique is the&nbsp;<a href="https:20px//en.wikipedia.org/wiki/Needleman%E2%80%93Wunsch_algorithm" title="Needleman–Wunsch algorithm">Needleman&ndash;Wunsch algorithm</a>, which is based on dynamic programming. Local alignments are more useful for dissimilar sequences that are suspected to contain regions of similarity or similar sequence motifs within their larger sequence context. The&nbsp;<a href="https://en.wikipedia.org/wiki/Smith%E2%80%93Waterman_algorithm" title="Smith–Waterman algorithm">Smith&ndash;Waterman algorithm</a>&nbsp;is a general local alignment method based on the same dynamic programming scheme but with additional choices to start and end at any place.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-Polyanovsky2011-4">[4]</a></sup></p> <p>Hybrid methods, known as semi-global or &quot;glocal&quot; (short for&nbsp;<strong>Referencesglo</strong>bal-lo<strong>cal</strong>) methods, search for the best possible partial alignment of the two sequences (a subset of one or both starts and one or both ends has to be chosen before aligning). This can be especially useful when the downstream part of one sequence overlaps with the upstream part of the other sequence. In this case, neither global nor local alignment is entirely appropriate: a global alignment would attempt to force the alignment to extend beyond the region of overlap, while a local alignment might not fully cover the region of overlap.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-brudno-7">[7]</a></sup>&nbsp;Another case where semi-global alignment is useful is when one sequence is short (for example a gene sequence) and the other is very long (for example a chromosome sequence). In that case, the short sequence should be globally (fully) aligned but only a local (partial) alignment is desired for the long sequence.</p> <h2>Pairwise alignment[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=5" title="Edit section: Pairwise alignment">edit</a>]</h2> <p>Pairwise sequence alignment methods are used to find the best-matching piecewise (local or global) alignments of two query sequences. Pairwise alignments can only be used between two sequences at a time, but they are efficient to calculate and are often used for methods that do not require extreme precision (such as searching a database for sequences with high similarity to a query). The three primary methods of producing pairwise alignments are dot-matrix methods, dynamic programming, and word methods;<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-mount-1">[1]</a></sup>&nbsp;however, multiple sequence alignment techniques can also align pairs of sequences. Although each method has its individual strengths and weaknesses, all three pairwise methods have difficulty with highly repetitive sequences of low&nbsp;<a href="https://en.wikipedia.org/wiki/Information_content" title="Information content">information content</a>&nbsp;- especially where the number of repetitions differ in the two sequences to be aligned. One way of quantifying the utility of a given pairwise alignment is the &#39;maximum unique match&#39; (MUM), or the longest subsequence that occurs in both query sequences. Longer MUM sequences typically reflect closer relatedness.</p> <h3>Dot-matrix methods[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=6" title="Edit section: Dot-matrix methods">edit</a>]</h3> <table> <tbody> <tr> <td><a href="https://en.wikipedia.org/wiki/File:Mup_locus_showing_DNA_repeats.jpg"><img alt="" src="https://upload.wikimedia.org/wikipedia/commons/thumb/7/72/Mup_locus_showing_DNA_repeats.jpg/200px-Mup_locus_showing_DNA_repeats.jpg" style="height:236px; width:200px" /></a>  <p>Self comparison of a part of a mouse strain genome. The dot-plot shows a patchwork of lines, demonstrating duplicated segments of DNA.</p> </td> </tr> </tbody></table> <p>See main article on&nbsp;<a href="https://en.wikipedia.org/wiki/Dot_plot_(bioinformatics)" title="Dot plot (bioinformatics)">dot plots (bioinformatics)</a>.</p> <table> <tbody> <tr> <td><a href="https://en.wikipedia.org/wiki/File:Zinc-finger-dot-plot.png"><img alt="" src="https://upload.wikimedia.org/wikipedia/commons/thumb/3/33/Zinc-finger-dot-plot.png/200px-Zinc-finger-dot-plot.png" style="height:200px; width:200px" /></a>  <p>A DNA&nbsp;<a href="https://en.wikipedia.org/wiki/Dot_plot_(bioinformatics)" title="Dot plot (bioinformatics)">dot plot</a>&nbsp;of a&nbsp;<a href="https://en.wikipedia.org/wiki/Human" title="Human">human</a><a href="https://en.wikipedia.org/wiki/Zinc_finger" title="Zinc finger">zinc finger</a>&nbsp;<a href="https://en.wikipedia.org/wiki/Transcription_factor" title="Transcription factor">transcription factor</a>(GenBank ID NM_002383), showing regional&nbsp;<a href="https://en.wikipedia.org/wiki/Self-similarity" title="Self-similarity">self-similarity</a>. The main diagonal represents the sequence&#39;s alignment with itself; lines off the main diagonal represent similar or repetitive patterns within the sequence. This is a typical example of a&nbsp;<a href="https://en.wikipedia.org/wiki/Recurrence_plot" title="Recurrence plot">recurrence plot</a>.</p> </td> </tr> </tbody></table> <p>The dot-matrix approach, which implicitly produces a family of alignments for individual sequence regions, is qualitative and conceptually simple, though time-consuming to analyze on a large scale. In the absence of noise, it can be easy to visually identify certain sequence features&mdash;such as insertions, deletions, repeats, or&nbsp;<a href="https://en.wikipedia.org/wiki/Inverted_repeat" title="Inverted repeat">inverted repeats</a>&mdash;from a dot-matrix plot. To construct a&nbsp;<a href="https://en.wikipedia.org/wiki/Dot_plot_(bioinformatics)" title="Dot plot (bioinformatics)">dot-matrix plot</a>, the two sequences are written along the top row and leftmost column of a two-dimensional&nbsp;<a href="https://en.wikipedia.org/wiki/Matrix_(mathematics)" title="Matrix (mathematics)">matrix</a>and a dot is placed at any point where the characters in the appropriate columns match&mdash;this is a typical&nbsp;<a href="https://en.wikipedia.org/wiki/Recurrence_plot" title="Recurrence plot">recurrence plot</a>. Some implementations vary the size or intensity of the dot depending on the degree of similarity of the two characters, to accommodate conservative substitutions. The dot plots of very closely related sequences will appear as a single line along the matrix&#39;s&nbsp;<a href="https://en.wikipedia.org/wiki/Main_diagonal" title="Main diagonal">main diagonal</a>.</p> <p>Problems with dot plots as an information display technique include: noise, lack of clarity, non-intuitiveness, difficulty extracting match summary statistics and match positions on the two sequences. There is also much wasted space where the match data is inherently duplicated across the diagonal and most of the actual area of the plot is taken up by either empty space or noise, and, finally, dot-plots are limited to two sequences. None of these limitations apply to Miropeats alignment diagrams but they have their own particular flaws.</p> <p>Dot plots can also be used to assess repetitiveness in a single sequence. A sequence can be plotted against itself and regions that share significant similarities will appear as lines off the main diagonal. This effect can occur when a protein consists of multiple similar&nbsp;<a href="https://en.wikipedia.org/wiki/Structural_domain" title="Structural domain">structural domains</a>.</p> <h3>Dynamic programming[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=7" title="Edit section: Dynamic programming">edit</a>]</h3> <p>The technique of&nbsp;<a href="https://en.wikipedia.org/wiki/Dynamic_programming" title="Dynamic programming">dynamic programming</a>&nbsp;can be applied to produce global alignments via the&nbsp;<a href="https://en.wikipedia.org/wiki/Needleman-Wunsch_algorithm" title="Needleman-Wunsch algorithm">Needleman-Wunsch algorithm</a>, and local alignments via the&nbsp;<a href="https://en.wikipedia.org/wiki/Smith-Waterman_algorithm" title="Smith-Waterman algorithm">Smith-Waterman algorithm</a>. In typical usage, protein alignments use a&nbsp;<a href="https://en.wikipedia.org/wiki/Substitution_matrix" title="Substitution matrix">substitution matrix</a>&nbsp;to assign scores to amino-acid matches or mismatches, and a&nbsp;<a href="https://en.wikipedia.org/wiki/Gap_penalty" title="Gap penalty">gap penalty</a>&nbsp;for matching an amino acid in one sequence to a gap in the other. DNA and RNA alignments may use a scoring matrix, but in practice often simply assign a positive match score, a negative mismatch score, and a negative gap penalty. (In standard dynamic programming, the score of each amino acid position is independent of the identity of its neighbors, and therefore&nbsp;<a href="https://en.wikipedia.org/wiki/Base_stacking" title="Base stacking">base stacking</a>&nbsp;effects are not taken into account. However, it is possible to account for such effects by modifying the algorithm.) A common extension to standard linear gap costs, is the usage of two different gap penalties for opening a gap and for extending a gap. Typically the former is much larger than the latter, e.g. -10 for gap open and -2 for gap extension. Thus, the number of gaps in an alignment is usually reduced and residues and gaps are kept together, which typically makes more biological sense. The Gotoh algorithm implements affine gap costs by using three matrices.</p> <p>Dynamic programming can be useful in aligning nucleotide to protein sequences, a task complicated by the need to take into account&nbsp;<a href="https://en.wikipedia.org/wiki/Frameshift" title="Frameshift">frameshift</a>&nbsp;mutations (usually insertions or deletions). The framesearch method produces a series of global or local pairwise alignments between a query nucleotide sequence and a search set of protein sequences, or vice versa. Its ability to evaluate frameshifts offset by an arbitrary number of nucleotides makes the method useful for sequences containing large numbers of indels, which can be very difficult to align with more efficient heuristic methods. In practice, the method requires large amounts of computing power or a system whose architecture is specialized for dynamic programming. The&nbsp;<a href="https://en.wikipedia.org/wiki/BLAST" title="BLAST">BLAST</a>&nbsp;and&nbsp;<a href="https://en.wikipedia.org/wiki/EMBOSS" title="EMBOSS">EMBOSS</a>&nbsp;suites provide basic tools for creating translated alignments (though some of these approaches take advantage of side-effects of sequence searching capabilities of the tools). More general methods are available from both commercial sources, such as&nbsp;<em>FrameSearch</em>, distributed as part of the&nbsp;<a href="https://en.wikipedia.org/wiki/Accelrys" title="Accelrys">Accelrys</a>&nbsp;<a href="https://en.wikipedia.org/wiki/GCG_(software)" title="GCG (software)">GCG package</a>, and&nbsp;<a href="https://en.wikipedia.org/wiki/Open_Source" title="Open Source">Open Source</a>&nbsp;software such as&nbsp;<a href="http://www.ebi.ac.uk/Tools/psa/genewise/" rel="nofollow">Genewise</a>.</p> <p>The dynamic programming method is guaranteed to find an optimal alignment given a particular scoring function; however, identifying a good scoring function is often an empirical rather than a theoretical matter. Although dynamic programming is extensible to more than two sequences, it is prohibitively slow for large numbers of sequences or extremely long sequences.</p> <h3>Word methods[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=8" title="Edit section: Word methods">edit</a>]</h3> <p>Word methods, also known as&nbsp;<em>k</em>-tuple methods, are&nbsp;<a href="https://en.wikipedia.org/wiki/Heuristic" title="Heuristic">heuristic</a>&nbsp;methods that are not guaranteed to find an optimal alignment solution, but are significantly more efficient than dynamic programming. These methods are especially useful in large-scale database searches where it is understood that a large proportion of the candidate sequences will have essentially no significant match with the query sequence. Word methods are best known for their implementation in the database search tools&nbsp;<a href="https://en.wikipedia.org/wiki/FASTA" title="FASTA">FASTA</a>&nbsp;and the&nbsp;<a href="https://en.wikipedia.org/wiki/BLAST" title="BLAST">BLAST</a>&nbsp;family.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-mount-1">[1]</a></sup>&nbsp;Word methods identify a series of short, nonoverlapping subsequences (&quot;words&quot;) in the query sequence that are then matched to candidate database sequences. The relative positions of the word in the two sequences being compared are subtracted to obtain an offset; this will indicate a region of alignment if multiple distinct words produce the same offset. Only if this region is detected do these methods apply more sensitive alignment criteria; thus, many unnecessary comparisons with sequences of no appreciable similarity are eliminated.</p> <p>In the FASTA method, the user defines a value&nbsp;<em>k</em>&nbsp;to use as the word length with which to search the database. The method is slower but more sensitive at lower values of&nbsp;<em>k</em>, which are also preferred for searches involving a very short query sequence. The BLAST family of search methods provides a number of algorithms optimized for particular types of queries, such as searching for distantly related sequence matches. BLAST was developed to provide a faster alternative to FASTA without sacrificing much accuracy; like FASTA, BLAST uses a word search of length&nbsp;<em>k</em>, but evaluates only the most significant word matches, rather than every word match as does FASTA. Most BLAST implementations use a fixed default word length that is optimized for the query and database type, and that is changed only under special circumstances, such as when searching with repetitive or very short query sequences. Implementations can be found via a number of web portals, such as&nbsp;<a href="http://www.ebi.ac.uk/fasta33/" rel="nofollow">EMBL FASTA</a>and&nbsp;<a href="https://www.ncbi.nlm.nih.gov/BLAST/" rel="nofollow">NCBI BLAST</a>.</p> <h2>Multiple sequence alignment[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=9" title="Edit section: Multiple sequence alignment">edit</a>]</h2> <p>Main article:&nbsp;<a href="https://en.wikipedia.org/wiki/Multiple_sequence_alignment" title="Multiple sequence alignment">Multiple sequence alignment</a></p> <p><a href="https://en.wikipedia.org/wiki/File:Hemagglutinin-alignments.png"><img alt="" src="https://upload.wikimedia.org/wikipedia/commons/thumb/b/b4/Hemagglutinin-alignments.png/300px-Hemagglutinin-alignments.png" style="height:322px; width:300px" /></a></p> <p>Alignment of 27&nbsp;<a href="https://en.wikipedia.org/wiki/Avian_influenza" title="Avian influenza">avian influenza</a>&nbsp;<a href="https://en.wikipedia.org/wiki/Hemagglutinin" title="Hemagglutinin">hemagglutinin</a>protein sequences colored by residue conservation (top) and residue properties (bottom)</p> <p><a href="https://en.wikipedia.org/wiki/Multiple_sequence_alignment" title="Multiple sequence alignment">Multiple sequence alignment</a>&nbsp;is an extension of pairwise alignment to incorporate more than two sequences at a time. Multiple alignment methods try to align all of the sequences in a given query set. Multiple alignments are often used in identifying&nbsp;<a href="https://en.wikipedia.org/wiki/Conservation_(genetics)" title="Conservation (genetics)">conserved</a>&nbsp;sequence regions across a group of sequences hypothesized to be evolutionarily related. Such conserved sequence motifs can be used in conjunction with structural and&nbsp;<a href="https://en.wikipedia.org/wiki/Reaction_mechanism" title="Reaction mechanism">mechanistic</a>information to locate the catalytic&nbsp;<a href="https://en.wikipedia.org/wiki/Active_site" title="Active site">active sites</a>&nbsp;of&nbsp;<a href="https://en.wikipedia.org/wiki/Enzyme" title="Enzyme">enzymes</a>. Alignments are also used to aid in establishing evolutionary relationships by constructing&nbsp;<a href="https://en.wikipedia.org/wiki/Phylogenetic_tree" title="Phylogenetic tree">phylogenetic trees</a>. Multiple sequence alignments are computationally difficult to produce and most formulations of the problem lead to&nbsp;<a href="https://en.wikipedia.org/wiki/NP-complete" title="NP-complete">NP-complete</a>&nbsp;combinatorial optimization problems.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-wang-8">[8]</a></sup><sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-elias-9">[9]</a></sup>&nbsp;Nevertheless, the utility of these alignments in bioinformatics has led to the development of a variety of methods suitable for aligning three or more sequences.</p> <h3>Dynamic programming[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=10" title="Edit section: Dynamic programming">edit</a>]</h3> <p>The technique of dynamic programming is theoretically applicable to any number of sequences; however, because it is computationally expensive in both time and&nbsp;<a href="https://en.wikipedia.org/wiki/Computer_memory" title="Computer memory">memory</a>, it is rarely used for more than three or four sequences in its most basic form. This method requires constructing the&nbsp;<em>n</em>-dimensional equivalent of the sequence matrix formed from two sequences, where&nbsp;<em>n</em>&nbsp;is the number of sequences in the query. Standard dynamic programming is first used on all pairs of query sequences and then the &quot;alignment space&quot; is filled in by considering possible matches or gaps at intermediate positions, eventually constructing an alignment essentially between each two-sequence alignment. Although this technique is computationally expensive, its guarantee of a global optimum solution is useful in cases where only a few sequences need to be aligned accurately. One method for reducing the computational demands of dynamic programming, which relies on the &quot;sum of pairs&quot;&nbsp;<a href="https://en.wikipedia.org/wiki/Objective_function" title="Objective function">objective function</a>, has been implemented in the&nbsp;<a href="https://www.ncbi.nlm.nih.gov/CBBresearch/Schaffer/msa.html" rel="nofollow">MSA</a>&nbsp;software package.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-lipman-10">[10]</a></sup></p> <h3>Progressive methods[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=11" title="Edit section: Progressive methods">edit</a>]</h3> <p>Progressive, hierarchical, or tree methods generate a multiple sequence alignment by first aligning the most similar sequences and then adding successively less related sequences or groups to the alignment until the entire query set has been incorporated into the solution. The initial tree describing the sequence relatedness is based on pairwise comparisons that may include heuristic pairwise alignment methods similar to&nbsp;<a href="https://en.wikipedia.org/wiki/FASTA" title="FASTA">FASTA</a>. Progressive alignment results are dependent on the choice of &quot;most related&quot; sequences and thus can be sensitive to inaccuracies in the initial pairwise alignments. Most progressive multiple sequence alignment methods additionally weight the sequences in the query set according to their relatedness, which reduces the likelihood of making a poor choice of initial sequences and thus improves alignment accuracy.</p> <p>Many variations of the&nbsp;<a href="https://en.wikipedia.org/wiki/Clustal" title="Clustal">Clustal</a>&nbsp;progressive implementation<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-higgins-11">[11]</a></sup><sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-thompson-12">[12]</a></sup><sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-chenna-13">[13]</a></sup>&nbsp;are used for multiple sequence alignment, phylogenetic tree construction, and as input for&nbsp;<a href="https://en.wikipedia.org/wiki/Protein_structure_prediction" title="Protein structure prediction">protein structure prediction</a>. A slower but more accurate variant of the progressive method is known as&nbsp;<a href="https://en.wikipedia.org/wiki/T-Coffee" title="T-Coffee">T-Coffee</a>.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-notredame-14">[14]</a></sup></p> <h3>Iterative methods[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=12" title="Edit section: Iterative methods">edit</a>]</h3> <p>Iterative methods attempt to improve on the heavy dependence on the accuracy of the initial pairwise alignments, which is the weak point of the progressive methods. Iterative methods optimize an&nbsp;<a href="https://en.wikipedia.org/wiki/Objective_function" title="Objective function">objective function</a>&nbsp;based on a selected alignment scoring method by assigning an initial global alignment and then realigning sequence subsets. The realigned subsets are then themselves aligned to produce the next iteration&#39;s multiple sequence alignment. Various ways of selecting the sequence subgroups and objective function are reviewed in.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-hirosawa-15">[15]</a></sup></p> <h3>Motif finding[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=13" title="Edit section: Motif finding">edit</a>]</h3> <p>Motif finding, also known as profile analysis, constructs global multiple sequence alignments that attempt to align short conserved&nbsp;<a href="https://en.wikipedia.org/wiki/Sequence_motif" title="Sequence motif">sequence motifs</a>&nbsp;among the sequences in the query set. This is usually done by first constructing a general global multiple sequence alignment, after which the highly&nbsp;<a href="https://en.wikipedia.org/wiki/Conservation_(genetics)" title="Conservation (genetics)">conserved</a>&nbsp;regions are isolated and used to construct a set of profile matrices. The profile matrix for each conserved region is arranged like a scoring matrix but its frequency counts for each amino acid or nucleotide at each position are derived from the conserved region&#39;s character distribution rather than from a more general empirical distribution. The profile matrices are then used to search other sequences for occurrences of the motif they characterize. In cases where the original&nbsp;<a href="https://en.wikipedia.org/wiki/Data_set" title="Data set">data set</a>&nbsp;contained a small number of sequences, or only highly related sequences,&nbsp;<a href="https://en.wikipedia.org/wiki/Pseudocount" title="Pseudocount">pseudocounts</a>&nbsp;are added to normalize the character distributions represented in the motif.</p> <h3>Techniques inspired by computer science[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=14" title="Edit section: Techniques inspired by computer science">edit</a>]</h3> <p>A variety of general&nbsp;<a href="https://en.wikipedia.org/wiki/Optimization_(mathematics)" title="Optimization (mathematics)">optimization</a>&nbsp;algorithms commonly used in computer science have also been applied to the multiple sequence alignment problem.&nbsp;<a href="https://en.wikipedia.org/wiki/Hidden_Markov_model" title="Hidden Markov model">Hidden Markov models</a>&nbsp;have been used to produce probability scores for a family of possible multiple sequence alignments for a given query set; although early HMM-based methods produced underwhelming performance, later applications have found them especially effective in detecting remotely related sequences because they are less susceptible to noise created by conservative or semiconservative substitutions.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-karplus-16">[16]</a></sup>&nbsp;<a href="https://en.wikipedia.org/wiki/Genetic_algorithm" title="Genetic algorithm">Genetic algorithms</a>and&nbsp;<a href="https://en.wikipedia.org/wiki/Simulated_annealing" title="Simulated annealing">simulated annealing</a>&nbsp;have also been used in optimizing multiple sequence alignment scores as judged by a scoring function like the sum-of-pairs method. More complete details and software packages can be found in the main article&nbsp;<a href="https://en.wikipedia.org/wiki/Multiple_sequence_alignment" title="Multiple sequence alignment">multiple sequence alignment</a>.</p> <p>The&nbsp;<a href="https://en.wikipedia.org/wiki/Burrows%E2%80%93Wheeler_transform" title="Burrows–Wheeler transform">Burrows&ndash;Wheeler transform</a>&nbsp;has been successfully applied to fast short read alignment in popular tools such as&nbsp;<a href="https://en.wikipedia.org/wiki/Bowtie_(sequence_analysis)" title="Bowtie (sequence analysis)">Bowtie</a>&nbsp;and BWA. See&nbsp;<a href="https://en.wikipedia.org/wiki/FM-index" title="FM-index">FM-index</a>.</p> <h2>Structural alignment[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=15" title="Edit section: Structural alignment">edit</a>]</h2> <p>Main article:&nbsp;<a href="https://en.wikipedia.org/wiki/Structural_alignment" title="Structural alignment">Structural alignment</a></p> <p>Structural alignments, which are usually specific to protein and sometimes RNA sequences, use information about the&nbsp;<a href="https://en.wikipedia.org/wiki/Secondary_structure" title="Secondary structure">secondary</a>&nbsp;and&nbsp;<a href="https://en.wikipedia.org/wiki/Tertiary_structure" title="Tertiary structure">tertiary structure</a>&nbsp;of the protein or RNA molecule to aid in aligning the sequences. These methods can be used for two or more sequences and typically produce local alignments; however, because they depend on the availability of structural information, they can only be used for sequences whose corresponding structures are known (usually through&nbsp;<a href="https://en.wikipedia.org/wiki/X-ray_crystallography" title="X-ray crystallography">X-ray crystallography</a>&nbsp;or&nbsp;<a href="https://en.wikipedia.org/wiki/NMR_spectroscopy" title="NMR spectroscopy">NMR spectroscopy</a>). Because both protein and RNA structure is more evolutionarily conserved than sequence,<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-chothia-17">[17]</a></sup>&nbsp;structural alignments can be more reliable between sequences that are very distantly related and that have diverged so extensively that sequence comparison cannot reliably detect their similarity.</p> <p>Structural alignments are used as the &quot;gold standard&quot; in evaluating alignments for homology-based&nbsp;<a href="https://en.wikipedia.org/wiki/Protein_structure_prediction" title="Protein structure prediction">protein structure prediction</a><sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-skolnick-18">[18]</a></sup>&nbsp;because they explicitly align regions of the protein sequence that are structurally similar rather than relying exclusively on sequence information. However, clearly structural alignments cannot be used in structure prediction because at least one sequence in the query set is the target to be modeled, for which the structure is not known. It has been shown that, given the structural alignment between a target and a template sequence, highly accurate models of the target protein sequence can be produced; a major stumbling block in homology-based structure prediction is the production of structurally accurate alignments given only sequence information.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-skolnick-18">[18]</a></sup></p> <h3>DALI[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=16" title="Edit section: DALI">edit</a>]</h3> <p>The DALI method, or&nbsp;<a href="https://en.wikipedia.org/wiki/Distance_matrix" title="Distance matrix">distance matrix</a>&nbsp;alignment, is a fragment-based method for constructing structural alignments based on contact similarity patterns between successive hexapeptides in the query sequences.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-holm-19">[19]</a></sup>&nbsp;It can generate pairwise or multiple alignments and identify a query sequence&#39;s structural neighbors in the&nbsp;<a href="https://en.wikipedia.org/wiki/Protein_Data_Bank" title="Protein Data Bank">Protein Data Bank</a>&nbsp;(PDB). It has been used to construct the&nbsp;<a href="https://en.wikipedia.org/wiki/Families_of_structurally_similar_proteins" title="Families of structurally similar proteins">FSSP</a>&nbsp;structural alignment database (Fold classification based on Structure-Structure alignment of Proteins, or Families of Structurally Similar Proteins). A DALI webserver can be accessed at&nbsp;<a href="https://web.archive.org/web/20090301064750/http://ekhidna.biocenter.helsinki.fi/dali_server/start" rel="nofollow">DALI</a>&nbsp;and the FSSP is located at&nbsp;<a href="https://web.archive.org/web/20051125045348/http://ekhidna.biocenter.helsinki.fi/dali/start" rel="nofollow">The Dali Database</a>.</p> <h3>SSAP[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=17" title="Edit section: SSAP">edit</a>]</h3> <p>SSAP (sequential structure alignment program) is a dynamic programming-based method of structural alignment that uses atom-to-atom vectors in structure space as comparison points. It has been extended since its original description to include multiple as well as pairwise alignments,<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-taylor-20">[20]</a></sup>&nbsp;and has been used in the construction of the&nbsp;<a href="https://en.wikipedia.org/wiki/CATH" title="CATH">CATH</a>&nbsp;(Class, Architecture, Topology, Homology) hierarchical database classification of protein folds.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-orengo-21">[21]</a></sup>&nbsp;The CATH database can be accessed at&nbsp;<a href="http://www.cathdb.info/" rel="nofollow">CATH Protein Structure Classification</a>.</p> <h3>Combinatorial extension[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=18" title="Edit section: Combinatorial extension">edit</a>]</h3> <p>The combinatorial extension method of structural alignment generates a pairwise structural alignment by using local geometry to align short fragments of the two proteins being analyzed and then assembles these fragments into a larger alignment.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-shindyalov-22">[22]</a></sup>&nbsp;Based on measures such as rigid-body&nbsp;<a href="https://en.wikipedia.org/wiki/Root_mean_square_deviation_(bioinformatics)" title="Root mean square deviation (bioinformatics)">root mean square distance</a>, residue distances, local secondary structure, and surrounding environmental features such as residue neighbor&nbsp;<a href="https://en.wikipedia.org/wiki/Hydrophobic" title="Hydrophobic">hydrophobicity</a>, local alignments called &quot;aligned fragment pairs&quot; are generated and used to build a similarity matrix representing all possible structural alignments within predefined cutoff criteria. A path from one protein structure state to the other is then traced through the matrix by extending the growing alignment one fragment at a time. The optimal such path defines the combinatorial-extension alignment. A web-based server implementing the method and providing a database of pairwise alignments of structures in the Protein Data Bank is located at the&nbsp;<a href="https://web.archive.org/web/19981203071023/http://cl.sdsc.edu/" rel="nofollow">Combinatorial Extension</a>&nbsp;website.</p> <h2>Phylogenetic analysis[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=19" title="Edit section: Phylogenetic analysis">edit</a>]</h2> <p>Main article:&nbsp;<a href="https://en.wikipedia.org/wiki/Computational_phylogenetics" title="Computational phylogenetics">Computational phylogenetics</a></p> <p>Phylogenetics and sequence alignment are closely related fields due to the shared necessity of evaluating sequence relatedness.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-ortet-23">[23]</a></sup>&nbsp;The field of&nbsp;<a href="https://en.wikipedia.org/wiki/Phylogenetics" title="Phylogenetics">phylogenetics</a>&nbsp;makes extensive use of sequence alignments in the construction and interpretation of&nbsp;<a href="https://en.wikipedia.org/wiki/Phylogenetic_tree" title="Phylogenetic tree">phylogenetic trees</a>, which are used to classify the evolutionary relationships between homologous&nbsp;<a href="https://en.wikipedia.org/wiki/Gene" title="Gene">genes</a>&nbsp;represented in the&nbsp;<a href="https://en.wikipedia.org/wiki/Genome" title="Genome">genomes</a>&nbsp;of divergent species. The degree to which sequences in a query set differ is qualitatively related to the sequences&#39; evolutionary distance from one another. Roughly speaking, high sequence identity suggests that the sequences in question have a comparatively young&nbsp;<a href="https://en.wikipedia.org/wiki/Most_recent_common_ancestor" title="Most recent common ancestor">most recent common ancestor</a>, while low identity suggests that the divergence is more ancient. This approximation, which reflects the &quot;<a href="https://en.wikipedia.org/wiki/Molecular_clock" title="Molecular clock">molecular clock</a>&quot; hypothesis that a roughly constant rate of evolutionary change can be used to extrapolate the elapsed time since two genes first diverged (that is, the&nbsp;<a href="https://en.wikipedia.org/wiki/Coalescence_(genetics)" title="Coalescence (genetics)">coalescence</a>&nbsp;time), assumes that the effects of mutation and&nbsp;<a href="https://en.wikipedia.org/wiki/Natural_selection" title="Natural selection">selection</a>&nbsp;are constant across sequence lineages. Therefore, it does not account for possible difference among organisms or species in the rates of&nbsp;<a href="https://en.wikipedia.org/wiki/DNA_repair" title="DNA repair">DNA repair</a>&nbsp;or the possible functional conservation of specific regions in a sequence. (In the case of nucleotide sequences, the molecular clock hypothesis in its most basic form also discounts the difference in acceptance rates between&nbsp;<a href="https://en.wikipedia.org/wiki/Silent_mutation" title="Silent mutation">silent mutations</a>&nbsp;that do not alter the meaning of a given&nbsp;<a href="https://en.wikipedia.org/wiki/Codon" title="Codon">codon</a>&nbsp;and other mutations that result in a different&nbsp;<a href="https://en.wikipedia.org/wiki/Amino_acid" title="Amino acid">amino acid</a>&nbsp;being incorporated into the protein). More statistically accurate methods allow the evolutionary rate on each branch of the phylogenetic tree to vary, thus producing better estimates of coalescence times for genes.</p> <p>Progressive multiple alignment techniques produce a phylogenetic tree by necessity because they incorporate sequences into the growing alignment in order of relatedness. Other techniques that assemble multiple sequence alignments and phylogenetic trees score and sort trees first and calculate a multiple sequence alignment from the highest-scoring tree. Commonly used methods of phylogenetic tree construction are mainly&nbsp;<a href="https://en.wikipedia.org/wiki/Heuristic" title="Heuristic">heuristic</a>&nbsp;because the problem of selecting the optimal tree, like the problem of selecting the optimal multiple sequence alignment, is&nbsp;<a href="https://en.wikipedia.org/wiki/NP-hard" title="NP-hard">NP-hard</a>.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-felsenstein-24">[24]</a></sup></p> <h3>Assessment of significance[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=20" title="Edit section: Assessment of significance">edit</a>]</h3> <p>Sequence alignments are useful in bioinformatics for identifying sequence similarity, producing phylogenetic trees, and developing homology models of protein structures. However, the biological relevance of sequence alignments is not always clear. Alignments are often assumed to reflect a degree of evolutionary change between sequences descended from a common ancestor; however, it is formally possible that&nbsp;<a href="https://en.wikipedia.org/wiki/Convergent_evolution" title="Convergent evolution">convergent evolution</a>&nbsp;can occur to produce apparent similarity between proteins that are evolutionarily unrelated but perform similar functions and have similar structures.</p> <p>In database searches such as BLAST, statistical methods can determine the likelihood of a particular alignment between sequences or sequence regions arising by chance given the size and composition of the database being searched. These values can vary significantly depending on the search space. In particular, the likelihood of finding a given alignment by chance increases if the database consists only of sequences from the same organism as the query sequence. Repetitive sequences in the database or query can also distort both the search results and the assessment of statistical significance; BLAST automatically filters such repetitive sequences in the query to avoid apparent hits that are statistical artifacts.</p> <p>Methods of statistical significance estimation for gapped sequence alignments are available in the literature.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-ortet-23">[23]</a></sup><sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-altschul-25">[25]</a></sup><sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-hartmann-26">[26]</a></sup><sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-newberg-27">[27]</a></sup><sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-eddy-28">[28]</a></sup><sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-bastien-29">[29]</a></sup><sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-agrawal11-30">[30]</a></sup><sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-agrawal08-31">[31]</a></sup></p> <h3>Assessment of credibility[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=21" title="Edit section: Assessment of credibility">edit</a>]</h3> <p>Statistical significance indicates the probability that an alignment of a given quality could arise by chance, but does not indicate how much superior a given alignment is to alternative alignments of the same sequences. Measures of alignment credibility indicate the extent to which the best scoring alignments for a given pair of sequences are substantially similar. Methods of alignment credibility estimation for gapped sequence alignments are available in the literature.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-NewbergLawrence2009-32">[32]</a></sup></p> <h3>Scoring functions[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=22" title="Edit section: Scoring functions">edit</a>]</h3> <p>The choice of a scoring function that reflects biological or statistical observations about known sequences is important to producing good alignments. Protein sequences are frequently aligned using&nbsp;<a href="https://en.wikipedia.org/wiki/Substitution_matrix" title="Substitution matrix">substitution matrices</a>&nbsp;that reflect the probabilities of given character-to-character substitutions. A series of matrices called&nbsp;<a href="https://en.wikipedia.org/wiki/Point_accepted_mutation" title="Point accepted mutation">PAM matrices</a>&nbsp;(Point Accepted Mutation matrices, originally defined by&nbsp;<a href="https://en.wikipedia.org/wiki/Margaret_Dayhoff" title="Margaret Dayhoff">Margaret Dayhoff</a>&nbsp;and sometimes referred to as &quot;Dayhoff matrices&quot;) explicitly encode evolutionary approximations regarding the rates and probabilities of particular amino acid mutations. Another common series of scoring matrices, known as&nbsp;<a href="https://en.wikipedia.org/wiki/BLOSUM" title="BLOSUM">BLOSUM</a>&nbsp;(Blocks Substitution Matrix), encodes empirically derived substitution probabilities. Variants of both types of matrices are used to detect sequences with differing levels of divergence, thus allowing users of BLAST or FASTA to restrict searches to more closely related matches or expand to detect more divergent sequences.&nbsp;<a href="https://en.wikipedia.org/wiki/Gap_penalty" title="Gap penalty">Gap penalties</a>account for the introduction of a gap - on the evolutionary model, an insertion or deletion mutation - in both nucleotide and protein sequences, and therefore the penalty values should be proportional to the expected rate of such mutations. The quality of the alignments produced therefore depends on the quality of the scoring function.</p> <p>It can be very useful and instructive to try the same alignment several times with different choices for scoring matrix and/or gap penalty values and compare the results. Regions where the solution is weak or non-unique can often be identified by observing which regions of the alignment are robust to variations in alignment parameters.</p> <h2>Other biological uses[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=23" title="Edit section: Other biological uses">edit</a>]</h2> <p>Sequenced RNA, such as&nbsp;<a href="https://en.wikipedia.org/wiki/Expressed_sequence_tags" title="Expressed sequence tags">expressed sequence tags</a>&nbsp;and full-length mRNAs, can be aligned to a sequenced genome to find where there are genes and get information about&nbsp;<a href="https://en.wikipedia.org/wiki/Alternative_splicing" title="Alternative splicing">alternative splicing</a><sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-33">[33]</a></sup>&nbsp;and&nbsp;<a href="https://en.wikipedia.org/wiki/RNA_editing" title="RNA editing">RNA editing</a>.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-34">[34]</a></sup>&nbsp;Sequence alignment is also a part of&nbsp;<a href="https://en.wikipedia.org/wiki/Genome_assembly" title="Genome assembly">genome assembly</a>, where sequences are aligned to find overlap so that&nbsp;<em><a href="https://en.wikipedia.org/wiki/Contig" title="Contig">contigs</a></em>&nbsp;(long stretches of sequence) can be formed.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-35">[35]</a></sup>&nbsp;Another use is&nbsp;<a href="https://en.wikipedia.org/wiki/Single_nucleotide_polymorphism" title="Single nucleotide polymorphism">SNP</a>analysis, where sequences from different individuals are aligned to find single basepairs that are often different in a population.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-36">[36]</a></sup></p> <h2>Non-biological uses[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=24" title="Edit section: Non-biological uses">edit</a>]</h2> <p>The methods used for biological sequence alignment have also found applications in other fields, most notably in&nbsp;<a href="https://en.wikipedia.org/wiki/Natural_language_processing" title="Natural language processing">natural language processing</a>&nbsp;and in social sciences, where the&nbsp;<a href="https://en.wikipedia.org/wiki/Needleman-Wunsch_algorithm" title="Needleman-Wunsch algorithm">Needleman-Wunsch algorithm</a>&nbsp;is usually referred to as&nbsp;<a href="https://en.wikipedia.org/wiki/Optimal_matching" title="Optimal matching">Optimal matching</a>.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-37">[37]</a></sup>&nbsp;Techniques that generate the set of elements from which words will be selected in natural-language generation algorithms have borrowed multiple sequence alignment techniques from bioinformatics to produce linguistic versions of computer-generated mathematical proofs.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-Barzilay-38">[38]</a></sup>&nbsp;In the field of historical and comparative&nbsp;<a href="https://en.wikipedia.org/wiki/Linguistics" title="Linguistics">linguistics</a>, sequence alignment has been used to partially automate the&nbsp;<a href="https://en.wikipedia.org/wiki/Comparative_method_(linguistics)" title="Comparative method (linguistics)">comparative method</a>&nbsp;by which linguists traditionally reconstruct languages.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-39">[39]</a></sup>&nbsp;Business and marketing research has also applied multiple sequence alignment techniques in analyzing series of purchases over time.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-prinzie-40">[40]</a></sup></p> <h2>Software[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=25" title="Edit section: Software">edit</a>]</h2> <p>Main article:&nbsp;<a href="https://en.wikipedia.org/wiki/Sequence_alignment_software" title="Sequence alignment software">Sequence alignment software</a></p> <p>A more complete list of available software categorized by algorithm and alignment type is available at&nbsp;<a href="https://en.wikipedia.org/wiki/Sequence_alignment_software" title="Sequence alignment software">sequence alignment software</a>, but common software tools used for general sequence alignment tasks include ClustalW2<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-41">[41]</a></sup>&nbsp;and T-coffee<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-42">[42]</a></sup>&nbsp;for alignment, and BLAST<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-43">[43]</a></sup>&nbsp;and FASTA3x<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-44">[44]</a></spansup>&nbsp;for database searching. Commercial tools such as&nbsp;<a href="https://en.wikipedia.org/wiki/DNASTAR" title="DNASTAR">DNASTAR Lasergene</a>,&nbsp;<a href="https://en.wikipedia.org/w/index.php?title=Geneious&amp;action=edit&amp;redlink=1" title="Geneious (page does not exist)">Geneious</a>, and&nbsp;<a href="https://en.wikipedia.org/wiki/PatternHunter" title="PatternHunter">PatternHunter</a>&nbsp;are also available. Tools annotated as performing&nbsp;<a href="http://edamontology.org/operation_0292" rel="nofollow">sequence alignment</a>&nbsp;are listed in the&nbsp;<a href="https://bio.tools/?page=1&amp;function=%22Sequence%20alignment%22&amp;sort=score" rel="nofollow">bio.tools</a>&nbsp;registry.</p> <p>Alignment algorithms and software can be directly compared to one another using a standardized set of&nbsp;<a href="https://en.wikipedia.org/wiki/Benchmark_(computing)" title="Benchmark (computing)">benchmark</a>&nbsp;reference multiple sequence alignments known as BAliBASE.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-thompson2-45">[45]</a></sup>&nbsp;The data set consists of structural alignments, which can be considered a standard against which purely sequence-based methods are compared. The relative performance of many common alignment methods on frequently encountered alignment problems has been tabulated and selected results published online at BAliBASE.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-46">[46]</a></sup><sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-thompson3-47">[47]</a></sup>&nbsp;A comprehensive list of BAliBASE scores for many (currently 12) different alignment tools can be computed within the protein workbench STRAP.<sup><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_note-48">[48]</a></sup></p> <h2>See also[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=26" title="Edit section: See also">edit</a>]</h2> <ul> <li><a href="https://en.wikipedia.org/wiki/Sequence_homology" title="Sequence homology">Sequence homology</a></li> <li><a href="https://en.wikipedia.org/wiki/Sequence_mining" title="Sequence mining">Sequence mining</a></li> <li><a href="https://en.wikipedia.org/wiki/BLAST" title="BLAST">BLAST</a></li> <li><a href="https://en.wikipedia.org/wiki/String_searching_algorithm" title="String searching algorithm">String searching algorithm</a></li> <li><a href="https://en.wikipedia.org/wiki/Alignment-free_sequence_analysis" title="Alignment-free sequence analysis">Alignment-free sequence analysis</a></li> <li><a href="https://en.wikipedia.org/wiki/UGENE" title="UGENE">UGENE</a></li> <li><a href="https://en.wikipedia.org/wiki/Needleman%E2%80%93Wunsch_algorithm" title="Needleman–Wunsch algorithm">Needleman&ndash;Wunsch algorithm</a></li></ul> <h2>References[<a href="https://en.wikipedia.org/w/index.php?title=Sequence_alignment&amp;action=edit&amp;section=27" title="Edit section: References">edit</a>]</h2>
<ol>
<li>^&nbsp;<span stylea href="font-size:14px">https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-mount_1-0">Jump up to:<sup><em><strong>a</strong></spanem></lisup> <li/a>&nbsp;<span stylea href="font-sizehttps:14px">Valery, O//en. P., Mikhail A, R., &amp; Vladimir G. T. (2011). Comparative analysis of the quality of a global algorithm and a local algorithm for alignment of two sequenceswikipedia. org/wiki/Sequence_alignment#cite_ref-mount_1-1"><sup><em>Algorithms Mol Biol</emstrong>. doi:&nbsp; 10.1186b</1748-7188-6-25.strong></spanem></lisup> <li/a>&nbsp;<span stylea href="font-sizehttps:14px">Philippe O//en., &amp; Olivier Bwikipedia. (2010). Where Does the Alignment Score Distribution Shape Come from? org/wiki/Sequence_alignment#cite_ref-mount_1-2"><sup><em>Evolutionary Bioinformatics<strong>c</strong></em></sup></a>&nbsp;<cite>Mount DM. (2004). &nbsp;<em>6Bioinformatics: Sequence and Genome Analysis</em>&nbsp;(2nd ed.). Cold Spring Harbor Laboratory Press: 159Cold Spring Harbor, NY.&nbsp;<a href="https://en.wikipedia.org/wiki/International_Standard_Book_Number" title="International Standard Book Number">ISBN</a>&ndashnbsp;<a href="https://en.wikipedia.org/wiki/Special:BookSources/978-0-87969-608-5" title="Special:BookSources/978-0-87969-608-5">978-0-87969-608-5</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-2" title="Jump up">^</a></strong>&nbsp;<cite><a href="http://www.ebi.ac.uk/Tools/msa/clustalw2/help/faq.html#23" rel="nofollow">&quot;Clustal FAQ #Symbols&quot;</a>.&nbsp;<em>Clustal</em>. Retrieved&nbsp;8 December&nbsp;2014.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-predict_3-0" title="Jump up">^</a></strong>&nbsp;<cite>Ng PC; Henikoff S (May 2001).&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC311071" rel="nofollow">&quot;Predicting deleterious amino acid substitutions&quot;</a>.&nbsp;<em>Genome Res</em>.&nbsp;<strong>11</strong>&nbsp;(5): 863&ndash;74.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1101%2Fgr.176601" rel="nofollow">10.1101/gr.176601</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Central" title="PubMed Central">PMC</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC311071" rel="nofollow">311071</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/11337480" rel="nofollow">11337480</a>.</cite></li> <li>^&nbsp;<a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-Polyanovsky2011_4-0">Jump up to:<sup><em><strong>a</strong></em></sup></a>&nbsp;<a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-Polyanovsky2011_4-1"><sup><em><strong>b</strong></em></sup></a>&nbsp;<cite>Polyanovsky, V. O.; Roytberg, M. A.; Tumanyan, V. G. (2011).&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3223492" rel="nofollow">&quot;Comparative analysis of the quality of a global algorithm and a local algorithm for alignment of two sequences&quot;</a>.&nbsp;<em>Algorithms for Molecular Biology</em>.&nbsp;<strong>6</strong>&nbsp;(1): 25.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1186%2F1748-7188-6-25" rel="nofollow">10.1186/1748-7188-6-25</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Central" title="PubMed Central">PMC</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3223492" rel="nofollow">3223492</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/22032267" rel="nofollow">22032267</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-Schneider_5-0" title="Jump up">^</a></strong>&nbsp;<cite>Schneider TD; Stephens RM (1990).&nbsp;<a href="http://nar.oxfordjournals.org/cgi/pmidlookup?view=long&amp;pmid=2172928" rel="nofollow">&quot;Sequence logos: a new way to display consensus sequences&quot;</a>.&nbsp;<em>Nucleic Acids Res</em>.&nbsp;<strong>18</strong>(20): 6097&ndash;6100.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1093%2Fnar%2F18.20.6097" rel="nofollow">10.1093/nar/18.20.6097</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Central" title="PubMed Central">PMC</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC332411" rel="nofollow">332411</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/2172928" rel="nofollow">2172928</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-6" title="Jump up">^</a></strong>&nbsp;<cite><a href="https://samtools.github.io/hts-specs/SAMv1.pdf" rel="nofollow">&quot;Sequence Alignment/Map Format Specification&quot;</a>&nbsp;(PDF).</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-brudno_7-0" title="Jump up">^</a></strong>&nbsp;<cite>Brudno M; Malde S; Poliakov A; Do CB; Couronne O; Dubchak I; Batzoglou S (2003).&nbsp;<a href="http://bioinformatics.oxfordjournals.org/cgi/pmidlookup?view=long&amp;pmid=12855437" rel="nofollow">&quot;Glocal alignment: finding rearrangements during alignment&quot;</a>.&nbsp;<em>Bioinformatics</em>. 19. Suppl 1 (90001): i54&ndash;62.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1093%2Fbioinformatics%2Fbtg1005" rel="nofollow">10.1093/bioinformatics/btg1005</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/12855437" rel="nofollow">12855437</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-wang_8-0" title="Jump up">^</a></strong>&nbsp;<cite>Wang L; Jiang T. (1994). &quot;On the complexity of multiple sequence alignment&quot;.&nbsp;<em>J Comput Biol</em>.&nbsp;<strong>1</strong>&nbsp;(4): 337&ndash;48.&nbsp;<a href="https://en.wikipedia.org/wiki/CiteSeerX" title="CiteSeerX">CiteSeerX</a>&nbsp;<a href="https://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.408.894" rel="nofollow">10.1.1.408.894</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1089%2Fcmb.1994.1.337" rel="nofollow">10.1089/cmb.1994.1.337</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/8790475" rel="nofollow">8790475</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-elias_9-0" title="Jump up">^</a></strong>&nbsp;<cite>Elias, Isaac (2006). &quot;Settling the intractability of multiple alignment&quot;.&nbsp;<em>J Comput Biol</em>.&nbsp;<strong>13</strong>&nbsp;(7): 1323&ndash;1339.&nbsp;<a href="https://en.wikipedia.org/wiki/CiteSeerX" title="CiteSeerX">CiteSeerX</a>&nbsp;<a href="https://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.6.256" rel="nofollow">10.1.1.6.256</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1089%2Fcmb.2006.13.1323" rel="nofollow">10.1089/cmb.2006.13.1323</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/17037961" rel="nofollow">17037961</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-lipman_10-0" title="Jump up">^</a></strong>&nbsp;<cite>Lipman DJ; Altschul SF; Kececioglu JD (1989).&nbsp;<a href="http://www.pnas.org/cgi/pmidlookup?view=long&amp;pmid=2734293" rel="nofollow">&quot;A tool for multiple sequence alignment&quot;</a>.&nbsp;<em>Proc Natl Acad Sci USA</em>.&nbsp;<strong>86</strong>&nbsp;(12): 4412&ndash;5.&nbsp;<a href="https://en.wikipedia.org/wiki/Bibcode" title="Bibcode">Bibcode</a>:<a href="http://adsabs.harvard.edu/abs/1989PNAS...86.4412L" rel="nofollow">1989PNAS...86.4412L</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1073%2Fpnas.86.12.4412" rel="nofollow">10.1073/pnas.86.12.4412</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Central" title="PubMed Central">PMC</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC287279" rel="nofollow">287279</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/2734293" rel="nofollow">2734293</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-higgins_11-0" title="Jump up">^</a></strong>&nbsp;<cite><a href="https://en.wikipedia.org/wiki/Desmond_G._Higgins" title="Desmond G. Higgins">Higgins DG</a>, Sharp PM (1988).&nbsp;<a href="http://linkinghub.elsevier.com/retrieve/pii/0378-1119(88)90330-7" rel="nofollow">&quot;CLUSTAL: a package for performing multiple sequence alignment on a microcomputer&quot;</a>.&nbsp;<em>Gene</em>.&nbsp;<strong>73</strong>&nbsp;(1): 237&ndash;44.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1016%2F0378-1119%2888%2990330-7" rel="nofollow">10.1016/0378-1119(88)90330-7</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/3243435" rel="nofollow">3243435</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-thompson_12-0" title="Jump up">^</a></strong>&nbsp;<cite>Thompson JD;&nbsp;<a href="https://en.wikipedia.org/wiki/Desmond_G._Higgins" title="Desmond G. Higgins">Higgins DG</a>; Gibson TJ. (1994).&nbsp;<a href="http://nar.oxfordjournals.org/content/22/22/4673" rel="nofollow">&quot;CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice&quot;</a>.&nbsp;<em>Nucleic Acids Res</em>.&nbsp;<strong>22</strong>&nbsp;(22): 4673&ndash;80.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1093%2Fnar%2F22.22.4673" rel="nofollow">10.1093/nar/22.22.4673</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Central" title="PubMed Central">PMC</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC308517" rel="nofollow">308517</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/7984417" rel="nofollow">7984417</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-chenna_13-0" title="Jump up">^</a></strong>&nbsp;<cite>Chenna R; Sugawara H; Koike T; Lopez R; Gibson TJ; Higgins DG; Thompson JD. (2003).&nbsp;<a href="http://nar.oxfordjournals.org/cgi/pmidlookup?view=long&amp;pmid=12824352" rel="nofollow">&quot;Multiple sequence alignment with the Clustal series of programs&quot;</a>.&nbsp;<em>Nucleic Acids Res</em>.&nbsp;<strong>31</strong>&nbsp;(13): 3497&ndash;500.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1093%2Fnar%2Fgkg500" rel="nofollow">10.1093/nar/gkg500</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Central" title="PubMed Central">PMC</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC168907" rel="nofollow">168907</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/12824352" rel="nofollow">12824352</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-notredame_14-0" title="Jump up">^</a></strong>&nbsp;<cite>Notredame C;&nbsp;<a href="https://en.wikipedia.org/wiki/Desmond_G._Higgins" title="Desmond G. Higgins">Higgins DG</a>; Heringa J. (2000).&nbsp;<a href="http://linkinghub.elsevier.com/retrieve/pii/S0022-2836(00)94042-7" rel="nofollow">&quot;T-Coffee: A novel method for fast and accurate multiple sequence alignment&quot;</a>.&nbsp;<em>J Mol Biol</em>.&nbsp;<strong>302</strong>&nbsp;(1): 205&ndash;17.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1006%2Fjmbi.2000.4042" rel="nofollow">10.1006/jmbi.2000.4042</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/10964570" rel="nofollow">10964570</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-hirosawa_15-0" title="Jump up">^</a></strong>&nbsp;<cite>Hirosawa M; Totoki Y; Hoshida M; Ishikawa M. (1995).&nbsp;<a href="http://bioinformatics.oxfordjournals.org/cgi/content/abstract/11/1/13" rel="nofollow">&quot;Comprehensive study on iterative algorithms of multiple sequence alignment&quot;</a>.&nbsp;<em>Comput Appl Biosci</em>.&nbsp;<strong>11</strong>&nbsp;(1): 13&ndash;8.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1093%2Fbioinformatics%2F11.1.13" rel="nofollow">10.1093/bioinformatics/11.1.13</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/7796270" rel="nofollow">7796270</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-karplus_16-0" title="Jump up">^</a></strong>&nbsp;<cite>Karplus K; Barrett C; Hughey R. (1998).&nbsp;<a href="http://bioinformatics.oxfordjournals.org/cgi/pmidlookup?view=long&amp;pmid=9927713" rel="nofollow">&quot;Hidden Markov models for detecting remote protein homologies&quot;</a>.&nbsp;<em>Bioinformatics</em>.&nbsp;<strong>14</strong>(10): 846&ndash;856.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1093%2Fbioinformatics%2F14.10.846" rel="nofollow">10.1093/bioinformatics/14.10.846</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/9927713" rel="nofollow">9927713</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-chothia_17-0" title="Jump up">^</a></strong>&nbsp;<cite>Chothia C; Lesk AM. (April 1986).&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1166865" rel="nofollow">&quot;The relation between the divergence of sequence and structure in proteins&quot;</a>.&nbsp;<em>EMBO J</em>.&nbsp;<strong>5</strong>(4): 823&ndash;6.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Central" title="PubMed Central">PMC</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1166865" rel="nofollow">1166865</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/3709526" rel="nofollow">3709526</a>.</cite></li> <li>^&nbsp;<a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-skolnick_18-0">Jump up to:<sup><em><strong>a</strong></em></sup></a>&nbsp;<a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-skolnick_18-1"><sup><em><strong>b</strong></em></sup></a>&nbsp;<cite>Zhang Y; Skolnick J. (2005).&nbsp;<a href="http://www.pnas.org/cgi/pmidlookup?view=long&amp;pmid=15653774" rel="nofollow">&quot;The protein structure prediction problem could be solved using the current PDB library&quot;</a>.&nbsp;<em>Proc Natl Acad Sci USA</em>.&nbsp;<strong>102</strong>&nbsp;(4): 1029&ndash;34.&nbsp;<a href="https://en.wikipedia.org/wiki/Bibcode" title="Bibcode">Bibcode</a>:<a href="http://adsabs.harvard.edu/abs/2005PNAS..102.1029Z" rel="nofollow">2005PNAS..102.1029Z</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1073%2Fpnas.0407152101" rel="nofollow">10.1073/pnas.0407152101</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Central" title="PubMed Central">PMC</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC545829" rel="nofollow">545829</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/15653774" rel="nofollow">15653774</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-holm_19-0" title="Jump up">^</a></strong>&nbsp;<cite>Holm L; Sander C (1996).&nbsp;<a href="http://www.sciencemag.org/cgi/pmidlookup?view=long&amp;pmid=8662544" rel="nofollow">&quot;Mapping the protein universe&quot;</a>.&nbsp;<em>Science</em>.&nbsp;<strong>273</strong>&nbsp;(5275): 595&ndash;603.&nbsp;<a href="https://en.wikipedia.org/wiki/Bibcode" title="Bibcode">Bibcode</a>:<a href="http://adsabs.harvard.edu/abs/1996Sci...273..595H" rel="nofollow">1996Sci...273..595H</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1126%2Fscience.273.5275.595" rel="nofollow">10.1126/science.273.5275.595</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/8662544" rel="nofollow">8662544</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-taylor_20-0" title="Jump up">^</a></strong>&nbsp;<cite>Taylor WR; Flores TP; Orengo CA. (1994).&nbsp;<a href="http://www.proteinscience.org/cgi/pmidlookup?view=long&amp;pmid=7849601" rel="nofollow">&quot;Multiple protein structure alignment&quot;</a>.&nbsp;<em>Protein Sci</em>.&nbsp;<strong>3</strong>&nbsp;(10): 1858&ndash;70.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1002%2Fpro.5560031025" rel="nofollow">10.1002/pro.5560031025</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Central" title="PubMed Central">PMC</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2142613" rel="nofollow">2142613</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/7849601" rel="nofollow">7849601</a>.</cite><sup>[<em><a href="https://en.wikipedia.org/wiki/Wikipedia:Link_rot" title="Wikipedia:Link rot">permanent dead link</a></em>]</sup></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-orengo_21-0" title="Jump up">^</a></strong>&nbsp;<cite>Orengo CA; Michie AD; Jones S; Jones DT; Swindells MB; Thornton JM (1997). &quot;CATH--a hierarchic classification of protein domain structures&quot;.&nbsp;<em>Structure</em>.&nbsp;<strong>5</strong>&nbsp;(8): 1093&ndash;108.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1016%2FS0969-2126%2897%2900260-8" rel="nofollow">10.1016/S0969-2126(97)00260-8</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/9309224" rel="nofollow">9309224</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-shindyalov_22-0" title="Jump up">^</a></strong>&nbsp;<cite>Shindyalov IN; Bourne PE. (1998).&nbsp;<a href="http://peds.oxfordjournals.org/cgi/pmidlookup?view=long&amp;pmid=9796821" rel="nofollow">&quot;Protein structure alignment by incremental combinatorial extension (CE) of the optimal path&quot;</a>.&nbsp;<em>Protein Eng</em>.&nbsp;<strong>11</strong>&nbsp;(9): 739&ndash;47.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1093%2Fprotein%2F11.9.739" rel="nofollow">10.1093/protein/11.9.739</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/9796821" rel="nofollow">9796821</a>.</cite></li> <li>^&nbsp;<a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-ortet_23-0">Jump up to:<sup><em><strong>a</strong></em></sup></a>&nbsp;<a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-ortet_23-1"><sup><em><strong>b</strong></em></sup></a>&nbsp;<cite>Ortet P; Bastien O (2010).&nbsp;<a href="http://www.la-press.com/where-does-the-alignment-score-distribution-shape-come-from-article-a2393" rel="nofollow">&quot;Where Does the Alignment Score Distribution Shape Come from?&quot;</a>.&nbsp;<em>Evolutionary Bioinformatics</em>.&nbsp;<strong>6</strong>: 159&ndash;187.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.4137%2FEBO.S5875" rel="nofollow">10.4137/EBO.S5875</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Central" title="PubMed Central">PMC</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3023300" rel="nofollow">3023300</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/21258650" rel="nofollow">21258650</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-felsenstein_24-0" title="Jump up">^</a></strong>&nbsp;<cite>Felsenstein J. (2004).&nbsp;<em>Inferring Phylogenies</em>. Sinauer Associates: Sunderland, MA.&nbsp;<a href="https://en.wikipedia.org/wiki/International_Standard_Book_Number" title="International Standard Book Number">ISBN</a>&nbsp;<a href="https://en.wikipedia.org/wiki/Special:BookSources/978-0-87893-177-4" title="Special:BookSources/978-0-87893-177-4">978-0-87893-177-4</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-altschul_25-0" title="Jump up">^</a></strong>&nbsp;<cite>Altschul SF; Gish W (1996).&nbsp;<em>Local Alignment Statistics</em>.&nbsp;<em>Meth.Enz</em>. Methods in Enzymology.&nbsp;<strong>266</strong>. pp.&nbsp;460&ndash;480.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1016%2FS0076-6879%2896%2966029-7" rel="nofollow">10.1016/S0076-6879(96)66029-7</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/International_Standard_Book_Number" title="International Standard Book Number">ISBN</a>&nbsp;<a href="https://en.wikipedia.org/wiki/Special:BookSources/9780121821678" title="Special:BookSources/9780121821678">9780121821678</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-hartmann_26-0" title="Jump up">^</a></strong>&nbsp;<cite>Hartmann AK (2002). &quot;Sampling rare events: statistics of local sequence alignments&quot;.&nbsp;<em>Phys. Rev. E</em>.&nbsp;<strong>65</strong>&nbsp;(5): 056102.&nbsp;<a href="https://en.wikipedia.org/wiki/ArXiv" title="ArXiv">arXiv</a>:<a href="https://arxiv.org/abs/cond-mat/0108201" rel="nofollow">cond-mat/0108201</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/Bibcode" title="Bibcode">Bibcode</a>:<a href="http://adsabs.harvard.edu/abs/2002PhRvE..65e6102H" rel="nofollow">2002PhRvE..65e6102H</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1103%2FPhysRevE.65.056102" rel="nofollow">10.1103/PhysRevE.65.056102</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/12059642" rel="nofollow">12059642</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-newberg_27-0" title="Jump up">^</a></strong>&nbsp;<cite>Newberg LA (2008).&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2737730" rel="nofollow">&quot;Significance of gapped sequence alignments&quot;</a>.&nbsp;<em>J Comput Biolo</em>.&nbsp;<strong>15</strong>&nbsp;(9): 1187&ndash;1194.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1089%2Fcmb.2008.0125" rel="nofollow">10.1089/cmb.2008.0125</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Central" title="PubMed Central">PMC</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2737730" rel="nofollow">2737730</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/18973434" rel="nofollow">18973434</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-eddy_28-0" title="Jump up">^</a></strong>&nbsp;<cite>Eddy SR; Rost, Burkhard (2008). Rost, Burkhard, ed.&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2396288" rel="nofollow">&quot;A probabilistic model of local sequence alignment that simplifies statistical significance estimation&quot;</a>.&nbsp;<em>PLoS Comput Biol</em>.&nbsp;<strong>4</strong>&nbsp;(5): e1000069.&nbsp;<a href="https://en.wikipedia.org/wiki/Bibcode" title="Bibcode">Bibcode</a>:<a href="http://adsabs.harvard.edu/abs/2008PLSCB...4E0069E" rel="nofollow">2008PLSCB...4E0069E</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1371%2Fjournal.pcbi.1000069" rel="nofollow">10.1371/journal.pcbi.1000069</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Central" title="PubMed Central">PMC</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2396288" rel="nofollow">2396288</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/18516236" rel="nofollow">18516236</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-bastien_29-0" title="Jump up">^</a></strong>&nbsp;<cite>Bastien O; Aude JC; Roy S; Marechal E (2004).&nbsp;<a href="http://bioinformatics.oxfordjournals.org/content/20/4/534.long" rel="nofollow">&quot;Fundamentals of massive automatic pairwise alignments of protein sequences: theoretical significance of Z-value statistics&quot;</a>.&nbsp;<em>Bioinformatics</em>.&nbsp;<strong>20</strong>(4): 534&ndash;537.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1093%2Fbioinformatics%2Fbtg440" rel="nofollow">10.1093/bioinformatics/btg440</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/14990449" rel="nofollow">14990449</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-agrawal11_30-0" title="Jump up">^</a></strong>&nbsp;<cite>Agrawal A; Huang X (2011).&nbsp;<a href="https://archive.is/20130415004914/http://ieeexplore.ieee.org/xpls/abs_all.jsp?arnumber=5276793" rel="nofollow">&quot;Pairwise Statistical Significance of Local Sequence Alignment Using Sequence-Specific and Position-Specific Substitution Matrices&quot;</a>.&nbsp;<em>IEEE/ACM Transactions on Computational Biology and Bioinformatics</em>.&nbsp;<strong>8</strong>&nbsp;(1): 194&ndash;205.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1109%2FTCBB.2009.69" rel="nofollow">10.1109/TCBB.2009.69</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/21071807" rel="nofollow">21071807</a>. Archived from&nbsp;<a href="http://ieeexplore.ieee.org/xpls/abs_all.jsp?arnumber=5276793" rel="nofollow">the original</a>&nbsp;on 2013-04-15.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-agrawal08_31-0" title="Jump up">^</a></strong>&nbsp;<cite>Agrawal A; Brendel VP; Huang X (2008).&nbsp;<a href="https://archive.is/20130128163812/http://inderscience.metapress.com/content/1558538106522500/" rel="nofollow">&quot;Pairwise statistical significance and empirical determination of effective gap opening penalties for protein local sequence alignment&quot;</a>.&nbsp;<em>International Journal of Computational Biology and Drug Design</em>.&nbsp;<strong>1</strong>&nbsp;(4): 347&ndash;367.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1504%2FIJCBDD.2008.022207" rel="nofollow">10.1504/IJCBDD.2008.022207</a>. Archived from&nbsp;<a href="http://inderscience.metapress.com/content/1558538106522500/" rel="nofollow">the original</a>&nbsp;on 28 January 2013.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-NewbergLawrence2009_32-0" title="Jump up">^</a></strong>&nbsp;<cite>Newberg LA; Lawrence CE (2009).&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2858568" rel="nofollow">&quot;Exact Calculation of Distributions on Integers, with Application to Sequence Alignment&quot;</a>.&nbsp;<em>J Comput Biolo</em>.&nbsp;<strong>16</strong>&nbsp;(1): 1&ndash;18.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1089%2Fcmb.2008.0137" rel="nofollow">10.1089/cmb.2008.0137</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Central" title="PubMed Central">PMC</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2858568" rel="nofollow">2858568</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/19119992" rel="nofollow">19119992</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-33" title="Jump up">^</a></strong>&nbsp;<cite>Kim N; Lee C (2008).&nbsp;<em>Bioinformatics detection of alternative splicing</em>.&nbsp;<em>Methods Mol. Biol</em>. Methods in Molecular Biology&trade;.&nbsp;<strong>452</strong>. pp.&nbsp;179&ndash;97.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1007%2F978-1-60327-159-2_9" rel="nofollow">10.1007/978-1-60327-159-2_9</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/International_Standard_Book_Number" title="International Standard Book Number">ISBN</a>&nbsp;<a href="https://en.wikipedia.org/wiki/Special:BookSources/978-1-58829-707-5" title="Special:BookSources/978-1-58829-707-5">978-1-58829-707-5</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/18566765" rel="nofollow">18566765</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-34" title="Jump up">^</a></strong>&nbsp;<cite>Li JB, Levanon EY, Yoon JK, et al. (May 2009). &quot;Genome-wide identification of human RNA editing sites by parallel DNA capturing and sequencing&quot;.&nbsp;<em>Science</em>.&nbsp;<strong>324</strong>&nbsp;(5931): 1210&ndash;3.&nbsp;<a href="https://en.wikipedia.org/wiki/Bibcode" title="Bibcode">Bibcode</a>:<a href="http://adsabs.harvard.edu/abs/2009Sci...324.1210L" rel="nofollow">2009Sci...324.1210L</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1126%2Fscience.1170995" rel="nofollow">10.1126/science.1170995</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/19478186" rel="nofollow">19478186</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-35" title="Jump up">^</a></strong>&nbsp;<cite>Blazewicz J, Bryja M, Figlerowicz M, et al. (June 2009). &quot;Whole genome assembly from 454 sequencing output via modified DNA graph concept&quot;.&nbsp;<em>Comput Biol Chem</em>.&nbsp;<strong>33</strong>&nbsp;(3): 224&ndash;30.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1016%2Fj.compbiolchem.2009.04.005" rel="nofollow">10.1016/j.compbiolchem.2009.04.005</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/19477687" rel="nofollow">19477687</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-36" title="Jump up">^</a></strong>&nbsp;<cite>Duran C; Appleby N; Vardy M; Imelfort M; Edwards D; Batley J (May 2009). &quot;Single nucleotide polymorphism discovery in barley using autoSNPdb&quot;.&nbsp;<em>Plant Biotechnol. J</em>.&nbsp;<strong>7</strong>&nbsp;(4): 326&ndash;33.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1111%2Fj.1467-7652.2009.00407.x" rel="nofollow">10.1111/j.1467-7652.2009.00407.x</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/19386041" rel="nofollow">19386041</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-37" title="Jump up">^</a></strong>&nbsp;<cite>Abbott A.; Tsay A. (2000). &quot;Sequence Analysis and Optimal Matching Methods in Sociology, Review and Prospect&quot;.&nbsp;<em>Sociological Methods and Research</em>.&nbsp;<strong>29</strong>&nbsp;(1): 3&ndash;33.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1177%2F0049124100029001001" rel="nofollow">10.1177/0049124100029001001</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-Barzilay_38-0" title="Jump up">^</a></strong>&nbsp;<cite>Barzilay R; Lee L. (2002).&nbsp;<a href="http://www.cs.cornell.edu/home/llee/papers/gen-msa.pdf" rel="nofollow">&quot;Bootstrapping Lexical Choice via Multiple-Sequence Alignment&quot;</a>&nbsp;(PDF).&nbsp;<em>Proceedings of the Conference on Empirical Methods in Natural Language Processing (EMNLP)</em>.&nbsp;<strong>10</strong>: 164&ndash;171.&nbsp;<a href="https://en.wikipedia.org/wiki/ArXiv" title="ArXiv">arXiv</a>:<a href="https://arxiv.org/abs/cs/0205065" rel="nofollow">cs/0205065</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.3115%2F1118693.1118715" rel="nofollow">10.3115/1118693.1118715</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-39" title="Jump up">^</a></strong>&nbsp;<cite>Kondrak, Grzegorz (2002).&nbsp;<a href="http://www.cs.ualberta.ca/~kondrak/papers/thesis.pdf" rel="nofollow">&quot;Algorithms for Language Reconstruction&quot;</a>&nbsp;(PDF). University of Toronto, Ontario. Retrieved&nbsp;2007-01-21.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-prinzie_40-0" title="Jump up">^</a></strong>&nbsp;<cite>Prinzie A.; D. Van den Poel (2006).&nbsp;<a href="http://econpapers.repec.org/paper/rugrugwps/05_2F292.htm" rel="nofollow">&quot;Incorporating sequential information into traditional classification models by using an element/position-sensitive SAM&quot;</a>.&nbsp;<em>Decision Support Systems</em>.&nbsp;<strong>42</strong>(2): 508&ndash;526.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1016%2Fj.dss.2005.02.004" rel="nofollow">10.1016/j.dss.2005.02.004</a>.</cite>&nbsp;See also Prinzie and Van den Poel&#39;s paper&nbsp;<cite>Prinzie, A; Vandenpoel, D (2007).&nbsp;<a href="http://econpapers.repec.org/paper/rugrugwps/07_2F442.htm" rel="nofollow">&quot;Predicting home-appliance acquisition sequences: Markov/Markov for Discrimination and survival analysis for modeling sequential information in NPTB models&quot;</a>.&nbsp;<em>Decision Support Systems</em>.&nbsp;<strong>44</strong>&nbsp;(1): 28&ndash;45.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1016%2Fj.dss.2007.02.008" rel="nofollow">10.1016/j.dss.2007.02.008</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-41" title="Jump up">^</a></strong>&nbsp;<cite>EMBL-EBI.&nbsp;<a href="http://www.ebi.ac.uk/Tools/msa/clustalw2/" rel="nofollow">&quot;ClustalW2 &lt; Multiple Sequence Alignment &lt; EMBL-EBI&quot;</a>.&nbsp;<em>www.EBI.ac.uk</em>. Retrieved&nbsp;12 June&nbsp;2017.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-42" title="Jump up">^</a></strong>&nbsp;<a href="https://web.archive.org/web/20080918022531/http://tcoffee.vital-it.ch/cgi-bin/Tcoffee/tcoffee_cgi/index.cgi" rel="nofollow">T-coffee</a></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-43" title="Jump up">^</a></strong>&nbsp;<cite><a href="http://blast.ncbi.nlm.nih.gov/Blast.cgi" rel="nofollow">&quot;BLAST: Basic Local Alignment Search Tool&quot;</a>.&nbsp;<em>blast.ncbi.nlm.NIH.gov</em>. Retrieved&nbsp;12 June&nbsp;2017.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-44" title="Jump up">^</a></strong>&nbsp;<cite><a href="http://fasta.bioch.virginia.edu/fasta_www2/fasta_list2.shtml" rel="nofollow">&quot;UVA FASTA Server&quot;</a>.&nbsp;<em>fasta.bioch.Virginia.edu</em>. Retrieved&nbsp;12 June&nbsp;2017.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-thompson2_45-0" title="Jump up">^</a></strong>&nbsp;<cite>Thompson JD; Plewniak F; Poch O (1999).&nbsp;<a href="http://bioinformatics.oxfordjournals.org/cgi/pmidlookup?view=long&amp;pmid=10068696" rel="nofollow">&quot;BAliBASE: a benchmark alignment database for the evaluation of multiple alignment programs&quot;</a>.&nbsp;<em>Bioinformatics</em>.&nbsp;<strong>15</strong>&nbsp;(1): 87&ndash;8.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1093%2Fbioinformatics%2F15.1.87" rel="nofollow">10.1093/bioinformatics/15.1.87</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/10068696" rel="nofollow">10068696</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-46" title="Jump up">^</a></strong>&nbsp;<a href="https://web.archive.org/web/20121130084356/http://bips.u-strasbg.fr/fr/Products/Databases/BAliBASE/prog_scores.html" rel="nofollow">BAliBASE</a></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-thompson3_47-0" title="Jump up">^</a></strong>&nbsp;<cite>Thompson JD; Plewniak F; Poch O. (1999).&nbsp;<a href="http://nar.oxfordjournals.org/cgi/pmidlookup?view=long&amp;pmid=10373585" rel="nofollow">&quot;A comprehensive comparison of multiple sequence alignment programs&quot;</a>.&nbsp;<em>Nucleic Acids Res</em>.&nbsp;<strong>27</strong>&nbsp;(13): 2682&ndash;90.&nbsp;<a href="https://en.wikipedia.org/wiki/Digital_object_identifier" title="Digital object identifier">doi</a>:<a href="https://doi.org/10.1093%2Fnar%2F27.13.2682" rel="nofollow">10.1093/nar/27.13.2682</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Central" title="PubMed Central">PMC</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC148477" rel="nofollow">148477</a>.&nbsp;<a href="https://en.wikipedia.org/wiki/PubMed_Identifier" title="PubMed Identifier">PMID</a>&nbsp;<a href="https://www.ncbi.nlm.nih.gov/pubmed/10373585" rel="nofollow">10373585</a>.</cite></li> <li><strong><a href="https://en.wikipedia.org/wiki/Sequence_alignment#cite_ref-48" title="Jump up">^</a></strong>&nbsp;<cite><a href="http://3d-alignment.eu/" rel="nofollow">&quot;Multiple sequence alignment: Strap&quot;</a>.&nbsp;<em>3d-alignment.eu</em>. Retrieved&nbsp;12 June&nbsp;2017.</spancite></li>
</ol>