# 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/MT-A70.hmm.txt # target sequence database: proteomes/Macaca_mulatta.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: MT-A70 [M=172] Accession: PF05063.13 Description: MT-A70 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.7e-63 213.8 0.0 2.4e-63 213.3 0.0 1.2 1 tr|G7MRV0|G7MRV0_MACMU Methyltransferase like 14 OS=Macaca m 1.8e-39 135.7 0.4 5.4e-39 134.1 0.1 1.9 2 tr|F6YE47|F6YE47_MACMU Methyltransferase like 4 OS=Macaca mu Domain annotation for each sequence (and alignments): >> tr|G7MRV0|G7MRV0_MACMU Methyltransferase like 14 OS=Macaca mulatta OX=9544 GN=METTL14 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 213.3 0.0 2.3e-67 2.4e-63 1 171 [. 186 362 .. 186 363 .. 0.96 Alignments for each domain: == domain 1 score: 213.3 bits; conditional E-value: 2.3e-67 ECEEEE---BHHCC..HTT--SS--B-HHHHHTS-GGGCECCEEEEEEEEETTCHHHHHHHHHHHCT-EEEEEEEEEEE-TTSCBS- CS MT-A70 1 fdviiidPpweeklvarkkelpyktlneeelkaLpipelldkrslvflWvtnrellelkrellkkWglelveeltwlkinklgepia 87 fdvi+++Pp+ee++++ ++++k+++++++++L+i+e++++rs++flW+++ e+l+l+r++l+kWg++++e+++w+k+n +++p++ tr|G7MRV0|G7MRV0_MACMU 186 FDVILLEPPLEEYYRETGITANEKCWTWDDIMKLEIDEIAAPRSFIFLWCGSGEGLDLGRVCLRKWGYRRCEDICWIKTN-KNNPGK 271 9*************66555*************************************************************.556555 PP ..BS-TTSSSBB-EEEEEEEEES---TTT.T.TB-TTSSBSEEEE..E-CTTTSS-THHHHHHHHHCTCS-EEEET--G..--TTEE CS MT-A70 88 ..tlesgavlkkskEhllvgvkgnvkrstdk.iinknvdtdliis..perehsrKPpelfellehllesarkLelFaRs..lrpgwl 167 tl+++av++++kEh+l+g+kg+vkrstd+ +i++nvd+dlii+ pe+++++KP e+f+++eh++ ++r+L+lF+R+ +rpgwl tr|G7MRV0|G7MRV0_MACMU 272 tkTLDPKAVFQRTKEHCLMGIKGTVKRSTDGdFIHANVDIDLIITeePEIGNIEKPVEIFHIIEHFCLGRRRLHLFGRDstIRPGWL 358 569******************************************999*************************************** PP EEET CS MT-A70 168 tvGn 171 tvG+ tr|G7MRV0|G7MRV0_MACMU 359 TVGP 362 ***7 PP >> tr|F6YE47|F6YE47_MACMU Methyltransferase like 4 OS=Macaca mulatta OX=9544 GN=METTL4 PE=3 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -1.2 0.0 0.19 2e+03 91 113 .. 106 128 .. 94 147 .. 0.71 2 ! 134.1 0.1 5.1e-43 5.4e-39 1 172 [] 281 455 .. 281 455 .. 0.92 Alignments for each domain: == domain 1 score: -1.2 bits; conditional E-value: 0.19 TTSSSBB-EEEEEEEEES---TT CS MT-A70 91 sgavlkkskEhllvgvkgnvkrs 113 + + + kE l+ gvk++ + s tr|F6YE47|F6YE47_MACMU 106 KECQQSNEKEDLMNGVKKEISIS 128 44456678999999999998744 PP == domain 2 score: 134.1 bits; conditional E-value: 5.1e-43 ECEEEE---BHHCC..HTT--SS--B-HHHHHTS-GGGCECCEEEEEEEEETTCHHHH..HHHHHHHCT-EEEEEEEEEEE-TTSCB CS MT-A70 1 fdviiidPpweeklvarkkelpyktlneeelkaLpipelldkrslvflWvtnrellel..krellkkWglelveeltwlkinklgep 85 fdvi+idPpw++k+v+r++ +y+ l+ ++k++pip+l+++++l++ Wvtnr+++ k+el ++W +e+v+e+ w+ki++ ge tr|F6YE47|F6YE47_MACMU 281 FDVIVIDPPWQNKSVKRSN--RYSYLSPLQIKQIPIPKLAAPNCLLVTWVTNRQKHLRfiKEELYPSWSVEVVAEWHWVKITNSGEF 365 9*****************8..9***********************************99999************************* PP S-BS-TTSSSBB-EEEEEEEEES---TTT.TTB-TTS....SBSEEEE.E-CTTTSS-THHHHHHHHHCTCS-EEEET--G--TTEE CS MT-A70 86 iatlesgavlkkskEhllvgvkgnvkrstdkiinknv....dtdliis.perehsrKPpelfellehllesarkLelFaRslrpgwl 167 +++l+s +kk++E l++g ++ ++ + nv d +li+s p + hs+KPp+ l + ++ + LelFaR+l+pgw+ tr|F6YE47|F6YE47_MACMU 366 VFPLDSP--HKKPYEGLILGRVQENTALPLRNADVNVlpipDHKLIVSvPCTLHSHKPPLAEVLKDYIKPDGEYLELFARNLQPGWT 450 **98888..99*********77666654433344444778899*****99*********9999999999****************** PP EEETT CS MT-A70 168 tvGne 172 ++Gne tr|F6YE47|F6YE47_MACMU 451 SWGNE 455 ***97 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (172 nodes) Target sequences: 21211 (10888574 residues searched) Passed MSV filter: 545 (0.0256942); expected 424.2 (0.02) Passed bias filter: 383 (0.0180567); expected 424.2 (0.02) Passed Vit filter: 17 (0.000801471); expected 21.2 (0.001) Passed Fwd filter: 2 (9.42907e-05); expected 0.2 (1e-05) Initial search space (Z): 21211 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.14u 0.01s 00:00:00.15 Elapsed: 00:00:00.06 # Mc/sec: 31213.91 // [ok]