# 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/Pan_troglodytes.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.2e-90 302.7 0.0 1.6e-89 300.8 0.0 1.7 1 tr|H2Q3C1|H2Q3C1_PANTR METTL15 isoform 3 OS=Pan troglod 1.3e-52 179.6 0.0 7.1e-52 177.2 0.0 2.1 1 tr|A0A2I3T743|A0A2I3T743_PANTR Uncharacterized protein OS=Pan t Domain annotation for each sequence (and alignments): >> tr|H2Q3C1|H2Q3C1_PANTR METTL15 isoform 3 OS=Pan troglodytes OX=9598 GN=METTL15 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 300.8 0.0 1.4e-93 1.6e-89 3 307 .. 71 405 .. 69 407 .] 0.89 Alignments for each domain: == domain 1 score: 300.8 bits; conditional E-value: 1.4e-93 XX--TTTHHHHHHT---TT-EEEES--TTSHHHHHHHTT--TTEEEEEES-HHHHHHHTT.....---EEEEES-GGGHHHHHHHTT CS Methyltransf_5 3 eHvpVLlqevvellepkedgvyiDaTlGaggHsealLeqlselrligidrDksalarakerLkefekrltlvkgnfkelkealaean 89 H+pV+++evv+ l+p+++++++D+T+G+ggH++a+L++ s+ l+++drD++a a a++ + + k+++ + g+f++ ++ l +a+ tr|H2Q3C1|H2Q3C1_PANTR 71 LHIPVMVDEVVHCLSPQKGQIFLDMTFGSGGHTKAILQKESDIVLYALDRDPTAYALAEHLSELYPKQIRAMLGQFSQAEALLMKAG 157 6***********************************************************77788899******************9 PP TT..-EEEEEEE-S--HHHHH-GGG---SSS--BS---S-TTS..S-.-HHHHHHH--HHHHHHHHHHHH--TTHHHHHHHHHHHHH CS Methyltransf_5 90 vt..kvdgillDLGvSslqldkaerGFsfaedapLDmRmdqsa..el.tAaevlnrlseeeLakilkeygeekfakriakaivekrk 171 v+ ++dg+l+DLG Ss+qld++erGFs ++d+pLDmRmd + ++ tAa+v+n l++++La+il++ygeek+ak+ia aiv++r+ tr|H2Q3C1|H2Q3C1_PANTR 158 VQpgTFDGVLMDLGCSSMQLDTPERGFSLRKDGPLDMRMDGGRypDMpTAADVVNALDQQALASILRTYGEEKHAKKIASAIVQARS 244 76669***********************************87622346*************************************** PP HS---BHHHHHHHHHHHS-------.........-TTHHHHHHHHHHHHTHHHHHHHHHHHHHHHEEEEEEEEEEESSHHHHHHHHH CS Methyltransf_5 172 kkpiettkelaevikkvvpakekrk.......kihpatrvFQAlRIavNdEleeleelleaalellakggrlavisFhSledrivkk 251 pi++t++la +++ ++p+++ + ++h at++FQAlRI+vN+El+el + l++a + l++ggrl+ +sFhSledrivk+ tr|H2Q3C1|H2Q3C1_PANTR 245 IYPITRTQQLAGIVAGAFPPSAVYTrkdllqrSTHIATKTFQALRIFVNNELNELYTGLKTAQKFLRPGGRLVALSFHSLEDRIVKR 331 *******************99877667899999****************************************************** PP HHHHTT------------------..................-EEEEE.E.E---HHHHHH-GGGTT-EEEEEE CS Methyltransf_5 252 vfkeaskkeevPaglpvaseelka..................klklitkkpilpseeeikeNprsrsakLrvie 307 ++ s +e+ ++ + ++++ +li kk++ p ++++++Npr rsakLr++ tr|H2Q3C1|H2Q3C1_PANTR 332 FLLGISMTERFNLSVRQQVMKTSQlgsdhenteevsmrraplMWELIHKKVLSPQDQDVQDNPRGRSAKLRAAI 405 *97666544332222211111111222334455677887777789**************************975 PP >> tr|A0A2I3T743|A0A2I3T743_PANTR Uncharacterized protein OS=Pan troglodytes OX=9598 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 177.2 0.0 6.2e-56 7.1e-52 3 306 .. 71 340 .. 69 343 .] 0.85 Alignments for each domain: == domain 1 score: 177.2 bits; conditional E-value: 6.2e-56 XX--TTTHHHHHHT---TT-EEEES--TTSHHHHHHHTT--TTEEEEEES-HHHHHHHTT.....---EEEEES-GGGH CS Methyltransf_5 3 eHvpVLlqevvellepkedgvyiDaTlGaggHsealLeqlselrligidrDksalarakerLkefekrltlvkgnfkel 81 H+pV ++evv+ l +++D+T+G+ggH++a+L++ s+ l+++d ++a a a++ + + k+++ + g+f++ tr|A0A2I3T743|A0A2I3T743_PANTR 71 LHTPVTVDEVVHCL------FFLDMTFGSGGHTKAILQKESDIVLYALDTNPTAYALAEHLSELYPKQIRAMLGQFSQA 143 69999999999987......9***************************************77788899*********** PP HHHHHHTTTT..-EEEEEEE-S--HHHHH-GGG---SSS--BS---S-TTSS--HHHHHHH--HHHHHHHHHHHH--TT CS Methyltransf_5 82 kealaeanvt..kvdgillDLGvSslqldkaerGFsfaedapLDmRmdqsaeltAaevlnrlseeeLakilkeygeekf 158 ++ l +a+v+ ++dg+l+DLG Ss+qld++er F ++ Dm tAa+v+n l++++La+il++ygeek+ tr|A0A2I3T743|A0A2I3T743_PANTR 144 EALLMKAGVQpgTFDGVLMDLGCSSMQLDTPERFFPLERRRYPDM-------PTAADVVNALDQQALASILRTYGEEKH 215 *******976669********************998888777777.......59************************* PP HHHHHHHHHHHHHHS---BHHHHHHHHHHHS-------..-TTHHHHHHHHHHHHTHHHHHHHHHHHHHHHEEEEEEEE CS Methyltransf_5 159 akriakaivekrkkkpiettkelaevikkvvpakekrkkihpatrvFQAlRIavNdEleeleelleaalellakggrla 237 ak+ia aiv++r ++ + +p++++r +h at++FQAlRI+vN+El+el + l++a + l++gg l+ tr|A0A2I3T743|A0A2I3T743_PANTR 216 AKKIASAIVQAR------------SICLSSLPTSSQR-PTHIATKTFQALRIFVNNELNELYTGLKTAQKFLRPGGHLV 281 ***********9............4556678888875.59*************************************** PP EEESSHHHHHHHHHHHHHTT-------------------EEEEE.E.E---..HHHHHH-GGGTT-EEEEE CS Methyltransf_5 238 visFhSledrivkkvfkeaskkeevPaglpvaseelkaklklitkkpilps..eeeikeNprsrsakLrvi 306 +sFhSl+d i v+k ++ s++ +++ + + p++ ++ + +Npr rsakLr+ tr|A0A2I3T743|A0A2I3T743_PANTR 282 TLSFHSLQDCIQ--VMKTLQL----------SSDHENTEGGSMRRAPLMWEliYKKVLNNPRGRSAKLRAT 340 *********883..3333222..........2233333334444444443211344678999999999975 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (309 nodes) Target sequences: 23008 (11942255 residues searched) Passed MSV filter: 732 (0.031815); expected 460.2 (0.02) Passed bias filter: 393 (0.017081); expected 460.2 (0.02) Passed Vit filter: 39 (0.00169506); expected 23.0 (0.001) Passed Fwd filter: 2 (8.69263e-05); expected 0.2 (1e-05) Initial search space (Z): 23008 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.27u 0.01s 00:00:00.28 Elapsed: 00:00:00.09 # Mc/sec: 41001.74 // [ok]