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
|
# hmmscan :: search sequence(s) against a profile 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 sequence file: ../sco_all_prots.fa
# target HMM database: ../biodata/Pfam-A.hmm
# output directed to file: ../sco_hmmscan.hsr
# number of worker threads: 2
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: SCO3574 [L=354]
Scores for complete sequence (score includes all domains):
--- full sequence --- --- best 1 domain --- -#dom-
E-value score bias E-value score bias exp N Model Description
------- ------ ----- ------- ------ ----- ---- -- -------- -----------
1.4e-19 70.2 0.6 4.9e-18 65.1 0.4 2.2 1 Abhydrolase_1 alpha/beta hydrolase fold
4.1e-07 29.4 0.0 1e-06 28.1 0.0 1.7 2 Hydrolase_4 Putative lysophospholipase
------ inclusion threshold ------
0.017 13.3 0.0 0.032 12.4 0.0 1.4 1 AXE1 Acetyl xylan esterase (AXE1)
0.054 12.5 0.3 0.13 11.3 0.1 1.7 2 Esterase Putative esterase
0.071 12.3 0.2 0.1 11.8 0.2 1.2 1 Abhydrolase_3 alpha/beta hydrolase fold
Domain annotation for each model (and alignments):
>> Abhydrolase_1 alpha/beta hydrolase fold
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 65.1 0.4 2.1e-21 4.9e-18 1 230 [. 77 346 .. 77 347 .. 0.87
Alignments for each domain:
== domain 1 score: 65.1 bits; conditional E-value: 2.1e-21
SEEEEEESCTTT--S.C...SGGGGGSHHHHHHHHHHHHHHSTSSSEEEEE-CT--HHHHHHHHHSGGGEEEEEEESBSS................ CS
Abhydrolase_1 1 fdviglDlrGvgySs.paggpelphyttddlakdleavldalglkkltlvGhSmGGmialayaaeypdrvkamvllgtva................ 79
f+ ++lD rG+g+Ss p + y++ dl++d av++alg ++++vGh++G+ ia++ a+ pd ++a+ ll+ ++
SCO3574 77 FRAVALDVRGYGRSSrP---DAVEAYRMLDLVADNVAVVEALGESSAVVVGHDWGANIAAHSALLRPDVFRAVGLLSVPYtppggprpseafagms 169
7889***********67...999************************************************************************* PP
..............................HHHSCCSCCTTCHHHHHHHHHHHHHHHHHHHHHHHHCCCCHHHHHHHHHHCHHCCHCHHHHHHHHH CS
Abhydrolase_1 80 ..............................svlsdglssdvlnrgstltlladnvagrlskgikplqgdasvqaqsllrpqvtdflkrfdansyir 145
+ + +gl +++ + +t + d+ + ++ + ++ +++ s +r + l++ d++ y+
SCO3574 170 dpagpfagqefyvsyfqepgraeaeiepdvRGWLAGLYAALSA--GTMPGPQDPDPH-FVAPGGRMRDRFP----SAGRL--PSWLTEEDLDVYAG 256
**************99999999998865554444444444444..566777777777.7788888888888....77777..8999********** PP
TEHHHHHHHCEH..HCHHHHHHHHHCCHHHTTTSCEEEEEETTTSSSHH..HHHHH.HHHHSTTEEEEEETTHHCTHHHHSHHHHHHHHH CS
Abhydrolase_1 146 egetraldgkll..aalgrdlvwdvtetlstilvptlvisgtdDplvpy..kasek.laelipksqlvvidgpgHlafldkpdeinkliv 230
e+e+++l+g l+ ++ rd+ + + +i++p+l ++g D +++ +a+e+ l + s +++dg+gH+++ ++p+e n+l++
SCO3574 257 EFERTGLTGALNryRNMDRDWADLAAHEGAPITQPSLFLGGALDASTTWlsDAIEAyPVTLPGLSASHLLDGCGHWLQQERPEETNRLLT 346
******************************************99955542245555347777779**********************987 PP
>> Hydrolase_4 Putative lysophospholipase
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 28.1 0.0 4.3e-10 1e-06 17 75 .. 50 109 .. 40 113 .. 0.82
2 ? -3.5 0.0 3.3 7.8e+03 38 48 .. 245 255 .. 242 276 .. 0.78
Alignments for each domain:
== domain 1 score: 28.1 bits; conditional E-value: 4.3e-10
Hydrolase_4 17 aavvlvhglaEhsgrYqelaedLaaqGyavyayDhrGhGkseg.ergyvpsfddyvdDla 75
+v+lvhg+ E+ + + +Laa+G++ +a+D rG G+s + ++ ++ d+v D
SCO3574 50 PLVLLVHGFPESWYSWRHQLPALAAAGFRAVALDVRGYGRSSRpDAVEAYRMLDLVADNV 109
579**************************************8624556777777777755 PP
== domain 2 score: -3.5 bits; conditional E-value: 3.3
Hydrolase_4 38 dLaaqGyavya 48
+L++++++vya
SCO3574 245 WLTEEDLDVYA 255
57888888887 PP
>> AXE1 Acetyl xylan esterase (AXE1)
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ? 12.4 0.0 1.4e-05 0.032 80 124 .. 46 91 .. 38 103 .. 0.85
Alignments for each domain:
== domain 1 score: 12.4 bits; conditional E-value: 1.4e-05
---B-EEEB----------HHHHH.HHHH---EEE--------..S CS
AXE1 80 eeklPalvefhGYegssgdwhdkl.kyaaaGyavvamdvRGqgges 124
++ P+++ hG+ +s +w +l + aaaG+ va+dvRG g +s
SCO3574 46 QGSGPLVLLVHGFPESWYSWRHQLpALAAAGFRAVALDVRGYGRSS 91
6789************9999999956899************99554 PP
>> Esterase Putative esterase
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ? -2.3 0.0 0.78 1.9e+03 31 247 .. 56 71 .. 44 88 .. 0.66
2 ? 11.3 0.1 5.5e-05 0.13 103 150 .. 102 152 .. 98 190 .. 0.77
Alignments for each domain:
== domain 1 score: -2.3 bits; conditional E-value: 0.78
------TCHHHHHCHHHHHHHH...EEEEEESSE----S-TTSBTTS-....BB--SCCCCCCHBHHHHTHHHHHHHHHHSEECCE.EEEE--THHHHHHH CS
Esterase 31 DGltgwdqntrakelleslaaegeiaeviivvipaggessfysdwdlgarfyenategpgaeayetfltqeLlplldanfrtnpdgdravvGqSmGGlgAL 131
+G+
SCO3574 56 HGF-------------------------------------------------------------------------------------------------- 58
443.................................................................................................. PP
HHHHC-TTT-SEEEEES----CTCSH---TTCCHHHCT-HHHHHCGTCTTT-EEEEEEEE----CCCCTTHHHHHHHHHCHHCTCCCCEEEEE-TT--SSH CS
Esterase 132 ilalkhpdrFgsvssfSpivnpsnamWgpeekeawqegdpvllakelsasnkrlriyldvGtseddlgeqlsveilevlarnrelkeqlaargveghdhvy 232
SCO3574 - ----------------------------------------------------------------------------------------------------- -
..................................................................................................... PP
H---STTTHHHHHHH CS
Esterase 233 dGghdweyWraqLia 247
+ +w Wr+qL+a
SCO3574 59 P--ESWYSWRHQLPA 71
3..455556666554 PP
== domain 2 score: 11.3 bits; conditional E-value: 5.5e-05
HHHHHHHHSE.ECCE...EEEE--THHHHHHHHHHHC-TTT-SEEEEES-- CS
Esterase 103 lplldanfrt.npdgd..ravvGqSmGGlgALilalkhpdrFgsvssfSpi 150
l+l+++n + ++ g+ +vvG+ G+ +A + al +pd F v+ +S +
SCO3574 102 LDLVADNVAVvEALGEssAVVVGHDWGANIAAHSALLRPDVFRAVGLLSVP 152
55555555441333446689*****************************33 PP
>> Abhydrolase_3 alpha/beta hydrolase fold
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ? 11.8 0.2 4.3e-05 0.1 28 95 .. 75 142 .. 61 196 .. 0.78
Alignments for each domain:
== domain 1 score: 11.8 bits; conditional E-value: 4.3e-05
HTSEEEEEE---TTTT-TTHHHHHHHHHHHHHHHTCC.CHTCCCCEEEEE-------HHHHHHHHHHH- CS
Abhydrolase_3 28 lgavvvsvdYrlAPehrfPaavedalaalrwlaeqaw.elgadpsrivvaGdSaGGnLaaavalrardr 95
+g ++v+ d r + P ave a+++l +a++ + + s+ vv+G+ G+n+aa+ al d
SCO3574 75 AGFRAVALDVRGYGRSSRPDAVE-AYRMLDLVADNVAvVEALGESSAVVVGHDWGANIAAHSALLRPDV 142
55566666666666666677775.788888999998866667789*****************9877664 PP
Internal pipeline statistics summary:
-------------------------------------
Query sequence(s): 1 (354 residues)
Target model(s): 11912 (2158902 nodes)
Passed MSV filter: 170 (0.0142713); expected 238.2 (0.02)
Passed bias filter: 157 (0.01318); expected 238.2 (0.02)
Passed Vit filter: 31 (0.00260242); expected 11.9 (0.001)
Passed Fwd filter: 5 (0.000419745); expected 0.1 (1e-05)
Initial search space (Z): 11912 [actual number of targets]
Domain search space (domZ): 5 [number of targets reported over threshold]
# CPU time: 0.25u 0.18s 00:00:00.43 Elapsed: 00:00:00.24
# Mc/sec: 3184.38
//
|
