# 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/Pan_troglodytes.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.9e-63 213.8 0.0 2.7e-63 213.3 0.0 1.2 1 tr|H2QQ29|H2QQ29_PANTR METTL14 isoform 1 OS=Pan troglod 1.2e-60 204.7 0.0 2.3e-60 203.8 0.0 1.5 1 tr|G2HF97|G2HF97_PANTR METTL3 isoform 2 OS=Pan troglody 1.9e-39 135.7 0.6 5.1e-39 134.3 0.1 1.8 2 tr|A0A2I3T433|A0A2I3T433_PANTR Methyltransferase like 4 OS=Pan Domain annotation for each sequence (and alignments): >> tr|H2QQ29|H2QQ29_PANTR METTL14 isoform 1 OS=Pan troglodytes OX=9598 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 3.5e-67 2.7e-63 1 171 [. 186 362 .. 186 363 .. 0.96 Alignments for each domain: == domain 1 score: 213.3 bits; conditional E-value: 3.5e-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|H2QQ29|H2QQ29_PANTR 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|H2QQ29|H2QQ29_PANTR 272 tkTLDPKAVFQRTKEHCLMGIKGTVKRSTDGdFIHANVDIDLIITeePEIGNIEKPVEIFHIIEHFCLGRRRLHLFGRDstIRPGWL 358 569******************************************999*************************************** PP EEET CS MT-A70 168 tvGn 171 tvG+ tr|H2QQ29|H2QQ29_PANTR 359 TVGP 362 ***7 PP >> tr|G2HF97|G2HF97_PANTR METTL3 isoform 2 OS=Pan troglodytes OX=9598 GN=METTL3 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 203.8 0.0 3e-64 2.3e-60 1 172 [] 389 550 .. 389 550 .. 0.98 Alignments for each domain: == domain 1 score: 203.8 bits; conditional E-value: 3e-64 ECEEEE---BHHCC..HTT--SS--B-HHHHHTS-GGGCECCEEEEEEEEETTCHHHHHHHHHHHCT-EEEEEEEEEEE-TTSCBS- CS MT-A70 1 fdviiidPpweeklvarkkelpyktlneeelkaLpipelldkrslvflWvtnrellelkrellkkWglelveeltwlkinklgepia 87 f+v+++dPpw++++ elpy+tl+++e+++L+ip+l+d++ ++flWvt+r+++ + re+l+ Wg+e+v+e++w+k+n+l+++i+ tr|G2HF97|G2HF97_PANTR 389 FAVVMADPPWDIHM-----ELPYGTLTDDEMRRLNIPVLQDDG-FLFLWVTGRAMELG-RECLNLWGYERVDEIIWVKTNQLQRIIR 468 89************.....9*********************99.*************9.**************************** PP BS-TTSSSBB-EEEEEEEEES---TTT.TTB-TTSSBSEEEEE-CTTTSS-THHHHHHHHHCTCS-EEEET--G..--TTEEEEETT CS MT-A70 88 tlesgavlkkskEhllvgvkgnvkrstdkiinknvdtdliisperehsrKPpelfellehllesarkLelFaRs..lrpgwltvGne 172 t+++g++l++ kEh+lvgvkgn++ +n++ d+d+i++++r++s+KP+e++ ++e+l++++rk+elF+R ++p+w+t+Gn+ tr|G2HF97|G2HF97_PANTR 469 TGRTGHWLNHGKEHCLVGVKGNPQ-----GFNQGLDCDVIVAEVRSTSHKPDEIYGMIERLSPGTRKIELFGRPhnVQPNWITLGNQ 550 ***********************8.....67899***************************************999*********97 PP >> tr|A0A2I3T433|A0A2I3T433_PANTR Methyltransferase like 4 OS=Pan troglodytes OX=9598 GN=METTL4 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -1.3 0.0 0.29 2.2e+03 92 112 .. 58 78 .. 46 96 .. 0.72 2 ! 134.3 0.1 6.6e-43 5.1e-39 1 172 [] 232 406 .. 232 406 .. 0.92 Alignments for each domain: == domain 1 score: -1.3 bits; conditional E-value: 0.29 TSSSBB-EEEEEEEEES---T CS MT-A70 92 gavlkkskEhllvgvkgnvkr 112 + + + kE l+ gvk++ + tr|A0A2I3T433|A0A2I3T433_PANTR 58 ECQQSNEKEDLMNGVKKEISI 78 445667899999999999864 PP == domain 2 score: 134.3 bits; conditional E-value: 6.6e-43 ECEEEE---BHHCC..HTT--SS--B-HHHHHTS-GGGCECCEEEEEEEEETTCHHHH..HHHHHHHCT-EEEEEEEEE CS MT-A70 1 fdviiidPpweeklvarkkelpyktlneeelkaLpipelldkrslvflWvtnrellel..krellkkWglelveeltwl 77 fdvi+idPpw++k+v+r++ +y+ l+ ++k++pip+l+++++l++ Wvtnr+++ k+el ++W +e+v+e+ w+ tr|A0A2I3T433|A0A2I3T433_PANTR 232 FDVIVIDPPWQNKSVKRSN--RYSYLSPLQIKQIPIPKLAAPNCLLVTWVTNRQKHLRfiKEELYPSWSVEVVAEWHWV 308 9*****************8..9***********************************99999***************** PP EE-TTSCBS-BS-TTSSSBB-EEEEEEEEES---TTT.TTB-TTS....SBSEEEE.E-CTTTSS-THHHHHHHHHCTC CS MT-A70 78 kinklgepiatlesgavlkkskEhllvgvkgnvkrstdkiinknv....dtdliis.perehsrKPpelfellehlles 151 ki++ ge +++l+s +kk++E l++g ++ ++ + nv d +li+s p + hs+KPp+ l + ++ tr|A0A2I3T433|A0A2I3T433_PANTR 309 KITNSGEFVFPLDSP--HKKPYEGLILGRVQEKTALPLRNADVNVlpipDHKLIVSvPCTLHSHKPPLAEVLKDYIKPD 385 **********98888..99*********88877765544444444778999*****99*********9999999999** PP S-EEEET--G--TTEEEEETT CS MT-A70 152 arkLelFaRslrpgwltvGne 172 + LelFaR+l+pgw+++Gne tr|A0A2I3T433|A0A2I3T433_PANTR 386 GEYLELFARNLQPGWTSWGNE 406 *******************97 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (172 nodes) Target sequences: 23008 (11942255 residues searched) Passed MSV filter: 576 (0.0250348); expected 460.2 (0.02) Passed bias filter: 389 (0.0169072); expected 460.2 (0.02) Passed Vit filter: 21 (0.000912726); expected 23.0 (0.001) Passed Fwd filter: 3 (0.000130389); expected 0.2 (1e-05) Initial search space (Z): 23008 [actual number of targets] Domain search space (domZ): 3 [number of targets reported over threshold] # CPU time: 0.14u 0.01s 00:00:00.15 Elapsed: 00:00:00.08 # Mc/sec: 25675.85 // [ok]