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