# 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/Monosiga_brevicollis.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 3.8e-74 247.6 0.0 4.5e-74 247.4 0.0 1.0 1 tr|A9UTL3|A9UTL3_MONBE Predicted protein OS=Monosiga brevico Domain annotation for each sequence (and alignments): >> tr|A9UTL3|A9UTL3_MONBE Predicted protein OS=Monosiga brevicollis OX=81824 GN=15167 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 247.4 0.0 4.9e-78 4.5e-74 16 194 .] 7 189 .. 2 189 .. 0.97 Alignments for each domain: == domain 1 score: 247.4 bits; conditional E-value: 4.9e-78 CGGSB--HHHHCGG.TTSTTGGGSS-CHHHHHHHHHHHCSCHHCTT-EEEEEEETTSEEEEE-TT----SSHHHHHHHHHHHHCCSE CS EMG1 16 kseqLLnsddhkkilkklkreakeaRpDIvHqcLlalldspLnkagllqvyihttknvlIevnpqtRiPrtfkrFagLmvqLLkkls 102 + +qLLns+dh++ilkk++r+ ++RpDI+HqcLl+lldspLnkagllqvyiht+knvlIev+pq+R+Pr ++rF+gLmvqLL++ls tr|A9UTL3|A9UTL3_MONBE 7 NGFQLLNSTDHAHILKKHGRDSLNIRPDITHQCLLMLLDSPLNKAGLLQVYIHTAKNVLIEVHPQCRVPRIYSRFSGLMVQLLHDLS 93 789************************************************************************************ PP EEETTEE.SCEEEECS-HHHHCCC-SEEEEE-TTSEE-..................EEEE-SSSS....---HHEEESSSS---HHH CS EMG1 103 irategeekLlkviknpitdhlppgskiillsesaekklvdleefv.....vkpvvfvigamskgsvkidykeeeisisnypLsaav 184 i a+++++kLlkviknpitdhlp+g+ +i+++ +a+ klv+ ++f+ +kpv++viga+++g+ ++ + +eei++s y +saa+ tr|A9UTL3|A9UTL3_MONBE 94 IAAKGSSQKLLKVIKNPITDHLPTGCLKIGTTYQAK-KLVRARDFAtknamDKPVCVVIGAIAHGQLDVPWITEEIAVSGYAMSAAG 179 ************************************.9889999999988899********************************** PP HHHHHHHHCH CS EMG1 185 vcskiltafE 194 vc k++ fE tr|A9UTL3|A9UTL3_MONBE 180 VCGKLTDGFE 189 ******9998 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (194 nodes) Target sequences: 9188 (5511963 residues searched) Passed MSV filter: 221 (0.0240531); expected 183.8 (0.02) Passed bias filter: 185 (0.020135); expected 183.8 (0.02) Passed Vit filter: 20 (0.00217675); expected 9.2 (0.001) Passed Fwd filter: 1 (0.000108838); expected 0.1 (1e-05) Initial search space (Z): 9188 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.10u 0.01s 00:00:00.11 Elapsed: 00:00:00.05 # Mc/sec: 21386.42 // [ok]