# 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/Homo_sapiens.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1.4e-41 144.8 0.0 3.6e-41 143.4 0.0 1.6 1 sp|Q9UIC8|LCMT1_HUMAN Leucine carboxyl methyltransferase 1 6.4e-40 139.4 0.0 1.6e-39 138.0 0.0 1.6 1 sp|Q9UIC8-2|LCMT1_HUMAN Isoform 2 of Leucine carboxyl methyl 3.3e-30 107.7 0.0 7.1e-23 83.8 0.0 2.7 2 sp|Q9UIC8-3|LCMT1_HUMAN Isoform 3 of Leucine carboxyl methyl 5.6e-27 97.2 0.0 1.6e-26 95.7 0.0 1.7 1 sp|O60294|LCMT2_HUMAN Leucine carboxyl methyltransferase 2 1.7e-24 89.0 0.0 7.1e-24 87.0 0.0 1.7 1 tr|I3L2Q8|I3L2Q8_HUMAN Leucine carboxyl methyltransferase 1 3.8e-19 71.6 0.0 4.7e-19 71.3 0.0 1.1 1 tr|I3L4A2|I3L4A2_HUMAN Leucine carboxyl methyltransferase 1 0.00017 23.8 0.1 0.00026 23.3 0.1 1.3 1 tr|I3L2E3|I3L2E3_HUMAN Leucine carboxyl methyltransferase 1 ------ inclusion threshold ------ 0.18 14.0 0.0 0.19 13.9 0.0 1.1 1 tr|H3BUP7|H3BUP7_HUMAN Leucine carboxyl methyltransferase 1 Domain annotation for each sequence (and alignments): >> sp|Q9UIC8|LCMT1_HUMAN Leucine carboxyl methyltransferase 1 OS=Homo sapiens GN=LCMT1 PE=1 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 143.4 0.0 3.2e-45 3.6e-41 3 187 .. 27 216 .. 25 217 .. 0.95 Alignments for each domain: == domain 1 score: 143.4 bits; conditional E-value: 3.2e-45 HHHHHHHHHHHH.H.HHTSSHTCCCCCCCHHHHHHHHCHCCCHCHHCHGCHCCCCT--HHHHHHHHHHHHHHHHHHHHHH.HCTTESE CS LCM 3 gvaaaralesrr.a.adglidDplaeslvsaaglellerkadkelekekeeekkekkapelnrgiavRtrffddllkaal.saagirq 87 + ++ +a ++r a + g+++Dp+++++v+ ++ e+kape+nrg+++R++ + +l+ka+l ++++++q sp|Q9UIC8|LCMT1_HUMAN 27 RGTCEDASLCKRfAvSIGYWHDPYIQHFVRLSK---------------------ERKAPEINRGYFARVHGVSQLIKAFLrKTECHCQ 93 678999******669***************993.....................3469****************************** PP EEEET-TT--HHHHHHHTT..TEEEEEEE-HHHHHHHHHHHHHSHHHHH..................EECSE..EE--TTS.HHHHHH CS LCM 88 vViLGAGlDtrayRlpwla..glkvfEvDqpevlefKretlkelgaepk..................aelry..vpvDlrd.ddwpea 152 +V+LGAG+Dt+++Rl++++ +k+fEvD+p+++++K + +k ++ +++ ++ ry +++Dlrd ++++e+ sp|Q9UIC8|LCMT1_HUMAN 94 IVNLGAGMDTTFWRLKDEDllPSKYFEVDFPMIVTRKLHSIKCKPPLSSpilelhsedtlqmdghilDSKRYavIGADLRDlSELEEK 181 *********************************************99999****************96667667************** PP HHHHTT-TTSEEEEEEESSGGGS-HHHHHHHHHHH CS LCM 153 LaaagfdpekptawvaeGvlvYLtpeasdalleai 187 L++++++++ pt+++ae+vlvY+tpe+s++ll++ sp|Q9UIC8|LCMT1_HUMAN 182 LKKCNMNTQLPTLLIAECVLVYMTPEQSANLLKWA 216 ********************************996 PP >> sp|Q9UIC8-2|LCMT1_HUMAN Isoform 2 of Leucine carboxyl methyltransferase 1 OS=Homo sapiens GN=LCMT1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 138.0 0.0 1.5e-43 1.6e-39 16 187 .. 65 239 .. 57 240 .. 0.93 Alignments for each domain: == domain 1 score: 138.0 bits; conditional E-value: 1.5e-43 HHTSSHTCCCCCCCHHHHHHHHCHCCCHCHHCHGCHCCCCT--HHHHHHHHHHHHHHHHHHHHHH.HCTTESEEEEET-TT--HHH CS LCM 16 adglidDplaeslvsaaglellerkadkelekekeeekkekkapelnrgiavRtrffddllkaal.saagirqvViLGAGlDtray 100 + g+++Dp+++++v+ + ke+kape+nrg+++R++ + +l+ka+l ++++++q+V+LGAG+Dt+++ sp|Q9UIC8-2|LCMT1_HUMAN 65 SIGYWHDPYIQHFVRLS---------------------KERKAPEINRGYFARVHGVSQLIKAFLrKTECHCQIVNLGAGMDTTFW 129 56889999999999988.....................34579******************************************* PP HHHHTT..TEEEEEEE-HHHHHHHHHHHHHSHHHHH..................EECSE..EE--TTS.HHHHHHHHHHTT-TTSE CS LCM 101 Rlpwla..glkvfEvDqpevlefKretlkelgaepk..................aelry..vpvDlrd.ddwpeaLaaagfdpekp 163 Rl++++ +k+fEvD+p+++++K + +k ++ +++ ++ ry +++Dlrd ++++e+L++++++++ p sp|Q9UIC8-2|LCMT1_HUMAN 130 RLKDEDllPSKYFEVDFPMIVTRKLHSIKCKPPLSSpilelhsedtlqmdghilDSKRYavIGADLRDlSELEEKLKKCNMNTQLP 215 ********************************99999****************96667667************************* PP EEEEEESSGGGS-HHHHHHHHHHH CS LCM 164 tawvaeGvlvYLtpeasdalleai 187 t+++ae+vlvY+tpe+s++ll++ sp|Q9UIC8-2|LCMT1_HUMAN 216 TLLIAECVLVYMTPEQSANLLKWA 239 *********************996 PP >> sp|Q9UIC8-3|LCMT1_HUMAN Isoform 3 of Leucine carboxyl methyltransferase 1 OS=Homo sapiens GN=LCMT1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 83.8 0.0 6.4e-27 7.1e-23 3 127 .. 27 135 .. 25 136 .. 0.93 2 ! 21.6 0.1 7.4e-08 0.00082 161 187 .. 135 161 .. 132 162 .. 0.92 Alignments for each domain: == domain 1 score: 83.8 bits; conditional E-value: 6.4e-27 HHHHHHHHHHHH.H.HHTSSHTCCCCCCCHHHHHHHHCHCCCHCHHCHGCHCCCCT--HHHHHHHHHHHHHHHHHHHHHH.HCTTE CS LCM 3 gvaaaralesrr.a.adglidDplaeslvsaaglellerkadkelekekeeekkekkapelnrgiavRtrffddllkaal.saagi 85 + ++ +a ++r a + g+++Dp+++++v+ ++ e+kape+nrg+++R++ + +l+ka+l +++++ sp|Q9UIC8-3|LCMT1_HUMAN 27 RGTCEDASLCKRfAvSIGYWHDPYIQHFVRLSK---------------------ERKAPEINRGYFARVHGVSQLIKAFLrKTECH 91 678999******669***************993.....................3469**************************** PP SEEEEET-TT--HHHHHHHTT..TEEEEEEE-HHHHHHHHHHHH CS LCM 86 rqvViLGAGlDtrayRlpwla..glkvfEvDqpevlefKretlk 127 +q+V+LGAG+Dt+++Rl++++ +k+fEvD+p+++++K + +k sp|Q9UIC8-3|LCMT1_HUMAN 92 CQIVNLGAGMDTTFWRLKDEDllPSKYFEVDFPMIVTRKLHSIK 135 ***************************************99886 PP == domain 2 score: 21.6 bits; conditional E-value: 7.4e-08 TSEEEEEEESSGGGS-HHHHHHHHHHH CS LCM 161 ekptawvaeGvlvYLtpeasdalleai 187 + pt+++ae+vlvY+tpe+s++ll++ sp|Q9UIC8-3|LCMT1_HUMAN 135 KLPTLLIAECVLVYMTPEQSANLLKWA 161 579*********************996 PP >> sp|O60294|LCMT2_HUMAN Leucine carboxyl methyltransferase 2 OS=Homo sapiens GN=LCMT2 PE=1 SV=3 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 95.7 0.0 1.4e-30 1.6e-26 9 187 .. 20 206 .. 12 207 .. 0.88 Alignments for each domain: == domain 1 score: 95.7 bits; conditional E-value: 1.4e-30 HHHHHHH..HHTSSHTCCCCCCCHHHHHHHHCHCCCHCHHCHGCHCCCCT--HHHHHHHHHHHHHHHHHHHHHHHCTT......ESEE CS LCM 9 alesrra..adglidDplaeslvsaaglellerkadkelekekeeekkekkapelnrgiavRtrffddllkaalsaag......irqv 88 + s+r+ a+g+++Dp+a+ lv a+ ++ap ++rg++vR+r++++ ++a+l++ g +q+ sp|O60294|LCMT2_HUMAN 20 SALSKRSlaARGYVQDPFAALLVPGAA----------------------RRAPLIHRGYYVRARAVRHCVRAFLEQIGapqaalRAQI 85 444444335788888888888888884......................469**********************657778888878** PP EEET-TT--HHHHHHHTT...TEEEEEEE-HHHHHHHHHHHHHSHHHHH..................EECSEEE--TTS.HHHHHHHH CS LCM 89 ViLGAGlDtrayRlpwla...glkvfEvDqpevlefKretlkelgaepk..................aelryvpvDlrd.ddwpeaLa 154 LGAG+D++ +Rl+ + v EvD+p+v ++K+e++ e +++ a+ ++Dlr+ ++++eaL sp|O60294|LCMT2_HUMAN 86 LSLGAGFDSLYFRLKTAGrlaRAAVWEVDFPDVARRKAERIGETPELCAltgpfergepasalcfesADYCILGLDLRQlQRVEEALG 173 ***************999889***********************99999**************996556666**************** PP HHTT-TTSEEEEEEESSGGGS-HHHHHHHHHHH CS LCM 155 aagfdpekptawvaeGvlvYLtpeasdalleai 187 aag+d++ pt+++ae vl YL+pe+++al+++ sp|O60294|LCMT2_HUMAN 174 AAGLDAASPTLLLAEAVLTYLEPESAAALIAWA 206 ******************************996 PP >> tr|I3L2Q8|I3L2Q8_HUMAN Leucine carboxyl methyltransferase 1 (Fragment) OS=Homo sapiens GN=LCMT1 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 87.0 0.0 6.4e-28 7.1e-24 3 127 .. 11 119 .] 9 119 .] 0.93 Alignments for each domain: == domain 1 score: 87.0 bits; conditional E-value: 6.4e-28 HHHHHHHHHHHH.H.HHTSSHTCCCCCCCHHHHHHHHCHCCCHCHHCHGCHCCCCT--HHHHHHHHHHHHHHHHHHHHHH.HCTTES CS LCM 3 gvaaaralesrr.a.adglidDplaeslvsaaglellerkadkelekekeeekkekkapelnrgiavRtrffddllkaal.saagir 86 + ++ +a ++r a + g+++Dp+++++v+ ++ e+kape+nrg+++R++ + +l+ka+l ++++++ tr|I3L2Q8|I3L2Q8_HUMAN 11 RGTCEDASLCKRfAvSIGYWHDPYIQHFVRLSK---------------------ERKAPEINRGYFARVHGVSQLIKAFLrKTECHC 76 679999******669***************993.....................3469***************************** PP EEEEET-TT--HHHHHHHTT..TEEEEEEE-HHHHHHHHHHHH CS LCM 87 qvViLGAGlDtrayRlpwla..glkvfEvDqpevlefKretlk 127 q+V+LGAG+Dt+++Rl++++ +k+fEvD+p+++++K + +k tr|I3L2Q8|I3L2Q8_HUMAN 77 QIVNLGAGMDTTFWRLKDEDllPSKYFEVDFPMIVTRKLHSIK 119 **************************************98875 PP >> tr|I3L4A2|I3L4A2_HUMAN Leucine carboxyl methyltransferase 1 (Fragment) OS=Homo sapiens GN=LCMT1 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 71.3 0.0 4.3e-23 4.7e-19 107 187 .. 5 106 .. 1 107 [. 0.94 Alignments for each domain: == domain 1 score: 71.3 bits; conditional E-value: 4.3e-23 TEEEEEEE-HHHHHHHHHHHHHSHHHHH..................EECSE..EE--TTS.HHHHHHHHHHTT-TTSEEEEEEESSG CS LCM 107 glkvfEvDqpevlefKretlkelgaepk..................aelry..vpvDlrd.ddwpeaLaaagfdpekptawvaeGvl 172 +k+fEvD+p+++++K + +k ++ +++ ++ ry +++Dlrd ++++e+L++++++++ pt+++ae+vl tr|I3L4A2|I3L4A2_HUMAN 5 PSKYFEVDFPMIVTRKLHSIKCKPPLSSpilelhsedtlqmdghilDSKRYavIGADLRDlSELEEKLKKCNMNTQLPTLLIAECVL 91 579*********************99999****************96667667********************************** PP GGS-HHHHHHHHHHH CS LCM 173 vYLtpeasdalleai 187 vY+tpe+s++ll++ tr|I3L4A2|I3L4A2_HUMAN 92 VYMTPEQSANLLKWA 106 ************996 PP >> tr|I3L2E3|I3L2E3_HUMAN Leucine carboxyl methyltransferase 1 (Fragment) OS=Homo sapiens GN=LCMT1 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 23.3 0.1 2.3e-08 0.00026 162 187 .. 2 27 .. 1 28 [. 0.93 Alignments for each domain: == domain 1 score: 23.3 bits; conditional E-value: 2.3e-08 SEEEEEEESSGGGS-HHHHHHHHHHH CS LCM 162 kptawvaeGvlvYLtpeasdalleai 187 pt+++ae+vlvY+tpe+s++ll++ tr|I3L2E3|I3L2E3_HUMAN 2 LPTLLIAECVLVYMTPEQSANLLKWA 27 59**********************96 PP >> tr|H3BUP7|H3BUP7_HUMAN Leucine carboxyl methyltransferase 1 OS=Homo sapiens GN=LCMT1 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.9 0.0 1.7e-05 0.19 3 37 .. 27 63 .. 25 92 .. 0.82 Alignments for each domain: == domain 1 score: 13.9 bits; conditional E-value: 1.7e-05 HHHHHHHHHHHH.H.HHTSSHTCCCCCCCHHHHHHHH CS LCM 3 gvaaaralesrr.a.adglidDplaeslvsaaglell 37 + ++ +a ++r a + g+++Dp+++++v+ ++ + + tr|H3BUP7|H3BUP7_HUMAN 27 RGTCEDASLCKRfAvSIGYWHDPYIQHFVRLSKERKA 63 67999*******669***************9876654 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (188 nodes) Target sequences: 88473 (35070517 residues searched) Passed MSV filter: 2488 (0.0281216); expected 1769.5 (0.02) Passed bias filter: 1843 (0.0208312); expected 1769.5 (0.02) Passed Vit filter: 120 (0.00135635); expected 88.5 (0.001) Passed Fwd filter: 8 (9.04231e-05); expected 0.9 (1e-05) Initial search space (Z): 88473 [actual number of targets] Domain search space (domZ): 8 [number of targets reported over threshold] # CPU time: 0.59u 0.03s 00:00:00.62 Elapsed: 00:00:00.27 # Mc/sec: 24419.47 // [ok]