# 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/Anolis_carolinensis.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 8.7e-142 472.8 0.0 1.1e-141 472.5 0.0 1.0 1 tr|H9G883|H9G883_ANOCA tRNA (guanine(26)-N(2))-dimethyltrans 3.8e-67 227.1 5.1 2.8e-57 194.7 1.2 3.1 3 tr|H9G676|H9G676_ANOCA tRNA methyltransferase 1 like OS=Anol 5.4e-05 22.6 0.1 8.2e-05 22.0 0.1 1.5 1 tr|H9GBK0|H9GBK0_ANOCA tRNA (guanine(37)-N1)-methyltransfera ------ inclusion threshold ------ 0.046 12.9 0.0 0.073 12.3 0.0 1.2 1 tr|G1KMS3|G1KMS3_ANOCA Uncharacterized protein OS=Anolis car Domain annotation for each sequence (and alignments): >> tr|H9G883|H9G883_ANOCA tRNA (guanine(26)-N(2))-dimethyltransferase OS=Anolis carolinensis OX=28377 GN=TRMT1 PE=3 SV= # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 472.5 0.0 2.3e-145 1.1e-141 4 376 .] 41 485 .. 38 485 .. 0.92 Alignments for each domain: == domain 1 score: 472.5 bits; conditional E-value: 2.3e-145 EEEE----SS--TT--SS--GGGHHHHHHHHHH........HHHHHHHH-S-................................... CS TRM 4 kikvedkevtvskkeeVFYNPvqefnRDlsvlvi....r.qlellkkklgek................................... 50 i++ ++ ++++ ++eVFYNPvqefnRDl+++vi r ql+ ++ +l+ + tr|H9G883|H9G883_ANOCA 41 LITEGKARIIFPCANEVFYNPVQEFNRDLTCAVItefaRlQLSSKGIRLAVPgeensskvvadlskeekddnleqhpsddqlqeegh 127 4555556699***********************9555444455555555555*********************************** PP ............EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHHHHC...--SS-EE CS TRM 51 ............ikvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneveeivvineddAnvlla...ekkkrfd 122 ++vl+aL+asGlR+iRfa+Evpg ++v+aND+s+kAveli +Nv+ln+v ++v+++ +dA+++++ + k fd tr|H9G883|H9G883_ANOCA 128 kvaavgqvcqdgLRVLEALAASGLRSIRFAKEVPGLRSVVANDFSSKAVELITQNVQLNDVGHLVTPSFSDARMVMYqckADKALFD 214 *****************************************************************************88777889** PP EEEE--SS--GGGHHHHHHHEEEEEEEEEEE--HHHHTTSSHHHHHHHHSSB---STTHHHHHHHHHHHHHHHHHHTTTEEEEEEEE CS TRM 123 viDlDPfGspapfldsavqsvkekglLlvTaTDtavLaGaspekclrkYgavalkselcHEvalRillkavareaakYekaveplls 209 viDlDP+Gsp++fld+avq+v+e+glL+vT+TD+av+aG+s e+c++kYga+++k ++cHE+alRi+l++++ +a++Y+++++plls tr|H9G883|H9G883_ANOCA 215 VIDLDPYGSPSTFLDAAVQAVSEGGLLCVTCTDMAVMAGNSGETCYSKYGAMSIKARFCHEMALRIILHSLDLRANCYQRYIVPLLS 301 *************************************************************************************** PP EEETTEEEEEEEEEE-HHHHHHHHTTEEEEEEETTT..........--EEEE--GGG--........SB-TTT-SBEEEEEEEE-S- CS TRM 210 lsadfYvrvfvkvktgaakvkkvleklglvlhCsgC..........lseekvkglakfk........aecehcgaklhlaGPlwlgp 278 +s+dfY+rv+v+v t + kvk++++k++lv++C+gC + ++ ++ ++k cehc ++++laGP+w++p tr|H9G883|H9G883_ANOCA 302 VSVDFYIRVYVRVFTAQSKVKASASKQALVYNCVGCgsfhiqrlgkMAHNGNN--INYKaacgppvgPSCEHCHQRHQLAGPFWAEP 386 ************************************99887776552222222..222255789999******************** PP S--HHHHHHHHHHHCT...-TTS-HHHHHHHHHHHHHHS--SS-EEHHHHHHHHT-SS---HHHHHHH....TT-EE-TTSTTEEE- CS TRM 279 lkdkefveevleeaekek.eeelkkrvkklLklikeElldaplyydlsqlasvlklsvpplkkvvsalkelgfevsrtHakpkaikT 364 l+d +fv++vl ++ ++ + ++ +r++++L+++ eEl+d+plyy+l++l++++++s+p+l +++sal+++g++vs++Ha+++aikT tr|H9G883|H9G883_ANOCA 387 LHDMAFVQQVLSAVGSNPgRFQTGQRIQGILSMALEELSDVPLYYTLDNLSKTVHCSSPSLLQFRSALLHAGYRVSLSHACKSAIKT 473 ***************9754445599************************************************************** PP SS-HHHHHHHHH CS TRM 365 napleavvevlr 376 +ap +v++++r tr|H9G883|H9G883_ANOCA 474 DAPPPVVWDIMR 485 **********97 PP >> tr|H9G676|H9G676_ANOCA tRNA methyltransferase 1 like OS=Anolis carolinensis OX=28377 GN=TRMT1L PE=3 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -2.3 0.1 0.44 2e+03 228 250 .. 172 207 .. 168 237 .. 0.64 2 ! 194.7 1.2 6e-61 2.8e-57 12 311 .. 240 561 .. 222 585 .. 0.85 3 ! 32.5 0.0 1.1e-11 5.3e-08 316 374 .. 626 684 .. 603 686 .. 0.82 Alignments for each domain: == domain 1 score: -2.3 bits; conditional E-value: 0.44 HHHHHHTTEEEEEEETTT.............--EEE CS TRM 228 kvkkvleklglvlhCsgC.............lseek 250 ++ + k+g +hC C +++++ tr|H9G676|H9G676_ANOCA 172 TADPISAKMGSHYHCLICsatigrrtdmighVKRHV 207 567788899999******766655544444322222 PP == domain 2 score: 194.7 bits; conditional E-value: 6e-61 SS--TT--SS--GGGHHHHHHHHHH...HHHHHHHH-S-EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT CS TRM 12 vtvskkeeVFYNPvqefnRDlsvlvirqlellkkklgekikvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklne 98 t ++k++ ++NP ++ nR l+v + l++l ++ ++++ lda+ a+G+ ++ +a+ + + +v +ND +ek v +i+eN +ln+ tr|H9G676|H9G676_ANOCA 240 YTTPQKTDSYFNPKMKLNRQLTVCA---LAVLAEE-RKPLECLDAFGATGIMGLQWAKHLGNSVKVTINDYNEKSVTMIQENCHLNK 322 457899999****************...7777777.689*******************************************88888 PP --......................GGEEEE-S-HHHHHC--SS-EEEEEE--SS--GGGHHHHHHHEEEEEEEEEEE--HHHHTTSS CS TRM 99 ve......................eivvineddAnvllaekkkrfdviDlDPfGspapfldsavqsvkekglLlvTaTDtavLaGas 163 ++ ++++ ++dAnv++ + ++fd+i lDPfGs++++ld+a + v++ g+++vT+TD ++L ++ tr|H9G676|H9G676_ANOCA 323 MKvvmnlreeedcdeaveeveenlPVIEVTKMDANVIM--HLRAFDFIHLDPFGSSVNYLDAAFRNVRNLGIVSVTSTDISSLYSKA 407 77777888888898888888875448999*********..889******************************************** PP HHHHHHHHSSB---STTHHHHHHHHHHHHHHHHHHTTTEEEEEEEEEEETTEEEEEEEEEE-HHHHHHHHTTEEEEEEETTT..... CS TRM 164 pekclrkYgavalkselcHEvalRillkavareaakYekavepllslsadfYvrvfvkvktgaakvkkvleklglvlhCsgC..... 245 ++ r+Yg++ +++e++ E a R+++ avar+aa+ +k +e+ll+++ +++v v v+v +g+ ++++++k+ +++hC C tr|H9G676|H9G676_ANOCA 408 QHVAQRHYGSNIVRTEYYKELAARMIIAAVARAAARCNKGIEVLLAVALEHFVLVEVRVLRGPSPADDTAKKIRYLIHCQWCeervf 494 **********************************************************************************55544 PP .....--EEEE--GGG--........SB-TTT-SBEEEEEEEE-S-S--HHHHHHHHHHHCT..-TTS.-HHHHHHHHHH CS TRM 246 .....lseekvkglakfk........aecehcgaklhlaGPlwlgplkdkefveevleeaekekeeel.kkrvkklLkli 311 + + + + g++ GPlw g l ++ f+ ++l ea + + ++++ lLk++ tr|H9G676|H9G676_ANOCA 495 kkegnM-----IEE---NpyrqlpcdCHSRMTGKTAIELGPLWSGSLFNTGFLRRMLFEAVQ-----YgLEDIQPLLKTL 561 222221.....222...2344555334445556666678*****************998887.....3256666666654 PP == domain 3 score: 32.5 bits; conditional E-value: 1.1e-11 S--SS-EEHHHHHHHHT-SS---HHHHHHH....TT-EE-TTSTTEEE-SS-HHHHHHH CS TRM 316 ldaplyydlsqlasvlklsvpplkkvvsalkelgfevsrtHakpkaikTnapleavvev 374 ++p +y + s+ +++p+l+k+ + l e+g++vsrtH+ p +++Tnapl ++ + tr|H9G676|H9G676_ANOCA 626 AEHPSFYYNIHRHSIKGMNMPKLNKFLHYLSEAGYRVSRTHFDPMGVRTNAPLVQFKSI 684 44555555555578899***********************************9887666 PP >> tr|H9GBK0|H9GBK0_ANOCA tRNA (guanine(37)-N1)-methyltransferase OS=Anolis carolinensis OX=28377 GN=TRMT5 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 22.0 0.1 1.8e-08 8.2e-05 53 115 .. 303 362 .. 92 371 .. 0.86 Alignments for each domain: == domain 1 score: 22.0 bits; conditional E-value: 1.8e-08 EEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHHHHC CS TRM 53 vldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneveeivvineddAnvlla 115 ++d +++ G ai a+ v+aND++++ ++ +++N kln+v+ v++ ++d + +++ tr|H9GBK0|H9GBK0_ANOCA 303 LFDVFAGVGPFAIPAAKK---NCIVFANDLNPESYKWLQRNCKLNKVDTKVQVFNMDGREFMK 362 567777777777777766...5579***********************99999999*****98 PP >> tr|G1KMS3|G1KMS3_ANOCA Uncharacterized protein OS=Anolis carolinensis OX=28377 GN=etfbkmt PE=4 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.3 0.0 1.6e-05 0.073 44 100 .. 112 166 .. 104 183 .. 0.83 Alignments for each domain: == domain 1 score: 12.3 bits; conditional E-value: 1.6e-05 HHHH-S-EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT-- CS TRM 44 kkklgekikvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneve 100 + k++++ kvld s+ G ai ++g+++v+aNDi + A + N +ln+v+ tr|G1KMS3|G1KMS3_ANOCA 112 NPKVVRRKKVLDLGSGCGATAIAAM--MSGASQVVANDIDPVAAAATVLNCELNNVN 166 5567777789998888888888655..56**************************98 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (376 nodes) Target sequences: 18528 (9796552 residues searched) Passed MSV filter: 670 (0.0361615); expected 370.6 (0.02) Passed bias filter: 462 (0.0249352); expected 370.6 (0.02) Passed Vit filter: 45 (0.00242876); expected 18.5 (0.001) Passed Fwd filter: 9 (0.000485751); expected 0.2 (1e-05) Initial search space (Z): 18528 [actual number of targets] Domain search space (domZ): 4 [number of targets reported over threshold] # CPU time: 0.38u 0.02s 00:00:00.40 Elapsed: 00:00:00.13 # Mc/sec: 28334.64 // [ok]