# 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/Anolis_carolinensis.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1.2e-117 392.8 7.4 1.7e-117 392.3 7.4 1.1 1 tr|R4GB35|R4GB35_ANOCA mRNA cap guanine-N7 methyltransferase ------ inclusion threshold ------ 0.14 11.3 0.0 0.27 10.3 0.0 1.4 1 tr|G1KEE2|G1KEE2_ANOCA Glycine N-methyltransferase OS=Anolis Domain annotation for each sequence (and alignments): >> tr|R4GB35|R4GB35_ANOCA mRNA cap guanine-N7 methyltransferase OS=Anolis carolinensis OX=28377 GN=RNMT PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 392.3 7.4 1.8e-121 1.7e-117 4 333 .] 179 515 .. 176 515 .. 0.97 Alignments for each domain: == domain 1 score: 392.3 bits; conditional E-value: 1.8e-121 Pox_MCEL 4 aetesevtdlvathynerpetgellkkRiastiiylrnlnNyiKslLiseyasktkldesnklkvLdldcgkGgDLeKyekgeIskl 90 +++e+++t++va+hyn+++e g l++R+ s+i+ylrn+nN+iKs+Li+e++ k++ ++++++vLdl+cgkGgDL K++kg+Is+l tr|R4GB35|R4GB35_ANOCA 179 ESEEKAHTEAVAAHYNKLEEVG--LQQRSMSRIFYLRNFNNWIKSVLIGEFLDKVRR-RKQNITVLDLGCGKGGDLLKWRKGRISRL 262 567899****************..********************************9.8999************************* PP Pox_MCEL 91 vatDiaeesieeckeRYnklrsekkskyykFdaelirsdsfvssvrevfedyelkfdivslQFaiHYsFeseekartvlrNvaella 177 v+tDia+ s+++c++RY+ +r++ +++ F+ae+i++ds+++ +++++++++++fdi+s+QF++HYsFe+ e+a+++l+N++e+l+ tr|R4GB35|R4GB35_ANOCA 263 VCTDIAAISVQQCEQRYQDMRNRGRRDERIFSAEFITADSTKELLSQKYKNQDMYFDICSCQFVYHYSFETYEQADMMLQNACERLC 349 *****************8888888888899********************************************************* PP Pox_MCEL 178 sGGkvlgTtpdgdelskltvkklkviekdlpsfgNsiyvvkfeeekievplfGakYeynledvvkdvpEyivkfealvklaeeYgle 264 +GG+++gTtpd++el vk+l++++ + sfgN+iy+vkfe +k ++ +fG+kY++nle+vv +vpE++v+f+ l ++a++Yg++ tr|R4GB35|R4GB35_ANOCA 350 PGGYFIGTTPDSYEL----VKRLEASQ--TDSFGNEIYTVKFE-KKGQYDMFGCKYDFNLEGVV-NVPEFLVYFPLLEEMAKKYGMK 428 ***************....*******6..99************.777778**************.********************** PP Pox_MCEL 265 Lvdkkdfaeiieesikke.lkklikrvskleerkstknf.................eelnrsalkstlgkdeeealsyYlvfvfeKr 333 L++kk+f+e++ee+ik+ k l+++++ +e++++ +n+ ++ +++ + +t++++e+ea+s+Yl+f+feK+ tr|R4GB35|R4GB35_ANOCA 429 LIYKKTFQEFYEEKIKNGeHKLLLQQMQGMEKYPADENStlasdkkedyehaevlrKDKEVKLPLGTMSTSEWEATSIYLIFAFEKQ 515 *****************989999*****************************9999889999999*********************6 PP >> tr|G1KEE2|G1KEE2_ANOCA Glycine N-methyltransferase OS=Anolis carolinensis OX=28377 GN=GNMT PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.3 0.0 2.9e-05 0.27 62 192 .. 55 182 .. 45 186 .. 0.73 Alignments for each domain: == domain 1 score: 10.3 bits; conditional E-value: 2.9e-05 Pox_MCEL 62 esnklkvLdldcgkGgDLeKyekgeIsklvatDiaeesieec.keRYnklrsekkskyykFdaelirsdsfvssvrevfedyelkfd 147 +++ ++vLd++cg G D + + + ++D +++ ++ + keR+ +++k++ Fd +i ++++ ++++ fd tr|G1KEE2|G1KEE2_ANOCA 55 SQRCRTVLDVACGTGVDSIMLVEEGF-QVTSVDASDKMLKYAlKERW---NRRKEEA---FDDWVIEEANW-LTLSQDLQKPGEGFD 133 567899***********988776665.78899*****9998757888...4444443...77766654444.445555556666777 PP Pox_MCEL 148 ivsl...QFaiHYsFeseeka.rtvlrNvaellasGGkvlgTtpdgdel 192 v + Fa F+ + + + +l+N+a ++++GG +++ + d + tr|G1KEE2|G1KEE2_ANOCA 134 AVIClgnSFAHLPDFKGDQGVhKQALKNIASMVRPGGIFVIDHRNYDYI 182 6655333799888999988773579**************9977776665 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (333 nodes) Target sequences: 18528 (9796552 residues searched) Passed MSV filter: 897 (0.0484132); expected 370.6 (0.02) Passed bias filter: 397 (0.021427); expected 370.6 (0.02) Passed Vit filter: 39 (0.00210492); expected 18.5 (0.001) Passed Fwd filter: 2 (0.000107945); expected 0.2 (1e-05) Initial search space (Z): 18528 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.24u 0.01s 00:00:00.25 Elapsed: 00:00:00.08 # Mc/sec: 40778.15 // [ok]