# 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/EMG1.hmm.txt
# target sequence database: proteomes/Trichomonas_vaginalis.fasta
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Query: EMG1 [M=194]
Accession: PF03587.13
Description: EMG1/NEP1 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.2e-27 98.3 0.1 1.4e-27 98.1 0.1 1.0 1 tr|A2EHM0|A2EHM0_TRIVA Uncharacterized protein OS=Trichomona
1.1e-11 46.2 0.2 1.5e-11 45.8 0.2 1.1 1 tr|A2E5Z3|A2E5Z3_TRIVA Uncharacterized protein OS=Trichomona
------ inclusion threshold ------
0.088 13.9 0.6 0.87 10.7 0.1 2.0 2 tr|A2D7S2|A2D7S2_TRIVA Uncharacterized protein OS=Trichomona
Domain annotation for each sequence (and alignments):
>> tr|A2EHM0|A2EHM0_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_312150 PE=4 SV=1
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 98.1 0.1 8.5e-32 1.4e-27 26 194 .] 44 211 .. 25 211 .. 0.90
Alignments for each domain:
== domain 1 score: 98.1 bits; conditional E-value: 8.5e-32
HCGG.TTSTTGGGSS-CHHH....HHHHHHHHCSCHHCTT-EEEEEEETTSEEEEE-TT----SSHHHHHHHHHHHHCCSEEEETTE CS
EMG1 26 hkkilkklkreakeaRpDIv....HqcLlalldspLnkagllqvyihttknvlIevnpqtRiPrtfkrFagLmvqLLkklsirateg 108
+ +++k +++DI+ +++L+ l+++ Ln +g+ +v ++t+++++I v+p++++P f+ + +++ L++++++ + e+
tr|A2EHM0|A2EHM0_TRIVA 44 E------EFESQKITKKDIIteqlYHVLYFLQNTALNLSGQQRVLLRTEDKQVIFVDPRLKVPDVFDDYLAMLHLLMEDRRVCSVEK 124
4......45677788899977777*************************************************************** PP
E.SCEEEECS-HHHHCCC-SEEEEE-TTSEE-..............EEEE-SSSS....---HHEEESSSS---HHHHHHHHHHHCH CS
EMG1 109 eekLlkviknpitdhlppgskiillsesaekklvdleefv.vkpvvfvigamskgsvkidykeeeisisnypLsaavvcskiltafE 194
++++k++k+++t+ +p+g+k i+ s +++++++ ef+ +p +++i++ ++ + ++ ++++sn Lsa++++s+il+++E
tr|A2EHM0|A2EHM0_TRIVA 125 GSEIMKFVKSDLTHSIPTGAKRIAAYLSKTDSVTSPDEFLkAEPTAVYIHLDPTVQSYDEVLLGKYCVSNDILSAVTLASRILHSME 211
*****************************999999*****888***********9555566667********************998 PP
>> tr|A2E5Z3|A2E5Z3_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_399550 PE=4 SV=1
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 45.8 0.2 8.9e-16 1.5e-11 42 194 .] 29 168 .. 26 168 .. 0.87
Alignments for each domain:
== domain 1 score: 45.8 bits; conditional E-value: 8.9e-16
CHHHHHHHHHHHCSCHHCTT-EEEEEEETTSEEEEE-TT----SSHHHHHHHHHHHHCCSEEEETTEE.SCEEEECS-HHHHCCC.- CS
EMG1 42 pDIvHqcLlalldspLnkagllqvyihttknvlIevnpqtRiPrtfkrFagLmvqLLkklsirategeekLlkviknpitdhlpp.g 127
D H +L+ l + Ln+ ++l+++++ + +ev+p + +P + +F +L+ ++L+ l+++ ++ + l + ++ lp+
tr|A2E5Z3|A2E5Z3_TRIVA 29 LDFAHRILYFLNSTSLNQYNKLRIFLKCK-TTYVEVSPMLLVPASLAEFINLLTNCLNGLEVKVDNIR--VLTL-----KKELPQlP 107
5999**********************986.5679******************************9976..4443.....34455534 PP
SEEEEE-TTSEE-.............EEEE-SSSS....---HHEEESSSS---HHHHHHHHHHHCH CS
EMG1 128 skiillsesaekklvdleefvvkpvvfvigamskgsvkidykeeeisisnypLsaavvcskiltafE 194
+ +++++ ++e k+vd++ +v +p+v+ ++ + + +e++++s+ypL+ + c +i++afE
tr|A2E5Z3|A2E5Z3_TRIVA 108 KIKYAVHPTSE-KVVDPSTIVFNPCVIELSRQPDM-----LDSEKVCFSKYPLDPITQCGSIVSAFE 168
56899999999.888******99988777777766.....7889*********************98 PP
>> tr|A2D7S2|A2D7S2_TRIVA Uncharacterized protein OS=Trichomonas vaginalis OX=5722 GN=TVAG_070180 PE=4 SV=1
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ? 1.9 0.1 0.025 4.3e+02 171 193 .. 36 58 .. 7 59 .. 0.75
2 ? 10.7 0.1 5.2e-05 0.87 85 141 .. 63 134 .. 62 179 .. 0.64
Alignments for each domain:
== domain 1 score: 1.9 bits; conditional E-value: 0.025
HEEESSSS---HHHHHHHHHHHC CS
EMG1 171 eeisisnypLsaavvcskiltaf 193
+++ is+y++ + v+s+i ta+
tr|A2D7S2|A2D7S2_TRIVA 36 TNVIISKYKMPISKVISQIFTAI 58
38899*************99975 PP
== domain 2 score: 10.7 bits; conditional E-value: 5.2e-05
SSHHHHHHHHHHHHCCSEEEETTEE.SCEEEEC.S-HHHHCCC-SEEEEE-TTSEE-................. CS
EMG1 85 rtfkrFagLmvqLLkklsirategeekLlkvik.npitdhlppgskiillsesaekk................l 141
r+f+ + +L+ q+L k si+++ ++ Ll+ + ++ +l++ ++++ ++ + +k
tr|A2D7S2|A2D7S2_TRIVA 63 RSFNLYIKLINQILLKYSITSKPSK-DLLEEAEfIHLQFILND-NNVYQITYDENKLffpkreeihdifidddI 134
89********************998.66654441555555555.566666666553355555555555554431 PP
Internal pipeline statistics summary:
-------------------------------------
Query model(s): 1 (194 nodes)
Target sequences: 50190 (17055200 residues searched)
Passed MSV filter: 2376 (0.0473401); expected 1003.8 (0.02)
Passed bias filter: 1077 (0.0214585); expected 1003.8 (0.02)
Passed Vit filter: 50 (0.000996214); expected 50.2 (0.001)
Passed Fwd filter: 3 (5.97729e-05); expected 0.5 (1e-05)
Initial search space (Z): 50190 [actual number of targets]
Domain search space (domZ): 3 [number of targets reported over threshold]
# CPU time: 0.30u 0.02s 00:00:00.32 Elapsed: 00:00:00.15
# Mc/sec: 22058.06
//
[ok]