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
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
|
# hmmsearch :: search profile(s) against a sequence database
# HMMER 3.0 (March 2010); http://hmmer.org/
# Copyright (C) 2010 Howard Hughes Medical Institute.
# Freely distributed under the GNU General Public License (GPLv3).
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# query HMM file: mult.hmm
# target sequence database: /home/bow/db/hmmer/uniprot_sprot.fasta
# output directed to file: hmmer_cases/text_hmmsearch_mult.out
# per-seq hits tabular output: hmmer_cases/tab_hmmsearch_mult.out
# per-dom hits tabular output: hmmer_cases/domtab_hmmsearch_mult.out
# sequence reporting threshold: score >= 480
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: globins4 [M=149]
Scores for complete sequences (score includes all domains):
--- full sequence --- --- best 1 domain --- -#dom-
E-value score bias E-value score bias exp N Sequence Description
------- ------ ----- ------- ------ ----- ---- -- -------- -----------
[No hits detected that satisfy reporting thresholds]
Domain annotation for each sequence (and alignments):
[No targets detected that satisfy reporting thresholds]
Internal pipeline statistics summary:
-------------------------------------
Query model(s): 1 (149 nodes)
Target sequences: 535698 (190107059 residues)
Passed MSV filter: 20698 (0.0386374); expected 10714.0 (0.02)
Passed bias filter: 16985 (0.0317063); expected 10714.0 (0.02)
Passed Vit filter: 2311 (0.004314); expected 535.7 (0.001)
Passed Fwd filter: 1104 (0.00206086); expected 5.4 (1e-05)
Initial search space (Z): 535698 [actual number of targets]
Domain search space (domZ): 0 [number of targets reported over threshold]
# CPU time: 14.07u 0.15s 00:00:14.22 Elapsed: 00:00:03.97
# Mc/sec: 7135.00
//
Query: Pkinase [M=260]
Accession: PF00069.17
Description: Protein kinase domain
Scores for complete sequences (score includes all domains):
--- full sequence --- --- best 1 domain --- -#dom-
E-value score bias E-value score bias exp N Sequence Description
------- ------ ----- ------- ------ ----- ---- -- -------- -----------
8.4e-147 492.3 0.0 1.1e-72 249.3 0.0 2.1 2 sp|Q9WUT3|KS6A2_MOUSE Ribosomal protein S6 kinase alpha-2 OS
2.7e-146 490.7 0.0 8.7e-73 249.7 0.0 2.1 2 sp|Q63531|KS6A1_RAT Ribosomal protein S6 kinase alpha-1 OS
2.2e-145 487.7 0.0 4.6e-72 247.4 0.0 2.1 2 sp|Q15418|KS6A1_HUMAN Ribosomal protein S6 kinase alpha-1 OS
2.5e-145 487.5 0.0 7.5e-72 246.7 0.0 2.1 2 sp|P18652|KS6AA_CHICK Ribosomal protein S6 kinase 2 alpha OS
2.4e-144 484.2 0.0 2.4e-70 241.7 0.0 2.1 2 sp|Q15349|KS6A2_HUMAN Ribosomal protein S6 kinase alpha-2 OS
5e-144 483.2 0.0 3.9e-70 241.0 0.0 2.2 2 sp|P51812|KS6A3_HUMAN Ribosomal protein S6 kinase alpha-3 OS
5e-144 483.2 0.0 3.9e-70 241.0 0.0 2.2 2 sp|P18654|KS6A3_MOUSE Ribosomal protein S6 kinase alpha-3 OS
Domain annotation for each sequence (and alignments):
>> sp|Q9WUT3|KS6A2_MOUSE Ribosomal protein S6 kinase alpha-2 OS=Mus musculus GN=Rps6ka2 PE=2 SV=1
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 241.2 0.0 4.6e-75 3.5e-70 1 260 [] 59 318 .. 59 318 .. 0.95
2 ! 249.3 0.0 1.5e-77 1.1e-72 1 260 [] 415 672 .. 415 672 .. 0.97
Alignments for each domain:
== domain 1 score: 241.2 bits; conditional E-value: 4.6e-75
EEEEEEEEEETTEEEEEEEE...TTTTEEEEEEEEEHHHCCCCCCHHHHHHHHHHHHHSSSSB--EEEEEEETTEEEEEEE--TS-BH CS
Pkinase 1 yelleklGsGsfGkVykakk...kktgkkvAvKilkkeeekskkektavrElkilkklsHpnivkllevfetkdelylvleyveggdl 85
+ell+ lG+Gs+GkV+ ++k ++ g+ +A+K+lkk + k + + +++ E il++++Hp+ivkl+ +f+t+ +lyl+l++++ggdl
sp|Q9WUT3|KS6A2_MOUSE 59 FELLKVLGQGSYGKVFLVRKvtgSDAGQLYAMKVLKKATLKVRDRVRSKMERDILAEVNHPFIVKLHYAFQTEGKLYLILDFLRGGDL 146
67899**********7666611155667*****************99999************************************** PP
HHHHHHHHST-HHHHHHHHHHHHHHHHHHHHTTEE-S--SGGGEEEETTTEEEE--GTT.E..EECSS-C-S--S-GGGS-HHHHCCS CS
Pkinase 86 fdllkkegklseeeikkialqilegleylHsngiiHrDLKpeNiLldkkgevkiaDFGlakkleksseklttlvgtreYmAPEvllka 173
f+ l+ke ++ee++k+++ +++ +l++lH gii+rDLKpeNiLld++g++ki+DFGl+k+++ +++++++++gt eYmAPEv++ +
sp|Q9WUT3|KS6A2_MOUSE 147 FTRLSKEVMFTEEDVKFYLAELALALDHLHGLGIIYRDLKPENILLDEEGHIKITDFGLSKEATDHDKRAYSFCGTIEYMAPEVVN-R 233
**************************************************************************************.* PP
-CTHHHHHHHHHHHHHHHHHHSS-TTSSSHHCCTHHHHSSHHH......TTS.....HHHHHHHHHHT-SSGGGSTT.....HHHHHT CS
Pkinase 174 keytkkvDvWslGvilyelltgklpfsgeseedqleliekilkkkleedepkssskseelkdlikkllekdpakRlt.....aeeilk 256
+++t+++D+Ws+Gv+++e+ltg+lpf+g+ d++e+++ ilk kl ++ s +e+++l++ l++++p +Rl +eei++
sp|Q9WUT3|KS6A2_MOUSE 234 RGHTQSADWWSFGVLMFEMLTGSLPFQGK---DRKETMALILKAKLGMPQFLS----AEAQSLLRALFKRNPCNRLGagvdgVEEIKR 314
*****************************...999999999999999998866....99******************9999999**** PP
SGGG CS
Pkinase 257 hpwl 260
hp++
sp|Q9WUT3|KS6A2_MOUSE 315 HPFF 318
*997 PP
== domain 2 score: 249.3 bits; conditional E-value: 1.5e-77
EEEEEEEEEETTEEEEEEEETTTTEEEEEEEEEHHHCCCCCCHHHHHHHHHHHHH.SSSSB--EEEEEEETTEEEEEEE--TS-BHHH CS
Pkinase 1 yelleklGsGsfGkVykakkkktgkkvAvKilkkeeekskkektavrElkilkkl.sHpnivkllevfetkdelylvleyveggdlfd 87
ye++e +G Gs+++++++++k t ++AvKi++k++++ ++ E++il + +Hpni++l +v+ + +++ylv+e+++gg+l d
sp|Q9WUT3|KS6A2_MOUSE 415 YEIKEDIGVGSYSVCKRCVHKATDAEYAVKIIDKSKRDPSE------EIEILLRYgQHPNIITLKDVYDDGKYVYLVMELMRGGELLD 496
7899***********************************98......9*******99******************************* PP
HHHHHHST-HHHHHHHHHHHHHHHHHHHHTTEE-S--SGGGEEEETTTEE....EE--GTT.E..EECSS-C-S--S-GGGS-HHHHC CS
Pkinase 88 llkkegklseeeikkialqilegleylHsngiiHrDLKpeNiLldkkgev....kiaDFGlakkleksseklttlvgtreYmAPEvll 171
+ +++++se+e+ +++++i++ ++ylHs+g++HrDLKp+NiL +++ +i+DFG+ak+l +++ l t + t +++APEvl+
sp|Q9WUT3|KS6A2_MOUSE 497 RILRQRCFSEREASDVLYTIARTMDYLHSQGVVHRDLKPSNILYMDESGNpesiRICDFGFAKQLRAENGLLMTPCYTANFVAPEVLK 584
*******************************************98544444888********************************** PP
CS-CTHHHHHHHHHHHHHHHHHHSS-TTSSSHHCCTHHHHSSHHH......TTS.....HHHHHHHHHHT-SSGGGSTTHHHHHTSGG CS
Pkinase 172 kakeytkkvDvWslGvilyelltgklpfsgeseedqleliekilkkkleedepkssskseelkdlikkllekdpakRltaeeilkhpw 259
+++y++++DvWslG++ly++l g +pf + +++ +e++++i ++k+ + +++s s+ +kd+++k+l+ dp++Rlta ++lkhpw
sp|Q9WUT3|KS6A2_MOUSE 585 -RQGYDAACDVWSLGILLYTMLAGFTPFANGPDDTPEEILARIGSGKYALSGGNWDSISDAAKDVVSKMLHVDPQQRLTAVQVLKHPW 671
.*************************************************************************************** PP
G CS
Pkinase 260 l 260
+
sp|Q9WUT3|KS6A2_MOUSE 672 I 672
7 PP
>> sp|Q63531|KS6A1_RAT Ribosomal protein S6 kinase alpha-1 OS=Rattus norvegicus GN=Rps6ka1 PE=1 SV=1
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 239.2 0.0 1.9e-74 1.4e-69 1 258 [. 62 319 .. 62 321 .. 0.94
2 ! 249.7 0.0 1.1e-77 8.7e-73 2 260 .] 419 675 .. 418 675 .. 0.97
Alignments for each domain:
== domain 1 score: 239.2 bits; conditional E-value: 1.9e-74
EEEEEEEEEETTEEE...EEEEETTTTEEEEEEEEEHHHCCCCCCHHHHHHHHHHHHHSSSSB--EEEEEEETTEEEEEEE--TS-BHHH CS
Pkinase 1 yelleklGsGsfGkV...ykakkkktgkkvAvKilkkeeekskkektavrElkilkklsHpnivkllevfetkdelylvleyveggdlfd 87
+ell+ lG+GsfGkV k+++ ++g+ +A+K+lkk + k + + +++ E il+ ++Hp++vkl+ +f+t+ +lyl+l++++ggdlf+
sp|Q63531|KS6A1_RAT 62 FELLKVLGQGSFGKVflvRKVTRPDNGHLYAMKVLKKATLKVRDRVRTKMERDILADVNHPFVVKLHYAFQTEGKLYLILDFLRGGDLFT 151
67899**********4444555666777****************999999**************************************** PP
HHHHHHST-HHHHHHHHHHHHHHHHHHHHTTEE-S--SGGGEEEETTTEEEE--GTT.E..EECSS-C-S--S-GGGS-HHHHCCS-CTH CS
Pkinase 88 llkkegklseeeikkialqilegleylHsngiiHrDLKpeNiLldkkgevkiaDFGlakkleksseklttlvgtreYmAPEvllkakeyt 177
l+ke ++ee++k+++ +++ gl++lHs gii+rDLKpeNiLld++g++k++DFGl+k++ +++k+++++gt eYmAPEv++ ++++t
sp|Q63531|KS6A1_RAT 152 RLSKEVMFTEEDVKFYLAELALGLDHLHSLGIIYRDLKPENILLDEEGHIKLTDFGLSKEAIDHEKKAYSFCGTVEYMAPEVVN-RQGHT 240
************************************************************************************.***** PP
HHHHHHHHHHHHHHHHHSS-TTSSSHHCCTHHHHSSHHH......TTS.....HHHHHHHHHHT-SSGGGSTT.....HHHHHTSG CS
Pkinase 178 kkvDvWslGvilyelltgklpfsgeseedqleliekilkkkleedepkssskseelkdlikkllekdpakRlt.....aeeilkhp 258
+++D+Ws+Gv+++e+ltg+lpf+g+ d++e++ ilk kl ++ s +e+++l++ l++++pa+Rl aeei++h
sp|Q63531|KS6A1_RAT 241 HSADWWSYGVLMFEMLTGSLPFQGK---DRKETMTLILKAKLGMPQFLS----TEAQSLLRALFKRNPANRLGsgpdgAEEIKRHI 319
*************************...999999999999999998866....89*****************98888888888885 PP
== domain 2 score: 249.7 bits; conditional E-value: 1.1e-77
EEEEEEEEETTEEEEEEEETTTTEEEEEEEEEHHHCCCCCCHHHHHHHHHHHHH.SSSSB--EEEEEEETTEEEEEEE--TS-BHHHHHH CS
Pkinase 2 elleklGsGsfGkVykakkkktgkkvAvKilkkeeekskkektavrElkilkkl.sHpnivkllevfetkdelylvleyveggdlfdllk 90
++e++G Gs+++++++++k t+ ++AvK+++k++++ ++ E++il + +Hpni++l +v+ +++++ylv+e+++gg+l d +
sp|Q63531|KS6A1_RAT 419 IVKETIGVGSYSVCKRCVHKATNMEYAVKVIDKSKRDPSE------EIEILLRYgQHPNIITLKDVYDDSKHVYLVTELMRGGELLDKIL 502
6789**********************************98......9*******99********************************** PP
HHHST-HHHHHHHHHHHHHHHHHHHHTTEE-S--SGGGEEEETTTEE....EE--GTT.E..EECSS-C-S--S-GGGS-HHHHCCS-CT CS
Pkinase 91 kegklseeeikkialqilegleylHsngiiHrDLKpeNiLldkkgev....kiaDFGlakkleksseklttlvgtreYmAPEvllkakey 176
+++ +se+e+ +++++i + +eylHs+g++HrDLKp+NiL +++ +i+DFG+ak+l +++ l t + t +++APEvl+ +++y
sp|Q63531|KS6A1_RAT 503 RQKFFSEREASFVLYTISKTVEYLHSQGVVHRDLKPSNILYVDESGNpeclRICDFGFAKQLRAENGLLMTPCYTANFVAPEVLK-RQGY 591
9***************************************98554445999**********************************.**** PP
HHHHHHHHHHHHHHHHHHSS-TTSSSHHCCTHHHHSSHHH......TTS.....HHHHHHHHHHT-SSGGGSTTHHHHHTSGGG CS
Pkinase 177 tkkvDvWslGvilyelltgklpfsgeseedqleliekilkkkleedepkssskseelkdlikkllekdpakRltaeeilkhpwl 260
++ +D+WslGv+ly++l g +pf + +++ +e++ +i+++k + ++++ se++kdl++k+l+ dp++Rlta+++l+hpw+
sp|Q63531|KS6A1_RAT 592 DEGCDIWSLGVLLYTMLAGYTPFANGPSDTPEEILTRISSGKFTLSGGNWNTVSETAKDLVSKMLHVDPHQRLTAKQVLQHPWI 675
***********************************************************************************7 PP
>> sp|Q15418|KS6A1_HUMAN Ribosomal protein S6 kinase alpha-1 OS=Homo sapiens GN=RPS6KA1 PE=1 SV=2
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 238.6 0.0 2.9e-74 2.2e-69 1 258 [. 62 319 .. 62 321 .. 0.94
2 ! 247.4 0.0 6e-77 4.6e-72 1 260 [] 418 675 .. 418 675 .. 0.97
Alignments for each domain:
== domain 1 score: 238.6 bits; conditional E-value: 2.9e-74
EEEEEEEEEETTEEE...EEEEETTTTEEEEEEEEEHHHCCCCCCHHHHHHHHHHHHHSSSSB--EEEEEEETTEEEEEEE--TS-BH CS
Pkinase 1 yelleklGsGsfGkV...ykakkkktgkkvAvKilkkeeekskkektavrElkilkklsHpnivkllevfetkdelylvleyveggdl 85
+ell+ lG+GsfGkV k+++ ++g+ +A+K+lkk + k + + +++ E il+ ++Hp++vkl+ +f+t+ +lyl+l++++ggdl
sp|Q15418|KS6A1_HUMAN 62 FELLKVLGQGSFGKVflvRKVTRPDSGHLYAMKVLKKATLKVRDRVRTKMERDILADVNHPFVVKLHYAFQTEGKLYLILDFLRGGDL 149
67899**********4444555566667****************999999************************************** PP
HHHHHHHHST-HHHHHHHHHHHHHHHHHHHHTTEE-S--SGGGEEEETTTEEEE--GTT.E..EECSS-C-S--S-GGGS-HHHHCCS CS
Pkinase 86 fdllkkegklseeeikkialqilegleylHsngiiHrDLKpeNiLldkkgevkiaDFGlakkleksseklttlvgtreYmAPEvllka 173
f+ l+ke ++ee++k+++ +++ gl++lHs gii+rDLKpeNiLld++g++k++DFGl+k++ +++k+++++gt eYmAPEv++ +
sp|Q15418|KS6A1_HUMAN 150 FTRLSKEVMFTEEDVKFYLAELALGLDHLHSLGIIYRDLKPENILLDEEGHIKLTDFGLSKEAIDHEKKAYSFCGTVEYMAPEVVN-R 236
**************************************************************************************.* PP
-CTHHHHHHHHHHHHHHHHHHSS-TTSSSHHCCTHHHHSSHHH......TTS.....HHHHHHHHHHT-SSGGGSTT.....HHHHHT CS
Pkinase 174 keytkkvDvWslGvilyelltgklpfsgeseedqleliekilkkkleedepkssskseelkdlikkllekdpakRlt.....aeeilk 256
+++++++D+Ws+Gv+++e+ltg+lpf+g+ d++e++ ilk kl ++ s +e+++l++ l++++pa+Rl aeei++
sp|Q15418|KS6A1_HUMAN 237 QGHSHSADWWSYGVLMFEMLTGSLPFQGK---DRKETMTLILKAKLGMPQFLS----TEAQSLLRALFKRNPANRLGsgpdgAEEIKR 317
*****************************...999999999999999998866....89*****************988888899888 PP
SG CS
Pkinase 257 hp 258
h
sp|Q15418|KS6A1_HUMAN 318 HV 319
85 PP
== domain 2 score: 247.4 bits; conditional E-value: 6e-77
EEEEEEEEEETTEEEEEEEETTTTEEEEEEEEEHHHCCCCCCHHHHHHHHHHHHH.SSSSB--EEEEEEETTEEEEEEE--TS-BHHH CS
Pkinase 1 yelleklGsGsfGkVykakkkktgkkvAvKilkkeeekskkektavrElkilkkl.sHpnivkllevfetkdelylvleyveggdlfd 87
y ++e++G Gs+++++++++k t+ ++AvK+++k++++ ++ E++il + +Hpni++l +v+ + +++ylv+e+++gg+l d
sp|Q15418|KS6A1_HUMAN 418 YVVKETIGVGSYSECKRCVHKATNMEYAVKVIDKSKRDPSE------EIEILLRYgQHPNIITLKDVYDDGKHVYLVTELMRGGELLD 499
67899**********************************98......9*******99******************************* PP
HHHHHHST-HHHHHHHHHHHHHHHHHHHHTTEE-S--SGGGEEEETTTEE....EE--GTT.E..EECSS-C-S--S-GGGS-HHHHC CS
Pkinase 88 llkkegklseeeikkialqilegleylHsngiiHrDLKpeNiLldkkgev....kiaDFGlakkleksseklttlvgtreYmAPEvll 171
+ +++ +se+e+ +++ +i + +eylHs+g++HrDLKp+NiL +++ +i+DFG+ak+l +++ l t + t +++APEvl+
sp|Q15418|KS6A1_HUMAN 500 KILRQKFFSEREASFVLHTIGKTVEYLHSQGVVHRDLKPSNILYVDESGNpeclRICDFGFAKQLRAENGLLMTPCYTANFVAPEVLK 587
***9***************************************98554445999********************************** PP
CS-CTHHHHHHHHHHHHHHHHHHSS-TTSSSHHCCTHHHHSSHHH......TTS.....HHHHHHHHHHT-SSGGGSTTHHHHHTSGG CS
Pkinase 172 kakeytkkvDvWslGvilyelltgklpfsgeseedqleliekilkkkleedepkssskseelkdlikkllekdpakRltaeeilkhpw 259
+++y++ +D+WslG++ly++l g +pf + +++ +e++ +i ++k + ++++ se++kdl++k+l+ dp++Rlta+++l+hpw
sp|Q15418|KS6A1_HUMAN 588 -RQGYDEGCDIWSLGILLYTMLAGYTPFANGPSDTPEEILTRIGSGKFTLSGGNWNTVSETAKDLVSKMLHVDPHQRLTAKQVLQHPW 674
.*************************************************************************************** PP
G CS
Pkinase 260 l 260
+
sp|Q15418|KS6A1_HUMAN 675 V 675
6 PP
>> sp|P18652|KS6AA_CHICK Ribosomal protein S6 kinase 2 alpha OS=Gallus gallus GN=RPS6KA PE=2 SV=1
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 239.0 0.0 2.1e-74 1.6e-69 1 259 [. 80 338 .. 80 339 .. 0.95
2 ! 246.7 0.0 9.8e-77 7.5e-72 1 260 [] 435 692 .. 435 692 .. 0.97
Alignments for each domain:
== domain 1 score: 239.0 bits; conditional E-value: 2.1e-74
EEEEEEEEEETTEEEEEEEET...TTTEEEEEEEEEHHHCCCCCCHHHHHHHHHHHHHSSSSB--EEEEEEETTEEEEEEE--TS-BH CS
Pkinase 1 yelleklGsGsfGkVykakkk...ktgkkvAvKilkkeeekskkektavrElkilkklsHpnivkllevfetkdelylvleyveggdl 85
+ell+ lG+GsfGkV+ ++k ++++ +A+K+lkk + k + + +++ E il+ ++Hp++vkl+ +f+t+ +lyl+l++++ggdl
sp|P18652|KS6AA_CHICK 80 FELLKVLGQGSFGKVFLVRKItppDSNHLYAMKVLKKATLKVRDRVRTKIERDILADVNHPFVVKLHYAFQTEGKLYLILDFLRGGDL 167
67899**********7777661115677**************99999888************************************** PP
HHHHHHHHST-HHHHHHHHHHHHHHHHHHHHTTEE-S--SGGGEEEETTTEEEE--GTT.E..EECSS-C-S--S-GGGS-HHHHCCS CS
Pkinase 86 fdllkkegklseeeikkialqilegleylHsngiiHrDLKpeNiLldkkgevkiaDFGlakkleksseklttlvgtreYmAPEvllka 173
f+ l+ke ++ee++k+++ +++ gl++lHs gii+rDLKpeNiLld++g++k++DFGl+k++ +++k+++++gt eYmAPEv++ +
sp|P18652|KS6AA_CHICK 168 FTRLSKEVMFTEEDVKFYLAELALGLDHLHSLGIIYRDLKPENILLDEEGHIKLTDFGLSKEAIDHEKKAYSFCGTVEYMAPEVVN-R 254
**************************************************************************************.* PP
-CTHHHHHHHHHHHHHHHHHHSS-TTSSSHHCCTHHHHSSHHH......TTS.....HHHHHHHHHHT-SSGGGSTT.....HHHHHT CS
Pkinase 174 keytkkvDvWslGvilyelltgklpfsgeseedqleliekilkkkleedepkssskseelkdlikkllekdpakRlt.....aeeilk 256
+++++++D+Ws+Gv+++e+ltg+lpf+g+ d++e++ ilk kl ++ s +e+++l++ l++++pa+Rl aeei++
sp|P18652|KS6AA_CHICK 255 QGHSHSADWWSYGVLMFEMLTGSLPFQGK---DRKETMTLILKAKLGMPQFLS----AEAQSLLRALFKRNPANRLGsgpdgAEEIKR 335
*****************************...999999999999999998866....99*****************999999****** PP
SGG CS
Pkinase 257 hpw 259
hp+
sp|P18652|KS6AA_CHICK 336 HPF 338
*98 PP
== domain 2 score: 246.7 bits; conditional E-value: 9.8e-77
EEEEEEEEEETTEEEEEEEETTTTEEEEEEEEEHHHCCCCCCHHHHHHHHHHHHH.SSSSB--EEEEEEETTEEEEEEE--TS-BHHH CS
Pkinase 1 yelleklGsGsfGkVykakkkktgkkvAvKilkkeeekskkektavrElkilkkl.sHpnivkllevfetkdelylvleyveggdlfd 87
y ++e +G Gs+++++++++k+t+ ++AvK+++k++++ ++ E++il + +Hpni++l +v+ + +++ylv+e+++gg+l d
sp|P18652|KS6AA_CHICK 435 YVVKEAIGVGSYSVCKRCIHKTTNMEYAVKVIDKSKRDPSE------EIEILLRYgQHPNIITLKDVYDDGKYVYLVTELMRGGELLD 516
568899*********************************98......9*******99******************************* PP
HHHHHHST-HHHHHHHHHHHHHHHHHHHHTTEE-S--SGGGEEEETTTEE....EE--GTT.E..EECSS-C-S--S-GGGS-HHHHC CS
Pkinase 88 llkkegklseeeikkialqilegleylHsngiiHrDLKpeNiLldkkgev....kiaDFGlakkleksseklttlvgtreYmAPEvll 171
+ +++ +se+e+ +++ +i++ +eylHs+g++HrDLKp+NiL +++ +i+DFG+ak+l +++ l t + t +++APEvl+
sp|P18652|KS6AA_CHICK 517 KILRQKFFSEREASSVLHTICKTVEYLHSQGVVHRDLKPSNILYVDESGNpesiRICDFGFAKQLRAENGLLMTPCYTANFVAPEVLK 604
***9***************************************98554444888********************************** PP
CS-CTHHHHHHHHHHHHHHHHHHSS-TTSSSHHCCTHHHHSSHHH......TTS.....HHHHHHHHHHT-SSGGGSTTHHHHHTSGG CS
Pkinase 172 kakeytkkvDvWslGvilyelltgklpfsgeseedqleliekilkkkleedepkssskseelkdlikkllekdpakRltaeeilkhpw 259
+++y++ +D+WslGv+ly++l g++pf + +++ +e++ +i ++k + + ++++ s+ +kdl++k+l+ dp++Rlta+++l+hpw
sp|P18652|KS6AA_CHICK 605 -RQGYDEGCDIWSLGVLLYTMLAGCTPFANGPSDTPEEILTRIGGGKFSVNGGNWDTISDVAKDLVSKMLHVDPHQRLTAKQVLQHPW 691
.*************************************************************************************** PP
G CS
Pkinase 260 l 260
+
sp|P18652|KS6AA_CHICK 692 I 692
7 PP
>> sp|Q15349|KS6A2_HUMAN Ribosomal protein S6 kinase alpha-2 OS=Homo sapiens GN=RPS6KA2 PE=1 SV=2
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 241.7 0.0 3.2e-75 2.4e-70 1 260 [] 59 318 .. 59 318 .. 0.95
2 ! 240.7 0.0 6.4e-75 4.9e-70 1 260 [] 415 672 .. 415 672 .. 0.97
Alignments for each domain:
== domain 1 score: 241.7 bits; conditional E-value: 3.2e-75
EEEEEEEEEETTEEEEEEEETTTTE...EEEEEEEEHHHCCCCCCHHHHHHHHHHHHHSSSSB--EEEEEEETTEEEEEEE--TS-BH CS
Pkinase 1 yelleklGsGsfGkVykakkkktgk...kvAvKilkkeeekskkektavrElkilkklsHpnivkllevfetkdelylvleyveggdl 85
+ell+ lG+Gs+GkV+ ++k k ++ +A+K+lkk + k + + +++ E il++++Hp+ivkl+ +f+t+ +lyl+l++++ggdl
sp|Q15349|KS6A2_HUMAN 59 FELLKVLGQGSYGKVFLVRKVKGSDagqLYAMKVLKKATLKVRDRVRSKMERDILAEVNHPFIVKLHYAFQTEGKLYLILDFLRGGDL 146
67899**********7777665544444*****************99999************************************** PP
HHHHHHHHST-HHHHHHHHHHHHHHHHHHHHTTEE-S--SGGGEEEETTTEEEE--GTT.E..EECSS-C-S--S-GGGS-HHHHCCS CS
Pkinase 86 fdllkkegklseeeikkialqilegleylHsngiiHrDLKpeNiLldkkgevkiaDFGlakkleksseklttlvgtreYmAPEvllka 173
f+ l+ke ++ee++k+++ +++ +l++lHs gii+rDLKpeNiLld++g++ki+DFGl+k++ +++++++++gt eYmAPEv++ +
sp|Q15349|KS6A2_HUMAN 147 FTRLSKEVMFTEEDVKFYLAELALALDHLHSLGIIYRDLKPENILLDEEGHIKITDFGLSKEAIDHDKRAYSFCGTIEYMAPEVVN-R 233
**************************************************************************************.* PP
-CTHHHHHHHHHHHHHHHHHHSS-TTSSSHHCCTHHHHSSHHH......TTS.....HHHHHHHHHHT-SSGGGSTT.....HHHHHT CS
Pkinase 174 keytkkvDvWslGvilyelltgklpfsgeseedqleliekilkkkleedepkssskseelkdlikkllekdpakRlt.....aeeilk 256
+++t+++D+Ws+Gv+++e+ltg+lpf+g+ d++e+++ ilk kl ++ s e+++l++ l++++p +Rl +eei++
sp|Q15349|KS6A2_HUMAN 234 RGHTQSADWWSFGVLMFEMLTGSLPFQGK---DRKETMALILKAKLGMPQFLS----GEAQSLLRALFKRNPCNRLGagidgVEEIKR 314
*****************************...999999999999999888866....89******************9999999**** PP
SGGG CS
Pkinase 257 hpwl 260
hp++
sp|Q15349|KS6A2_HUMAN 315 HPFF 318
*997 PP
== domain 2 score: 240.7 bits; conditional E-value: 6.4e-75
EEEEEEEEEETTEEEEEEEETTTTEEEEEEEEEHHHCCCCCCHHHHHHHHHHHHH.SSSSB--EEEEEEETTEEEEEEE--TS-BHHH CS
Pkinase 1 yelleklGsGsfGkVykakkkktgkkvAvKilkkeeekskkektavrElkilkkl.sHpnivkllevfetkdelylvleyveggdlfd 87
ye++e +G Gs+++++++++k t +++AvKi++k++++ ++ E++il + +Hpni++l +v+ + + +ylv+e+++gg+l d
sp|Q15349|KS6A2_HUMAN 415 YEIKEDIGVGSYSVCKRCVHKATDTEYAVKIIDKSKRDPSE------EIEILLRYgQHPNIITLKDVYDDGKFVYLVMELMRGGELLD 496
7899***********************************98......9*******99******************************* PP
HHHHHHST-HHHHHHHHHHHHHHHHHHHHTTEE-S--SGGGEEEETTTEE....EE--GTT.E..EECSS-C-S--S-GGGS-HHHHC CS
Pkinase 88 llkkegklseeeikkialqilegleylHsngiiHrDLKpeNiLldkkgev....kiaDFGlakkleksseklttlvgtreYmAPEvll 171
+ +++ +se+e+ +++ +i + ++ylHs+g++HrDLKp+NiL +++ +++DFG+ak+l ++ l t + t +++APEvl+
sp|Q15349|KS6A2_HUMAN 497 RILRQRYFSEREASDVLCTITKTMDYLHSQGVVHRDLKPSNILYRDESGSpesiRVCDFGFAKQLRAGNGLLMTPCYTANFVAPEVLK 584
********************************************9554445888********************************** PP
CS-CTHHHHHHHHHHHHHHHHHHSS-TTSSSHHCCTHHHHSSHHH......TTS.....HHHHHHHHHHT-SSGGGSTTHHHHHTSGG CS
Pkinase 172 kakeytkkvDvWslGvilyelltgklpfsgeseedqleliekilkkkleedepkssskseelkdlikkllekdpakRltaeeilkhpw 259
+++y++++D+WslG++ly++l g +pf + +++ +e++++i ++k+ + +++s s+ +kd+++k+l+ dp++Rlta+++lkhpw
sp|Q15349|KS6A2_HUMAN 585 -RQGYDAACDIWSLGILLYTMLAGFTPFANGPDDTPEEILARIGSGKYALSGGNWDSISDAAKDVVSKMLHVDPHQRLTAMQVLKHPW 671
.*************************************************************************************** PP
G CS
Pkinase 260 l 260
+
sp|Q15349|KS6A2_HUMAN 672 V 672
6 PP
>> sp|P51812|KS6A3_HUMAN Ribosomal protein S6 kinase alpha-3 OS=Homo sapiens GN=RPS6KA3 PE=1 SV=1
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 240.3 0.0 8.7e-75 6.6e-70 1 260 [] 68 327 .. 68 327 .. 0.95
2 ! 241.0 0.0 5.1e-75 3.9e-70 1 260 [] 422 679 .. 422 679 .. 0.97
Alignments for each domain:
== domain 1 score: 240.3 bits; conditional E-value: 8.7e-75
EEEEEEEEEETTEEEEEEEETTTTE...EEEEEEEEHHHCCCCCCHHHHHHHHHHHHHSSSSB--EEEEEEETTEEEEEEE--TS-BH CS
Pkinase 1 yelleklGsGsfGkVykakkkktgk...kvAvKilkkeeekskkektavrElkilkklsHpnivkllevfetkdelylvleyveggdl 85
+ell+ lG+GsfGkV+ +kk + ++ +A+K+lkk + k + + +++ E il +++Hp+ivkl+ +f+t+ +lyl+l++++ggdl
sp|P51812|KS6A3_HUMAN 68 FELLKVLGQGSFGKVFLVKKISGSDarqLYAMKVLKKATLKVRDRVRTKMERDILVEVNHPFIVKLHYAFQTEGKLYLILDFLRGGDL 155
67899**********8888876655455****************999999************************************** PP
HHHHHHHHST-HHHHHHHHHHHHHHHHHHHHTTEE-S--SGGGEEEETTTEEEE--GTT.E..EECSS-C-S--S-GGGS-HHHHCCS CS
Pkinase 86 fdllkkegklseeeikkialqilegleylHsngiiHrDLKpeNiLldkkgevkiaDFGlakkleksseklttlvgtreYmAPEvllka 173
f+ l+ke ++ee++k+++ +++ +l++lHs gii+rDLKpeNiLld++g++k++DFGl+k+ +++k+++++gt eYmAPEv++ +
sp|P51812|KS6A3_HUMAN 156 FTRLSKEVMFTEEDVKFYLAELALALDHLHSLGIIYRDLKPENILLDEEGHIKLTDFGLSKESIDHEKKAYSFCGTVEYMAPEVVN-R 242
**************************************************************9***********************.* PP
-CTHHHHHHHHHHHHHHHHHHSS-TTSSSHHCCTHHHHSSHHH......TTS.....HHHHHHHHHHT-SSGGGSTT.....HHHHHT CS
Pkinase 174 keytkkvDvWslGvilyelltgklpfsgeseedqleliekilkkkleedepkssskseelkdlikkllekdpakRlt.....aeeilk 256
+++t+++D+Ws+Gv+++e+ltg+lpf+g+ d++e++ +ilk kl ++ s +e+++l++ l++++pa+Rl +eei++
sp|P51812|KS6A3_HUMAN 243 RGHTQSADWWSFGVLMFEMLTGTLPFQGK---DRKETMTMILKAKLGMPQFLS----PEAQSLLRMLFKRNPANRLGagpdgVEEIKR 323
*****************************...999999999999988888866....9******************988888899999 PP
SGGG CS
Pkinase 257 hpwl 260
h+++
sp|P51812|KS6A3_HUMAN 324 HSFF 327
9886 PP
== domain 2 score: 241.0 bits; conditional E-value: 5.1e-75
EEEEEEEEEETTEEEEEEEETTTTEEEEEEEEEHHHCCCCCCHHHHHHHHHHHHH.SSSSB--EEEEEEETTEEEEEEE--TS-BHHH CS
Pkinase 1 yelleklGsGsfGkVykakkkktgkkvAvKilkkeeekskkektavrElkilkkl.sHpnivkllevfetkdelylvleyveggdlfd 87
ye++e +G Gs+++++++++k t+ ++AvKi++k++++ ++ E++il + +Hpni++l +v+ + +++y+v+e+++gg+l d
sp|P51812|KS6A3_HUMAN 422 YEVKEDIGVGSYSVCKRCIHKATNMEFAVKIIDKSKRDPTE------EIEILLRYgQHPNIITLKDVYDDGKYVYVVTELMKGGELLD 503
7899***********************************88......9*******99******************************* PP
HHHHHHST-HHHHHHHHHHHHHHHHHHHHTTEE-S--SGGGEEEETTTEE....EE--GTT.E..EECSS-C-S--S-GGGS-HHHHC CS
Pkinase 88 llkkegklseeeikkialqilegleylHsngiiHrDLKpeNiLldkkgev....kiaDFGlakkleksseklttlvgtreYmAPEvll 171
+ +++ +se+e+ ++ +i + +eylH +g++HrDLKp+NiL +++ +i+DFG+ak+l +++ l t + t +++APEvl+
sp|P51812|KS6A3_HUMAN 504 KILRQKFFSEREASAVLFTITKTVEYLHAQGVVHRDLKPSNILYVDESGNpesiRICDFGFAKQLRAENGLLMTPCYTANFVAPEVLK 591
***9***************************************98554444888********************************** PP
CS-CTHHHHHHHHHHHHHHHHHHSS-TTSSSHHCCTHHHHSSHHH......TTS.....HHHHHHHHHHT-SSGGGSTTHHHHHTSGG CS
Pkinase 172 kakeytkkvDvWslGvilyelltgklpfsgeseedqleliekilkkkleedepkssskseelkdlikkllekdpakRltaeeilkhpw 259
+++y++++D+WslGv+ly++ltg +pf + +++ +e++++i ++k + + ++s s+++kdl++k+l+ dp++Rlta+ +l+hpw
sp|P51812|KS6A3_HUMAN 592 -RQGYDAACDIWSLGVLLYTMLTGYTPFANGPDDTPEEILARIGSGKFSLSGGYWNSVSDTAKDLVSKMLHVDPHQRLTAALVLRHPW 678
.*************************************************************************************** PP
G CS
Pkinase 260 l 260
+
sp|P51812|KS6A3_HUMAN 679 I 679
7 PP
>> sp|P18654|KS6A3_MOUSE Ribosomal protein S6 kinase alpha-3 OS=Mus musculus GN=Rps6ka3 PE=1 SV=2
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 240.3 0.0 8.7e-75 6.6e-70 1 260 [] 68 327 .. 68 327 .. 0.95
2 ! 241.0 0.0 5.1e-75 3.9e-70 1 260 [] 422 679 .. 422 679 .. 0.97
Alignments for each domain:
== domain 1 score: 240.3 bits; conditional E-value: 8.7e-75
EEEEEEEEEETTEEEEEEEETTTTE...EEEEEEEEHHHCCCCCCHHHHHHHHHHHHHSSSSB--EEEEEEETTEEEEEEE--TS-BH CS
Pkinase 1 yelleklGsGsfGkVykakkkktgk...kvAvKilkkeeekskkektavrElkilkklsHpnivkllevfetkdelylvleyveggdl 85
+ell+ lG+GsfGkV+ +kk + ++ +A+K+lkk + k + + +++ E il +++Hp+ivkl+ +f+t+ +lyl+l++++ggdl
sp|P18654|KS6A3_MOUSE 68 FELLKVLGQGSFGKVFLVKKISGSDarqLYAMKVLKKATLKVRDRVRTKMERDILVEVNHPFIVKLHYAFQTEGKLYLILDFLRGGDL 155
67899**********8888876655455****************999999************************************** PP
HHHHHHHHST-HHHHHHHHHHHHHHHHHHHHTTEE-S--SGGGEEEETTTEEEE--GTT.E..EECSS-C-S--S-GGGS-HHHHCCS CS
Pkinase 86 fdllkkegklseeeikkialqilegleylHsngiiHrDLKpeNiLldkkgevkiaDFGlakkleksseklttlvgtreYmAPEvllka 173
f+ l+ke ++ee++k+++ +++ +l++lHs gii+rDLKpeNiLld++g++k++DFGl+k+ +++k+++++gt eYmAPEv++ +
sp|P18654|KS6A3_MOUSE 156 FTRLSKEVMFTEEDVKFYLAELALALDHLHSLGIIYRDLKPENILLDEEGHIKLTDFGLSKESIDHEKKAYSFCGTVEYMAPEVVN-R 242
**************************************************************9***********************.* PP
-CTHHHHHHHHHHHHHHHHHHSS-TTSSSHHCCTHHHHSSHHH......TTS.....HHHHHHHHHHT-SSGGGSTT.....HHHHHT CS
Pkinase 174 keytkkvDvWslGvilyelltgklpfsgeseedqleliekilkkkleedepkssskseelkdlikkllekdpakRlt.....aeeilk 256
+++t+++D+Ws+Gv+++e+ltg+lpf+g+ d++e++ +ilk kl ++ s +e+++l++ l++++pa+Rl +eei++
sp|P18654|KS6A3_MOUSE 243 RGHTQSADWWSFGVLMFEMLTGTLPFQGK---DRKETMTMILKAKLGMPQFLS----PEAQSLLRMLFKRNPANRLGagpdgVEEIKR 323
*****************************...999999999999988888866....9******************988888899999 PP
SGGG CS
Pkinase 257 hpwl 260
h+++
sp|P18654|KS6A3_MOUSE 324 HSFF 327
9886 PP
== domain 2 score: 241.0 bits; conditional E-value: 5.1e-75
EEEEEEEEEETTEEEEEEEETTTTEEEEEEEEEHHHCCCCCCHHHHHHHHHHHHH.SSSSB--EEEEEEETTEEEEEEE--TS-BHHH CS
Pkinase 1 yelleklGsGsfGkVykakkkktgkkvAvKilkkeeekskkektavrElkilkkl.sHpnivkllevfetkdelylvleyveggdlfd 87
ye++e +G Gs+++++++++k t+ ++AvKi++k++++ ++ E++il + +Hpni++l +v+ + +++y+v+e+++gg+l d
sp|P18654|KS6A3_MOUSE 422 YEVKEDIGVGSYSVCKRCIHKATNMEFAVKIIDKSKRDPTE------EIEILLRYgQHPNIITLKDVYDDGKYVYVVTELMKGGELLD 503
7899***********************************88......9*******99******************************* PP
HHHHHHST-HHHHHHHHHHHHHHHHHHHHTTEE-S--SGGGEEEETTTEE....EE--GTT.E..EECSS-C-S--S-GGGS-HHHHC CS
Pkinase 88 llkkegklseeeikkialqilegleylHsngiiHrDLKpeNiLldkkgev....kiaDFGlakkleksseklttlvgtreYmAPEvll 171
+ +++ +se+e+ ++ +i + +eylH +g++HrDLKp+NiL +++ +i+DFG+ak+l +++ l t + t +++APEvl+
sp|P18654|KS6A3_MOUSE 504 KILRQKFFSEREASAVLFTITKTVEYLHAQGVVHRDLKPSNILYVDESGNpesiRICDFGFAKQLRAENGLLMTPCYTANFVAPEVLK 591
***9***************************************98554444888********************************** PP
CS-CTHHHHHHHHHHHHHHHHHHSS-TTSSSHHCCTHHHHSSHHH......TTS.....HHHHHHHHHHT-SSGGGSTTHHHHHTSGG CS
Pkinase 172 kakeytkkvDvWslGvilyelltgklpfsgeseedqleliekilkkkleedepkssskseelkdlikkllekdpakRltaeeilkhpw 259
+++y++++D+WslGv+ly++ltg +pf + +++ +e++++i ++k + + ++s s+++kdl++k+l+ dp++Rlta+ +l+hpw
sp|P18654|KS6A3_MOUSE 592 -RQGYDAACDIWSLGVLLYTMLTGYTPFANGPDDTPEEILARIGSGKFSLSGGYWNSVSDTAKDLVSKMLHVDPHQRLTAALVLRHPW 678
.*************************************************************************************** PP
G CS
Pkinase 260 l 260
+
sp|P18654|KS6A3_MOUSE 679 I 679
7 PP
Internal pipeline statistics summary:
-------------------------------------
Query model(s): 1 (260 nodes)
Target sequences: 535698 (190107059 residues)
Passed MSV filter: 21115 (0.0394159); expected 10714.0 (0.02)
Passed bias filter: 15859 (0.0296044); expected 10714.0 (0.02)
Passed Vit filter: 4968 (0.00927388); expected 535.7 (0.001)
Passed Fwd filter: 3953 (0.00737916); expected 5.4 (1e-05)
Initial search space (Z): 535698 [actual number of targets]
Domain search space (domZ): 7 [number of targets reported over threshold]
# CPU time: 39.03u 0.22s 00:00:39.25 Elapsed: 00:00:10.55
# Mc/sec: 4685.10
//
|
