# 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/TRM.hmm.txt # target sequence database: proteomes/Natronomonas_moolapensis.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: TRM [M=376] Accession: PF02005.15 Description: N2,N2-dimethylguanosine tRNA 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.7e-86 288.1 0.0 1.9e-86 287.9 0.0 1.0 1 tr|M1XK62|M1XK62_NATM8 tRNA (guanine(26)-N(2))-dimethyltrans ------ inclusion threshold ------ 0.011 12.3 0.0 0.016 11.7 0.0 1.2 1 tr|M1Y3X8|M1Y3X8_NATM8 tRNA (Guanine(37)-N(1))-methyltransfe Domain annotation for each sequence (and alignments): >> tr|M1XK62|M1XK62_NATM8 tRNA (guanine(26)-N(2))-dimethyltransferase OS=Natronomonas moolapensis (strain DSM 18674 / J # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 287.9 0.0 1.4e-89 1.9e-86 1 373 [. 6 363 .. 6 366 .] 0.90 Alignments for each domain: == domain 1 score: 287.9 bits; conditional E-value: 1.4e-89 TTEEEEE----SS..--TT--SS--GGGHHHHHHHHHH...HHHHHHHH-S-EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HH CS TRM 1 Geakikvedkevt..vskkeeVFYNPvqefnRDlsvlvirqlellkkklgekikvldaLsasGlRaiRfalEvpgvkevlaNDisek 85 G+++i+v+++ + + VFYNPvqe nRD++ +v+r ++ + ++lda++asG+Ra+R a ++ Di e tr|M1XK62|M1XK62_NATM8 6 GAVTIEVPEERHGasEGSGSGVFYNPVQELNRDVTAAVLRVVDDECD------SYLDAMTASGVRAVRAAEA---GYAATGCDIDEA 83 566777777653322456688*****************665554444......48*************8855...4556678***** PP HHHHHHHHHHHTT--GGEEEE-S-HHHHHC--SS-EEEEEE--SS--GGGHHHHHHHEEEEEEEEEEE--HHHHTTSSHHHHHHHHS CS TRM 86 AvelikeNvklneveeivvineddAnvllaekkkrfdviDlDPfGspapfldsavqsvkekglLlvTaTDtavLaGaspekclrkYg 172 Av+l + N + n + +i++ + n+ + +++ +dv+DlDPfG+p pf d+a +s++e +++vTaTDta+L+Ga+ e+++r Yg tr|M1XK62|M1XK62_NATM8 84 AVDLARSNLERNGLSG--TIHHRNVNAHM--HEHGHDVVDLDPFGTPIPFADAAFRSARE--YVCVTATDTAPLCGAHFESGVRSYG 164 ***************9..55555677778..67789**********************99..************************* PP SB---STTHHHHHHHHHHHHHHHHHHTTTEEEEEEEEEEETTEEEEEEEEEE-HHHHHHHHTTEEEEEEETTT--E.EEE--GGG-- CS TRM 173 avalkselcHEvalRillkavareaakYekavepllslsadfYvrvfvkvktgaakvkkvleklglvlhCsgClse.ekvkglakfk 258 av+ ++e+ E++lR+ll+a++r+aa+Y+ a +p+ls+++ +Yvr +++++ ga+ +++++ +lg+v hC +Cl + ++ + +a tr|M1XK62|M1XK62_NATM8 165 AVPRNTEFHPEMGLRVLLSALVRTAARYDLAATPVLSHVSSHYVRTYLRLDGGARAADELFGELGYVDHCQRCLWRaHERTLVADPT 251 *************************************************************************77624444455566 PP SB-TTT-SBEEEEEEEE-S-S--HHHHHHHHHHHCT..-TTS-HHHHHHHHHHHHHHS--SS-EEHHHHHHHHT-SS---HHHHHHH CS TRM 259 aecehcgaklhlaGPlwlgplkdkefveevleeaekekeeelkkrvkklLklikeElldaplyydlsqlasvlklsvpplkkvvsal 345 a c++cg+++ aGP+wlgp +d +fve+v e++ + + ++ + lL ++ E l+ p+ yd ++l ++ ++++ al tr|M1XK62|M1XK62_NATM8 252 AACPECGQSTWTAGPIWLGPAHDVAFVERVGEAVPE--SFGTASAARDLLGTVAGE-LHRPTHYDQHRLYKRWGEPNTAMDEFLGAL 335 9*********************************99..4566888999********.9***************************** PP ....TT-EE-TTSTTEEE-SS-HHHHHH CS TRM 346 kelgfevsrtHakpkaikTnapleavve 373 +e+g+e+srtH+ ++++kT+a + ++ e tr|M1XK62|M1XK62_NATM8 336 REAGHEASRTHYGGTTFKTDAGVAEIRE 363 **********************999876 PP >> tr|M1Y3X8|M1Y3X8_NATM8 tRNA (Guanine(37)-N(1))-methyltransferase OS=Natronomonas moolapensis (strain DSM 18674 / JCM # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.7 0.0 1.2e-05 0.016 79 149 .. 204 273 .. 180 275 .. 0.79 Alignments for each domain: == domain 1 score: 11.7 bits; conditional E-value: 1.2e-05 EE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHHHHC--SS-EEEEEE--SS--GGGHHHHHHHEEEEEEE CS TRM 79 aNDisekAvelikeNvklneveeivvineddAnvllaekkkrfdviDlDPfGspapfldsavqsvkekglL 149 D++e Ave +++N++ n v+++v + ++d + + + d i + s+ fld+av+ + e+ lL tr|M1Y3X8|M1Y3X8_NATM8 204 GVDLNEAAVEYLRRNAERNGVSDRVIAINADVRDVT-GYDDWADRIVMNLPHSADGFLDTAVELAGEDCLL 273 569********************8877777888776.5555556666665567888888888777777666 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (376 nodes) Target sequences: 2723 (809592 residues searched) Passed MSV filter: 70 (0.0257069); expected 54.5 (0.02) Passed bias filter: 67 (0.0246052); expected 54.5 (0.02) Passed Vit filter: 7 (0.00257069); expected 2.7 (0.001) Passed Fwd filter: 2 (0.000734484); expected 0.0 (1e-05) Initial search space (Z): 2723 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.03u 0.00s 00:00:00.03 Elapsed: 00:00:00.02 # Mc/sec: 15220.33 // [ok]