# 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/Tetrahymena_thermophila.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 5e-69 232.4 0.2 7.8e-69 231.8 0.2 1.3 1 tr|Q237J4|Q237J4_TETTS EMG1/NEP1 methyltransferase OS=Tetrah Domain annotation for each sequence (and alignments): >> tr|Q237J4|Q237J4_TETTS EMG1/NEP1 methyltransferase OS=Tetrahymena thermophila (strain SB210) OX=312017 GN=TTHERM_003 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 231.8 0.2 2.9e-73 7.8e-69 1 194 [] 105 303 .. 105 303 .. 0.97 Alignments for each domain: == domain 1 score: 231.8 bits; conditional E-value: 2.9e-73 EEEEEES--B--GGSCGGSB--HHHHCGG..TTSTTGGGSS-CHHHHHHHHHHHCSCHHCTT-EEEEEEETTSEEEEE-TT----SS CS EMG1 1 ivLeeasLelvkvkkkseqLLnsddhkki.lkklkreakeaRpDIvHqcLlalldspLnkagllqvyihttknvlIevnpqtRiPrt 86 ++Le+++Le++k + k++ LLnsd+hk + kklk++ +++RpDIvH++Ll+l+dspLnkagllqvyiht++nvlI+++p+t+iPrt tr|Q237J4|Q237J4_TETTS 105 VILEKCPLETAKLG-KDFVLLNSDEHKGYiSKKLKKDFSTYRPDIVHHSLLSLMDSPLNKAGLLQVYIHTENNVLIYISPNTKIPRT 190 89***********9.************986788****************************************************** PP HHHHHHHHHHHHCCSEEEETTEE.SCEEEECS-HHHHCCC-SEEEEE-TTSEE-..................EEEE-SSSS...... CS EMG1 87 fkrFagLmvqLLkklsirategeekLlkviknpitdhlppgskiillsesaekklvdleefv.....vkpvvfvigamskg..svki 166 +krFa L++qLL kl++ra +++e+Llk++knp+t+hlp+++ +i++s++a+ lv+++e++ +kpvv ++g++skg ++++ tr|Q237J4|Q237J4_TETTS 191 YKRFAALFAQLLTKLKVRAVQSSETLLKIVKNPVTQHLPSDAMKIGMSTQAR--LVSFKEYIeklpkNKPVVYIVGGVSKGnpAMEV 275 ****************************************************..889*****9999999************888999 PP ---HHEEESSSS---HHHHHHHHHHHCH CS EMG1 167 dykeeeisisnypLsaavvcskiltafE 194 dy +e+i+is+y Lsa +++s++++ +E tr|Q237J4|Q237J4_TETTS 276 DYADEHICISKYGLSAGYCISRMMNCYE 303 ***********************98776 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (194 nodes) Target sequences: 26971 (16912921 residues searched) Passed MSV filter: 2592 (0.0961032); expected 539.4 (0.02) Passed bias filter: 369 (0.0136814); expected 539.4 (0.02) Passed Vit filter: 31 (0.00114938); expected 27.0 (0.001) Passed Fwd filter: 1 (3.70769e-05); expected 0.3 (1e-05) Initial search space (Z): 26971 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.25u 0.01s 00:00:00.26 Elapsed: 00:00:00.11 # Mc/sec: 29828.24 // [ok]