# 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]