# 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/Trichoplax_adhaerens.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 8.9e-30 103.4 0.0 1.9e-29 102.4 0.0 1.5 1 tr|B3RN30|B3RN30_TRIAD Leucine carboxyl methyltransferase 1 5.4e-09 35.6 0.0 8.7e-09 35.0 0.0 1.3 1 tr|B3RQR9|B3RQR9_TRIAD Uncharacterized protein OS=Trichoplax Domain annotation for each sequence (and alignments): >> tr|B3RN30|B3RN30_TRIAD Leucine carboxyl methyltransferase 1 OS=Trichoplax adhaerens OX=10228 GN=TRIADDRAFT_20408 PE= # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 102.4 0.0 3.2e-33 1.9e-29 16 187 .. 23 198 .. 6 199 .. 0.86 Alignments for each domain: == domain 1 score: 102.4 bits; conditional E-value: 3.2e-33 HHTSSHTCCCCCCCHHHHHHHHCHCCCHCHHCHGCHCCCCT--HHHHHHHHHHHHHHHHHHHHHH.HCTTESEEEEET-TT--HHHH CS LCM 16 adglidDplaeslvsaaglellerkadkelekekeeekkekkapelnrgiavRtrffddllkaal.saagirqvViLGAGlDtrayR 101 + g+++D ++ + ++ +++k+pe+nrg+++R++++ dll+ +l +a +qv+ LGAG+Dt ++ tr|B3RN30|B3RN30_TRIAD 23 QLGYWKDRFVNLMSRN----------------------RSRKTPEINRGYYARVKAVGDLLRMFLkTTANNCQVISLGAGFDTNFWN 87 4444444444444444......................45679**********************88899***************** PP HHHTT..TEEEEEEE-HHHHHHHHHHHHHSHHHHH......................EECSEEE--TTS.HHHHHHHHHHTT-TTSE CS LCM 102 lpwla..glkvfEvDqpevlefKretlkelgaepk......................aelryvpvDlrd.ddwpeaLaaagfdpekp 163 l++++ l+++EvD+ ev + K + +k+ + + + ++ +++p+Dlr + ++++L + g+d + p tr|B3RN30|B3RN30_TRIAD 88 LQEEDnlPLNYIEVDFAEVTAAKCAYIKARKPLHQalsedseidfifglsgndavrsKNYNLIPCDLRYiNILEQRLNEIGIDYSLP 174 *****999********************99999999****************9887666777*******9***************** PP EEEEEESSGGGS-HHHHHHHHHHH CS LCM 164 tawvaeGvlvYLtpeasdalleai 187 t+++ e+vlvY+ p++s+++++++ tr|B3RN30|B3RN30_TRIAD 175 TICITECVLVYMAPDKSEEIINWM 198 *********************996 PP >> tr|B3RQR9|B3RQR9_TRIAD Uncharacterized protein OS=Trichoplax adhaerens OX=10228 GN=TRIADDRAFT_55092 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 35.0 0.0 1.5e-12 8.7e-09 15 112 .. 7 113 .. 1 123 [. 0.76 Alignments for each domain: == domain 1 score: 35.0 bits; conditional E-value: 1.5e-12 H.HHTSSHTCCCCCCCHHHHHHHHCHCCCHCHHCHGCHCCCCT--HHHHHHHHHHHHHHHHHHHHHH.HCTT.......ESEEEEET CS LCM 15 a.adglidDplaeslvsaaglellerkadkelekekeeekkekkapelnrgiavRtrffddllkaal.saag.......irqvViLG 92 + +d ++++++ s+vs+ + + ++d+ l++ ++k ++++p +nrg+++R++++d +++++l + +q+ LG tr|B3RQR9|B3RQR9_TRIAD 7 SrSDTAVQGTNDNSIVSKLSSAKHGYYQDDFLQYV--VQKPARRSPLINRGYYIRAKAIDFVIRKFLsY-SNndeaankNKQILSLG 90 34567777777777777776654333333333333..6777899**********************953.33367888899****** PP -TT--HHHHHHHTT...TEEEEE CS LCM 93 AGlDtrayRlpwla...glkvfE 112 AG+Dt+ yR+++l+ ++ v+E tr|B3RQR9|B3RQR9_TRIAD 91 AGFDTLYYRYHSLNmldSTIVYE 113 *************9887777777 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (188 nodes) Target sequences: 11518 (5222075 residues searched) Passed MSV filter: 327 (0.0283903); expected 230.4 (0.02) Passed bias filter: 229 (0.0198819); expected 230.4 (0.02) Passed Vit filter: 19 (0.00164959); expected 11.5 (0.001) Passed Fwd filter: 2 (0.000173641); expected 0.1 (1e-05) Initial search space (Z): 11518 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.09u 0.00s 00:00:00.09 Elapsed: 00:00:00.04 # Mc/sec: 24543.75 // [ok]