# 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/TRM13.hmm.txt # target sequence database: proteomes/Petromyzon_marinus.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: TRM13 [M=262] Accession: PF05206.13 Description: Methyltransferase TRM13 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.2e-73 246.6 0.0 3.8e-73 245.8 0.0 1.4 1 tr|S4R6T0|S4R6T0_PETMA tRNA methyltransferase 13 homolog (S. ------ inclusion threshold ------ 0.096 11.8 0.1 0.12 11.5 0.1 1.2 1 tr|S4RUU2|S4RUU2_PETMA Uncharacterized protein OS=Petromyzon Domain annotation for each sequence (and alignments): >> tr|S4R6T0|S4R6T0_PETMA tRNA methyltransferase 13 homolog (S. cerevisiae) 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 ! 245.8 0.0 5.7e-77 3.8e-73 1 261 [. 172 475 .] 172 475 .] 0.93 Alignments for each domain: == domain 1 score: 245.8 bits; conditional E-value: 5.7e-77 TRM13 1 lqqsslignleelkllekssayvEfGaGkgelsryvnqaleedkkenskyvLidrasnrlkaDrklrkdeselaikRlkidIkdLnl 87 +qq+s++g+l+ +ll ++s+ E GaG+g+ls+++ +al + + ++++L++r s+r+k+D k+++ + + ++Rl idI++L+l tr|S4R6T0|S4R6T0_PETMA 172 QQQASIVGHLSTYGLLAHKSVTHECGAGRGKLSHWLWRALPGA-TRGVSFLLVERCSTRFKVDGKHQEGD--VDFRRLLIDIQHLHL 255 59**************************************777.8899********************99..7************** PP TRM13 88 saleele.kkkkvvavsKHLCGaATDLtLrcllnsekask........................kakleglliAlCChhvCswkeyv 149 ++++l+ k+++vva++KHLCGaATDL+Lrcl+++ + + g++i lCChh+Csw yv tr|S4R6T0|S4R6T0_PETMA 256 GKVPQLQdKGTSVVAIGKHLCGAATDLALRCLMETLPRHGtraaapglaacgeeapgqapcetgTGPTLGIAITLCCHHRCSWPAYV 342 *******999************************99876689********************9988999****************** PP TRM13 150 nkeyleelgitkeefqiltkmvSWavsgkrkse.....................aeeedekeekekeekkeeeeseeekksklssee 215 +key+ ++g+ +ef+++++m+SWa++g r+ + e +++ee e ++ + ++e++ + s++e tr|S4R6T0|S4R6T0_PETMA 343 GKEYFCSQGLGPKEFDLIQRMSSWATCGWRRPTdqvpdngnedatreerasrqgGVEWHAEEEMEAIDHAAPCQAEDQISYGRSRAE 429 ****************************9977799**********99998776623333455566666677888899999999**** PP TRM13 216 reeiGlkakrlidegRllalkekgfeaelvkYvekevslEnvlLla 261 re++G+ +kr++degR+ +l+e g+ a+l+ Yve+ ++ EnvlL+a tr|S4R6T0|S4R6T0_PETMA 430 REQLGRLCKRVLDEGRVRYLREWGLTARLQRYVEAGTTVENVLLTA 475 *******************************************986 PP >> tr|S4RUU2|S4RUU2_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 ? 11.5 0.1 1.8e-05 0.12 126 172 .. 164 224 .. 129 277 .. 0.80 Alignments for each domain: == domain 1 score: 11.5 bits; conditional E-value: 1.8e-05 TRM13 126 k.........kakleglliAlCChhvCswkeyvn.....keyleelgitkeefqiltkmvS 172 + ++++gl +A CC h ++++ + e++e++g+++ + +lt S tr|S4RUU2|S4RUU2_PETMA 164 QtvggcfmslMNQMHGLYLAYCCNHPWAVRTLQEhsevlGEFMESKGVSSPGLLMLTTSLS 224 24555555557899*****************999998999999****99998888876655 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (262 nodes) Target sequences: 13539 (5048892 residues searched) Passed MSV filter: 661 (0.0488219); expected 270.8 (0.02) Passed bias filter: 319 (0.0235616); expected 270.8 (0.02) Passed Vit filter: 36 (0.00265899); expected 13.5 (0.001) Passed Fwd filter: 2 (0.000147721); expected 0.1 (1e-05) Initial search space (Z): 13539 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.12u 0.00s 00:00:00.12 Elapsed: 00:00:00.04 # Mc/sec: 33070.24 // [ok]