# 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/Trichoplax_adhaerens.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1.9e-96 322.6 0.2 1.3e-95 319.8 0.2 1.9 1 tr|B3RLV1|B3RLV1_TRIAD Uncharacterized protein OS=Trichoplax ------ inclusion threshold ------ 0.034 13.0 0.0 0.06 12.2 0.0 1.4 1 tr|B3RUJ8|B3RUJ8_TRIAD Uncharacterized protein OS=Trichoplax Domain annotation for each sequence (and alignments): >> tr|B3RLV1|B3RLV1_TRIAD Uncharacterized protein OS=Trichoplax adhaerens OX=10228 GN=TRIADDRAFT_19892 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 319.8 0.2 2.3e-99 1.3e-95 2 306 .. 29 367 .. 28 369 .. 0.90 Alignments for each domain: == domain 1 score: 319.8 bits; conditional E-value: 2.3e-99 XXX--TTTHHHHHHT---TT-EEEES--TTSHHHHHHHTT--TTEEEEEES-HHHHHHHTT........---EEEEES-GGGHHHHH CS Methyltransf_5 2 leHvpVLlqevvellepkedgvyiDaTlGaggHsealLeqlselrligidrDksalarakerLk...efekrltlvkgnfkelkeal 85 l+H+pVLl+ vv+ll+p +++ yiDaT+G+ggH+ a+L++ ++ ++++id D+ a+arak+ + ++++rl++v+ +f +++++l tr|B3RLV1|B3RLV1_TRIAD 29 LKHQPVLLDTVVKLLAPANGQSYIDATFGGGGHTSAILKAATDCQIYAIDTDPLAIARAKQLASttsNWIQRLKVVQCRFGQMENKL 115 58**********************************************************965411155689********9998877 PP HHTT...TT-EEEEEEE-S--HHHHH-GGG---SSS--BS---S-TTS........S--HHHHHHH--HHHHHHHHHHHH--TTHHH CS Methyltransf_5 86 aean...vtkvdgillDLGvSslqldkaerGFsfaedapLDmRmdqsa........eltAaevlnrlseeeLakilkeygeekfakr 161 ++ +++dgi++DLG+Ss+q+d++ rGFsf+ ++pLDmRmd ++ + tAa v+n++++e+La++l +ygee+ a+ tr|B3RLV1|B3RLV1_TRIAD 116 LAQEklpPNSMDGIIFDLGMSSMQIDDPLRGFSFQSSGPLDMRMDATSdeigkhqhRITAAHVINNIDKEDLADVLWNYGEERKARA 202 655546678***********************************97542222222157***************************** PP HHHHHHHHHHHS---BHHHHHHHHHHHS-------...-TTHHHHHHHHHHHHTHHHHHHHHHHHHHHHEEEEEEEEEEESSHHHHH CS Methyltransf_5 162 iakaivekrkkkpiettkelaevikkvvpakekrk.kihpatrvFQAlRIavNdEleeleelleaalellakggrlavisFhSledr 247 ia+ai ++rkk+pi+tt +l+++i + +p k+++ +ihpatr+FQA+RI +NdEl +l + l+a ll++ggrl visFhSledr tr|B3RLV1|B3RLV1_TRIAD 203 IATAIIKARKKQPITTTDHLVDIIDRQFPYKYRKDgQIHPATRTFQAIRILINDELWQLVQGLNACKALLKNGGRLIVISFHSLEDR 289 *****************************88765439************************************************** PP HHHHHHHHTT------------------...................-EEEEE.E.E---HHHHHH-GGGTT-EEEEE CS Methyltransf_5 248 ivkkvfkeaskkeevPaglpvaseelka...................klklitkkpilpseeeikeNprsrsakLrvi 306 ivkk+f+ + + ++ p s+++ + kli kkpi ps eei++N r+rsakLr++ tr|B3RLV1|B3RLV1_TRIAD 290 IVKKFFNGEDRPDFTSRSSPTSSKQYARlseegkskiesksiikrnsPWKLICKKPIYPSPEEIAHNSRARSAKLRCA 367 *******9999777788888888877766777777888888888887899**************************96 PP >> tr|B3RUJ8|B3RUJ8_TRIAD Uncharacterized protein OS=Trichoplax adhaerens OX=10228 GN=TRIADDRAFT_24037 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.2 0.0 1e-05 0.06 18 62 .. 123 175 .. 114 212 .. 0.76 Alignments for each domain: == domain 1 score: 12.2 bits; conditional E-value: 1e-05 -........-TT-EEEES--TTSHHHHHHHTT--TTEEEEEES-HHHHHHHTT CS Methyltransf_5 18 p........kedgvyiDaTlGaggHsealLeqlselrligidrDksalarake 62 +e +++D G+g lL + ++ + i+id+D++a+ ++ tr|B3RUJ8|B3RUJ8_TRIAD 123 LhefhhkreNESITFLDICCGSGAIGLSLLCENPQATCIAIDKDPNAISLTEL 175 44445555478889*********************************987654 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (309 nodes) Target sequences: 11518 (5222075 residues searched) Passed MSV filter: 467 (0.0405452); expected 230.4 (0.02) Passed bias filter: 220 (0.0191005); expected 230.4 (0.02) Passed Vit filter: 20 (0.00173641); expected 11.5 (0.001) Passed Fwd filter: 2 (0.000173641); expected 0.1 (1e-05) Initial search space (Z): 11518 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.13u 0.00s 00:00:00.13 Elapsed: 00:00:00.05 # Mc/sec: 32272.42 // [ok]