# 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]