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