# 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/Danio_rerio.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 7.4e-92 308.1 0.1 1.1e-91 307.5 0.1 1.3 1 tr|Q58EK9|Q58EK9_DANRE Coiled-coil domain containing 76 OS=D 8.2e-08 32.6 2.7 0.0021 18.1 0.3 2.6 3 tr|B8A5G3|B8A5G3_DANRE Glutathione S-transferase, C-terminal Domain annotation for each sequence (and alignments): >> tr|Q58EK9|Q58EK9_DANRE Coiled-coil domain containing 76 OS=Danio rerio OX=7955 GN=trmt13 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 307.5 0.1 8.8e-96 1.1e-91 1 261 [. 160 452 .. 160 453 .. 0.97 Alignments for each domain: == domain 1 score: 307.5 bits; conditional E-value: 8.8e-96 TRM13 1 lqqsslignleelkllekssayvEfGaGkgelsryvnqaleedkkenskyvLidrasnrlkaDrklrkdeselaikRlkidIkdLnl 87 +qq+s++gn+e l+ll+++++y+EfGaGkg+ls++++ al+++ en++++L++r+s+r+k+D k+++ +s +++Rl++dI++L+l tr|Q58EK9|Q58EK9_DANRE 160 KQQASILGNMEALELLSPNRCYIEFGAGKGKLSHWIHIALKSA--ENVHFLLVERSSTRFKVDGKHKNADS--TFDRLQVDIQHLDL 242 69**************************************777..************************96..************** PP TRM13 88 saleele.kkkkvvavsKHLCGaATDLtLrcllnsekask...................................kakleglliAlC 138 ++++ l+ k+ +v++v+KHLCGaATDL+Lrcl++++ +sk + ++gl+iAlC tr|Q58EK9|Q58EK9_DANRE 243 RKVPLLReKGLPVIGVGKHLCGAATDLALRCLFEHNCTSKddeppfkrvklnqaegaggkdmkpvpstegretgkEIVVSGLAIALC 329 *******99*************************99998899*******************************988899******** PP TRM13 139 ChhvCswkeyvnkeyleelgitkeefqiltkmvSWavsgkrkseaeeedekeekekeekkeeeeseeekksklsseereeiGlkakr 225 Chh+C+w++yv+ke+++e+g+ eef+++++m+SWa++g+rk ++++++k+++eke+++++e+ ++ +s+e+re+iG+ +kr tr|Q58EK9|Q58EK9_DANRE 330 CHHRCDWRHYVGKEFFRERGLGPEEFAAFQRMSSWATCGMRKAAETTSEHKRDEEKEHEMSDEAVLDTLNGFMSVEDREQIGRLCKR 416 *****************************************999888888999999999999999999999**************** PP TRM13 226 lidegRllalkekgfeaelvkYvekevslEnvlLla 261 lid+gR+++l++kgf++++++Y+ ++vslEnvlL a tr|Q58EK9|Q58EK9_DANRE 417 LIDHGRVHYLQKKGFKTSMKYYTSRDVSLENVLLAA 452 *********************************986 PP >> tr|B8A5G3|B8A5G3_DANRE Glutathione S-transferase, C-terminal domain-containing OS=Danio rerio OX=7955 GN=gstcd PE=4 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -4.0 0.3 0.96 1.2e+04 58 70 .. 199 211 .. 187 231 .. 0.50 2 ! 18.1 0.3 1.7e-07 0.0021 17 146 .. 420 527 .. 400 558 .. 0.71 3 ! 14.1 0.0 2.9e-06 0.037 215 260 .. 564 609 .. 554 611 .. 0.93 Alignments for each domain: == domain 1 score: -4.0 bits; conditional E-value: 0.96 TRM13 58 nrlkaDrklrkde 70 ++ k D+++++ e tr|B8A5G3|B8A5G3_DANRE 199 QQKKSDSQVKDAE 211 2222233222222 PP == domain 2 score: 18.1 bits; conditional E-value: 1.7e-07 TRM13 17 ekssayvEfGaGkgelsryvnqaleedkkenskyvLidrasnrlkaDrklrkdeselaikRlkidIkdLnlsaleelekkkkvvavs 103 ++ +++v+f++G g++ +++ +l + +++L++ k + +r +e + +++ + i + nl+ ++ + ++v+ tr|B8A5G3|B8A5G3_DANRE 420 QPGHTVVDFCSGGGHVGIVLAYML-----PKCQVILVEN-----KEESLIRARERSAQLTLTNIGFIQTNLDYFTG----NFNIGVA 492 556678888888888888888877.....4568888875.....55677888877777777777777778887776....5689*** PP TRM13 104 KHLCGaATDLtLrcllnsekaskkakleglliAlCChhvCswk 146 H CG+ATD+ L+ +l+ ++ g++i CC + tr|B8A5G3|B8A5G3_DANRE 493 LHACGVATDMVLDRCLQARA--------GFVISPCCYGFIQNT 527 ************99987765........899999997665544 PP == domain 3 score: 14.1 bits; conditional E-value: 2.9e-06 TRM13 215 ereeiGlkakrlidegRllalkekgfeaelvkYvekevslEnvlLl 260 er+ iG+++ l+d R +a + +g+++++++ ++k s+ n +L+ tr|B8A5G3|B8A5G3_DANRE 564 ERRSIGKRCMGLVDLDRSWAAEAHGYSVRVMTMIPKGCSPKNNMLV 609 899****************************************997 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (262 nodes) Target sequences: 25687 (13798380 residues searched) Passed MSV filter: 2204 (0.0858022); expected 513.7 (0.02) Passed bias filter: 791 (0.0307938); expected 513.7 (0.02) Passed Vit filter: 62 (0.00241367); expected 25.7 (0.001) Passed Fwd filter: 3 (0.000116791); expected 0.3 (1e-05) Initial search space (Z): 25687 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.36u 0.01s 00:00:00.37 Elapsed: 00:00:00.12 # Mc/sec: 30126.46 // [ok]