# 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/Tetrahymena_thermophila.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 7.6e-47 160.6 10.5 7.6e-47 160.6 10.5 1.9 2 tr|Q231D4|Q231D4_TETTS mRNA cap guanine-N7 methyltransferase 3.4e-37 128.9 14.6 1.2e-20 74.5 4.7 2.3 2 tr|Q234L5|Q234L5_TETTS mRNA capping enzyme, large subunit fa Domain annotation for each sequence (and alignments): >> tr|Q231D4|Q231D4_TETTS mRNA cap guanine-N7 methyltransferase OS=Tetrahymena thermophila (strain SB210) OX=312017 GN= # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -5.2 4.1 1.5 2.1e+04 282 312 .. 104 134 .. 89 164 .. 0.49 2 ! 160.6 10.5 5.7e-51 7.6e-47 6 332 .. 179 523 .. 172 524 .. 0.78 Alignments for each domain: == domain 1 score: -5.2 bits; conditional E-value: 1.5 Pox_MCEL 282 elkklikrvsklee.rkstknfeelnrsalks 312 ++kk+ +r + + + r+ + n+ e +rs++++ tr|Q231D4|Q231D4_TETTS 104 DVKKIQERKT-MQQvRERHDNSGEKKRSTMQQ 134 3444333332.222123333443333333332 PP == domain 2 score: 160.6 bits; conditional E-value: 5.7e-51 Pox_MCEL 6 tesevtdlvathynerpetgellkkRiastiiylrnlnNyiKslLiseyasktkl.........desnklkvLdldcgkGgDLeKye 83 + +e ++ +++ y++++ + +R+ + lrn+ N+iKs++i++y k+ e+ +l v++++cg+GgDL+K+ tr|Q231D4|Q231D4_TETTS 179 RVTEQQKTIQSGYEKKEALT--VLTRSVLETAELRNFHNWIKSVIIKKYSDDMKQiikkkedlrTETGQLFVMEIACGQGGDLKKWL 263 55667777888887777655..889999999******************8776665544433444566788**************** PP Pox_MCEL 84 kgeIsklvatDiaeesieeckeRY..nklrsekkskyykFdaelirsd......sfvssvrevfed.......yelkfdivslQFai 155 ++I +v++Di+ +s++e+++R + + +y kF++ +++ d f +++++ + +++ fdivs+Q a+ tr|Q231D4|Q231D4_TETTS 264 HADIGLYVGVDISFNSLKEASRRTkdIMEKLPPHWNYKKFKYGFYQKDgsastdEFWKHIHD---KdkdqdksKRFFFDIVSCQMAM 347 ************************65444445566788899977765411100034444444...155599888999********** PP Pox_MCEL 156 HYsFeseekartvlrNvaellasGGkvlgTtpdgdelskltvkklkvi.ekdlpsfgNsiyvvkfeeekiev.plfGakYeynledv 240 HY F se++ar + N++++l+ G++l+T d++ + k+++ + k +++ fgN+ +++ fe+ + v ++G kYe++l+d+ tr|Q231D4|Q231D4_TETTS 348 HYMFGSEQHARNFFSNATQRLNDQGYLLVTCSDSNAIVKKMRSRGKLDsTNNKYTFGNKYFSMAFENLNFPVgKPYGLKYEFYLQDA 434 *************************************999997666541567889***********99988889************9 PP Pox_MCEL 241 v..........kdvpEyivkfealvklaeeYgleLvdkkdfaeiieesikkelkklikrvskleerkstknfeelnrsalkstlgkd 317 v k pEy+v+ ++l kla eY l+ ++ +f e+++++ ++++ l k++ l++ +++ ++ + tr|Q231D4|Q231D4_TETTS 435 VgekdeatgqiKYTPEYLVELNNLNKLAMEYSLVVKENLNFIEFYKNY-RNKYGYLFKSMG-LDK---------FTNDH--PSIDPE 508 844333332224579*******************************99.677999999998.333.........22222..235556 PP Pox_MCEL 318 eeealsyYlvfvfeK 332 +e+ + Y v vf+K tr|Q231D4|Q231D4_TETTS 509 LWEISHCYRVIVFQK 523 666666666666665 PP >> tr|Q234L5|Q234L5_TETTS mRNA capping enzyme, large subunit family protein OS=Tetrahymena thermophila (strain SB210) O # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 59.5 2.3 3.1e-20 4.2e-16 4 132 .. 255 393 .. 252 404 .. 0.74 2 ! 74.5 4.7 8.7e-25 1.2e-20 142 293 .. 479 641 .. 466 668 .. 0.79 Alignments for each domain: == domain 1 score: 59.5 bits; conditional E-value: 3.1e-20 Pox_MCEL 4 aetesevtdlvathynerpetgellkkRiastiiylrnlnNyiKslLiseyasktkl........desnklkvLdldcgkGgDLeKy 82 ++++s ++ v+++yn+ + +++R++ + +++r + N++Ks+Li++y + ++ n l +Ldl+cgkGgD +K+ tr|Q234L5|Q234L5_TETTS 255 HHHQSSSNNQVQNFYNNNKGFN--QQERKNLETFHIRVFHNWVKSVLIQKYGEIQRSivkaenakLPENYLYILDLACGKGGDHKKW 339 5667788899******988766..99*************************765544333211112346688*************** PP Pox_MCEL 83 ekg.eIsklvatDiaeesieeckeRY.nklrsekkskyykFda..elirsdsfv 132 + + ++++Dia e+++++ +R +l++ k++++y+ + +++ d+ tr|Q234L5|Q234L5_TETTS 340 LMHsHACFYIGVDIAMEALNQAYQRAiITLNEIKQNSRYHSNIkfGFFQKDCSC 393 875156679****************98555555555555444433555555554 PP == domain 2 score: 74.5 bits; conditional E-value: 8.7e-25 Pox_MCEL 142 yelkfdivslQFaiHYsFeseekartvlrNvaellasGGkvlgTtpdgdelskltvkklkviekdlpsfgNsiyvvkfeeekiev.. 226 + fd+vs+Q a+HY +eseek+r l N +++l+ G++l T dg+ + ++ k + + ++ + + +++kf++ + tr|Q234L5|Q234L5_TETTS 479 GKYMFDVVSCQMAMHYMHESEEKVRNFLDNCTKRLNDQGFLLLTFTDGNAVLDIMKSKGQPTPEGGTIYSSKHFSMKFDNPVEQIdl 565 4678**********************************************9777777777777777778888888999965444477 PP Pox_MCEL 227 ..plfGakYeynl........edvvkdvpEyivkfealvklaeeYgleLvdkkdfaeiieesikkelkklikrvskl 293 ++G+kY ++l e+v+k vpEy+v+ + + + a+eY le ++ ++++e++ k + l+k+++ tr|Q234L5|Q234L5_TETTS 566 qqNPYGNKYGFYLqesvgsqvEGVIKYVPEYLVNVDLFIETAKEYSLEVKENNLLTDFYENN-KTPFADLLKQMKVT 641 889********9833333333457899**********************************9.66699999999822 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (333 nodes) Target sequences: 26971 (16912921 residues searched) Passed MSV filter: 5396 (0.200067); expected 539.4 (0.02) Passed bias filter: 150 (0.00556153); expected 539.4 (0.02) Passed Vit filter: 15 (0.000556153); expected 27.0 (0.001) Passed Fwd filter: 2 (7.41537e-05); expected 0.3 (1e-05) Initial search space (Z): 26971 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.42u 0.01s 00:00:00.43 Elapsed: 00:00:00.13 # Mc/sec: 43323.10 // [ok]