# 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/Pseudomonas_aeruginosa.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
    ------- ------ -----    ------- ------ -----   ---- --  --------               -----------
   2.6e-155  514.6   0.0     3e-155  514.4   0.0    1.0  1  sp|Q9HXG4|RLMF_PSEAE    Ribosomal RNA large subunit methyltra
  ------ inclusion threshold ------
       0.03   12.0   0.0      0.043   11.5   0.0    1.1  1  tr|Q9I281|Q9I281_PSEAE  Uncharacterized protein OS=Pseudomona


Domain annotation for each sequence (and alignments):
>> sp|Q9HXG4|RLMF_PSEAE  Ribosomal RNA large subunit methyltransferase F OS=Pseudomonas aeruginosa (strain ATCC 15692 / 
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 !  514.4   0.0  1.1e-158    3e-155       2     299 .]      14     314 ..      13     314 .. 0.99

  Alignments for each domain:
  == domain 1  score: 514.4 bits;  conditional E-value: 1.1e-158
                           XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXS-TTSHHHHHHHHHHHHHHHH----.---TTS----HHHHHHHHHHHH CS
       Methyltransf_10   2 alksglhprnrhkdkydfaalikvvPelkqfvilnpkGrvsidfadPlavkalnkallrefygvsewdiPdgfliPpvPgradyihyva 90 
                           +++s+lhprnrh ++ydf++li+ +Pel++fvilnp+Gr+sidfadP+avka+n+all++fy+v+ewdiPdg+l+Pp+Pgrady+hy+a
  sp|Q9HXG4|RLMF_PSEAE  14 KSASPLHPRNRHLGRYDFPRLIAGSPELERFVILNPYGRQSIDFADPAAVKAFNRALLQQFYDVREWDIPDGYLCPPIPGRADYLHYLA 102
                           6899************************************************************************************* PP

                           HHHT--TCG.CS-..--EEEEES-TTTTHHHHHHHHHH--EEEEEES-HHHHHHHHHHHHHT.T-TTTEEEEE--ST-SSTTTSTT--S CS
       Methyltransf_10  91 dllgkadsd.isa..qvraldiGvGanCiyPllGvteyGWrfvgsevddvslasakaivesnsklsdaiklrkqkekkaifdglieene 176
                           dllg+++++ i++   +rald+G+GanCiyPllG++eyGWrfvg+++d++slasa+ai+++n+++++ai+lr+q+++++if+gli+++e
  sp|Q9HXG4|RLMF_PSEAE 103 DLLGASHDGlIPRgpGLRALDVGTGANCIYPLLGHHEYGWRFVGADIDPQSLASAAAILAANPRFAAAIELRRQPDRRQIFQGLIGVDE 191
                           *******998888999************************************************************************* PP

                           -EEEEEE-----SS------------...---.....---------TTTTHHHHTHHHHHHHHHHHHHHHGGGEEEEEEEESSGGGHHH CS
       Methyltransf_10 177 rydvtlCnPPfhaslaeakagserklknlgkanktralpslnfgGqkaelyceGGevafikklieeskqlakqvlWfttlvskkenlap 265
                           r+d+tlCnPPfhasl+ea++gs+rk+knlgk+++tr+lp lnfgGq aelyceGGe+af++ ++ees+++a+qv+Wfttlvsk++nl++
  sp|Q9HXG4|RLMF_PSEAE 192 RFDMTLCNPPFHASLDEATRGSRRKWKNLGKLDPTRTLPLLNFGGQGAELYCEGGEAAFLAGMAEESRAFATQVFWFTTLVSKASNLPN 280
                           ***************************************************************************************** PP

                           HHHHHHHTT-SEEEEEEEEETTEEEEEEEEES-- CS
       Methyltransf_10 266 lkeelkklGaskvtvvelaqGqklsrfiaWsfld 299
                           l+e+lk+lGas+++vv++aqGqk+srf+aW++ld
  sp|Q9HXG4|RLMF_PSEAE 281 LQERLKTLGASDIRVVDMAQGQKQSRFVAWTYLD 314
                           ********************************87 PP

>> tr|Q9I281|Q9I281_PSEAE  Uncharacterized protein OS=Pseudomonas aeruginosa (strain ATCC 15692 / DSM 22644 / CIP 104116
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 ?   11.5   0.0   1.5e-05     0.043      99     148 ..      37      86 ..      22      97 .. 0.80

  Alignments for each domain:
  == domain 1  score: 11.5 bits;  conditional E-value: 1.5e-05
                             GCS---EEEEES-TTTTHHHHHHHHHH--EEEEEES-HHHHHHHHHHHHH CS
         Methyltransf_10  99 disaqvraldiGvGanCiyPllGvteyGWrfvgsevddvslasakaives 148
                              +++++r l +G Ga      l     GWrfv+++++   la a++  + 
  tr|Q9I281|Q9I281_PSEAE  37 RLPERARLLSVGCGAGGELAELAGLRPGWRFVAADLSADMLALARQRCAR 86 
                             478899*********998888888899***********999999887655 PP



Internal pipeline statistics summary:
-------------------------------------
Query model(s):                              1  (299 nodes)
Target sequences:                         5563  (1857172 residues searched)
Passed MSV filter:                       123  (0.0221104); expected 111.3 (0.02)
Passed bias filter:                      115  (0.0206723); expected 111.3 (0.02)
Passed Vit filter:                        15  (0.00269639); expected 5.6 (0.001)
Passed Fwd filter:                         2  (0.000359518); expected 0.1 (1e-05)
Initial search space (Z):               5563  [actual number of targets]
Domain search space  (domZ):               2  [number of targets reported over threshold]
# CPU time: 0.04u 0.00s 00:00:00.04 Elapsed: 00:00:00.02
# Mc/sec: 27764.72
//
[ok]