# 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/Methyltransf_10.hmm.txt
# target sequence database:        proteomes/Trichomonas_vaginalis.fasta
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

Query:       Methyltransf_10  [M=299]
Accession:   PF05971.11
Description: Protein of unknown function (DUF890)
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
    ------- ------ -----    ------- ------ -----   ---- --  --------               -----------
  ------ inclusion threshold ------
       0.37   11.5   1.2        2.3    8.9   0.1    2.0  2  tr|A2F938|A2F938_TRIVA  Uncharacterized protein OS=Trichomona


Domain annotation for each sequence (and alignments):
>> tr|A2F938|A2F938_TRIVA  Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_429190 PE=4 SV=1
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 ?    1.7   0.1    0.0073   3.7e+02     148     189 ..      32      76 ..      18     132 .. 0.79
   2 ?    8.9   0.1   4.6e-05       2.3     143     211 ..     156     230 ..     143     288 .. 0.74

  Alignments for each domain:
  == domain 1  score: 1.7 bits;  conditional E-value: 0.0073
                             HT.T-TTTEEEEE--ST-SSTTTSTT--S..-EEE.EEE-----S CS
         Methyltransf_10 148 snsklsdaiklrkqkekkaifdglieene..rydv.tlCnPPfha 189
                             + sk   +i+++  k+   i+++++++ +   +++ + C PP + 
  tr|A2F938|A2F938_TRIVA  32 QSSKNILTIRCKELKDLIDIYNSIVDKYQelMFKMpSNCFPPSFP 76 
                             4445556799999***********998652667764569999764 PP

  == domain 2  score: 8.9 bits;  conditional E-value: 4.6e-05
                             HHHHHHT.T-TTTEEEEE--ST-SSTTTSTT--S-EE.EEEE-----SS------------...---......... CS
         Methyltransf_10 143 kaivesnsklsdaiklrkqkekkaifdglieeneryd.vtlCnPPfhaslaeakagserklknlgkank......t 211
                             ++  ++n+ +++aik  k  +++++ +gl    + y+ + lC   f +++++ ++ s r+  n+  ++       t
  tr|A2F938|A2F938_TRIVA 156 QKLFKTNTSINTAIKNEKIVNQEMVSKGLDLYRNTYKaFKLCKDFFSSQKKDNNENSAREF-NIEDLKIalslvnT 230
                             678899************************9999997256787666666666666666766.77554223333330 PP



Internal pipeline statistics summary:
-------------------------------------
Query model(s):                              1  (299 nodes)
Target sequences:                        50190  (17055200 residues searched)
Passed MSV filter:                      1005  (0.0200239); expected 1003.8 (0.02)
Passed bias filter:                      637  (0.0126918); expected 1003.8 (0.02)
Passed Vit filter:                        46  (0.000916517); expected 50.2 (0.001)
Passed Fwd filter:                         1  (1.99243e-05); expected 0.5 (1e-05)
Initial search space (Z):              50190  [actual number of targets]
Domain search space  (domZ):               1  [number of targets reported over threshold]
# CPU time: 0.32u 0.01s 00:00:00.33 Elapsed: 00:00:00.12
# Mc/sec: 42495.87
//
[ok]