# 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/Trichomonas_vaginalis.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 3.1e-37 129.8 0.0 4.3e-37 129.3 0.0 1.2 1 tr|A2G151|A2G151_TRIVA Methyltransferase, putative, TIGR0002 1.2e-26 95.3 0.0 2.7e-26 94.1 0.0 1.6 1 tr|A2E9W0|A2E9W0_TRIVA Leucine carboxyl methyltransferase 1 1.3e-26 95.2 0.0 2.7e-26 94.1 0.0 1.5 1 tr|A2FN76|A2FN76_TRIVA Leucine carboxyl methyltransferase 1 7.3e-26 92.7 0.0 1.7e-25 91.6 0.0 1.6 1 tr|A2EQ50|A2EQ50_TRIVA Leucine carboxyl methyltransferase 1 1.2e-19 72.4 0.0 2.3e-19 71.5 0.0 1.4 1 tr|A2EY25|A2EY25_TRIVA Leucine carboxyl methyltransferase 1 7.4e-17 63.4 0.0 3.2e-16 61.3 0.0 2.0 1 tr|A2F5C5|A2F5C5_TRIVA O-methyltransferase N-terminus family Domain annotation for each sequence (and alignments): >> tr|A2G151|A2G151_TRIVA Methyltransferase, putative, TIGR00027 family protein OS=Trichomonas vaginalis OX=5722 GN=TVA # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 129.3 0.0 5.1e-41 4.3e-37 1 186 [. 13 183 .. 13 185 .. 0.93 Alignments for each domain: == domain 1 score: 129.3 bits; conditional E-value: 5.1e-41 HHHHHHHHHHHHHHH.HHTSSHTCCCCCCCHHHHHHHHCHCCCHCHHCHGCHCCCCT--HHHHHHHHHHHHHHHHHHHHHHHCTTES CS LCM 1 algvaaaralesrra.adglidDplaeslvsaaglellerkadkelekekeeekkekkapelnrgiavRtrffddllkaalsaagir 86 a+++a++r l+ + ++ +i+D++ +l+++++ ++ ++ d+ k ++ ++ ++ R+r dd+ ++ ++gi+ tr|A2G151|A2G151_TRIVA 13 AVRTALWRGLHGLIDaKPLVIEDSVGLQLANPQEGWQ--ERPDM------------KFTKPVRGAMVGRARLLDDYAIECS-KKGIH 84 799***********99*****************9984..33333............236789999****************.6**** PP EEEEET-TT--HHHHHHHTTTEEEEEEE-HHHHHHHHHHHHHSHHHHHEECSEEE--TTSHHHHHHHHHHTT-TTSEEEEEEESSGG CS LCM 87 qvViLGAGlDtrayRlpwlaglkvfEvDqpevlefKretlkelgaepkaelryvpvDlrdddwpeaLaaagfdpekptawvaeGvlv 173 q+V LGAGlDt+a Rl + ++++fE+D + +l++K+ +lkel+ +ka +++vpv++ +++w++ L+++ f++++p+++ + Gv+ tr|A2G151|A2G151_TRIVA 85 QYVSLGAGLDTLAQRLGDSLNMQIFEIDEEPTLTWKEGRLKELNL-SKAYQHFVPVNFIKSSWFDILKQSSFKSDQPSLLTCVGVTL 170 *****************************************9996.6889************************************* PP GS-HHHHHHHHHH CS LCM 174 YLtpeasdallea 186 YL+++a+++ll++ tr|A2G151|A2G151_TRIVA 171 YLDNNAVHELLKN 183 ***********97 PP >> tr|A2E9W0|A2E9W0_TRIVA Leucine carboxyl methyltransferase 1 OS=Trichomonas vaginalis OX=5722 GN=TVAG_184240 PE=3 SV= # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 94.1 0.0 3.3e-30 2.7e-26 16 186 .. 26 177 .. 9 179 .. 0.86 Alignments for each domain: == domain 1 score: 94.1 bits; conditional E-value: 3.3e-30 HHTSSHTCCCCCCCHHHHHHHHCHCCCHCHHCHGCHCCCCT--HHHHHHHHHHHHHHHHHHHHHHHCTT.ESEEEEET-TT--HHHH CS LCM 16 adglidDplaeslvsaaglellerkadkelekekeeekkekkapelnrgiavRtrffddllkaalsaag.irqvViLGAGlDtrayR 101 ++g+++D ++ + ++++ k+ p +n g++vRt ++ ++k++ + +g qvV+LG+G+Dt+++R tr|A2E9W0|A2E9W0_TRIVA 26 QRGYYHDRYIGCF----------------------VPQELKQLPPMNLGYYVRTQAIYVAVKKFHEIHGsNMQVVVLGCGYDTLFWR 90 3444444444444......................55566779**********************6677788*************** PP HHHTT..TEEEEEEE-HHHHHHHHHHHHHSHHHHHEECSEEE--TTS.HHHHHHHHHHTT-TTSEEEEEEESSGGGS-HHHHHHHHH CS LCM 102 lpwla..glkvfEvDqpevlefKretlkelgaepkaelryvpvDlrd.ddwpeaLaaagfdpekptawvaeGvlvYLtpeasdalle 185 l +++ ++f++D++ v++ K ++++ ep ++ ++ +Dl + + ++++L ++gf+ + pt+++ e+ l+Y +p a+d++l+ tr|A2E9W0|A2E9W0_TRIVA 91 LRDEKvtVKNWFDLDMQHVVKKKSQVINSNSFEPMDNYHLFYADLSKpEMVKQELNSHGFE-DIPTIFIDECTLIYVDPVAVDEILK 176 *****977789************************************99999********9.************************9 PP H CS LCM 186 a 186 + tr|A2E9W0|A2E9W0_TRIVA 177 F 177 7 PP >> tr|A2FN76|A2FN76_TRIVA Leucine carboxyl methyltransferase 1 OS=Trichomonas vaginalis OX=5722 GN=TVAG_003240 PE=3 SV= # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 94.1 0.0 3.2e-30 2.7e-26 20 186 .. 30 177 .. 9 179 .. 0.85 Alignments for each domain: == domain 1 score: 94.1 bits; conditional E-value: 3.2e-30 SHTCCCCCCCHHHHHHHHCHCCCHCHHCHGCHCCCCT--HHHHHHHHHHHHHHHHHHHHHHHCTT.ESEEEEET-TT--HHHHHHHT CS LCM 20 idDplaeslvsaaglellerkadkelekekeeekkekkapelnrgiavRtrffddllkaalsaag.irqvViLGAGlDtrayRlpwl 105 ++D +++ + +++++k+ p +n g++vRt ++ ++k++ + +g qvV+LG+G+Dt+++Rl ++ tr|A2FN76|A2FN76_TRIVA 30 YKDRFIACF----------------------VTHEQKQLPPMNLGYYVRTQAIYVAVKKFYELHGkDMQVVVLGCGYDTLFWRLRDE 94 333333333......................45566779**********************6677677******************* PP T..TEEEEEEE-HHHHHHHHHHHHHSHHHHHEECSEEE--TTS.HHHHHHHHHHTT-TTSEEEEEEESSGGGS-HHHHHHHHHH CS LCM 106 a..glkvfEvDqpevlefKretlkelgaepkaelryvpvDlrd.ddwpeaLaaagfdpekptawvaeGvlvYLtpeasdallea 186 ++f++D+p +++ K ++lk +p ++ + +Dl + ++++++L+a+gf+ + pt+++ e+ l+Y +p+a+d++l++ tr|A2FN76|A2FN76_TRIVA 95 HitVKNWFDLDMPHIVNKKSSVLKSEAFHPLDNYYLFVADLSKpEEVKKELTAHGFE-DLPTIFIDECTLIYVDPNAVDEILKF 177 9977789************************99999********************9.************************97 PP >> tr|A2EQ50|A2EQ50_TRIVA Leucine carboxyl methyltransferase 1 OS=Trichomonas vaginalis OX=5722 GN=TVAG_474040 PE=3 SV= # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 91.6 0.0 2e-29 1.7e-25 16 186 .. 26 177 .. 9 179 .. 0.86 Alignments for each domain: == domain 1 score: 91.6 bits; conditional E-value: 2e-29 HHTSSHTCCCCCCCHHHHHHHHCHCCCHCHHCHGCHCCCCT--HHHHHHHHHHHHHHHHHHHHHHHCTT.ESEEEEET-TT--HHHH CS LCM 16 adglidDplaeslvsaaglellerkadkelekekeeekkekkapelnrgiavRtrffddllkaalsaag.irqvViLGAGlDtrayR 101 ++g+++D ++ + ++++ k+ p +n g++vRt ++ ++k++ + +g qvV+LG+G+Dt+++R tr|A2EQ50|A2EQ50_TRIVA 26 QRGYYHDRYIGCF----------------------VPQELKQLPPMNLGYYVRTQAIYVAVKKFHEIHGsNMQVVVLGCGYDTLFWR 90 3444444444444......................55566779**********************6677788*************** PP HHHTT..TEEEEEEE-HHHHHHHHHHHHHSHHHHHEECSEEE--TTS.HHHHHHHHHHTT-TTSEEEEEEESSGGGS-HHHHHHHHH CS LCM 102 lpwla..glkvfEvDqpevlefKretlkelgaepkaelryvpvDlrd.ddwpeaLaaagfdpekptawvaeGvlvYLtpeasdalle 185 l +++ ++f++D++ v++ K ++++ ep ++ ++ +Dl + + ++++L ++gf+ + pt+++ e+ l+Y +p ++d++l+ tr|A2EQ50|A2EQ50_TRIVA 91 LRDEKitVKNWFDLDMQHVVNKKSKVINSNSFEPMDNYHLFYADLSKpEMVKQELNSHGFE-DIPTIFIDECTLIYVDPVSVDEILK 176 *****977789************************************99999********9.***********************99 PP H CS LCM 186 a 186 + tr|A2EQ50|A2EQ50_TRIVA 177 F 177 7 PP >> tr|A2EY25|A2EY25_TRIVA Leucine carboxyl methyltransferase 1 OS=Trichomonas vaginalis OX=5722 GN=TVAG_044530 PE=3 SV= # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 71.5 0.0 2.7e-23 2.3e-19 52 187 .. 43 180 .. 8 181 .. 0.84 Alignments for each domain: == domain 1 score: 71.5 bits; conditional E-value: 2.7e-23 CCCCT--HHHHHHHHHHHHHHHHHHHHHH.HCTTESEEEEET-TT--HHHHHHHTT..TEEEEEEE-HHHHHHHHHHHHHSHHHHHE CS LCM 52 ekkekkapelnrgiavRtrffddllkaal.saagirqvViLGAGlDtrayRlpwla..glkvfEvDqpevlefKretlkelgaepka 135 ++ p+ n g+++R +++ ++k++ ++ qvV LG+G+Dt+ +Rl +++ +k++EvD++ v++ K++ l+++ +p + tr|A2EY25|A2EY25_TRIVA 43 ILYGSQLPHQNFGTFIRSYMIRYVVKSFFdKFGPKSQVVSLGSGFDTLYWRLRDENiqFSKFIEVDKQFVVAKKQKYLENQVFQPLD 129 444567788889999999999999999996655599*********************999*************************** PP ECSEEE--TTSHHHHHHHHHHTT-TTSEEEEEEESSGGGS-HHHHHHHHHHH CS LCM 136 elryvpvDlrdddwpeaLaaagfdpekptawvaeGvlvYLtpeasdalleai 187 + + +Dl+ + +++ + + p+ pt+++ e l+YL+p+ +++ll+++ tr|A2EY25|A2EY25_TRIVA 130 NYILCDLDLN-LEESTNVLNTLLIPDLPTLFIDEFSLIYLQPQVIQRLLKFF 180 *999******.5666666677899*************************986 PP >> tr|A2F5C5|A2F5C5_TRIVA O-methyltransferase N-terminus family protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_222490 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 61.3 0.0 3.9e-20 3.2e-16 6 186 .. 18 187 .. 14 189 .. 0.85 Alignments for each domain: == domain 1 score: 61.3 bits; conditional E-value: 3.9e-20 HHHHHHHHHHHHTSSHTCCCCCCCHHHHHHHHCHCCCHCHHCHGCHCCCCT--HHHHHHHHHHHHHHHHHHHHHHHCTT...ESEEE CS LCM 6 aaralesrraadglidDplaeslvsaaglellerkadkelekekeeekkekkapelnrgiavRtrffddllkaalsaag...irqvV 89 ra es+r +d +++Dp a l+++ ++ + + + + g+ vR ++fd+ + +++++++ ++ vV tr|A2F5C5|A2F5C5_TRIVA 18 YYRATESKRGEDAILHDPKAVELIGKIDYDF-------SVYPM---------DIKPSLGCPVRSWYFDQNIIEFIKSHQneeTCIVV 88 579*****988****************7775.......22222.........223356889**************5455578899** PP EET-TT--HHHHHHHTT.........TEEEEEEE-HHHHHHHHHHHHHSHHHHHEECSEEE--TTSHHHHHHHHHHTT-TTSEEEEE CS LCM 90 iLGAGlDtrayRlpwla.........glkvfEvDqpevlefKretlkelgaepkaelryvpvDlrdddwpeaLaaagfdpekptawv 167 +LG+GlDtr+ R++ la ++k++E+D p+v+e++++++ pk+++ +++ + + +w + L ++ +++ ++v tr|A2F5C5|A2F5C5_TRIVA 89 NLGCGLDTRMDRIEYLAnknndfslqNVKFVEIDLPDVIEYRETLIP-----PKENQIHISMSIFEVKWMRDLLEK--YGKVNYIFV 168 ****************99************************99994.....666666799998889999999999..44799**** PP EESSGGGS-HHHHHHHHHH CS LCM 168 aeGvlvYLtpeasdallea 186 aeGv +Y ++e+ ++ ++ tr|A2F5C5|A2F5C5_TRIVA 169 AEGVFPYFENETLKKFFDD 187 ***********99988875 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (188 nodes) Target sequences: 50190 (17055200 residues searched) Passed MSV filter: 1462 (0.0291293); expected 1003.8 (0.02) Passed bias filter: 978 (0.019486); expected 1003.8 (0.02) Passed Vit filter: 65 (0.00129508); expected 50.2 (0.001) Passed Fwd filter: 6 (0.000119546); expected 0.5 (1e-05) Initial search space (Z): 50190 [actual number of targets] Domain search space (domZ): 6 [number of targets reported over threshold] # CPU time: 0.27u 0.01s 00:00:00.28 Elapsed: 00:00:00.12 # Mc/sec: 26719.81 // [ok]