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
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
|
!!SEQUENCE_LIST 1.0
(Peptide) FASTA of: test.gcg from: 1 to: 146 August 25, 2003 13:25
REFORMAT of: b124_sp.pep check: -1 from: 1 to: 146 January 28, 1999 16:22
(No documentation)
TO: PIR:* Sequences: 283,308 Symbols: 96,168,669 Word Size: 2
Databases searched:
NBRF, Release 76.0, Released on 31Mar2003, Formatted on 7Apr2003
Scoring matrix: GenRunData:blosum50.cmp
Variable pamfactor used
Gap creation penalty: 12 Gap extension penalty: 2
Histogram Key:
Each histogram symbol represents 474 search set sequences
Each inset symbol represents 4 search set sequences
z-scores computed from opt scores
z-score obs exp
(=) (*)
< 20 789 0:==
22 0 0:
24 4 0:=
26 8 6:*
28 9 64:*
30 101 390:*
32 407 1509:= *
34 2185 4092:===== *
36 7555 8404:================ *
38 16600 13889:=============================*======
40 25000 19373:========================================*============
42 27813 23681:=================================================*=========
44 28394 26123:=======================================================*====
46 26152 26607:========================================================*
48 23191 25473:================================================= *
50 20419 23244:============================================ *
52 18108 20435:======================================= *
54 15701 17455:================================== *
56 13874 14581:==============================*
58 11026 11970:======================== *
60 9392 9697:====================*
62 7678 7774:================*
64 6295 6183:=============*
66 4986 4887:==========*
68 3909 3844:========*
70 3131 3012:======*
72 2497 2354:====*=
74 1858 1835:===*
76 1469 1428:===*
78 1160 1110:==*
80 845 862:=*
82 665 659:=*
84 515 522:=*
86 376 404:*
88 261 313:*
90 225 242:*
92 157 187:* :=======================================*
94 132 145:* :================================= *
96 93 112:* :======================== *
98 63 87:* :================ *
100 73 67:* :================*==
102 44 52:* :=========== *
104 32 40:* :======== *
106 27 31:* :=======*
108 18 24:* :=====*
110 18 19:* :====*
112 11 14:* :===*
114 11 11:* :==*
116 10 9:* :==*
118 8 7:* :=*
>120 13 5:* :=*==
Joining threshold: 36, opt. threshold: 24, opt. width: 16, reg.-scaled
The best scores are: init1 initn opt z-sc E(283250)..
PIR2:S44629 Begin: 342 End: 470
! F22B7.10 protein - Caenorhabditis e... 108 143 241 304.1 1.1e-09
PIR1:WMBELM Begin: 307 End: 385
! membrane protein LMP-2A - human her... 59 91 99 130.6 5.1
PIR2:AG0762 Begin: 63 End: 144
! probable membrane protein STY2265 [... 65 65 96 128.9 6.4
PIR2:B83179 Begin: 9 End: 86
! hypothetical protein PA3730 [import... 40 40 92 127.0 8.2
\\End of List
test.gcg
PIR2:S44629
P1;S44629 - F22B7.10 protein - Caenorhabditis elegans
C;Species: Caenorhabditis elegans
C;Date: 20-Feb-1995 #sequence_revision 20-Feb-1995 #text_change 04-Mar-2000
C;Accession: S44629
R;Anderson, K.
submitted to the EMBL Data Library, March 1993 . . .
SCORES Init1: 108 Initn: 143 Opt: 241 z-score: 304.1 E(): 1.1e-09
>>PIR2:S44629 (628 aa)
initn: 143 init1: 108 opt: 241 Z-score: 304.1 expect(): 1.1e-09
Smith-Waterman score: 241; 32.6% identity in 135 aa overlap
(3-135:342-470)
10 20 30
test.gcg VXCAAEFDFMEKETPLRYTKTLLLPVVLVVFV
|:|||||:: | : |||:|::|: :|
S44629 GLGIEDDAHIFDILRSKFTSFANFHTRLYTCSAEFDFIQYSTIEKLCGTLLIPLALISLV
320 330 340 350 360 370
40 50 60 70 80 90
test.gcg AIVRKIISDMWGVLAKQQTHVRKHQFDHGELVYHALQLLAYTALGILIMRLKLFLTPYMC
::| ::::: ::| ::: :: ::||::|:::|| |::::||||||||:||::|
S44629 TFVFNFVKNT-NLLWRNSEEIG----ENGEILYNVVQLCCSTVMAFLIMRLKLFMTPHLC
380 390 400 410 420
100 110 120 130 140
test.gcg VMASLICSRQLFG--WLFCKVHPGAIVFVILAAMSIQGSANLQTQWKSTASLALET
::|:|: : :|:| : :: :|:| || | : :| |:: |
S44629 IVAALFANSKLLGGDRISKTIRVSALVGVI-AILFYRGIPNIRQQLNVKGEYSNPDQEML
430 440 450 460 470 480
S44629 FDWIQHNTKQDAVFAGTMPVMANVKLTTLRPIVNHPHYEHVGIRERTLKVYSMFSKKPIA
490 500 510 520 530 540
test.gcg
PIR1:WMBELM
P1;WMBELM - membrane protein LMP-2A - human herpesvirus 4
N;Contains: membrane protein LMP-2B
C;Species: human herpesvirus 4, Epstein-Barr virus
A;Note: host Homo sapiens (man)
C;Date: 31-Dec-1989 #sequence_revision 31-Dec-1989 #text_change 16-Jul-1999
C;Accession: A30178; B30178; S00392 . . .
SCORES Init1: 59 Initn: 91 Opt: 99 z-score: 130.6 E(): 5.1
>>PIR1:WMBELM (497 aa)
initn: 91 init1: 59 opt: 99 Z-score: 130.6 expect(): 5.1
Smith-Waterman score: 99; 32.9% identity in 79 aa overlap
(67-141:307-385)
40 50 60 70 80 90
test.gcg KIISDMWGVLAKQQTHVRKHQFDHGELVYHALQLLAYTALGILIMRLKLFLTPYMCVMAS
|| ||| || | : ::| :::
WMBELM MTLLLLAFVLWLSSPGGLGTLGAALLTLAAALALLASLILGTLNLTTMFLLMLLWTLVVL
280 290 300 310 320 330
100 110 120 130 140
test.gcg LICSR----QLFGWLFCKVHPGAIVFVILAAMSIQGSANLQTQWKSTASLALET
|||| | |: :: |:::::||: | |:: |||::|| :|
WMBELM LICSSCSSCPLSKILLARLFLYALALLLLASALIAGGSILQTNFKSLSSTEFIPNLFCML
340 350 360 370 380 390
WMBELM LLIVAGILFILAILTEWGSGNRTYGPVFMCLGGLLTMVAGAVWLTVMSNTLLSAWILTAG
400 410 420 430 440 450
test.gcg
PIR2:AG0762
P1;AG0762 - probable membrane protein STY2265 [imported] - Salmonella enterica
subsp. enterica serovar Typhi (strain CT18)
C;Species: Salmonella enterica subsp. enterica serovar Typhi
A;Note: this species has also been called Salmonella typhi
C;Date: 09-Nov-2001 #sequence_revision 09-Nov-2001 #text_change 18-Nov-2002
C;Accession: AG0762
R;Parkhill, J.; Dougan, G.; James, K.D.; Thomson, N.R.; Pickard, D.; Wain, J.;
Churcher, C.; Mungall, K.L.; Bentley, S.D.; Holden, M.T.G.; Sebaihia, M.;
Baker, S.; Basham, D.; Brooks, K.; Chillingworth, T.; Connerton, P.; Cronin,
A.; Davis, P.; Davies, R.M.; Dowd, L.; White, N.; Farrar, J.; Feltwell, T.;
Hamlin, N.; Haque, A.; Hien, T.T.; Holroyd, S.; Jagels, K.; Krogh, A.; Larsen,
T.S.; Leather, S.; Moule, S.; O'Gaora, P
SCORES Init1: 65 Initn: 65 Opt: 96 z-score: 128.9 E(): 6.4
>>PIR2:AG0762 (352 aa)
initn: 65 init1: 65 opt: 96 Z-score: 128.9 expect(): 6.4
Smith-Waterman score: 96; 27.6% identity in 87 aa overlap
(61-137:63-144)
40 50 60 70 80
test.gcg FVAIVRKIISDMWGVLAKQQTHVRKHQFDHGELVYHALQLLAYT----ALGILIMRLKLF
|::| :|:: :: | |||:: :||:||
AG0762 TFLLVRLFSIPEGTWPLITLVVIMGPISFWGNVVPRAFERIGGTILGAALGLVALRLELF
40 50 60 70 80 90
90 100 110 120 130 140
test.gcg LTPYM---CVMASLICSRQLFGWLFCKVHP--GAIVFVILAAMSIQGSANLQTQ-WKSTA
| | |::| ::| ||| :| : :: : ||:: :::::| |::
AG0762 SLPLMLVWCAIAMFLC-----GWLALGKKPYQALLIGITLAVVVGAPAGDMNTALWRGGD
100 110 120 130 140
test.gcg SLALET
AG0762 VILGALLAMLFTGIWPQRAFLHWRIQLAHCVTAYNRVYQAALSPNLLERPRLDKYLQRLL
150 160 170 180 190 200
test.gcg
PIR2:B83179
P1;B83179 - hypothetical protein PA3730 [imported] - Pseudomonas aeruginosa
(strain PAO1)
C;Species: Pseudomonas aeruginosa
C;Date: 15-Sep-2000 #sequence_revision 15-Sep-2000 #text_change 31-Dec-2000
C;Accession: B83179
R;Stover, C.K.; Pham, X.Q.; Erwin, A.L.; Mizoguchi, S.D.; Warrener, P.; Hickey,
M.J.; Brinkman, F.S.L.; Hufnagle, W.O.; Kowalik, D.J.; Lagrou, M.; Garber,
R.L.; Goltry, L.; Tolentino, E.; Westbrook-Wadman, S.; Yuan, Y.; Brody, L.L.;
Coulter, S.N.; Folger, K.R.; Kas, A.; Larbig, K.; Lim, R.M.; Smith, K.A.;
Spencer, D.H.; Wong, G.K.S.; Wu, Z.; Paulsen, I.T.; Reizer, J.; Saier, M.H.;
Hancock, R.E.W.; Lory, S.; Olson, M.V.
Nature 406, 959-964, 2000 . . .
SCORES Init1: 40 Initn: 40 Opt: 92 z-score: 127.0 E(): 8.2
>>PIR2:B83179 (213 aa)
initn: 40 init1: 40 opt: 92 Z-score: 127.0 expect(): 8.2
Smith-Waterman score: 92; 28.4% identity in 88 aa overlap
(22-109:9-86)
10 20 30 40 50 60
test.gcg VXCAAEFDFMEKETPLRYTKTLLLPVVLVVFVAIVRKIISDMWGVLAKQQTHVRKHQFDH
| :|:|| |: |: | :||::| :::: ::|
B83179 MEGFLQTALSFPTVLFSFLLILAII---YWGIVALGMVEIDVLDLDA
10 20 30 40
70 80 90 100 110 120
test.gcg GELVYHALQLLAYTALGILIMRLKLFLTPYMCVMASLICSRQLFGWLFCKVHPGAIVFVI
:| | | :|: |: :||| :| |:: | ::|:|::|
B83179 ESVVDGAGQA---EGLAALLAKLKLNGVPVTLVLTLL----SFFAWFLCYFVQLWLLSAL
50 60 70 80 90
130 140
test.gcg LAAMSIQGSANLQTQWKSTASLALET
B83179 PLGWLRYPLGAVVAVGALFLAAPLAATLCRPLRPLFRKLESTSSKSVLGQVAVVRSGRVT
100 110 120 130 140 150
! Distributed over 1 thread.
! Start time: Mon Aug 25 13:23:54 2003
! Completion time: Mon Aug 25 13:25:12 2003
! CPU time used:
! Database scan: 0:01:34.1
! Post-scan processing: 0:00:00.6
! Total CPU time: 0:01:34.7
! Output File: test.fasta
|
