# 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/LCM.hmm.txt # target sequence database: proteomes/Drosophila_melanogaster.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: LCM [M=188] Accession: PF04072.14 Description: Leucine carboxyl 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.2e-33 115.4 0.0 5.6e-33 114.1 0.0 1.6 1 tr|M9NF32|M9NF32_DROME Leucine carboxyl methyltransferase 1 Domain annotation for each sequence (and alignments): >> tr|M9NF32|M9NF32_DROME Leucine carboxyl methyltransferase 1 OS=Drosophila melanogaster OX=7227 GN=Dmel\CG3793 PE=1 S # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 114.1 0.0 4.1e-37 5.6e-33 1 188 [] 22 211 .. 22 211 .. 0.95 Alignments for each domain: == domain 1 score: 114.1 bits; conditional E-value: 4.1e-37 HHHHHHHHHHHHHH.H.HHTSSHTCCCCCCCHHHHHHHHCHCCCHCHHCHGCHCCCCT--HHHHHHHHHHHHHHHHHHHHHH.HCTT CS LCM 1 algvaaaralesrr.a.adglidDplaeslvsaaglellerkadkelekekeeekkekkapelnrgiavRtrffddllkaal.saag 84 a+++++ +a ++r a + g+++D+++ ++v+++++ kape+nrg+++R++ ++ ++++l +++g tr|M9NF32|M9NF32_DROME 22 AVIATNDDASDCKRcAvRLGYWKDDYIGYFVRNQER----------------------KAPEINRGYFARVKGVEMCVEKFLkKTSG 86 6889999******9769999************9955......................48*************************** PP ESEEEEET-TT--HHHHHHHTT..TEEEEEEE-HHHHHHHHHHHHHSHHHHH..................EECSEEE--TTS.HHHH CS LCM 85 irqvViLGAGlDtrayRlpwla..glkvfEvDqpevlefKretlkelgaepk..................aelryvpvDlrd.ddwp 150 +q+++LG+G+Dt+ +Rl ++a +++E+D+p v ++K t+k+ +a+ + ++++vDlr+ d++ tr|M9NF32|M9NF32_DROME 87 NCQIINLGCGFDTLYFRLRDTAhqVKNFIELDFPTVTARKCYTIKRNKALLArihdedgevrlsptdlhgPSYHLMGVDLRNlDEVD 173 **********************9667799**********************9*************9987667888************ PP HHHHHHTT-TTSEEEEEEESSGGGS-HHHHHHHHHHHH CS LCM 151 eaLaaagfdpekptawvaeGvlvYLtpeasdalleaia 188 ++L++a d + pt+++ae+vlvY++++++++ll++ia tr|M9NF32|M9NF32_DROME 174 SKLQQAEVDYSLPTIFLAECVLVYIEAQNCRNLLKWIA 211 ************************************96 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (188 nodes) Target sequences: 13786 (7394964 residues searched) Passed MSV filter: 537 (0.0389526); expected 275.7 (0.02) Passed bias filter: 317 (0.0229943); expected 275.7 (0.02) Passed Vit filter: 16 (0.0011606); 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.01s 00:00:00.13 Elapsed: 00:00:00.05 # Mc/sec: 27805.06 // [ok]