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