# 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/Monodelphis_domestica.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 6.1e-146 486.7 0.0 7.7e-146 486.3 0.0 1.1 1 tr|F6UGF2|F6UGF2_MONDO tRNA (guanine(26)-N(2))-dimethyltrans 1e-66 225.9 2.0 5.5e-55 187.3 0.2 3.0 3 tr|F7DB20|F7DB20_MONDO tRNA methyltransferase 1 like OS=Mono 0.0074 15.7 0.0 0.013 15.0 0.0 1.3 1 tr|F6PG51|F6PG51_MONDO tRNA (guanine(37)-N1)-methyltransfera ------ inclusion threshold ------ 0.054 12.9 5.4 0.021 14.3 0.9 2.1 2 tr|F7E330|F7E330_MONDO Adenylate kinase 9 OS=Monodelphis dom Domain annotation for each sequence (and alignments): >> tr|F6UGF2|F6UGF2_MONDO tRNA (guanine(26)-N(2))-dimethyltransferase OS=Monodelphis domestica OX=13616 GN=TRMT1 PE=3 S # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 486.3 0.0 1.4e-149 7.7e-146 3 376 .] 34 472 .. 32 472 .. 0.92 Alignments for each domain: == domain 1 score: 486.3 bits; conditional E-value: 1.4e-149 EEEEE----SS--TT--SS--GGGHHHHHHHHHH.......HHHHHHHH-S-................................... CS TRM 3 akikvedkevtvskkeeVFYNPvqefnRDlsvlvi....rqlellkkklgek................................... 50 + i + +++++++eVFYNPvqefnRDl+++vi r ++ll+k + + tr|F6UGF2|F6UGF2_MONDO 34 TVISEGAVRIVFPNSNEVFYNPVQEFNRDLTCAVItefaR-MQLLSKGI--Qievpgeeklqkvvvnlaeseegsmepcprdqpqta 117 55555555689***********************943322.23333333..14477******************************* PP .........EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHHHHC...--SS-EEEEE CS TRM 51 .........ikvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneveeivvineddAnvlla...ekkkrfdviD 125 ++vl++L+asGlR+iRfa+Evpg + v+aND+s +Av+li+ Nv+ln+++++v++n +dA++l++ + ++rfdviD tr|F6UGF2|F6UGF2_MONDO 118 avgnvceegLRVLEGLAASGLRSIRFAKEVPGLQVVIANDASARAVDLIRCNVQLNNISHLVQPNLADARLLMYehqRGSRRFDVID 204 **************************************************************************66556789***** PP E--SS--GGGHHHHHHHEEEEEEEEEEE--HHHHTTSSHHHHHHHHSSB---STTHHHHHHHHHHHHHHHHHHTTTEEEEEEEEEEE CS TRM 126 lDPfGspapfldsavqsvkekglLlvTaTDtavLaGaspekclrkYgavalkselcHEvalRillkavareaakYekavepllslsa 212 lDP+Gspa fld+avq+v+e+glL+vT+TD+av+aG+s e+c++kYga+++ks++cHE+alRi+l++++ +a++Y+++++plls+sa tr|F6UGF2|F6UGF2_MONDO 205 LDPYGSPAIFLDAAVQAVSEGGLLCVTCTDMAVMAGNSSETCYSKYGAMSIKSKACHEMALRIVLHSLDLRANCYQRYIVPLLSISA 291 *************************************************************************************** PP TTEEEEEEEEEE-HHHHHHHHTTEEEEEEETTT........--E.EEE--GGG--........SB-TTT-SBEEEEEEEE-S-S--H CS TRM 213 dfYvrvfvkvktgaakvkkvleklglvlhCsgC........lse.ekvkglakfk........aecehcgaklhlaGPlwlgplkdk 282 dfY+rvfv+v tg+akvk++++k+g+v++C+gC + + g+ kf+ ecehcg++++l+GPlw++p++d tr|F6UGF2|F6UGF2_MONDO 292 DFYIRVFVRVYTGQAKVKASASKQGFVFQCVGCgsfhlqrlG-KaVNSGGRDKFSaatgppvgSECEHCGQRHQLGGPLWAEPIHDL 377 *********************************999977742.22445567777788899999************************ PP HHHHHHHHHHCT...-TTS-HHHHHHHHHHHHHHS--SS-EEHHHHHHHHT-SS---HHHHHHH....TT-EE-TTSTTEEE-SS-H CS TRM 283 efveevleeaekek.eeelkkrvkklLklikeElldaplyydlsqlasvlklsvpplkkvvsalkelgfevsrtHakpkaikTnapl 368 +fv++vl++++ + + +++ r++++L++++eEl+d+plyy+l+ql+s++++++p+l +++sal+++gf+vs++Ha+++a+kT+ap tr|F6UGF2|F6UGF2_MONDO 378 DFVGQVLAAVSRNPgRFRTSARIRGVLSMVTEELPDVPLYYTLDQLSSTIRCNTPRLLQLRSALLHAGFRVSLSHACKNAVKTDAPP 464 ***********97645555******************************************************************** PP HHHHHHHH CS TRM 369 eavvevlr 376 +++++++r tr|F6UGF2|F6UGF2_MONDO 465 SVLWDIMR 472 *******7 PP >> tr|F7DB20|F7DB20_MONDO tRNA methyltransferase 1 like OS=Monodelphis domestica OX=13616 GN=TRMT1L PE=3 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -1.1 0.0 0.18 9.5e+02 229 251 .. 286 312 .. 281 332 .. 0.64 2 ! 187.3 0.2 1e-58 5.5e-55 14 299 .. 357 668 .. 339 702 .. 0.85 3 ! 36.4 0.0 7.2e-13 3.8e-09 302 375 .. 730 800 .. 716 801 .. 0.82 Alignments for each domain: == domain 1 score: -1.1 bits; conditional E-value: 0.18 HHHHHTTEEEEEEETTT....--EEEE CS TRM 229 vkkvleklglvlhCsgC....lseekv 251 ++++ +k+g +hC+ C +++ + tr|F7DB20|F7DB20_MONDO 286 GDQTSSKMGAHYHCIICsatiTRRTDM 312 58999************6554222222 PP == domain 2 score: 187.3 bits; conditional E-value: 1e-58 --TT--SS--GGGHHHHHHHHHH...HHHHHHHH-S-EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT-- CS TRM 14 vskkeeVFYNPvqefnRDlsvlvirqlellkkklgekikvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneve 100 ++k++ ++NP ++ nR l + l++l ++ ++++ lda+ a+G+ ++ +a+ + + +v +ND++e+ v +i+eN +ln+++ tr|F7DB20|F7DB20_MONDO 357 TPQKTDSYFNPKMKLNRQLIFST---LAVLAEE-RKPLECLDAFGATGIMGLQWAKHLGNSVKVTINDFNENSVTMIQENCHLNKMK 439 6788999************9999...7777777.689*******************************************9999999 PP ......................GGEEEE-S-HHHHHC--SS-EEEEEE--SS--GGGHHHHHHHEEEEEEEEEEE--HHHHTTSSHH CS TRM 101 ......................eivvineddAnvllaekkkrfdviDlDPfGspapfldsavqsvkekglLlvTaTDtavLaGaspe 165 +++ ++dAnvl+ + + fd+i lDPfG+++++ldsa + v++ g+++vT+TD ++L ++ + tr|F7DB20|F7DB20_MONDO 440 vmldskeedecdetleegeenvGAIEVTKMDANVLM--HLRSFDFIHLDPFGTSVNYLDSAFRNVRNLGIVSVTSTDISSLYAKALH 524 999999999999999999875547999*********..7899********************************************* PP HHHHHHSSB---STTHHHHHHHHHHHHHHHHHHTTTEEEEEEEEEEETTEEEEEEEEEE-HHHHHHHHTTEEEEEEETTT....... CS TRM 166 kclrkYgavalkselcHEvalRillkavareaakYekavepllslsadfYvrvfvkvktgaakvkkvleklglvlhCsgC....... 245 r+Yg + +++e++ E a Ri++ avar+a++ +k +e+l++++ +++v v v+v +g++ ++++++k+ +++hC C tr|F7DB20|F7DB20_MONDO 525 VARRHYGCNIVRTEYYKELAARIVVAAVARAATRCNKGIEVLFAVALEHFVLVVVRVLRGPTPADETAKKIRYLIHCQWCeerifqk 611 ********************************************************************************6554432 PP ...--EEEE--GGG--........SB-TTT-SBEEEEEEEE-S-S--HHHHHHHHHHHCT..-TT CS TRM 246 ...lseekvkglakfk........aecehcgaklhlaGPlwlgplkdkefveevleeaekekeee 299 + + + + + g++ GPlw g l ++ f+ ++l ea ++ + tr|F7DB20|F7DB20_MONDO 612 dgnM-----IEE---NpyrqlpcdCHGSMPGKTAIELGPLWSGSLFNTGFLRRMLFEAVQHGLDD 668 2221.....222...2333544224444556667778*****************98887744444 PP == domain 3 score: 36.4 bits; conditional E-value: 7.2e-13 HHHHHHHHHHHHHHS--SS-EEHHHHHHHHT-SS---HHHHHHH....TT-EE-TTSTTEEE-SS-HHHHHHHH CS TRM 302 krvkklLklikeElldaplyydlsqlasvlklsvpplkkvvsalkelgfevsrtHakpkaikTnapleavvevl 375 + +kk + +++E + p+yy++++ s+ +++p+l+k+ + l ++gf+vsrtH+ p +++T+apl+++ +l tr|F7DB20|F7DB20_MONDO 730 NATKKQKADVNAE--HPPFYYNIHRH-SIKGMNMPKLNKFLQYLSQAGFRVSRTHFDPMGVRTDAPLSQFKAIL 800 5555555566667..88999999985.6789************************************9998776 PP >> tr|F6PG51|F6PG51_MONDO tRNA (guanine(37)-N1)-methyltransferase OS=Monodelphis domestica OX=13616 GN=TRMT5 PE=3 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 15.0 0.0 2.4e-06 0.013 63 115 .. 314 363 .. 262 372 .. 0.82 Alignments for each domain: == domain 1 score: 15.0 bits; conditional E-value: 2.4e-06 HHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHHHHC CS TRM 63 RaiRfalEvpgvkevlaNDisekAvelikeNvklneveeivvineddAnvlla 115 ai a+ +v+aND++++ ++ + +N kln+v++ v+i + d + +l+ tr|F6PG51|F6PG51_MONDO 314 FAIPVAKK---DCTVFANDLNPESYKWLLHNCKLNKVDQKVKIFNLDGKDFLQ 363 55555544...5689***********************999998888887775 PP >> tr|F7E330|F7E330_MONDO Adenylate kinase 9 OS=Monodelphis domestica OX=13616 GN=AK9 PE=4 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 14.3 0.9 3.9e-06 0.021 262 347 .. 157 243 .. 125 264 .. 0.86 2 ? -0.3 0.2 0.1 5.4e+02 259 315 .. 984 1042 .. 963 1053 .. 0.74 Alignments for each domain: == domain 1 score: 14.3 bits; conditional E-value: 3.9e-06 TTT-SB.EEEEEEEE-S-S--HHHHHHHHHHHCT..-TTS-HHHHHHHHHHHHHHS--SS-EEHHHHHHHHT-SS---HHHHHHH.. CS TRM 262 ehcgak.lhlaGPlwlgplkdkefveevleeaekekeeelkkrvkklLklikeElldaplyydlsqlasvlklsvpplkkvvsalke 347 cg++ ++l+G +++ + +d ef+e+ ++ ++ ++e+ kk +++ + ++eE +++++ +++ +a++ + v++ +++++++++ tr|F7E330|F7E330_MONDO 157 RLCGQRqHSLTGFVYMREQWDPEFLESRRKKKKEFHKEKDKKAEEEVEEEVEEEEEEQVFIAEMQMMAEIFRHLVQKPEDFMENIEA 243 56776516799********************999888888999999999*****9***************************99975 PP == domain 2 score: -0.3 bits; conditional E-value: 0.1 SB-TTT-SBEEEEEEEE-S-S--HHHHHHHHHHHCT..-TTS...-HHHHHHHHHHHHHH CS TRM 259 aecehcgaklhlaGPlwlgplkdkefveevleeaekekeeel...kkrvkklLklikeEl 315 +e++ +++ +GP + ++++d+ + lee+ ++ + +l kkr +k L+l +eE tr|F7E330|F7E330_MONDO 984 QEKIMVKT-ERKVGPEYDEEIEDETAARQELEEIAAQANVTLdeeKKRPRKELALTDEEE 1042 44444444.4567************99998888888766655666889999999999885 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (376 nodes) Target sequences: 21272 (11266248 residues searched) Passed MSV filter: 728 (0.0342234); expected 425.4 (0.02) Passed bias filter: 503 (0.0236461); expected 425.4 (0.02) Passed Vit filter: 44 (0.00206845); expected 21.3 (0.001) Passed Fwd filter: 4 (0.000188041); expected 0.2 (1e-05) Initial search space (Z): 21272 [actual number of targets] Domain search space (domZ): 4 [number of targets reported over threshold] # CPU time: 0.31u 0.01s 00:00:00.32 Elapsed: 00:00:00.09 # Mc/sec: 47067.88 // [ok]