# 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/Ciona_savignyi.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.1e-130 433.4 0.0 6.7e-130 433.0 0.0 1.0 1 tr|H2ZLY0|H2ZLY0_CIOSA tRNA (guanine(26)-N(2))-dimethyltrans ------ inclusion threshold ------ 0.048 12.2 0.0 0.065 11.8 0.0 1.2 1 tr|H2Y5Z4|H2Y5Z4_CIOSA tRNA (guanine(37)-N1)-methyltransfera Domain annotation for each sequence (and alignments): >> tr|H2ZLY0|H2ZLY0_CIOSA tRNA (guanine(26)-N(2))-dimethyltransferase OS=Ciona savignyi OX=51511 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 433.0 0.0 1.1e-133 6.7e-130 7 376 .] 13 436 .. 7 436 .. 0.94 Alignments for each domain: == domain 1 score: 433.0 bits; conditional E-value: 1.1e-133 E----SS--TT--SS--GGGHHHHHHHHHH....HHHHHHHH-S-.................................EEEEESS-T CS TRM 7 vedkevtvskkeeVFYNPvqefnRDlsvlvir.qlellkkklgek.................................ikvldaLsa 59 + ++++ +++++ VFYNPvqefnRD++++vi +++++ kk + + +++l+aLsa tr|H2ZLY0|H2ZLY0_CIOSA 13 EGKASILFPSENSVFYNPVQEFNRDVTIAVITeFIKAKLKKDKINavefprndssvislnddncsdvtiypgeklengVRILEALSA 99 3444578***********************98788888777666599**************************************** PP TSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHHHHC...--SS-EEEEEE--SS--GGGHHHHHHHE CS TRM 60 sGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneveeivvineddAnvlla...ekkkrfdviDlDPfGspapfldsavqsv 143 +GlR++Rfa+E+pgvk+++aND+se A +lik+N++ n+v ++v++n++dA++ ++ k+k+fdviD+DP+Gs+a fldsavq+v tr|H2ZLY0|H2ZLY0_CIOSA 100 TGLRSVRFAKELPGVKQIVANDLSEDAFHLIKQNIDHNNVCDVVEANQGDATTYMYinrGKSKEFDVIDIDPYGSAACFLDSAVQAV 186 ********************************************************998899************************* PP EEEEEEEEEE--HHHHTTSSHHHHHHHHSSB---STTHHHHHHHHHHHHHHHHHHTTTEEEEEEEEEEETTEEEEEEEEEE-HHHHH CS TRM 144 kekglLlvTaTDtavLaGaspekclrkYgavalkselcHEvalRillkavareaakYekavepllslsadfYvrvfvkvktgaakvk 230 k +glL+vT+TD+av+aG+ p +c +kY++v+lk+++cHE+alR++l+a+ ++aa+Y+++++pllsls+dfYvr+fv+v+ ++ ++ tr|H2ZLY0|H2ZLY0_CIOSA 187 KPGGLLCVTCTDMAVMAGNFPGTCWAKYQSVPLKTKACHEQALRMVLSAISSHAARYSRYIVPLLSLSVDFYVRIFVQVHYKPSIAQ 273 *************************************************************************************** PP HHHTTEEEEEEETTT...........--EEEE--GGG--..............SB-TTT-SBEEEEEEEE-S-S--HHHHHHHHHHH CS TRM 231 kvleklglvlhCsgC...........lseekvkglakfk..............aecehcgaklhlaGPlwlgplkdkefveevleea 292 ++ +k+ ++ C+gC + ++c +c+++++++GP w ++d +fv++vl + tr|H2ZLY0|H2ZLY0_CIOSA 274 QTSSKQIKLYRCVGCssfhsqhlvdrK---------GTDknptyfasvgppvgERCLECNHRHKIGGPYWGSAIHDLDFVNKVLTSV 351 ***************876554443221.........2113345555678***9********************************** PP CT..-TTS..-HHHHHHHHHHHHHHS--SS-EEHHHHHHHHT-SS---HHHHHHH....TT-EE-TTSTTEEE-SS-HHHHHHHHH CS TRM 293 ekekeeel..kkrvkklLklikeElldaplyydlsqlasvlklsvpplkkvvsalkelgfevsrtHakpkaikTnapleavvevlr 376 ++ ++++ ++r+ ++L+++ eEl+d p+yy++++l+ +l++ +pp+k ++sa++++g++vs+ H + ++ikT+ap ++++++r tr|H2ZLY0|H2ZLY0_CIOSA 352 KS-FPTRFgtAERIIGMLSMVVEELPDIPFYYVTDELSMTLRCGSPPSKLIRSAILNAGYRVSGFHGRVNSIKTDAPPTVIWDIMR 436 **.6788778***************************************************************************7 PP >> tr|H2Y5Z4|H2Y5Z4_CIOSA tRNA (guanine(37)-N1)-methyltransferase OS=Ciona savignyi OX=51511 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.8 0.0 1.1e-05 0.065 50 113 .. 254 312 .. 187 333 .. 0.75 Alignments for each domain: == domain 1 score: 11.8 bits; conditional E-value: 1.1e-05 -EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--.GGEEEE-S-HHHH CS TRM 50 kikvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneve.eivvineddAnvl 113 + +v+d +++ G ai a+ v+aND+++ ++ + +Nvkln+ e e + d + + tr|H2Y5Z4|H2Y5Z4_CIOSA 254 DATVFDVFAGVGPFAIPIAKK---GCLVYANDLNPDSYKWLVHNVKLNKTEvECF---NCDGRDF 312 378888899999888888877...5679********************9985432...2344444 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (376 nodes) Target sequences: 11592 (4231397 residues searched) Passed MSV filter: 375 (0.0323499); expected 231.8 (0.02) Passed bias filter: 219 (0.0188923); expected 231.8 (0.02) Passed Vit filter: 16 (0.00138026); expected 11.6 (0.001) Passed Fwd filter: 2 (0.000172533); expected 0.1 (1e-05) Initial search space (Z): 11592 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.11u 0.00s 00:00:00.11 Elapsed: 00:00:00.04 # Mc/sec: 39775.13 // [ok]