# 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/Caenorhabditis_elegans.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1.6e-72 243.4 0.0 1.8e-72 243.2 0.0 1.0 1 sp|Q9XX15|NEP1_CAEEL Ribosomal RNA small subunit methyltrans Domain annotation for each sequence (and alignments): >> sp|Q9XX15|NEP1_CAEEL Ribosomal RNA small subunit methyltransferase nep-1 OS=Caenorhabditis elegans OX=6239 GN=Y39A1A # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 243.2 0.0 9.2e-77 1.8e-72 1 194 [] 30 225 .. 30 225 .. 0.97 Alignments for each domain: == domain 1 score: 243.2 bits; conditional E-value: 9.2e-77 EEEEEES--B--GGSCGGSB--HHHHCGG.TTSTTGGGSS-CHHHHHHHHHHHCSCHHCTT-EEEEEEETTSEEEEE-TT----SSHHH CS EMG1 1 ivLeeasLelvkvkkkseqLLnsddhkkilkklkreakeaRpDIvHqcLlalldspLnkagllqvyihttknvlIevnpqtRiPrtfkr 89 +vLe +sLe++kv+ ++ +L+sd+h+++l+k k++++++RpDI+HqcLl+lldspLn+ag+l+v+ +t+knvl+ v+pq+RiPrtf+r sp|Q9XX15|NEP1_CAEEL 30 VVLEGCSLETAKVG-GEYAILSSDKHANFLRKQKKDPADYRPDILHQCLLNLLDSPLNRAGKLRVFFRTSKNVLVDVSPQCRIPRTFDR 117 8************9.************************************************************************** PP HHHHHHHHHCCSEEEETTEE.SCEEEECS-HHHHCCC-SEEEEE-TTSEE-.................EEEE-SSSS....---HHEEE CS EMG1 90 FagLmvqLLkklsirategeekLlkviknpitdhlppgskiillsesaekklvdleefv....vkpvvfvigamskgsvkidykeeeis 174 F+gLmvqLL+klsira+e ++kL++v+knp+++hlp gs+++l+s++ + +l+ +++v ++p+v++ig++++g++ +dy+++e + sp|Q9XX15|NEP1_CAEEL 118 FCGLMVQLLHKLSIRAAETTQKLMSVVKNPVSNHLPVGSRKMLMSFNVP-ELTMANKLVapetDEPLVLIIGGIARGKIVVDYNDSETK 205 *************************************************.555667777899999************************ PP SSSS---HHHHHHHHHHHCH CS EMG1 175 isnypLsaavvcskiltafE 194 isnypLsaa++c+k+++ +E sp|Q9XX15|NEP1_CAEEL 206 ISNYPLSAALTCAKVTSGLE 225 ****************9887 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (194 nodes) Target sequences: 20008 (8229220 residues searched) Passed MSV filter: 681 (0.0340364); expected 400.2 (0.02) Passed bias filter: 459 (0.0229408); expected 400.2 (0.02) Passed Vit filter: 25 (0.0012495); expected 20.0 (0.001) Passed Fwd filter: 1 (4.998e-05); expected 0.2 (1e-05) Initial search space (Z): 20008 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.12u 0.01s 00:00:00.13 Elapsed: 00:00:00.06 # Mc/sec: 26607.81 // [ok]