# 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/Monodelphis_domestica.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 4.1e-60 202.9 0.0 7.7e-60 202.0 0.0 1.5 1 tr|F7FEI2|F7FEI2_MONDO Methyltransferase like 3 OS=Monodelph 7.1e-41 140.2 0.0 1.1e-40 139.7 0.0 1.2 1 tr|F7G0D7|F7G0D7_MONDO Methyltransferase like 14 OS=Monodelp 2.8e-40 138.3 0.8 4.5e-40 137.6 0.8 1.3 1 tr|F7GDI0|F7GDI0_MONDO Methyltransferase like 4 OS=Monodelph Domain annotation for each sequence (and alignments): >> tr|F7FEI2|F7FEI2_MONDO Methyltransferase like 3 OS=Monodelphis domestica OX=13616 GN=METTL3 PE=3 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 202.0 0.0 1.1e-63 7.7e-60 1 172 [] 389 550 .. 389 550 .. 0.98 Alignments for each domain: == domain 1 score: 202.0 bits; conditional E-value: 1.1e-63 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|F7FEI2|F7FEI2_MONDO 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+++i++++r++s+KP+e++ ++e+l++++rk+elF+R ++p+w+t+Gn+ tr|F7FEI2|F7FEI2_MONDO 469 TGRTGHWLNHGKEHCLVGVKGNPQ-----GFNRGLDCEVIVAEVRSTSHKPDEIYGMIERLSPGTRKIELFGRPhnVQPNWITLGNQ 550 ***********************8.....67899***************************************999*********97 PP >> tr|F7G0D7|F7G0D7_MONDO Methyltransferase like 14 OS=Monodelphis domestica OX=13616 GN=METTL14 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 139.7 0.0 1.5e-44 1.1e-40 61 171 .. 182 298 .. 172 299 .. 0.94 Alignments for each domain: == domain 1 score: 139.7 bits; conditional E-value: 1.5e-44 HHHHHCT-EEEEEEEEEEE-TTSCBS-..BS-TTSSSBB-EEEEEEEEES---TTT.T.TB-TTSSBSEEEE..E-CTTTSS-THHH CS MT-A70 61 ellkkWglelveeltwlkinklgepia..tlesgavlkkskEhllvgvkgnvkrstdk.iinknvdtdliis..perehsrKPpelf 142 +l+kWg++++e+++w+k+n +++p++ tl+++av++++kEh+l+g+kg+v+rstd+ +i++nvd+dlii+ pe+++++KP e+f tr|F7G0D7|F7G0D7_MONDO 182 ICLRKWGYRRCEDICWIKTN-KNNPGKtkTLDPKAVFQRTKEHCLMGIKGTVRRSTDGdFIHANVDIDLIITeePEIGNIEKPVEIF 267 59******************.556555569******************************************999************ PP HHHHHHCTCS-EEEET--G..--TTEEEEET CS MT-A70 143 ellehllesarkLelFaRs..lrpgwltvGn 171 +++eh++ ++r+L+lF+R+ +rpgwltvG+ tr|F7G0D7|F7G0D7_MONDO 268 HIIEHFCLGRRRLHLFGRDstIRPGWLTVGP 298 ******************************7 PP >> tr|F7GDI0|F7GDI0_MONDO Methyltransferase like 4 OS=Monodelphis domestica OX=13616 GN=METTL4 PE=3 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 137.6 0.8 6.4e-44 4.5e-40 1 172 [] 293 467 .. 293 467 .. 0.91 Alignments for each domain: == domain 1 score: 137.6 bits; conditional E-value: 6.4e-44 ECEEEE---BHHCC..HTT--SS--B-HHHHHTS-GGGCECCEEEEEEEEETTCHHHH..HHHHHHHCT-EEEEEEEEEEE-TTSCB CS MT-A70 1 fdviiidPpweeklvarkkelpyktlneeelkaLpipelldkrslvflWvtnrellel..krellkkWglelveeltwlkinklgep 85 +dvi+idPpw++k+v+r+k +y+ l+ ++k++pip+l+++++lv Wvtnr+++ + k+el + W +e ++e+ w+ki++ ge tr|F7GDI0|F7GDI0_MONDO 293 YDVIVIDPPWQNKSVKRSK--RYNYLSSWQIKQIPIPKLAAPNCLVITWVTNRQKHLQfvKEELYPYWSVETLTEWHWVKITRSGEF 377 8*****************9..9999*******************************999999************************* PP S-BS-TTSSSBB-EEEEEEE.EES---TTT.TTB-TTS....SBSEEEE.E-CTTTSS-THHHHHHHHHCTCS-EEEET--G--TTE CS MT-A70 86 iatlesgavlkkskEhllvg.vkgnvkrstdkiinknv....dtdliis.perehsrKPpelfellehllesarkLelFaRslrpgw 166 +++l+s +kk++E l++g v+ ++ ++ + ++ v d +li+s p + hs+KPp+ l ++++a++LelFaR+l+pgw tr|F7GDI0|F7GDI0_MONDO 378 VFPLDSP--HKKPYEGLVLGrVRRKMVSTLRE-PEMTVlpipDHKLIVSvPCTLHSHKPPLTEVLKDYIVPDAKCLELFARNLQPGW 461 **98888..99*********445554444444.44444567999*****99*********9999999999***************** PP EEEETT CS MT-A70 167 ltvGne 172 +++Gne tr|F7GDI0|F7GDI0_MONDO 462 TSWGNE 467 ****97 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (172 nodes) Target sequences: 21272 (11266248 residues searched) Passed MSV filter: 551 (0.0259026); expected 425.4 (0.02) Passed bias filter: 352 (0.0165476); expected 425.4 (0.02) Passed Vit filter: 27 (0.00126927); expected 21.3 (0.001) Passed Fwd filter: 3 (0.00014103); expected 0.2 (1e-05) Initial search space (Z): 21272 [actual number of targets] Domain search space (domZ): 3 [number of targets reported over threshold] # CPU time: 0.15u 0.01s 00:00:00.16 Elapsed: 00:00:00.08 # Mc/sec: 24222.43 // [ok]