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

Query:       EMG1  [M=194]
Accession:   PF03587.13
Description: EMG1/NEP1 methyltransferase
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.5e-84  281.4   0.0      3e-84  281.2   0.0    1.0  1  sp|Q9W4J5|NEP1_DROME  Ribosomal RNA small subunit methyltrans


Domain annotation for each sequence (and alignments):
>> sp|Q9W4J5|NEP1_DROME  Ribosomal RNA small subunit methyltransferase NEP1 OS=Drosophila melanogaster OX=7227 GN=CG3527
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 !  281.2   0.0   2.2e-88     3e-84       1     194 []      45     246 ..      45     246 .. 0.98

  Alignments for each domain:
  == domain 1  score: 281.2 bits;  conditional E-value: 2.2e-88
                           EEEEEES--B--GGSCGGSB--HHHHCGG.TTSTTGGGSS-CHHHHHHHHHHHCSCHHCTT-EEEEEEETTSEEEEE-TT----SSHHH CS
                  EMG1   1 ivLeeasLelvkvkkkseqLLnsddhkkilkklkreakeaRpDIvHqcLlalldspLnkagllqvyihttknvlIevnpqtRiPrtfkr 89 
                           ivLe a+Le+vkv+ ++++LLn+ddh+ i++k++r++ ++RpDI+HqcLl+l+dspLn+agllqv+++t++nvlIe+npqtRiPrtfkr
  sp|Q9W4J5|NEP1_DROME  45 IVLEGAQLETVKVH-NTFELLNCDDHAGIMRKNQRDPGSCRPDITHQCLLMLFDSPLNRAGLLQVFVRTEHNVLIEINPQTRIPRTFKR 132
                           8************7.************************************************************************** PP

                           HHHHHHHHHCCSEEEETTEE.SCEEEECS-HHHHCCC-SEEEEE-TTSEE-.......................EEEE-SSSS....-- CS
                  EMG1  90 FagLmvqLLkklsirategeekLlkviknpitdhlppgskiillsesaekklvdleefv..........vkpvvfvigamskgsvkidy 168
                           FagLmvqLL+k +ira+++++ L++viknpitdh+p g+k++++s+s++  l +++++v          ++pvv+viga+++g  k+dy
  sp|Q9W4J5|NEP1_DROME 133 FAGLMVQLLHKFQIRANDSSRRLMSVIKNPITDHVPVGCKKYAMSFSGK-LLPNCRDLVphgdetsasyDEPVVIVIGAFAHGVLKTDY 220
                           *************************************************.666889999999999999999****************** PP

                           -HHEEESSSS---HHHHHHHHHHHCH CS
                  EMG1 169 keeeisisnypLsaavvcskiltafE 194
                           +ee +sisnypLsaa++cski++afE
  sp|Q9W4J5|NEP1_DROME 221 TEELFSISNYPLSAAIACSKICSAFE 246
                           *************************8 PP



Internal pipeline statistics summary:
-------------------------------------
Query model(s):                              1  (194 nodes)
Target sequences:                        13786  (7394964 residues searched)
Passed MSV filter:                       426  (0.0309009); expected 275.7 (0.02)
Passed bias filter:                      291  (0.0211084); expected 275.7 (0.02)
Passed Vit filter:                        21  (0.00152328); expected 13.8 (0.001)
Passed Fwd filter:                         1  (7.25374e-05); expected 0.1 (1e-05)
Initial search space (Z):              13786  [actual number of targets]
Domain search space  (domZ):               1  [number of targets reported over threshold]
# CPU time: 0.12u 0.00s 00:00:00.12 Elapsed: 00:00:00.05
# Mc/sec: 28692.46
//
[ok]