# 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/Pox_MCEL.hmm.txt # target sequence database: proteomes/Homo_sapiens.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: Pox_MCEL [M=333] Accession: PF03291.15 Description: mRNA capping enzyme 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2.3e-121 407.3 7.2 3.6e-121 406.6 7.2 1.3 1 sp|O43148|MCES_HUMAN mRNA cap guanine-N7 methyltransferase 5.5e-116 389.6 5.5 8.9e-116 388.9 5.5 1.3 1 sp|O43148-2|MCES_HUMAN Isoform 2 of mRNA cap guanine-N7 meth 2.4e-105 354.6 5.1 2.7e-105 354.4 5.1 1.0 1 tr|K7EP06|K7EP06_HUMAN mRNA cap guanine-N7 methyltransferase 2.8e-73 249.2 5.1 3.1e-73 249.0 5.1 1.0 1 tr|F5H2D9|F5H2D9_HUMAN mRNA cap guanine-N7 methyltransferase ------ inclusion threshold ------ 0.042 15.2 0.0 0.57 11.5 0.0 2.0 1 sp|Q9NZJ6|COQ3_HUMAN Hexaprenyldihydroxybenzoate methyltra 0.14 13.5 0.2 0.25 12.7 0.2 1.4 1 tr|K7ERH6|K7ERH6_HUMAN mRNA cap guanine-N7 methyltransferase Domain annotation for each sequence (and alignments): >> sp|O43148|MCES_HUMAN mRNA cap guanine-N7 methyltransferase OS=Homo sapiens GN=RNMT PE=1 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 406.6 7.2 2.4e-125 3.6e-121 2 333 .] 137 475 .. 136 475 .. 0.98 Alignments for each domain: == domain 1 score: 406.6 bits; conditional E-value: 2.4e-125 Pox_MCEL 2 geaetesevtdlvathynerpetgellkkRiastiiylrnlnNyiKslLiseyasktkldesnklkvLdldcgkGgDLeKyekgeIskl 90 ++++ e+++++ va+hyne++e g l+kR++s+i+ylrn+nN++Ks+Li+e+++k++++++++++vLdl+cgkGgDL K++kg+I+kl sp|O43148|MCES_HUMAN 137 KQKNLEEGHSSTVAAHYNELQEVG--LEKRSQSRIFYLRNFNNWMKSVLIGEFLEKVRQKKKRDITVLDLGCGKGGDLLKWKKGRINKL 223 67788999****************..*************************************************************** PP Pox_MCEL 91 vatDiaeesieeckeRYnklrsekkskyykFdaelirsdsfvssvrevfedyelkfdivslQFaiHYsFeseekartvlrNvaellasG 179 v+tDia++s+++c++RY+ +++++ s+ y F+ae+i++ds ++ + ++f+d+++ fdi+s+QF++HYsFes e+a+++lrN++e+l +G sp|O43148|MCES_HUMAN 224 VCTDIADVSVKQCQQRYEDMKNRRDSE-YIFSAEFITADSSKELLIDKFRDPQMCFDICSCQFVCHYSFESYEQADMMLRNACERLSPG 311 *****************9999999996.9************************************************************ PP Pox_MCEL 180 GkvlgTtpdgdelskltvkklkviekdlpsfgNsiyvvkfeeekievplfGakYeynledvvkdvpEyivkfealvklaeeYgleLvdk 268 G+++gTtp+++el +++l+++e ++sfgN+iy+vkf+ +k +plfG+kY++nle+vv dvpE++v+f+ l ++a++Y+++Lv+k sp|O43148|MCES_HUMAN 312 GYFIGTTPNSFEL----IRRLEASE--TESFGNEIYTVKFQ-KKGDYPLFGCKYDFNLEGVV-DVPEFLVYFPLLNEMAKKYNMKLVYK 392 *************....*******6..**************.89999***************.************************** PP Pox_MCEL 269 kdfaeiieesikke.lkklikrvskleerkstknf.................eelnrsalkstlgkdeeealsyYlvfvfeKr 333 k+f e++ee+ik++ +k+l+kr+++le +++++ + ++++++ + +tl+k+e+ea+s+Ylvf+feK+ sp|O43148|MCES_HUMAN 393 KTFLEFYEEKIKNNeNKMLLKRMQALEPYPANESSklvsekvddyehaakymKNSQVRLPLGTLSKSEWEATSIYLVFAFEKQ 475 *************99******************999***************999****************************6 PP >> sp|O43148-2|MCES_HUMAN Isoform 2 of mRNA cap guanine-N7 methyltransferase OS=Homo sapiens GN=RNMT # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 388.9 5.5 6.1e-120 8.9e-116 2 326 .. 137 468 .. 136 472 .. 0.97 Alignments for each domain: == domain 1 score: 388.9 bits; conditional E-value: 6.1e-120 Pox_MCEL 2 geaetesevtdlvathynerpetgellkkRiastiiylrnlnNyiKslLiseyasktkldesnklkvLdldcgkGgDLeKyekgeIs 88 ++++ e+++++ va+hyne++e g l+kR++s+i+ylrn+nN++Ks+Li+e+++k++++++++++vLdl+cgkGgDL K++kg+I+ sp|O43148-2|MCES_HUMAN 137 KQKNLEEGHSSTVAAHYNELQEVG--LEKRSQSRIFYLRNFNNWMKSVLIGEFLEKVRQKKKRDITVLDLGCGKGGDLLKWKKGRIN 221 67788999****************..************************************************************* PP Pox_MCEL 89 klvatDiaeesieeckeRYnklrsekkskyykFdaelirsdsfvssvrevfedyelkfdivslQFaiHYsFeseekartvlrNvael 175 klv+tDia++s+++c++RY+ +++++ s+ y F+ae+i++ds ++ + ++f+d+++ fdi+s+QF++HYsFes e+a+++lrN++e+ sp|O43148-2|MCES_HUMAN 222 KLVCTDIADVSVKQCQQRYEDMKNRRDSE-YIFSAEFITADSSKELLIDKFRDPQMCFDICSCQFVCHYSFESYEQADMMLRNACER 307 *******************9999999996.9******************************************************** PP Pox_MCEL 176 lasGGkvlgTtpdgdelskltvkklkviekdlpsfgNsiyvvkfeeekievplfGakYeynledvvkdvpEyivkfealvklaeeYg 262 l +GG+++gTtp+++el +++l+++e ++sfgN+iy+vkf+ +k +plfG+kY++nle+vv dvpE++v+f+ l ++a++Y+ sp|O43148-2|MCES_HUMAN 308 LSPGGYFIGTTPNSFEL----IRRLEASE--TESFGNEIYTVKFQ-KKGDYPLFGCKYDFNLEGVV-DVPEFLVYFPLLNEMAKKYN 386 *****************....*******6..**************.89999***************.******************** PP Pox_MCEL 263 leLvdkkdfaeiieesikke.lkklikrvskleerkstknf.................eelnrsalkstlgkdeeealsyYl 326 ++Lv+kk+f e++ee+ik++ +k+l+kr+++le +++++ + ++++++ + +tl+k+e+ea+++ + sp|O43148-2|MCES_HUMAN 387 MKLVYKKTFLEFYEEKIKNNeNKMLLKRMQALEPYPANESSklvsekvddyehaakymKNSQVRLPLGTLSKSEWEATRLTV 468 *******************99******************999***************9999****************99755 PP >> tr|K7EP06|K7EP06_HUMAN mRNA cap guanine-N7 methyltransferase OS=Homo sapiens GN=RNMT PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 354.4 5.1 1.8e-109 2.7e-105 46 333 .] 1 297 [. 1 297 [. 0.98 Alignments for each domain: == domain 1 score: 354.4 bits; conditional E-value: 1.8e-109 Pox_MCEL 46 iKslLiseyasktkldesnklkvLdldcgkGgDLeKyekgeIsklvatDiaeesieeckeRYnklrsekkskyykFdaelirsdsfv 132 +Ks+Li+e+++k++++++++++vLdl+cgkGgDL K++kg+I+klv+tDia++s+++c++RY+ +++++ s+ y F+ae+i++ds + tr|K7EP06|K7EP06_HUMAN 1 MKSVLIGEFLEKVRQKKKRDITVLDLGCGKGGDLLKWKKGRINKLVCTDIADVSVKQCQQRYEDMKNRRDSE-YIFSAEFITADSSK 86 7*************************************************************9999999996.9************* PP Pox_MCEL 133 ssvrevfedyelkfdivslQFaiHYsFeseekartvlrNvaellasGGkvlgTtpdgdelskltvkklkviekdlpsfgNsiyvvkf 219 + + ++f+d+++ fdi+s+QF++HYsFes e+a+++lrN++e+l +GG+++gTtp+++el +++l+++e ++sfgN+iy+vkf tr|K7EP06|K7EP06_HUMAN 87 ELLIDKFRDPQMCFDICSCQFVCHYSFESYEQADMMLRNACERLSPGGYFIGTTPNSFEL----IRRLEASE--TESFGNEIYTVKF 167 ************************************************************....*******6..************* PP Pox_MCEL 220 eeekievplfGakYeynledvvkdvpEyivkfealvklaeeYgleLvdkkdfaeiieesikke.lkklikrvskleerkstknf... 302 + +k +plfG+kY++nle+vv dvpE++v+f+ l ++a++Y+++Lv+kk+f e++ee+ik++ +k+l+kr+++le +++++ + tr|K7EP06|K7EP06_HUMAN 168 Q-KKGDYPLFGCKYDFNLEGVV-DVPEFLVYFPLLNEMAKKYNMKLVYKKTFLEFYEEKIKNNeNKMLLKRMQALEPYPANESSklv 252 *.89999***************.***************************************99******************999** PP Pox_MCEL 303 ..............eelnrsalkstlgkdeeealsyYlvfvfeKr 333 ++++++ + +tl+k+e+ea+s+Ylvf+feK+ tr|K7EP06|K7EP06_HUMAN 253 sekvddyehaakymKNSQVRLPLGTLSKSEWEATSIYLVFAFEKQ 297 *************999****************************6 PP >> tr|F5H2D9|F5H2D9_HUMAN mRNA cap guanine-N7 methyltransferase OS=Homo sapiens GN=RNMT PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 249.0 5.1 2.1e-77 3.1e-73 115 333 .] 5 233 .. 1 233 [. 0.97 Alignments for each domain: == domain 1 score: 249.0 bits; conditional E-value: 2.1e-77 Pox_MCEL 115 kskyykFdaelirsdsfvssvrevfedyelkfdivslQFaiHYsFeseekartvlrNvaellasGGkvlgTtpdgdelskltvkklk 201 + + y F+ae+i++ds ++ + ++f+d+++ fdi+s+QF++HYsFes e+a+++lrN++e+l +GG+++gTtp+++el +++l+ tr|F5H2D9|F5H2D9_HUMAN 5 RDSEYIFSAEFITADSSKELLIDKFRDPQMCFDICSCQFVCHYSFESYEQADMMLRNACERLSPGGYFIGTTPNSFEL----IRRLE 87 344699************************************************************************....***** PP Pox_MCEL 202 viekdlpsfgNsiyvvkfeeekievplfGakYeynledvvkdvpEyivkfealvklaeeYgleLvdkkdfaeiieesikke.lkkli 287 ++e ++sfgN+iy+vkf+ +k +plfG+kY++nle+vv dvpE++v+f+ l ++a++Y+++Lv+kk+f e++ee+ik++ +k+l+ tr|F5H2D9|F5H2D9_HUMAN 88 ASE--TESFGNEIYTVKFQ-KKGDYPLFGCKYDFNLEGVV-DVPEFLVYFPLLNEMAKKYNMKLVYKKTFLEFYEEKIKNNeNKMLL 170 **6..**************.89999***************.***************************************99***** PP Pox_MCEL 288 krvskleerkstknf.................eelnrsalkstlgkdeeealsyYlvfvfeKr 333 kr+++le +++++ + ++++++ + +tl+k+e+ea+s+Ylvf+feK+ tr|F5H2D9|F5H2D9_HUMAN 171 KRMQALEPYPANESSklvsekvddyehaakymKNSQVRLPLGTLSKSEWEATSIYLVFAFEKQ 233 *************999***************999****************************6 PP >> sp|Q9NZJ6|COQ3_HUMAN Hexaprenyldihydroxybenzoate methyltransferase, mitochondrial OS=Homo sapiens GN=COQ3 PE=1 SV=3 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.5 0.0 3.9e-05 0.57 60 196 .. 142 262 .. 116 333 .. 0.66 Alignments for each domain: == domain 1 score: 11.5 bits; conditional E-value: 3.9e-05 Pox_MCEL 60 ldesnklkvLdldcgkGgDLeKyekgeI.sklvatDiaeesieeckeRYnklrsekkskyykFda.elirsdsfvssvrevfedyelkf 146 + +k+Ld++cg G L g++ ++ +++D+ +e i++++ s Fd r + v s++e++e+++ f sp|Q9NZJ6|COQ3_HUMAN 142 GKPLLGMKILDVGCGGG--LLTEPLGRLgASVIGIDPVDENIKTAQCHK--------S----FDPvLDKRIEYRVCSLEEIVEETAETF 216 55555677888888643..555555555256788888888887665322........2....332122233444456666666666777 PP Pox_MCEL 147 divslQFaiHYsFeseekartvlrNvaellasGGkvlgTtpdgdelsklt 196 d v + ++ e +t l+ +++l++GG +++Tt++ +ls+ + sp|Q9NZJ6|COQ3_HUMAN 217 DAVVASEVV----EHVIDLETFLQCCCQVLKPGGSLFITTINKTQLSYAL 262 765544443....23345668999******************99998754 PP >> tr|K7ERH6|K7ERH6_HUMAN mRNA cap guanine-N7 methyltransferase (Fragment) OS=Homo sapiens GN=RNMT PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.7 0.2 1.7e-05 0.25 275 326 .. 2 78 .. 1 82 [. 0.86 Alignments for each domain: == domain 1 score: 12.7 bits; conditional E-value: 1.7e-05 Pox_MCEL 275 ieesikke.lkklikrvskleerkstknf........................eelnrsalkstlgkdeeealsyYl 326 +ee+ik++ +k+l+kr+++le+ + ++++++ + +tl+k+e+ea+++ + tr|K7ERH6|K7ERH6_HUMAN 2 YEEKIKNNeNKMLLKRMQALEDFHFFPQPypanessklvsekvddyehaakymKNSQVRLPLGTLSKSEWEATRLTV 78 689999998************987655546778889999999999999999999999999999********998755 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (333 nodes) Target sequences: 88473 (35070517 residues searched) Passed MSV filter: 3258 (0.0368248); expected 1769.5 (0.02) Passed bias filter: 1555 (0.017576); expected 1769.5 (0.02) Passed Vit filter: 156 (0.00176325); expected 88.5 (0.001) Passed Fwd filter: 6 (6.78173e-05); expected 0.9 (1e-05) Initial search space (Z): 88473 [actual number of targets] Domain search space (domZ): 6 [number of targets reported over threshold] # CPU time: 0.81u 0.03s 00:00:00.84 Elapsed: 00:00:00.25 # Mc/sec: 46713.93 // [ok]