# 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/Dictyostelium_discoideum.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1e-92 308.6 0.7 7.1e-92 305.9 0.7 1.9 1 tr|Q7KWM6|Q7KWM6_DICDI Uncharacterized protein OS=Dictyostel Domain annotation for each sequence (and alignments): >> tr|Q7KWM6|Q7KWM6_DICDI Uncharacterized protein OS=Dictyostelium discoideum OX=44689 GN=DDB0168247 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 305.9 0.7 5.6e-96 7.1e-92 1 194 [] 18 332 .. 18 332 .. 0.98 Alignments for each domain: == domain 1 score: 305.9 bits; conditional E-value: 5.6e-96 EEEEEES--B--GGSCGGSB--HHHHCGG.TTSTTGGGSS-CHHHHHHHHHHHCSCHHCTT-EEEEEEETTSEEEEE-TT----SSH CS EMG1 1 ivLeeasLelvkvkkkseqLLnsddhkkilkklkreakeaRpDIvHqcLlalldspLnkagllqvyihttknvlIevnpqtRiPrtf 87 ivLe+a+L++vkvk +s+qLLn+ddh++ilkk+++ea+eaRpDI+HqcLlal+dspLnkagllqv+i+ttknvlIev+pqtRiPrtf tr|Q7KWM6|Q7KWM6_DICDI 18 IVLEHATLDTVKVK-ESFQLLNCDDHADILKKHGKEASEARPDILHQCLLALFDSPLNKAGLLQVFIKTTKNVLIEVHPQTRIPRTF 103 8************7.************************************************************************ PP HHHHHHHHHHHCCSEEEETTEE.SCEEEECS-HHHHCCC-SEEEEE-TTSEE-.................................. CS EMG1 88 krFagLmvqLLkklsirategeekLlkviknpitdhlppgskiillsesaekklvdleefv..........vkpvv........... 153 +rFagLmvqLLkklsirat+g++kL+kviknpitdhl+pg+ki+++s+sa+ k+vdl+efv v+pv+ tr|Q7KWM6|Q7KWM6_DICDI 104 NRFAGLMVQLLKKLSIRATNGPDKLFKVIKNPITDHLAPGTKIYATSFSAP-KCVDLFEFVpeifgydplpVIPVSttptgvalgss 189 ***************************************************.*********999998877777777999******** PP ....................................................................................... CS EMG1 154 ....................................................................................... 153 tr|Q7KWM6|Q7KWM6_DICDI 190 tnnnnnvnsiqslnslggaqvrsagesyyegmedeeaaeeeekatetnnknknkrkqdeksdndkkkqktsnttttttstattvqkk 276 *************************************************************************************** PP ...............EEEE-SSSS....---HHEEESSSS---HHHHHHHHHHHCH CS EMG1 154 ...............fvigamskgsvkidykeeeisisnypLsaavvcskiltafE 194 +vigams+g++ki+yk+eei++snypLsaa++c+ki+t+fE tr|Q7KWM6|Q7KWM6_DICDI 277 pekksdqlnglqqvcIVIGAMSTGAMKIEYKHEEIAFSNYPLSAAGACFKITTVFE 332 ******************************************************98 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (194 nodes) Target sequences: 12739 (6856946 residues searched) Passed MSV filter: 827 (0.0649188); expected 254.8 (0.02) Passed bias filter: 226 (0.0177408); expected 254.8 (0.02) Passed Vit filter: 17 (0.00133448); expected 12.7 (0.001) Passed Fwd filter: 1 (7.84991e-05); expected 0.1 (1e-05) Initial search space (Z): 12739 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.10u 0.00s 00:00:00.10 Elapsed: 00:00:00.05 # Mc/sec: 26604.95 // [ok]