# 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/Petromyzon_marinus.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 5.4e-137 456.6 0.1 6.9e-137 456.2 0.1 1.1 1 tr|S4RHT4|S4RHT4_PETMA tRNA (guanine(26)-N(2))-dimethyltrans 9.4e-67 225.4 0.2 5.8e-56 189.9 0.1 2.3 2 tr|S4RMN3|S4RMN3_PETMA tRNA (guanine(26)-N(2))-dimethyltrans Domain annotation for each sequence (and alignments): >> tr|S4RHT4|S4RHT4_PETMA tRNA (guanine(26)-N(2))-dimethyltransferase OS=Petromyzon marinus OX=7757 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 456.2 0.1 1e-140 6.9e-137 7 376 .] 52 488 .. 46 488 .. 0.93 Alignments for each domain: == domain 1 score: 456.2 bits; conditional E-value: 1e-140 E----.SS--TT--SS--GGGHHHHHHHHHH...HHHHHHHH-S-.......................................... CS TRM 7 vedke.vtvskkeeVFYNPvqefnRDlsvlvirqlellkkklgek.......................................... 50 +e+k+ + +++ +eVFYNPvqefnRD++ +v+ e +++l e+ tr|S4RHT4|S4RHT4_PETMA 52 IEGKArIFFPNLSEVFYNPVQEFNRDITSAVV--TEFIRERLHERrirvvvpgeesrvvvslsdgpespgsvsaaspqyaekqayvg 136 4555547799*******************995..445555554445667899*********************************** PP ......EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHHHHC...--SS-EEEEEE-- CS TRM 51 ......ikvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneveeivvineddAnvlla...ekkkrfdviDlDP 128 +kvl+aL+asGlR+iRfa+E+pg k+++ ND+s++Av+ i++N+++n+v+++v+++ +dA++l++ ++k+ fdv+DlDP tr|S4RHT4|S4RHT4_PETMA 137 vkceegVKVLEALAASGLRSIRFAKEIPGIKSIITNDFSKNAVAAINKNIEMNQVQHLVTASYSDASMLMYkskNHKESFDVVDLDP 223 ***********************************************************************98778899******** PP SS--GGGHHHHHHHEEEEEEEEEEE--HHHHTTSSHHHHHHHHSSB---STTHHHHHHHHHHHHHHHHHHTTTEEEEEEEEEEETTE CS TRM 129 fGspapfldsavqsvkekglLlvTaTDtavLaGaspekclrkYgavalkselcHEvalRillkavareaakYekavepllslsadfY 215 +Gsp++fld+avqsv+e+glL+vT+TD+av+aG+s e+c++kYg+++lk ++cHE+alRi+l+ ++ + ++Y+++++pll++sadfY tr|S4RHT4|S4RHT4_PETMA 224 YGSPVQFLDAAVQSVNEGGLLCVTCTDMAVMAGNSGETCYSKYGSMSLKAKFCHEMALRIILHCLDLHSNRYHRYIVPLLAISADFY 310 *************************************************************************************** PP EEEEEEEEE-HHHHHHHHTTEEEEEEETTT.........--EEEE--GGG--........SB-TTT-SBEEEEEEEE-S-S--HHHH CS TRM 216 vrvfvkvktgaakvkkvleklglvlhCsgC.........lseekvkglakfk........aecehcgaklhlaGPlwlgplkdkefv 285 +rvfv++ +g+akvk +++k++lv++C+gC s+ +++gl+ +ce+cg+++hl+GP+w g ++d ++v tr|S4RHT4|S4RHT4_PETMA 311 IRVFVRIFSGQAKVKMSASKQALVYNCVGCgsfhlqpmgKSTAHNSGLKF-TaatgppvgPNCENCGQRHHLGGPIWTGAIHDINYV 396 ******************************99998886533344444433.247899999*************************** PP HHHHHHHCT...-TTS-HHHHHHHHHHHHHHS--SS-EEHHHHHHHHT-SS---HHHHHHH....TT-EE-TTSTTEEE-SS-HHHH CS TRM 286 eevleeaekek.eeelkkrvkklLklikeElldaplyydlsqlasvlklsvpplkkvvsalkelgfevsrtHakpkaikTnapleav 371 ++++ ++e+++ + +++ r+ ++L++++eEl d plyy l+ql+sv+++++p+l +++sa++++g++vs++Ha+++aikT+ap +++ tr|S4RHT4|S4RHT4_PETMA 397 QKIISAVEANRtRFKTSDRIIGMLAVVTEELQDCPLYYGLDQLSSVIRCNTPSLLEFRSAVLHAGYRVSLSHACKNAIKTDAPPHVI 483 ********975344459********************************************************************** PP HHHHH CS TRM 372 vevlr 376 ++++r tr|S4RHT4|S4RHT4_PETMA 484 WDIMR 488 ****7 PP >> tr|S4RMN3|S4RMN3_PETMA tRNA (guanine(26)-N(2))-dimethyltransferase OS=Petromyzon marinus OX=7757 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 189.9 0.1 8.6e-60 5.8e-56 12 331 .. 163 510 .. 148 539 .. 0.88 2 ! 33.4 0.0 3e-12 2e-08 317 375 .. 558 615 .. 553 616 .. 0.89 Alignments for each domain: == domain 1 score: 189.9 bits; conditional E-value: 8.6e-60 SS--TT--SS--GGGH.HHHHHHHHH...HHHHHHHH-S-EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHT CS TRM 12 vtvskkeeVFYNPvqe.fnRDlsvlvirqlellkkklgekikvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvkln 97 v++k++ F+NP+++ nR l+ + l++l+++ +++ lda+ a+G+ ++ +a+ + ++ +v +ND+se v +++eN +ln tr|S4RMN3|S4RMN3_PETMA 163 HNVKRKSDSFFNPRMKkPNRQLVFCA---LAALTQE-RKPLACLDAFGATGIMGLQWAKHLGNTVRVTINDVSEASVRMMEENCRLN 245 3489999*******9637****9999...6666665.589********************************************999 PP T--....................G.GEEEE-S-HHHHHC--SS-EEEEEE--SS--GGGHHHHHHHEEEEEEEEEEE--HHHHTTSS CS TRM 98 eve....................e.ivvineddAnvllaekkkrfdviDlDPfGspapfldsavqsvkekglLlvTaTDtavLaGas 163 + + +v++ + dAnvl+ + ++fd+i lDPfG+++++ld+a + v+++g++++T+TD+++L + tr|S4RMN3|S4RMN3_PETMA 246 GMTvteswpdsdeggeevegeaeDrLVEVTQYDANVLM--HLRAFDYIHLDPFGTAVDYLDAAFRNVRNHGIVSITSTDLSSLYCKV 330 9988899999999**9999*998469************..889******************************************** PP HHHHHHHHSSB---STTHHHHHHHHHHHH.....HHHHHHTTTEEEEEEEEEEETTEEEEEEEEEE-HHHHHHHHTTEEEEEEETTT CS TRM 164 pekclrkYgavalkselcHEvalRillka.....vareaakYekavepllslsadfYvrvfvkvktgaakvkkvleklglvlhCsgC 245 p+ lr+Y+ + +++e++HE a R++l+a + r+a + +k ve+ll+++ +++v v v+v +g+ ++++ +++l ++hC C tr|S4RMN3|S4RMN3_PETMA 331 PNVALRHYHCNLVRTEYYHELAARMVLTAavsacLCRAAGRCNKGVEVLLAVALEHFVLVVVRVLRGPGHADECARRLRRLVHCQWC 417 **************************98622222567899*********************************************** PP ..........--EEEE--GGG--SB-TTT-SBEEEEEEEE-S-S--HHHHHHHHHHHCT..-TTS-HHHHHHHHHHHHHH.S--SS- CS TRM 246 ..........lseekvkglakfkaecehcgaklhlaGPlwlgplkdkefveevleeaekekeeelkkrvkklLklikeEl.ldaply 321 + +e+ +++ + + + g + GP+w g l ++ f+ ++l ea + + ++v+ lLk++ E+ ++ tr|S4RMN3|S4RMN3_PETMA 418 eervfapdshMLSENPYSVLPCDCHNSMPGRTAIELGPMWSGSLFNSGFLRRMLFEAVQ----HEVEDVQPLLKTLAGEAeCTNVRK 500 877776665533366666666667777778888899******************99988....346899*******99998899999 PP EEHHHHHHHH CS TRM 322 ydlsqlasvl 331 y l ++a tr|S4RMN3|S4RMN3_PETMA 501 YALTNVAAHT 510 9999998754 PP == domain 2 score: 33.4 bits; conditional E-value: 3e-12 --SS-EEHHHHHHHHT-SS---HHHHHHH....TT-EE-TTSTTEEE-SS-HHHHHHHH CS TRM 317 daplyydlsqlasvlklsvpplkkvvsalkelgfevsrtHakpkaikTnapleavvevl 375 + p+yy++++ s+ +++p+l+++ + l e+gf++srtH+ p++i+T+a+l ++ vl tr|S4RMN3|S4RMN3_PETMA 558 HPPFYYSIHRH-SIKGMNMPKLNHFLASLTEAGFRASRTHFDPTSIRTDATLPQINAVL 615 57899999985.67899**********************************99887765 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (376 nodes) Target sequences: 13539 (5048892 residues searched) Passed MSV filter: 345 (0.0254819); expected 270.8 (0.02) Passed bias filter: 280 (0.020681); expected 270.8 (0.02) Passed Vit filter: 21 (0.00155107); expected 13.5 (0.001) Passed Fwd filter: 2 (0.000147721); expected 0.1 (1e-05) Initial search space (Z): 13539 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.14u 0.01s 00:00:00.15 Elapsed: 00:00:00.04 # Mc/sec: 47459.58 // [ok]