# 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/Petromyzon_marinus.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2e-38 131.9 0.0 2.4e-38 131.7 0.0 1.0 1 tr|S4RXH9|S4RXH9_PETMA Uncharacterized protein OS=Petromyzon 0.00042 19.1 0.0 0.00054 18.8 0.0 1.1 1 tr|S4RC50|S4RC50_PETMA Phenylethanolamine N-methyltransferas ------ inclusion threshold ------ 0.039 12.7 0.0 0.19 10.4 0.0 1.8 2 tr|S4R684|S4R684_PETMA EEF1A lysine methyltransferase 2 OS=P 0.2 10.3 0.0 0.25 10.0 0.0 1.1 1 tr|S4RMF5|S4RMF5_PETMA HemK methyltransferase family member Domain annotation for each sequence (and alignments): >> tr|S4RXH9|S4RXH9_PETMA Uncharacterized protein OS=Petromyzon marinus OX=7757 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 131.7 0.0 7e-42 2.4e-38 11 133 .. 3 123 .. 1 127 [. 0.92 Alignments for each domain: == domain 1 score: 131.7 bits; conditional E-value: 7e-42 Pox_MCEL 11 tdlvathynerpetgellkkRiastiiylrnlnNyiKslLiseyasktkld.esnklkvLdldcgkGgDLeKyekgeIsklvatDia 96 ++vathyn+++e g + R++s+i+y+rn+nN++Ks+Li+e+++++++d + +l+vLdl+cgkGgD+ K++ + I +lv++Dia tr|S4RXH9|S4RXH9_PETMA 3 GSVVATHYNQLQEVG--RDARSHSRILYMRNFNNWLKSVLINEFVQRVRADgRGGALSVLDLGCGKGGDMLKWRRASIGHLVCADIA 87 589************..8*******************************9636789******************************* PP Pox_MCEL 97 eesieeckeRYnklrsekkskyykFdaelirsdsfvs 133 ++s+++c++RYn +r + F+ae+i++d+ ++ tr|S4RXH9|S4RXH9_PETMA 88 ATSVDQCRQRYNDMRGRGGG-GPIFSAEFIAADCSKV 123 ***********888766665.478**********875 PP >> tr|S4RC50|S4RC50_PETMA Phenylethanolamine N-methyltransferase OS=Petromyzon marinus OX=7757 GN=PNMT PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 18.8 0.0 1.6e-07 0.00054 140 184 .. 154 198 .. 115 216 .. 0.84 Alignments for each domain: == domain 1 score: 18.8 bits; conditional E-value: 1.6e-07 Pox_MCEL 140 edyelkfdivslQFaiHYsFeseekartvlrNvaellasGGkvlg 184 ++ ++fd+v+ F++ + ++e r +lrNva+ll sGG +l tr|S4RC50|S4RC50_PETMA 154 PEQPAPFDVVQSSFCLECACDTELSFRAALRNVAALLVSGGWFLL 198 346799***********************************9985 PP >> tr|S4R684|S4R684_PETMA EEF1A lysine methyltransferase 2 OS=Petromyzon marinus OX=7757 GN=EEF1AKMT2 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.4 0.0 5.6e-05 0.19 64 104 .. 78 118 .. 67 138 .. 0.89 2 ? -0.3 0.0 0.1 3.4e+02 164 187 .. 174 196 .. 170 203 .. 0.78 Alignments for each domain: == domain 1 score: 10.4 bits; conditional E-value: 5.6e-05 Pox_MCEL 64 nklkvLdldcgkGgDLeKyekgeIsklvatDiaeesieeck 104 ++ +vLd++ g+G L k + +l ++D ++++++ ++ tr|S4R684|S4R684_PETMA 78 KDSSVLDIGTGNGVLLLELAKVGFTRLTGIDYSAAAVRLAQ 118 5679********************************98776 PP == domain 2 score: -0.3 bits; conditional E-value: 0.1 Pox_MCEL 164 kartvlrNvaellasGGkvlgTtp 187 +a ++ + ++la+GG +++T tr|S4R684|S4R684_PETMA 174 RA-DYVEALRRMLAPGGLFVITSC 196 55.445778899*********976 PP >> tr|S4RMF5|S4RMF5_PETMA HemK methyltransferase family member 1 OS=Petromyzon marinus OX=7757 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.0 0.0 7.4e-05 0.25 62 105 .. 56 100 .. 36 119 .. 0.85 Alignments for each domain: == domain 1 score: 10.0 bits; conditional E-value: 7.4e-05 Pox_MCEL 62 esnklkvLdldcgkGg.DLeKyekgeIsklvatDiaeesieecke 105 + +vLd++cg G+ L ++ +++ va+D ++++++ ++e tr|S4RMF5|S4RMF5_PETMA 56 QGGPYRVLDIGCGSGAiSLAILSASPLASCVAVDKSADAVALARE 100 556789*********977999999****************99887 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (333 nodes) Target sequences: 13539 (5048892 residues searched) Passed MSV filter: 392 (0.0289534); expected 270.8 (0.02) Passed bias filter: 252 (0.0186129); expected 270.8 (0.02) Passed Vit filter: 16 (0.00118177); expected 13.5 (0.001) Passed Fwd filter: 4 (0.000295443); expected 0.1 (1e-05) Initial search space (Z): 13539 [actual number of targets] Domain search space (domZ): 4 [number of targets reported over threshold] # CPU time: 0.11u 0.00s 00:00:00.11 Elapsed: 00:00:00.04 # Mc/sec: 42032.03 // [ok]