# 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/Methyltransf_5.hmm.txt
# target sequence database: proteomes/Trypanosoma_brucei.fasta
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: Methyltransf_5 [M=309]
Accession: PF01795.18
Description: MraW 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
------- ------ ----- ------- ------ ----- ---- -- -------- -----------
4.1e-32 110.9 0.0 5.9e-32 110.3 0.0 1.2 1 tr|Q57W60|Q57W60_TRYB2 S-adenosyl-methyltransferase mraW-lik
0.0047 15.4 0.2 0.0087 14.5 0.2 1.4 1 tr|Q4FKE9|Q4FKE9_TRYB2 Variant surface glycoprotein (VSG), p
------ inclusion threshold ------
0.051 12.0 0.0 0.079 11.4 0.0 1.2 1 tr|Q388Y3|Q388Y3_TRYB2 Uncharacterized protein OS=Trypanosom
Domain annotation for each sequence (and alignments):
>> tr|Q57W60|Q57W60_TRYB2 S-adenosyl-methyltransferase mraW-like protein, putative OS=Trypanosoma brucei brucei (strain
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 110.3 0.0 2.1e-35 5.9e-32 25 255 .. 131 369 .. 116 404 .. 0.87
Alignments for each domain:
== domain 1 score: 110.3 bits; conditional E-value: 2.1e-35
EES--TTSHHHHHHHTT--TT.EEEEEES-HHHHHHHTT.....-.--EEEEES-GGGHHHHHHHTTTT-EEEEEEE-S--HHHHH- CS
Methyltransf_5 25 iDaTlGaggHsealLeqlsel.rligidrDksalarakerLkefe.krltlvkgnfkelkealaeanvtkvdgillDLGvSslqldk 109
+D+TlG+g H+ a+Le+ r++++d D++a+++a++ ++ef +rl+++ + +e k++ e+ ++d+i++D GvS ql++
tr|Q57W60|Q57W60_TRYB2 131 LDCTLGSGYHAGAVLENGGPYtRVVALDCDHDAMHAARDLVEEFGgDRLRFYCCKMSEAKAMFGER---SFDAIMIDGGVSDTQLED 214
8***************998877**********************9455****99999999999988...9***************** PP
GGG---SSS-..-BS---S-TTSS--HHHHHHH--HHHHHHHHHHHH..--TTHHHHHHHHHHHHHHS---BHHHHHHHHHHHS--- CS
Methyltransf_5 110 aerGFsfaed..apLDmRmdqsaeltAaevlnrlseeeLakilkeyg..eekfakriakaivekrkkkpiettkelaevikkvvpak 192
+erGF + + LDmR+ ++ + A e ln++s++ L++ l yg e a ++++ai++++ p+ +e+ i+++ +
tr|Q57W60|Q57W60_TRYB2 215 PERGFLLDDEggHRLDMRFGPQMGVGALEYLNTVSQHTLVSSLLAYGllEYGQAMKMSRAITRRK---PFVDSREVLTCIEQAGDEL 298
******98654489*********************************444445678899997655...99999***99999986554 PP
----.........-TTH.HHHHHHHHHHHTHHHHHHHHHHHHHHHEEEEEEEEEEESSHHHHHHHHHHHHH CS
Methyltransf_5 193 ekrk.......kihpat.rvFQAlRIavNdEleeleelleaalellakggrlavisFhSledrivkkvfke 255
+ + p + +lR +N+E+ el++ +e+al +l+ +grl+v s e+r+v+ ++ +
tr|Q57W60|Q57W60_TRYB2 299 PEGGwrsqgsrRKSPMSwKFLTSLRCIINNEMYELRQGIENALLMLRDDGRLVVFSRLPWEERLVRGTVDD 369
44334445555556665166689******************************************988766 PP
>> tr|Q4FKE9|Q4FKE9_TRYB2 Variant surface glycoprotein (VSG), putative OS=Trypanosoma brucei brucei (strain 927/4 GUTat
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 14.5 0.2 3.1e-06 0.0087 156 209 .. 107 159 .. 105 181 .. 0.89
Alignments for each domain:
== domain 1 score: 14.5 bits; conditional E-value: 3.1e-06
-TTHHHHHHHHHHHHHHS---BHHHHHHHHHHHS-------..-T.THHHHHHHH CS
Methyltransf_5 156 ekfakriakaivekrkkkpiettkelaevikkvvpakekrkkihp.atrvFQAlR 209
e+++k+ a+a++++ ++++i +t+ la+v+++ ++++++ +++p a ++FQAl
tr|Q4FKE9|Q4FKE9_TRYB2 107 EEYWKQAAAAVKDSDTRRQILQTAGLADVTEEQLANYRR--QLNPiAEHAFQALQ 159
78******************************9999988..45774789****95 PP
>> tr|Q388Y3|Q388Y3_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 ? 11.4 0.0 2.8e-05 0.079 11 78 .. 27 100 .. 23 115 .. 0.76
Alignments for each domain:
== domain 1 score: 11.4 bits; conditional E-value: 2.8e-05
HHHHHT---TT-EEEES--TTSHHHHHHHT.T--TT...EEEEEES-HHHHHHHTT.....---..EEEEES-G CS
Methyltransf_5 11 evvellepkedgvyiDaTlGaggHsealLe.qlsel...rligidrDksalarakerLkefekr..ltlvkgnf 78
+v +ll+++ +iDa +gg +++L+ ql+e+ r+++id Dk al++ake +++ + l++
tr|Q388Y3|Q388Y3_TRYB2 27 KVQTLLRENGCSSFIDAGCSSGGLLRHILTsQLQEHsfsRALAIDLDKVALHEAKEAITALGFSspVALLHPMH 100
5667888899999**************99624444334599***************998754332366666555 PP
Internal pipeline statistics summary:
-------------------------------------
Query model(s): 1 (309 nodes)
Target sequences: 8561 (4332363 residues searched)
Passed MSV filter: 259 (0.0302535); expected 171.2 (0.02)
Passed bias filter: 143 (0.0167037); expected 171.2 (0.02)
Passed Vit filter: 11 (0.0012849); expected 8.6 (0.001)
Passed Fwd filter: 3 (0.000350426); expected 0.1 (1e-05)
Initial search space (Z): 8561 [actual number of targets]
Domain search space (domZ): 3 [number of targets reported over threshold]
# CPU time: 0.10u 0.00s 00:00:00.10 Elapsed: 00:00:00.05
# Mc/sec: 26774.00
//
[ok]