# 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/Danio_rerio.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 5e-42 144.5 0.0 1.2e-41 143.3 0.0 1.6 1 tr|A0A2R8Q906|A0A2R8Q906_DANRE Leucine carboxyl methyltransfera 1.4e-27 97.4 0.0 2.3e-27 96.7 0.0 1.3 1 tr|A8WHX0|A8WHX0_DANRE Leucine carboxyl methyltransfera Domain annotation for each sequence (and alignments): >> tr|A0A2R8Q906|A0A2R8Q906_DANRE Leucine carboxyl methyltransferase 1 OS=Danio rerio OX=7955 GN=lcmt1 PE=1 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 143.3 0.0 9e-46 1.2e-41 1 188 [] 16 208 .. 16 208 .. 0.96 Alignments for each domain: == domain 1 score: 143.3 bits; conditional E-value: 9e-46 HHHHHHHHHHHHHH.H.HHTSSHTCCCCCCCHHHHHHHHCHCCCHCHHCHGCHCCCCT--HHHHHHHHHHHHHHHHHHH CS LCM 1 algvaaaralesrr.a.adglidDplaeslvsaaglellerkadkelekekeeekkekkapelnrgiavRtrffddllk 77 a+++++ +a ++r a ++g+++D +++++v++ g + kape+nrg++ R++ ++ ll+ tr|A0A2R8Q906|A0A2R8Q906_DANRE 16 AVRATCDDASICKRfAtSKGYWTDLYVQYFVRQIGER---------------------KAPEINRGYYGRVKGMNLLLD 73 7899**********659*****************665.....................59******************* PP HHH.HCTTESEEEEET-TT--HHHHHHHTT..TEEEEEEE-HHHHHHHHHHHHHSHHHHH..................E CS LCM 78 aal.saagirqvViLGAGlDtrayRlpwla..glkvfEvDqpevlefKretlkelgaepk..................a 135 ++l ++++ +qvV+LGAGlDt+++Rl++++ k+fEvD+p+++++K++++k+++ ++k + tr|A0A2R8Q906|A0A2R8Q906_DANRE 74 DFLkKTQCDCQVVNLGAGLDTTFWRLKDENimPKKFFEVDFPMIVARKIHNIKTKPPLSKpliethstdsllldghslD 152 ***99999********************************************************************996 PP ECSE..EE--TTS.HHHHHHHHHHTT-TTSEEEEEEESSGGGS-HHHHHHHHHHHH CS LCM 136 elry..vpvDlrd.ddwpeaLaaagfdpekptawvaeGvlvYLtpeasdalleaia 188 + ry +++Dlrd +++++L++ +++pe pt++++e+vlvY++pe+s++l+++ia tr|A0A2R8Q906|A0A2R8Q906_DANRE 153 SDRYciIGADLRDiPTLEDKLKKFQINPELPTLFLSECVLVYMSPEQSSRLVHWIA 208 777667************************************************96 PP >> tr|A8WHX0|A8WHX0_DANRE Leucine carboxyl methyltransferase 2 OS=Danio rerio OX=7955 GN=lcmt2 PE=1 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 96.7 0.0 1.8e-31 2.3e-27 15 187 .. 9 207 .. 3 208 .. 0.86 Alignments for each domain: == domain 1 score: 96.7 bits; conditional E-value: 1.8e-31 H.HHTSSHTCCCCCCCHHHHHHHHCHCCCHCHHCHGCHCCCCT--HHHHHHHHHHHHHHHHHHHHHH.HCTT..ESEEEEET-TT-- CS LCM 15 a.adglidDplaeslvsaaglellerkadkelekekeeekkekkapelnrgiavRtrffddllkaal.saag..irqvViLGAGlDt 97 + +d ++++++ +s vs+++ ++ + + ++++ ++++ k++++ap +nrg++vR +++d+ +k++l +++ rq+ LGAG+D+ tr|A8WHX0|A8WHX0_DANRE 9 QgKDAAVQGTNDSSVVSKVSAAA-QGYFHDDF-LKHFVCKMSRRAPLINRGYYVRWKAVDHCVKQFLhATKScsRRQILSLGAGFDS 93 35677788888888888887664.33333333.3336789999************************544447889*********** PP HHHHHHHTT...TEEEEEEE-HHHHHHHHHHHHHSHHHHH....................EECSEEE--TTS.HHHHHHHHHHTT-T CS LCM 98 rayRlpwla...glkvfEvDqpevlefKretlkelgaepk....................aelryvpvDlrd.ddwpeaLaaagfdp 160 +++Rl + g++vfEvD+pev ++K+++++ ++ + +++++D r+ +++++ L aag++ tr|A8WHX0|A8WHX0_DANRE 94 LFFRLRAEGalgGVTVFEVDFPEVARRKAALINSNACLKDllpdwetvlnkqtnavfissGHYNLIGADVRKeQEVEASLSAAGLQW 180 *******999**********************99888888**************99996556666********************** PP TSEEEEEEESSGGGS-HHHHHHHHHHH CS LCM 161 ekptawvaeGvlvYLtpeasdalleai 187 + pt++++e+vl Y++ + sda++ + tr|A8WHX0|A8WHX0_DANRE 181 DNPTLILSEVVLTYMETQWSDAVIGWA 207 **********************99875 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (188 nodes) Target sequences: 25687 (13798380 residues searched) Passed MSV filter: 723 (0.0281465); expected 513.7 (0.02) Passed bias filter: 518 (0.0201658); expected 513.7 (0.02) Passed Vit filter: 37 (0.00144042); expected 25.7 (0.001) Passed Fwd filter: 2 (7.78604e-05); expected 0.3 (1e-05) Initial search space (Z): 25687 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.23u 0.03s 00:00:00.26 Elapsed: 00:00:00.28 # Mc/sec: 8945.16 // [ok]