# 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]