Number of input and output field
| Number of input field |
: |
Two |
| Number of output field |
: |
One or More |
This function is used to align two deoxyribonucleic acid (DNA) or protein sequences in order to identify regions of similarity between them. These regions of similarity are useful for understanding the functional, structural, and evolutionary relationships between sequences.
Two commonly used sequence alignment algorithms are global alignment and local alignment:
- Global alignment:
Global alignment compares two sequences by aligning their entire lengths and maximizing the overall similarity. This method is typically used when the sequences are similar and of approximately equal length.
- Local alignment:
Local alignment focuses on aligning only the regions with the highest density of matches, rather than the entire sequence length. This approach is useful for identifying short, conserved regions within protein or nucleotide sequences.
Parameters of the SequenceAlignment Function:
- Alignment Algorithm
Specifies the type of alignment to use: Global Alignment or Local Alignment.
- Similarity Matrix
A similarity (substitution) matrix used to score residue similarity. Available matrices include:
DEFAULT, DNA_1, DNA_2, BLOSUM62, BLOSUM45, BLOSUM80, BLOSUM90, PAM10, PAM30, PAM70, PAM120, and PAM250
(see Substitution Matrices for more detail).
When using the DEFAULT similarity matrix, residue similarity scoring is based on the Match Score and the Mismatch Score parameters.
- Match Score
The score assigned when two characters in the sequences match.
- Mismatch Score
The score assigned when two characters in the sequences do not match.
- Opening Gap Penalty
The penalty applied to the first gap in a gap cluster, which is the gap that immediately follows an amino acid or nucleotide in the aligned sequence.
- Extension Gap Penalty
The penalty applied to each additional gap following the opening gap.
- Selected Outputs
Specifies which alignment results will be forwarded to the output. Available output options include:
- Identities(%)
- Optimum Score
- Total Matches
- Total Mismatches
- Total Opening Gaps
- Total Extension Gaps
- Aligned Seq 1
- Aligned Seq 2
- Pairwise Alignment Result
If more than one output item is selected, the results are returned in the order listed above.
Use Case : Sample of using SequenceAlignment function
Sequence 1 :
GACCATATGATTGACGCCTATGTCAATCTCTACACTACATTGCTGGAAAGCAAATCCTGAGAGATGCTACCCCCGCCGTTGCTGCGGGGGCCAACGCGTTAATGCCGATTCTTCAGATTATCAATCACTTCTCCGAGATCCAGCCCTTGATCCTGCAACAACACCAGCAGGTGATACATCAAATCAGATGCCTCGTTGGTCAGTTCAAAGCGGTCATGTACTGTCGCTGCCAGTGCGGTTT
Sequence 2 :
GGTKAAATGCGTAGAGATCTGGAGGAATACCGGTGGCGAAGGCGGCCCYCTGGACAGATACTGACACTCAGATGCGTAAGCGTGGGGAGCAAACAGGATTAGATACCCTNGTAGTCCACGCCGTAAACGATGTCTACTTGGAGGTTGTGGCCTTNAGCCGTGGCTTTCGGACGTAACGCGTTAAGTAGACCGCCTGGGGAGTACGGTYGCAAGA
Case 1 : Parameter
Alignment Algorithm : Global Alignment
Similarity Matrix : DEFAULT
Match Score : 1
Mismatch Score : -2
Opening Gap Penalty : -5
Extension Gap Penalty : -2
Results of SequneceAlignment Function :
Identities (%):
41.05691
Optimum Score :
-249
Total Matches :
101
Total Mismatches :
108
Total Opening Gaps :
12
Total Gaps :
37
Aligned Seq 1 :
GACCATATGATTGACGCCTATGTCAATCTCTACACTACATTGCTGGAAAGCAAATCCTGAGAGATGCT-ACCCCC----GCCGTTGCTGCGGGGGCCAACGCGTTAATGCCGATTCTTCAGATTATCAATCACTTCTCCGAGATCCAGCCCTTGATCCTGCAACAACACCAGCAGGTGATACATCAAATCAGATGCCTCGTTGGTCAGTTCAAAGCGGTCATGTACTGTCGCTGCCAGTGCGGTTT
Aligned Seq 2 :
GGTKAAATGCGTAGAG---ATCTGGAGGAATAC---CGGTGGCGAAGGCGGCCCYCTGGACAGATACTGACACTCAGATGCGTAAGCGTGGGGAGCAAACAGGATTA----GATAC-CCTNGTAGTCCACGCCGTAAACGATGTCTA--CTTGGAGGTTGTGGC--CTTNAGCCGTGGCT-------TTCGGACGTAACG-CGTTAAGT--AGACCGCCTGGGGAGTACGGTYGCAAGA-------
Pairwise Alignment Result:
S1 1 GACCATATGATTGACGCCTATGTCAATCTCTACACTACATTGCTGGAAAGCAAATCCTGA 60
| | ||| | | || | | ||| | || | | ||
S2 1 GGTKAAATGCGTAGAG---ATCTGGAGGAATAC---CGGTGGCGAAGGCGGCCCYCTGGA 54
S1 61 GAGATGCT-ACCCCC----GCCGTTGCTGCGGGGGCCAACGCGTTAATGCCGATTCTTCA 115
|||| || || | | || || ||| || ||| | | | ||| | |
S2 55 CAGATACTGACACTCAGATGCGTAAGCGTGGGGAGCAAACAGGATTA----GATAC-CCT 109
S1 116 GATTATCAATCACTTCTCCGAGATCCAGCCCTTGATCCTGCAACAACACCAGCAGGTGAT 175
| || | | | ||| || | | | || || | | ||| | | |
S2 110 NGTAGTCCACGCCGTAAACGATGTCTA--CTTGGAGGTTGTGGC--CTTNAGCCGTGGCT 165
S1 176 ACATCAAATCAGATGCCTCGTTGGTCAGTTCAAAGCGGTCATGTACTGTCGCTGCCAGTG 235
|| || | || | | ||| | | || | | | | || ||
S2 166 -------TTCGGACGTAACG-CGTTAAGT--AGACCGCCTGGGGAGTACGGTYGCAAGA- 214
S1 236 CGGTTT 241
S2 215 ------ 214
Case 2 : Parameter
Alignment Algorithm : Local Alignment
Similarity Matrix : DEFAULT
Match Score : 1
Mismatch Score : -2
Opening Gap Penalty : -5
Extension Gap Penalty : -2
Results of SequenceAlignment Function :
Identities (%):
100.0
Optimum Score :
10
Total Matches :
10
Total Mismatches :
0
Total Opening Gaps :
0
Total Gaps :
0
Aligned Seq 1 :
AACGCGTTAA
Aligned Seq 2 :
AACGCGTTAA
Pairwise Alignment Result:
s1 AACGCGTTAA
||||||||||
s2 AACGCGTTAA
Protein Residue | A | B | C | D | E | F | G | H | I | J | K | L | M | N | P | Q | R | S | T | V | W | X | Y | Z | * |
| A | 2 | 0 | -2 | 0 | 0 | -3 | 1 | -1 | -1 | -1 | -1 | -2 | -1 | 0 | 1 | 0 | -2 | 1 | 1 | 0 | -6 | -1 | -3 | 0 | -8 |
|---|
| B | 0 | 3 | -4 | 3 | 3 | -4 | 0 | 1 | -2 | -3 | 1 | -3 | -2 | 2 | -1 | 1 | -1 | 0 | 0 | -2 | -5 | -1 | -3 | 2 | -8 |
|---|
| C | -2 | -4 | 12 | -5 | -5 | -4 | -3 | -3 | -2 | -5 | -5 | -6 | -5 | -4 | -3 | -5 | -4 | 0 | -2 | -2 | -8 | -1 | 0 | -5 | -8 |
|---|
| D | 0 | 3 | -5 | 4 | 3 | -6 | 1 | 1 | -2 | -3 | 0 | -4 | -3 | 2 | -1 | 2 | -1 | 0 | 0 | -2 | -7 | -1 | -4 | 3 | -8 |
|---|
| E | 0 | 3 | -5 | 3 | 4 | -5 | 0 | 1 | -2 | -3 | 0 | -3 | -2 | 1 | -1 | 2 | -1 | 0 | 0 | -2 | -7 | -1 | -4 | 3 | -8 |
|---|
| F | -3 | -4 | -4 | -6 | -5 | 9 | -5 | -2 | 1 | 2 | -5 | 2 | 0 | -3 | -5 | -5 | -4 | -3 | -3 | -1 | 0 | -1 | 7 | -5 | -8 |
|---|
| G | 1 | 0 | -3 | 1 | 0 | -5 | 5 | -2 | -3 | -4 | -2 | -4 | -3 | 0 | 0 | -1 | -3 | 1 | 0 | -1 | -7 | -1 | -5 | 0 | -8 |
|---|
| H | -1 | 1 | -3 | 1 | 1 | -2 | -2 | 6 | -2 | -2 | 0 | -2 | -2 | 2 | 0 | 3 | 2 | -1 | -1 | -2 | -3 | -1 | 0 | 2 | -8 |
|---|
| I | -1 | -2 | -2 | -2 | -2 | 1 | -3 | -2 | 5 | 3 | -2 | 2 | 2 | -2 | -2 | -2 | -2 | -1 | 0 | 4 | -5 | -1 | -1 | -2 | -8 |
|---|
| J | -1 | -3 | -5 | -3 | -3 | 2 | -4 | -2 | 3 | 5 | -3 | 5 | 3 | -3 | -2 | -2 | -3 | -2 | -1 | 2 | -3 | -1 | -1 | -2 | -8 |
|---|
| K | -1 | 1 | -5 | 0 | 0 | -5 | -2 | 0 | -2 | -3 | 5 | -3 | 0 | 1 | -1 | 1 | 3 | 0 | 0 | -2 | -3 | -1 | -4 | 0 | -8 |
|---|
| L | -2 | -3 | -6 | -4 | -3 | 2 | -4 | -2 | 2 | 5 | -3 | 6 | 4 | -3 | -3 | -2 | -3 | -3 | -2 | 2 | -2 | -1 | -1 | -3 | -8 |
|---|
| M | -1 | -2 | -5 | -3 | -2 | 0 | -3 | -2 | 2 | 3 | 0 | 4 | 6 | -2 | -2 | -1 | 0 | -2 | -1 | 2 | -4 | -1 | -2 | -2 | -8 |
|---|
| N | 0 | 2 | -4 | 2 | 1 | -3 | 0 | 2 | -2 | -3 | 1 | -3 | -2 | 2 | 0 | 1 | 0 | 1 | 0 | -2 | -4 | -1 | -2 | 1 | -8 |
|---|
| P | 1 | -1 | -3 | -1 | -1 | -5 | 0 | 0 | -2 | -2 | -1 | -3 | -2 | 0 | 6 | 0 | 0 | 1 | 0 | -1 | -6 | -1 | -5 | 0 | -8 |
|---|
| Q | 0 | 1 | -5 | 2 | 2 | -5 | -1 | 3 | -2 | -2 | 1 | -2 | -1 | 1 | 0 | 4 | 1 | -1 | -1 | -2 | -5 | -1 | -4 | 3 | -8 |
|---|
| R | -2 | -1 | -4 | -1 | -1 | -4 | -3 | 2 | -2 | -3 | 3 | -3 | 0 | 0 | 0 | 1 | 6 | 0 | -1 | -2 | 2 | -1 | -4 | 0 | -8 |
|---|
| S | 1 | 0 | 0 | 0 | 0 | -3 | 1 | -1 | -1 | -2 | 0 | -3 | -2 | 1 | 1 | -1 | 0 | 2 | 1 | -1 | -2 | -1 | -3 | 0 | -8 |
|---|
| T | 1 | 0 | -2 | 0 | 0 | -3 | 0 | -1 | 0 | -1 | 0 | -2 | -1 | 0 | 0 | -1 | -1 | 1 | 3 | 0 | -5 | -1 | -3 | -1 | -8 |
|---|
| V | 0 | -2 | -2 | -2 | -2 | -1 | -1 | -2 | 4 | 2 | -2 | 2 | 2 | -2 | -1 | -2 | -2 | -1 | 0 | 4 | -6 | -1 | -2 | -2 | -8 |
|---|
| W | -6 | -5 | -8 | -7 | -7 | 0 | -7 | -3 | -5 | -3 | -3 | -2 | -4 | -4 | -6 | -5 | 2 | -2 | -5 | -6 | 17 | -1 | 0 | -6 | -8 |
|---|
| X | -1 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | -1 | -8 |
|---|
| Y | -3 | -3 | 0 | -4 | -4 | 7 | -5 | 0 | -1 | -1 | -4 | -1 | -2 | -2 | -5 | -4 | -4 | -3 | -3 | -2 | 0 | -1 | 10 | -4 | -8 |
|---|
| Z | 0 | 2 | -5 | 3 | 3 | -5 | 0 | 2 | -2 | -2 | 0 | -3 | -2 | 1 | 0 | 3 | 0 | 0 | -1 | -2 | -6 | -1 | -4 | 3 | -8 |
|---|
| * | -8 | -8 | -8 | -8 | -8 | -8 | -8 | -8 | -8 | -8 | -8 | -8 | -8 | -8 | -8 | -8 | -8 | -8 | -8 | -8 | -8 | -8 | -8 | -8 | 1 |
|---|