# hmmsearch :: search profile(s) against a sequence database
# HMMER 3.1b2 (February 2015); http://hmmer.org/
# Copyright (C) 2015 Howard Hughes Medical Institute.
# Freely distributed under the GNU General Public License (GPLv3).
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# query HMM file: pfam_list/SpoU_methylase.hmm.txt
# target sequence database: proteomes/Trypanosoma_brucei.fasta
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: SpoU_methylase [M=142]
Accession: PF00588.18
Description: SpoU rRNA Methylase family
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
------- ------ ----- ------- ------ ----- ---- -- -------- -----------
1.2e-25 89.3 0.0 3.8e-24 84.4 0.0 2.4 2 tr|Q381K3|Q381K3_TRYB2 Uncharacterized protein OS=Trypanosom
3.1e-17 62.0 0.0 2.9e-15 55.7 0.0 3.0 3 tr|Q57ZH7|Q57ZH7_TRYB2 Uncharacterized protein OS=Trypanosom
3.4e-16 58.7 0.0 1.6e-15 56.5 0.0 2.0 1 tr|Q57TV7|Q57TV7_TRYB2 Uncharacterized protein OS=Trypanosom
8.6e-09 34.7 0.0 4.9e-08 32.2 0.0 2.2 2 tr|Q38DC9|Q38DC9_TRYB2 Uncharacterized protein OS=Trypanosom
2.1e-08 33.4 1.8 0.002 17.2 0.0 4.0 4 tr|Q381U1|Q381U1_TRYB2 Uncharacterized protein OS=Trypanosom
Domain annotation for each sequence (and alignments):
>> tr|Q381K3|Q381K3_TRYB2 Uncharacterized protein OS=Trypanosoma brucei brucei (strain 927/4 GUTat10.1) OX=185431 GN=Tb
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ? 2.4 0.0 0.045 78 2 32 .. 1264 1293 .. 1263 1313 .. 0.80
2 ! 84.4 0.0 2.2e-27 3.8e-24 2 142 .] 1559 1701 .. 1558 1701 .. 0.94
Alignments for each domain:
== domain 1 score: 2.4 bits; conditional E-value: 0.045
EEEEEES-SSHHHHHHHHHHHHHTT-SEEEE CS
SpoU_methylase 2 livvlveiedpgNlGaiaRtakafgaselvl 32
l+vv+ + + +N+G++ R+++ f + +l+
tr|Q381K3|Q381K3_TRYB2 1264 LHVVVMKCLTLTNMGSVRRLMELFAL-QLLQ 1293
578888888999************99.5554 PP
== domain 2 score: 84.4 bits; conditional E-value: 2.2e-27
EEEEEES-SSHHHHHHHHHHHHHTT-SEEEEESGSSCTCGSHHHHHHCTTHHHHHH..EEEESSHHHHHHHHCCCTEEEEEEST. CS
SpoU_methylase 2 livvlveiedpgNlGaiaRtakafgaselvlvkprsadllnpeakraaagalevle..lkvvdsleealkelkeegfalvatsl. 83
+iv+ +++p N+ ++R+ ++f ++++v+++ ++ l +p++++aa ++ +l+ +++v+ l+ +l+ l++ g++lv++++
tr|Q381K3|Q381K3_TRYB2 1559 VIVIASLLQNPVNVAGLFRCGEVFAVEKVVVSD-AAV-LEHPHFVAAARSSDLWLPwsAVQVKALSGYLSSLRQDGYTLVGIEQt 1641
7899999*******************9999966.555.9*****************777888********************997 PP
.T-SEGCCSHHHTTEEEEEEBTTTBS-HHHHCTSSEEEE----TTSSS--HHHHHHHHHH CS
SpoU_methylase 84 .kgkslseidlkekvalvfGnEdkGlsnelleladrlvkIpmsgeveSLNvavAvailly 142
++ ++s +++++++++v+G E++G++ +ll+ +d +v+Ip++g ++SLNv+v+ +i++y
tr|Q381K3|Q381K3_TRYB2 1642 aGSVPMSSYQFPKRAVIVLGAEGHGIPAQLLTVLDVCVEIPQYGLIRSLNVHVTGSIVMY 1701
666677799*************************************************98 PP
>> tr|Q57ZH7|Q57ZH7_TRYB2 Uncharacterized protein OS=Trypanosoma brucei brucei (strain 927/4 GUTat10.1) OX=185431 GN=Tb
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ? -0.8 0.0 0.43 7.4e+02 9 36 .. 96 122 .. 93 128 .. 0.82
2 ! 55.7 0.0 1.7e-18 2.9e-15 3 127 .. 177 322 .. 175 326 .. 0.74
3 ? 2.2 0.1 0.05 86 126 141 .. 360 375 .. 354 376 .. 0.84
Alignments for each domain:
== domain 1 score: -0.8 bits; conditional E-value: 0.43
-SSHHHHHHHHHHHHHTT-SEEEEESGS CS
SpoU_methylase 9 iedpgNlGaiaRtakafgaselvlvkpr 36
+++ N+ i+R+a+a+ + +++ v++r
tr|Q57ZH7|Q57ZH7_TRYB2 96 PKQYPNMQRIVRLARALSV-AVHPVERR 122
667789999********99.88887765 PP
== domain 2 score: 55.7 bits; conditional E-value: 1.7e-18
EEEEES-SSHHHHHHHHHHHHHTT-SEEEEESGSSCTCGSHHHHHHCTTHHHHHHEEEE.....SSHHHHHHHHCCC.T...EEEEE CS
SpoU_methylase 3 ivvlveiedpgNlGaiaRtakafgaselvlvkprsadllnpeakraaagalevlelkvv.....dsleealkelkee.g...falva 80
+v l+++ dp N+Gai+R+a +g++++vl++ +sa + + r++ag e++ ++++ + le+++++ k + g ++++a
tr|Q57ZH7|Q57ZH7_TRYB2 177 VVFLDHVIDPANVGAIMRSAFFMGLRRIVLSS-DSA-GCTAAVCRSSAGFMEYMSVYRAvvptlTFLENTVAQQKARdGghlLEIYA 261
699***************************87.555.7999************9965544444444555555544443235566788 PP
EST...T-SEGCCS...........HHHTTEEEEEEBTTTBS-HHHHCTSSEEEE----TT CS
SpoU_methylase 81 tsl...kgkslsei...........dlkekvalvfGnEdkGlsnelleladrlvkIpmsge 127
ts+ +g +++ + + +++ +l++GnE++Gl+ e++ l++++ +Ip ++e
tr|Q57ZH7|Q57ZH7_TRYB2 262 TSPmprRGVTTCSArelpsvscgpvNKPRRRLLLLGNEGAGLPSEVMGLCTHTAHIPLAPE 322
8887443333332243344444433444*****************************9886 PP
== domain 3 score: 2.2 bits; conditional E-value: 0.05
TTSSS--HHHHHHHHH CS
SpoU_methylase 126 geveSLNvavAvaill 141
+e SLNv+ A a++l
tr|Q57ZH7|Q57ZH7_TRYB2 360 PEEVSLNVSAASALIL 375
5667*********987 PP
>> tr|Q57TV7|Q57TV7_TRYB2 Uncharacterized protein OS=Trypanosoma brucei brucei (strain 927/4 GUTat10.1) OX=185431 GN=Tb
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 56.5 0.0 9.4e-19 1.6e-15 3 141 .. 174 342 .. 172 343 .. 0.77
Alignments for each domain:
== domain 1 score: 56.5 bits; conditional E-value: 9.4e-19
EEEEES-SSHHHHHHHHHHHHHTT-SEEEEESGSSCTCGSHHHHHHCTTHHHHHHEEEE...SSHHHHHHHHCCC.T.EEEEEEST. CS
SpoU_methylase 3 ivvlveiedpgNlGaiaRtakafgaselvlvkprsadllnpeakraaagalevlelkvv...dsleealkelkee.g.falvatsl. 83
++vlv++++p N+G+++R+a +g++ ++l++ + d++ ++a+ra+ g+ + ++k++ ++ +al+ + e + + l+a+ +
tr|Q57TV7|Q57TV7_TRYB2 174 MLVLVGLRIPSNVGTLIRAAVEMGFESVLLIN--CLDPFGEKALRASEGTVFSPQFKIFepgSDPVSALNSIAVEhNlLPLLALPSq 258
68****************************99..999***************99999994455666677775555233455666554 PP
T-SEGCCS........HHH................TTEEEEEEBTTTBS-HHHHCTSS..EEEE---.-TTSSS--HHHHHHHHH CS
SpoU_methylase 84 kgkslsei........dlk................ekvalvfGnEdkGlsnellelad..rlvkIpm.sgeveSLNvavAvaill 141
k+++ e+ + + +l++G+E kGl++ + e ++v++p + veSLNv+vA +il+
tr|Q57TV7|Q57TV7_TRYB2 259 KAETAFEVaknlhkinA-MrrsqeaegdaaashnhIGPLLILGSEAKGLRDLAGEWSVprKFVSVPLpNSTVESLNVSVAGSILI 342
55544444344443440.155566677777888776678***********98776554449******88889***********97 PP
>> tr|Q38DC9|Q38DC9_TRYB2 Uncharacterized protein OS=Trypanosoma brucei brucei (strain 927/4 GUTat10.1) OX=185431 GN=Tb
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ? -0.8 0.0 0.45 7.7e+02 104 104 .. 123 123 .. 61 156 .. 0.54
2 ! 32.2 0.0 2.9e-11 4.9e-08 3 58 .. 239 293 .. 237 305 .. 0.90
Alignments for each domain:
== domain 1 score: -0.8 bits; conditional E-value: 0.45
T CS
SpoU_methylase 104 E 104
E
tr|Q38DC9|Q38DC9_TRYB2 123 E 123
3 PP
== domain 2 score: 32.2 bits; conditional E-value: 2.9e-11
EEEEES-SSHHHHHHHHHHHHHTT-SEEEEESGSSCTCGSHHHHHHCTTHHHHHHE CS
SpoU_methylase 3 ivvlveiedpgNlGaiaRtakafgaselvlvkprsadllnpeakraaagalevlel 58
++vl++++dpg lG+++Rta+ f ++++ ++ ++adl++ +++raa ga+++ +
tr|Q38DC9|Q38DC9_TRYB2 239 VLVLDNVDDPGLLGTVLRTAAGFHYDAVIATN-HCADLYDHRVIRAARGAHFQKAV 293
79***********************9999988.99***************987654 PP
>> tr|Q381U1|Q381U1_TRYB2 Uncharacterized protein OS=Trypanosoma brucei brucei (strain 927/4 GUTat10.1) OX=185431 GN=Tb
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ? -2.8 0.0 1.8 3.1e+03 2 27 .. 53 78 .. 52 80 .. 0.83
2 ! 17.2 0.0 1.2e-06 0.002 75 142 .] 136 219 .. 113 219 .. 0.73
3 ! 10.6 0.1 0.00014 0.23 4 38 .. 324 358 .. 321 367 .. 0.84
4 ? -1.3 0.0 0.62 1.1e+03 93 141 .. 473 525 .. 438 526 .. 0.61
Alignments for each domain:
== domain 1 score: -2.8 bits; conditional E-value: 1.8
EEEEEES-SSHHHHHHHHHHHHHTT- CS
SpoU_methylase 2 livvlveiedpgNlGaiaRtakafga 27
l v l++ +++ N + +Rt+ +g
tr|Q381U1|Q381U1_TRYB2 53 LSVGLEGCTSAFNAMNSIRTCMFYGT 78
56788889999999999****98885 PP
== domain 2 score: 17.2 bits; conditional E-value: 1.2e-06
TEEEEEEST...T-SEGCCSHHHTTEEEEEEBTTTBS-HHHHC............T.SSEEEE----TTSSS--HHHHHHHHHH CS
SpoU_methylase 75 gfalvatsl...kgkslseidlkekvalvfGnEdkGlsnelle............l.adrlvkIpmsgeveSLNvavAvailly 142
+f lva+++ +++sl + ++++ + l++G E++G+s + +e + ad +v++p++g ++SLNv ++++i+l+
tr|Q381U1|Q381U1_TRYB2 136 NFPLVALENfthRSQSLFTCRISSVARLLIGHENRGVSSKYIEcdnatisdesdgQcADCVVYVPQYGTISSLNVVTSMGIALF 219
46777776665577788888999999999999999999877753333333333333144699*******************985 PP
== domain 3 score: 10.6 bits; conditional E-value: 0.00014
EEEES-SSHHHHHHHHHHHHHTT-SEEEEESGSSC CS
SpoU_methylase 4 vvlveiedpgNlGaiaRtakafgaselvlvkprsa 38
v++ + d +NlG i+R+a+af ++++++ ++r+
tr|Q381U1|Q381U1_TRYB2 324 VLCQNEIDLRNLGGIVRSANAFLVDNIFYFGRRKI 358
5666667889******************9998765 PP
== domain 4 score: -1.3 bits; conditional E-value: 0.62
HHTTEEEEEEBTTTBS-HHHHCTSSEEEE-.-...--TTSSS--HHHHHHHHH CS
SpoU_methylase 93 lkekvalvfGnEdkGlsnelleladrlvkI.p...msgeveSLNvavAvaill 141
++e v+lv+ +E+ + +l+++++++++ p ++ L vavA++i+l
tr|Q381U1|Q381U1_TRYB2 473 IREGVLLVVPQEGILPPLHILQRCEKILTVlPhdfTDTVSAGLPVAVAAGIAL 525
44667777777776667777777777766522111233334567777777765 PP
Internal pipeline statistics summary:
-------------------------------------
Query model(s): 1 (142 nodes)
Target sequences: 8561 (4332363 residues searched)
Passed MSV filter: 227 (0.0265156); expected 171.2 (0.02)
Passed bias filter: 191 (0.0223105); expected 171.2 (0.02)
Passed Vit filter: 17 (0.00198575); expected 8.6 (0.001)
Passed Fwd filter: 5 (0.000584044); expected 0.1 (1e-05)
Initial search space (Z): 8561 [actual number of targets]
Domain search space (domZ): 5 [number of targets reported over threshold]
# CPU time: 0.08u 0.00s 00:00:00.08 Elapsed: 00:00:00.04
# Mc/sec: 15379.89
//
[ok]