# 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/DREV.hmm.txt # target sequence database: proteomes/Trypanosoma_brucei.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: DREV [M=265] Accession: PF05219.11 Description: DREV 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2.1e-56 189.9 0.0 2.6e-56 189.6 0.0 1.0 1 tr|Q584Q2|Q584Q2_TRYB2 DREV methyltransferase, putative OS=T 0.0026 15.9 0.0 0.52 8.3 0.0 2.3 2 tr|Q38FT0|Q38FT0_TRYB2 Uncharacterized protein OS=Trypanosom 0.0099 14.0 0.0 5.2 5.0 0.0 3.2 3 tr|Q582P2|Q582P2_TRYB2 Uncharacterized protein OS=Trypanosom ------ inclusion threshold ------ 0.019 13.0 0.0 0.032 12.3 0.0 1.3 1 tr|Q585P9|Q585P9_TRYB2 rRNA adenine N(6)-methyltransferase O 0.021 12.9 0.0 0.05 11.6 0.0 1.5 1 tr|Q38DV8|Q38DV8_TRYB2 Ubiquinone biosynthesis O-methyltrans Domain annotation for each sequence (and alignments): >> tr|Q584Q2|Q584Q2_TRYB2 DREV methyltransferase, putative OS=Trypanosoma brucei brucei (strain 927/4 GUTat10.1) OX=185 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 189.6 0.0 1.5e-59 2.6e-56 2 264 .. 24 298 .. 23 299 .. 0.88 Alignments for each domain: == domain 1 score: 189.6 bits; conditional E-value: 1.5e-59 DREV 2 yvlnre.lseslqalfvaslpdaetqsflenskklsanlfiqlyhalassilslfmsktdinGllgrGsmfvlseeqfrkllaie.. 85 y +n + l++++ +++ ++ d+et +l+++ ++sa + l +l f+s t nG++grG mfv+s eq + ll tr|Q584Q2|Q584Q2_TRYB2 24 YRANIDvLPQHIVTKYFILEKDKETTDWLKEAANVSA------LRILIANLLRHFVSLTTANGIVGRGGMFVVSTEQAANLLLGRti 104 66777779999999***************99999887......58899999*****************************9964411 PP DREV 86 kespasdqed.....kklldlGaGdGeitkrvaptfeevyatelsetmrdrlkkkkynvlsale.lqetd...vkldlilclnlldr 163 ++++ +++ ++lld+GaGdG +t r++p f++v ate s mr+rl +++y+vl + + + d +d+i+cln+ldr tr|Q584Q2|Q584Q2_TRYB2 105 AQQWTRYRDTmeppfETLLDVGAGDGGVTSRLQPLFKRVTATEFSIPMRWRLWRRGYEVLPYQDpFTNKDgsrRYYDVISCLNVLDR 191 35676544432234469***********************************************988876332469*********** PP DREV 164 cfdplklledirsalepngrvivalvlPvshyvetntgksekpsellevkG......rsyeeevarli.evlkkagfeveaftrlPy 243 + pl ll +r l p+g ++a+vlP+ ++ve + + +ps+++ ++G +++e+ + l+ +vl + gf++e++ +lPy tr|Q584Q2|Q584Q2_TRYB2 192 ADRPLDLLASMRDSLGPEGILLLAVVLPWCPFVEDGS-IKRRPSQMMPMDGgecckgATFEQSLQCLMdNVLVPFGFQLERWCKLPY 277 **********************************988.77888888887665666559**********9****************** PP DREV 244 lceGdlnnsfyylddaifvlk 264 lceG++ ++ +l+da++vl tr|Q584Q2|Q584Q2_TRYB2 278 LCEGNMIAEYAVLHDAVMVLT 298 *******************96 PP >> tr|Q38FT0|Q38FT0_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 ! 8.3 0.0 0.0003 0.52 85 129 .. 126 170 .. 119 183 .. 0.86 2 ! 4.9 0.0 0.0032 5.5 149 191 .. 260 302 .. 251 337 .. 0.81 Alignments for each domain: == domain 1 score: 8.3 bits; conditional E-value: 0.0003 DREV 85 ekespasdqedkklldlGaGdGeitkrvaptfeevyatelsetmr 129 +++ pa d ++++ ldlG G + ++++f++v + s + tr|Q38FT0|Q38FT0_TRYB2 126 HEFVPANDGKERRALDLGCAVGASSLEMSKYFSRVVGIDYSVAFI 170 56789999**************************99888887655 PP == domain 2 score: 4.9 bits; conditional E-value: 0.0032 DREV 149 dvkldlilclnlldrcfdplklledirsalepngrvivalvlP 191 + ++d +l+ nll r +p kll+ + + l + g +++ tr|Q38FT0|Q38FT0_TRYB2 260 GERFDAVLVANLLCRVPNPRKLLDMLPLLLVSGGILVISSPYS 302 5689********************9999999988888775555 PP >> tr|Q582P2|Q582P2_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 ? 0.4 0.0 0.08 1.4e+02 108 144 .. 102 138 .. 98 145 .. 0.84 2 ! 3.8 0.0 0.0074 13 83 127 .. 534 578 .. 524 584 .. 0.86 3 ! 5.0 0.0 0.003 5.2 151 170 .. 639 658 .. 627 676 .. 0.83 Alignments for each domain: == domain 1 score: 0.4 bits; conditional E-value: 0.08 DREV 108 itkrvaptfeevyatelsetmrdrlkkkkynvlsale 144 ++ r++p++ee +a + e d l+ +yn ++ e tr|Q582P2|Q582P2_TRYB2 102 LSMRIRPHLEEAFAVGVDEMSWDDLNEEHYNWPTVEE 138 6889***********9999999999999999766555 PP == domain 2 score: 3.8 bits; conditional E-value: 0.0074 DREV 83 aiekespasdqedkklldlGaGdGeitkrvaptfeevyatelset 127 ++++++ + +k+ ld+ G G i+ ++++f++v + s + tr|Q582P2|Q582P2_TRYB2 534 VMNQFRCCARTGPKRALDFACGAGRIAFELTSMFDQVVGSDFSAR 578 5677888888899************************99988876 PP == domain 3 score: 5.0 bits; conditional E-value: 0.003 DREV 151 kldlilclnlldrcfdplkl 170 ++ lilc n ldr ++p + tr|Q582P2|Q582P2_TRYB2 639 DFSLILCWNTLDRSYNPTAV 658 699*************9865 PP >> tr|Q585P9|Q585P9_TRYB2 rRNA adenine N(6)-methyltransferase OS=Trypanosoma brucei brucei (strain 927/4 GUTat10.1) OX= # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.3 0.0 1.9e-05 0.032 89 135 .. 83 129 .. 57 140 .. 0.84 Alignments for each domain: == domain 1 score: 12.3 bits; conditional E-value: 1.9e-05 DREV 89 pasdqedkklldlGaGdGeitkrvaptfeevyatelsetmrdrlkkk 135 a+ + + +l++G G G++t+++ + ++v a e+ +m l+k+ tr|Q585P9|Q585P9_TRYB2 83 KAAVKPTDIVLEIGPGTGNLTEKLLQAAKKVIAFEVDPRMVAELNKR 129 3444556789**********************************986 PP >> tr|Q38DV8|Q38DV8_TRYB2 Ubiquinone biosynthesis O-methyltransferase, mitochondrial OS=Trypanosoma brucei brucei (stra # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.6 0.0 2.9e-05 0.05 97 187 .. 69 174 .. 58 180 .. 0.77 Alignments for each domain: == domain 1 score: 11.6 bits; conditional E-value: 2.9e-05 DREV 97 klldlGaGdGeitkrvaptfee...........vyatelsetmrdrl...kkkkynvlsale.lqetdvkldlilclnlldrcfdpl 168 +ld+G G G +++ +a + a e +++ + + y+ + e ++ + +++l+++ ++++ d+ tr|Q38DV8|Q38DV8_TRYB2 69 SVLDVGCGGGILSESIARLGGDvlgidmckesiAVAEERRQRVLRDItpqASLAYRCVPLHEvVEGEGRQFELVVASEVIEHVSDAA 155 688999999999998888877766666666664445555554443331214578999999995566789****************** PP DREV 169 klledirsalepngrviva 187 +l+di a +p g +++ tr|Q38DV8|Q38DV8_TRYB2 156 AFLHDICAATKPGGVLVIS 174 ************9988887 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (265 nodes) Target sequences: 8561 (4332363 residues searched) Passed MSV filter: 129 (0.0150683); expected 171.2 (0.02) Passed bias filter: 106 (0.0123817); expected 171.2 (0.02) Passed Vit filter: 11 (0.0012849); expected 8.6 (0.001) Passed Fwd filter: 5 (0.000584044); expected 0.1 (1e-05) Initial search space (Z): 8561 [actual number of targets] Domain search space (domZ): 5 [number of targets reported over threshold] # CPU time: 0.07u 0.00s 00:00:00.07 Elapsed: 00:00:00.04 # Mc/sec: 28701.90 // [ok]