# 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/Trypanosoma_brucei.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
------- ------ ----- ------- ------ ----- ---- -- -------- -----------
2e-59 199.5 0.0 2.6e-59 199.1 0.0 1.1 1 tr|Q57UH6|Q57UH6_TRYB2 Uncharacterized protein OS=Trypanosom
2.7e-37 127.3 0.0 3.3e-37 127.0 0.0 1.1 1 tr|Q586T6|Q586T6_TRYB2 Uncharacterized protein OS=Trypanosom
Domain annotation for each sequence (and alignments):
>> tr|Q57UH6|Q57UH6_TRYB2 Uncharacterized protein OS=Trypanosoma brucei brucei (strain 927/4 GUTat10.1) OX=185431 GN=Tb
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 199.1 0.0 6.2e-63 2.6e-59 1 194 [] 25 267 .. 25 267 .. 0.87
Alignments for each domain:
== domain 1 score: 199.1 bits; conditional E-value: 6.2e-63
EEEEEES--B--GGSCGGSB--HHHHCGG.TTSTTGGGSS-CHHHHHHHHHHHCSCHHCTT-EEEEEEETTSEEEEE-TT----SSH CS
EMG1 1 ivLeeasLelvkvkkkseqLLnsddhkkilkklkreakeaRpDIvHqcLlalldspLnkagllqvyihttknvlIevnpqtRiPrtf 87
i+Le+++L++v+ k + ++LL +++h+ +++k++++e+RpD+vHq+Ll+l+dspLn+ag+lqv+++t+k+v+I+v+p++R+Pr++
tr|Q57UH6|Q57UH6_TRYB2 25 IILEHCPLAVVRGK-YGFELL-AERHRAHHARYKQDPAEWRPDVVHQALLHLMDSPLNRAGMLQVFLRTKKGVVIAVDPRLRVPRHM 109
8************7.******.7889899999******************************************************* PP
HHHHHHHHHHHCCSEEEETTEE.SCEEEECS-HHHHCCC-SEEEEE-TTSEE-.................................. CS
EMG1 88 krFagLmvqLLkklsirategeekLlkviknpitdhlppgskiillsesaekk.................................. 140
+ F+++mv++L ++++ra++g+ +Ll+v+knp+tdh+p+++++i ++++++ +
tr|Q57UH6|Q57UH6_TRYB2 110 RLFEKMMVSCLYRMKVRAANGSVSLLRVVKNPVTDHIPANCTLIRVEKDGDAVpdifaycatcgtnrhadgstsenvgtatacsssa 196
***************************************************778888888888887777777777777777777777 PP
..........................EEEE-SSSS....---HH....EEESSSS---HHHHHHHHHHHCH CS
EMG1 141 .........lvdleefv....vkpvvfvigamskgsvkidykee....eisisnypLsaavvcskiltafE 194
+++ ++++ +p++f+ig+mskg+v++ y + +i++ ++ +saa++cs i++ fE
tr|Q57UH6|Q57UH6_TRYB2 197 gsdakephfAQKHRDLAdakqFRPFAFIIGGMSKGDVDAPYAPRaqvkSIRLGDRGMSAAAACSVIIHGFE 267
76665555522333333554555888*********9998777645555********************998 PP
>> tr|Q586T6|Q586T6_TRYB2 Uncharacterized protein OS=Trypanosoma brucei brucei (strain 927/4 GUTat10.1) OX=185431 GN=Tb
# score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc
--- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ----
1 ! 127.0 0.0 7.7e-41 3.3e-37 1 193 [. 63 270 .. 63 271 .. 0.87
Alignments for each domain:
== domain 1 score: 127.0 bits; conditional E-value: 7.7e-41
EEEEEES--B--GGSCGGSB--HHHHCGG.TTSTTGGGSS-CHHHHHHHHHHHCSCHHC.TT-EEEEEE..ETTSEEEEE-TT---- CS
EMG1 1 ivLeeasLelvkvkkkseqLLnsddhkkilkklkreakeaRpDIvHqcLlalldspLnk.agllqvyih..ttknvlIevnpqtRiP 84
++Le+a L l +++ ++++ d + + + +++RpDIvHqcLlal+ds L + + +l+v+i+ + ++ +ev+p++R P
tr|Q586T6|Q586T6_TRYB2 63 VILEKAGLLLGREG-----VVDAYDRAGSTTIANASLRDVRPDIVHQCLLALFDSGLAQeQRQLRVFINlfSRGGRTVEVSPSLRPP 144
67888888877777.....8888888888999999**********************873679*****977899************* PP
SSHHHHHHHHHHHHCCSEEEETTEE.SCEEEECS-HHHHCCC-SEEEEE-TTSEE-..............................E CS
EMG1 85 rtfkrFagLmvqLLkklsirategeekLlkviknpitdhlppgskiillsesaekklvdleefv....vkpvv.............f 154
rt +rF+gL++ LL+++ +r+++ge +Ll+v+ +++ +lp g+++i+l++ + + ++++++ ++pv f
tr|Q586T6|Q586T6_TRYB2 145 RTLRRFKGLVAALLRDGVVRSQDGE-VLLQVLPGSLAPVLPYGAEVIGLTNAHSAPVRTPTQLAraalENPVDdklqgglkgvagfF 230
************************9.***********************99877444445555555565555577778888888889 PP
EEE-SSSS.....---HHEEESSSS---HHHHHHHHHHHC CS
EMG1 155 vigamskgsvki.dykeeeisisnypLsaavvcskiltaf 193
+i + gs + dy+++++++s yp++a+v+cs++++ f
tr|Q586T6|Q586T6_TRYB2 231 CICCSEDGSLDGvDYVTQQVCLSAYPMTAHVMCSRLCEGF 270
9999999976656************************987 PP
Internal pipeline statistics summary:
-------------------------------------
Query model(s): 1 (194 nodes)
Target sequences: 8561 (4332363 residues searched)
Passed MSV filter: 207 (0.0241794); expected 171.2 (0.02)
Passed bias filter: 173 (0.0202079); expected 171.2 (0.02)
Passed Vit filter: 10 (0.00116809); expected 8.6 (0.001)
Passed Fwd filter: 2 (0.000233618); expected 0.1 (1e-05)
Initial search space (Z): 8561 [actual number of targets]
Domain search space (domZ): 2 [number of targets reported over threshold]
# CPU time: 0.07u 0.01s 00:00:00.08 Elapsed: 00:00:00.05
# Mc/sec: 16809.57
//
[ok]