-
Notifications
You must be signed in to change notification settings - Fork 1
/
identity.html
183 lines (180 loc) · 6.63 KB
/
identity.html
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01//EN" "http://www.w3.org/TR/html4/strict.dtd">
<html lang="en">
<head>
<meta content="en-us" http-equiv="Content-Language"/>
<meta content="text/html; charset=utf-8" http-equiv="Content-Type"/>
<meta content="no-cache, no-store, must-revalidate" http-equiv="Cache-Control"/>
<meta content="no-cache" http-equiv="Pragma"/>
<meta content="0" http-equiv="Expires"/>
<title>
allpairs_global command
</title>
<link href="stylesx.css" rel="stylesheet" type="text/css"/>
<style type="text/css">
body.c4 {background-color:#c0c0c0;}
div.c3 {position:absolute; top:45px; left:20px; width:830px; background-color:#ffffff; border-width:10px; border-style:solid;border-color:white;}
span.c2 {font-weight: bold}
div.c1 {position:absolute; top:10px; left:20px; width:850px; height:60px;}
.TopButtonPara { color:white; background-color:rgb(50,100,150); border-color:rgb(50,100,150); font-family:Arial, Helvetica, sans-serif; font-weight:normal; font-size:9pt; text-align:center; border-width:4px; border-style:solid; }
.TopButton { color:white; }
a.TopButton:link { text-decoration:none; }
a.TopButton:visited { text-decoration:none; }
a.TopButton:hover { color:orange; }
.NewButtonPara { color:white; background-color:rgb(50,100,150); border-color:rgb(50,100,150); font-family:Arial, Helvetica, sans-serif; font-weight:normal; font-size:9pt; text-align:center; border-width:4px; border-style:solid; }
.NewButton { color:white; }
a.NewButton:link { text-decoration:none; }
a.NewButton:visited { text-decoration:none; }
a.NewButton:hover { color:orange; }
.SideButtonPara { color:white; font-family:Arial, Helvetica, sans-serif; font-size:9pt; font-weight:normal; text-align:center; line-height:18px; }
.SideButton { color:white; }
a.SideButton:link { text-decoration:none; }
a.SideButton:visited { text-decoration:none; }
a.SideButton:hover { color:orange; }
</style>
</head>
<body style="background-color:#c0c0c0;">
<div>
<a href="https://drive5.com/usearch">
<img alt="USEARCH v12" src="usearch12_banner.jpg" style="position:absolute; top:40px; left:10px; padding:0px; border:0px;"/>
</a>
</div>
<div style="position:absolute; top:115px; left:10px; width:850px; background-color:#ffffff; min-height:500px">
<div style="position:relative; float:left; background-color:#696969; width:125px; left: 0px; min-height:500px; padding:5px; height: 125px;">
<div class="SideButtonPara" style="text-align:center; padding-top:5px;">
<a class="SideButton" href="index.html">
Docs home
</a>
<br/>
<hr style="border:0; border-bottom: 1px solid white;"/>
<a class="SideButton" href="cmds.html">
Commands
</a>
<br/>
<a class="SideButton" href="topics.html">
Topics
</a>
<br/>
<a class="SideButton" href="citation.html">
Publications
</a>
<br/>
</div>
</div>
<div class="ManText" style="left:20px; position: absolute; left:135px; width:695px; background-color:white; padding:10px">
<h1>
Pair-wise sequence identity
</h1>
<span class="ManText">
<b>
See also
<br/>
</b>
<a href="opt_id.html">
-id option
</a>
<b>
<br/>
</b>
<a href="more_clusters.html">
Identity and clustering
</a>
<br/>
<a href="accept_options.html">
Accept options
</a>
<br/>
<a href="aln_params.html">
Alignment parameters
</a>
</span>
<p>
<span class="ManText">
USEARCH uses the BLAST definition of identity, which is the number of identities divided by the number of alignment columns. This is the definition used by most bioinformatics programs.
</span>
</p>
<p>
<span class="ManText">
See
<a href="opt_id.html">
-id option
</a>
for rules that determine when a pair of letters match.
</span>
</p>
<p>
<span class="ManText">
With a
<a href="local_global.html">
global alignment
</a>
, columns containing
<a href="terminal_gaps.html">
terminal gaps
</a>
are discarded before calculating identity. Internal gaps always count as differences.
</span>
</p>
<p class="ManText">
<b>
<span>
Different programs may report different identities for the same sequences
</span>
</b>
<span class="c2">
<b>
<br/>
</b>
The pair-wise identity between two sequences depends on the alignment and the definition of identity. Alignments vary due to the use of different parameters such as gap penalties and substitution scores. Definitions of identity used by different programs may vary depending on the treatment of gaps and detailed rules for deciding whether letters match (e.g., upper- and lower-case and wildcards).
</span>
</p>
<p class="ManText">
<b>
<span>
USEARCH, BLAST and CD-HIT definitions of %id
</span>
</b>
<span>
<br/>
An "identity" is a column with two matching letters.
</span>
</p>
<p class="ManText">
<span class="c2">
USEARCH and BLAST = Identities / Columns
</span>
</p>
<p class="ManText">
<span class="c2">
CD-HIT = Identities / (Length of shorter sequence)
</span>
</p>
<p class="ManText">
<b>
<span class="c2">
Problems with the CD-HIT definition
</span>
</b>
<span>
<br/>
For historical reasons, versions 5 and earlier of USEARCH used the CD-HIT definition of identity. In versions 6 and later, the BLAST definition is used. I made the change because I felt the CD-HIT definition had several important weaknesses. The CD-HIT definition is not symmetrical between the longer and shorter sequence. Gaps in the longer sequence reduce %id but gaps in the shorter sequence do not. Gappier alignments therefore tend to have higher identities according to CD-HIT compared to other methods, and the CD-HIT %id correlates less well with evolutionary distance. A measure of %id that counts gaps as differences is more robust against the choice of alignment parameters (gap penalties and substitution matrices).
</span>
</p>
<p>
<span>
<a href="/usearch/cdhitv4_alignment86.html">
Example where CD-HIT id is 97% and USEARCH id is 86%
</a>
</span>
</p>
<p>
<span>
<a href="/usearch/cdhitv4_alignment95.html">
Example where CD-HIT id is 97% and USEARCH id is 95%
</a>
<br/>
</span>
</p>
</div>
</div>
</body>
</html>