# 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/Saccharomyces_cerevisiae.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 7e-39 132.2 0.0 1.1e-38 131.5 0.0 1.3 1 sp|Q08282|TYW4_YEAST tRNA wybutosine-synthesizing protein 4 5.6e-31 106.4 0.0 7e-31 106.1 0.0 1.1 1 sp|Q04081|LCMT1_YEAST Leucine carboxyl methyltransferase 1 O Domain annotation for each sequence (and alignments): >> sp|Q08282|TYW4_YEAST tRNA wybutosine-synthesizing protein 4 OS=Saccharomyces cerevisiae (strain ATCC 204508 / S288c) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 131.5 0.0 3.7e-42 1.1e-38 2 187 .. 7 242 .. 6 243 .. 0.95 Alignments for each domain: == domain 1 score: 131.5 bits; conditional E-value: 3.7e-42 HHHHHHHHHHHHHH..HHTSSHTCCCCCCCHHHHHHHH............CHCCCHCHHCHG.CHCCCCT--HHHHHHHHHHHHHHHHH CS LCM 2 lgvaaaralesrra..adglidDplaeslvsaaglell............erkadkelekek.eeekkekkapelnrgiavRtrffddl 75 ++++++++ ++rr+ ad +i++++++s++s++++ell ++k++k le++k +++kk+k++p++nrg+++R ++++++ sp|Q08282|TYW4_YEAST 7 IKQTNKNVKQERRKkyADLAIQGTNNSSIASKRSVELLylpklssannfqMDKNNKLLEYFKfFVPKKIKRSPCINRGYWLRLFAIRSR 95 7899*********66**********************9************************99************************* PP HHHHH.HCTT..ESEEEEET-TT--HHHHHHHTT.........TEEEEEEE-HHHHHHHHHHHHHSHHHHH.................. CS LCM 76 lkaal.saag..irqvViLGAGlDtrayRlpwla.........glkvfEvDqpevlefKretlkelgaepk.................. 134 l++++ ++ + vV+LG+G+D+++++l +++ +++++++D++++l++K+e++k++++++k sp|Q08282|TYW4_YEAST 96 LNSIIeQTPQdkKIVVVNLGCGYDPLPFQLLDTNniqsqqyhdRVSFIDIDYSDLLKIKIELIKTIPELSKiiglsedkdyvddsnvdf 184 *****88777764489**********************************************************************999 PP .EECSE..EE--TTSHHHHHHHHHH..TT-TTSEEEEEEESSGGGS-HHHHHHHHHHH CS LCM 135 .aelry..vpvDlrdddwpeaLaaa..gfdpekptawvaeGvlvYLtpeasdalleai 187 ++ +y p+Dl+d++++++L ++ +dp++ +++vae+ l+Y++pe+sd+++ea+ sp|Q08282|TYW4_YEAST 185 lTTPKYlaRPCDLNDSKMFSTLLNEcqLYDPNVVKVFVAEVSLAYMKPERSDSIIEAT 242 98888788**************999999***************************986 PP >> sp|Q04081|LCMT1_YEAST Leucine carboxyl methyltransferase 1 OS=Saccharomyces cerevisiae (strain ATCC 204508 / S288c) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 106.1 0.0 2.3e-34 7e-31 2 187 .. 5 219 .. 4 220 .. 0.89 Alignments for each domain: == domain 1 score: 106.1 bits; conditional E-value: 2.3e-34 HHHHHHHHHHHHHH..HHTSSHTCCCCCCCHHHHHHHH..CHCCCHCHHCHG.C.....HCCCCT--HHHHHHHHHHHHHHHHHHHHH CS LCM 2 lgvaaaralesrra..adglidDplaeslvsaaglell..erkadkelekek.e.....eekkekkapelnrgiavRtrffddllkaa 79 ++++ +al+++ a + g++ + ++l + +++++ +r++ +l+k++ + ++++++++p +n g+++Rt +d+++ ++ sp|Q04081|LCMT1_YEAST 5 IQQTDYDALSCKLAaiSVGYLPSSGLQRLSVDLSKKYTewHRSYLITLKKFSrRafgkvDKAMRSSFPVMNYGTYLRTVGIDAAILEF 92 68999*******9669****************99999988566666666655436777788999************************ PP H.HCTTESEEEEET-TT--HHHHHHHTT.TEEEEEEE-HHHHHHHHHHHHHSHHHHH................EECSE..EE--TTSH CS LCM 80 l.saagirqvViLGAGlDtrayRlpwla.glkvfEvDqpevlefKretlkelgaepk................aelry..vpvDlrdd 147 l +++++ qvV+LG+G D r++ l ++ +l ++++D++e +e+K+++l+e+++++ ++ ry ++Dl+d sp|Q04081|LCMT1_YEAST 93 LvANEKV-QVVNLGCGSDLRMLPLLQMFpHLAYVDIDYNESVELKNSILRESEILRIslglskedtakspfliDQGRYklAACDLNDI 179 *876777.*********************************************99999**************96666666*******6 PP HHHHHHHHHTT-TTSEEEEEEESSGGGS-HHHHHHHHHHH CS LCM 148 dwpeaLaaagfdpekptawvaeGvlvYLtpeasdalleai 187 + +++L + ++e pt++++e++l+Y+++++s+ l+++i sp|Q04081|LCMT1_YEAST 180 TETTRLLDVCTKREIPTIVISECLLCYMHNNESQLLINTI 219 6666677777999************************998 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (188 nodes) Target sequences: 6049 (2936363 residues searched) Passed MSV filter: 184 (0.0304183); expected 121.0 (0.02) Passed bias filter: 133 (0.0219871); expected 121.0 (0.02) Passed Vit filter: 11 (0.00181848); expected 6.0 (0.001) Passed Fwd filter: 2 (0.000330633); expected 0.1 (1e-05) Initial search space (Z): 6049 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.06u 0.00s 00:00:00.06 Elapsed: 00:00:00.03 # Mc/sec: 18401.21 // [ok]