# 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/Monosiga_brevicollis.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
    ------- ------ -----    ------- ------ -----   ---- --  --------               -----------
    7.7e-60  201.8   0.1    2.6e-31  108.1   0.2    2.9  2  tr|A9V4W7|A9V4W7_MONBE  Predicted protein OS=Monosiga brevico
  ------ inclusion threshold ------
      0.018   13.5   0.0      0.022   13.2   0.0    1.1  1  tr|A9UWC2|A9UWC2_MONBE  Predicted protein OS=Monosiga brevico


Domain annotation for each sequence (and alignments):
>> tr|A9V4W7|A9V4W7_MONBE  Predicted protein OS=Monosiga brevicollis OX=81824 GN=27301 PE=4 SV=1
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 !  108.1   0.2   5.6e-35   2.6e-31       7     163 ..       1     172 [.       1     179 [. 0.85
   2 !   89.8   0.0   2.1e-29   9.8e-26     176     298 ..     212     324 ..     193     325 .. 0.89

  Alignments for each domain:
  == domain 1  score: 108.1 bits;  conditional E-value: 5.6e-35
                             XXXXXXXXXXX.XXXXXXXXXXXXXXXXX.......XXXXX.........XXXS-TTSHHHHHHHHHHHHHHHH----.---TTS-- CS
         Methyltransf_10   7 lhprnrhkdky.dfaalikvvPelkqfvi.......lnpkG.........rvsidfadPlavkalnkallrefygvsewdiPdgfli 76 
                             +hprnr+  ++ +fa+l+++  +l+ +v        ++ +G         r  idf+d  a++al+ a l+ ++++ + diPdg+li
  tr|A9V4W7|A9V4W7_MONBE   1 MHPRNRYARRPpNFAQLADKCAKLRPYVSypddseaKEAEGharsqhrprRAFIDFRDTNAIRALTEATLECDFNLLT-DIPDGHLI 86 
                             8*********99***************97222222234455333222222456************************9.******** PP

                             --HHHHHHHHHHHHHHHT--TCGCS---EEEEES-TTTTHHHHHHHHHH.--EEEEEES-HHHHHHHHHHHHHT.T-TTTEEEEE-- CS
         Methyltransf_10  77 PpvPgradyihyvadllgkadsdisaqvraldiGvGanCiyPllGvtey.GWrfvgsevddvslasakaivesnsklsdaiklrkqk 162
                             P+vP +++y+h+++dll       + +vr++diG+G   iy +lG+ ++  W+fvg+e+d+++  +a+  ++ n +l++ ++++  +
  tr|A9V4W7|A9V4W7_MONBE  87 PTVPQKLNYLHWLEDLLRTLVA-PPDKVRVVDIGTGPIAIYCVLGARLHpEWQFVGTELDPKACRHAQETIARN-QLTESVHIVSLP 171
                             ****************987654.5689*******************98747***********************.999999998877 PP

                             S CS
         Methyltransf_10 163 e 163
                             +
  tr|A9V4W7|A9V4W7_MONBE 172 H 172
                             6 PP

  == domain 2  score: 89.8 bits;  conditional E-value: 2.1e-29
                             S.-EEEEEE-----SS------------...---.....---------TTTTHHHHTHHHHHHHHHHHHHHHGGGEEEEEEEESSGG CS
         Methyltransf_10 176 e.rydvtlCnPPfhaslaeakagserklknlgkanktralpslnfgGqkaelyceGGevafikklieeskqlakqvlWfttlvskke 261
                             e ++ +++CnPPf+++ +        + + ++        p++  +G  +e++++GGevaf+++++ees  l+++v W+t +++kk 
  tr|A9V4W7|A9V4W7_MONBE 212 ElEFAISMCNPPFYDTHESP-----ANYTGIRP------APKVEATGSSSEMHVQGGEVAFVTRMMEESARLRTRVGWYTVMLGKKR 287
                             44799*********987654.....34444533......588999****************************************** PP

                             GHHHHHHHHHHTT-SEEEEEEEEETTEEEEEEEEES- CS
         Methyltransf_10 262 nlaplkeelkklGaskvtvvelaqGqklsrfiaWsfl 298
                             +l +l+ +l+ lG ++v+   ++qG++ ++++aW+fl
  tr|A9V4W7|A9V4W7_MONBE 288 SLLDLRPQLTALGCTHVVESSFVQGKTHRWAVAWTFL 324
                             ***********************************96 PP

>> tr|A9UWC2|A9UWC2_MONBE  Predicted protein OS=Monosiga brevicollis OX=81824 GN=36597 PE=4 SV=1
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 ?   13.2   0.0   4.7e-06     0.022      88     149 ..      36      93 ..      26     146 .. 0.77

  Alignments for each domain:
  == domain 1  score: 13.2 bits;  conditional E-value: 4.7e-06
                             HHHHHHT--TCGCS---EEEEES-TTTTHHHHHHHHHH--EEEEEES-HHHHHHHHHHHHHT CS
         Methyltransf_10  88 yvadllgkadsdisaqvraldiGvGanCiyPllGvteyGWrfvgsevddvslasakaivesn 149
                             ++a l  +  ++++ ++  l  G+Ga+C y      e G++ vg +v + +l++a++  e+ 
  tr|A9UWC2|A9UWC2_MONBE  36 RLATLCHSMRENKAPRALELGCGTGASCCY----LAEQGYNVVGIDVVQAALDQAEKQREQL 93 
                             56666666666676677778899*******....589******************9976664 PP



Internal pipeline statistics summary:
-------------------------------------
Query model(s):                              1  (299 nodes)
Target sequences:                         9188  (5511963 residues searched)
Passed MSV filter:                       219  (0.0238354); expected 183.8 (0.02)
Passed bias filter:                      173  (0.0188289); expected 183.8 (0.02)
Passed Vit filter:                        11  (0.00119721); expected 9.2 (0.001)
Passed Fwd filter:                         2  (0.000217675); expected 0.1 (1e-05)
Initial search space (Z):               9188  [actual number of targets]
Domain search space  (domZ):               2  [number of targets reported over threshold]
# CPU time: 0.12u 0.01s 00:00:00.13 Elapsed: 00:00:00.05
# Mc/sec: 32961.54
//
[ok]