# 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/Drosophila_melanogaster.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2.7e-132 441.1 0.0 3.5e-132 440.8 0.0 1.0 1 sp|Q9VK89|TRM1_DROME Probable tRNA (guanine(26)-N(2))-dime ------ inclusion threshold ------ 0.051 12.4 0.0 0.075 11.8 0.0 1.3 1 tr|Q9VEA4|Q9VEA4_DROME Trimethylguanosine synthase 1 OS=Dros Domain annotation for each sequence (and alignments): >> sp|Q9VK89|TRM1_DROME Probable tRNA (guanine(26)-N(2))-dimethyltransferase OS=Drosophila melanogaster OX=7227 GN=CG63 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 440.8 0.0 5e-136 3.5e-132 11 376 .] 24 452 .. 12 452 .. 0.92 Alignments for each domain: == domain 1 score: 440.8 bits; conditional E-value: 5e-136 -SS--TT--SS--GGGHHHHHHHHHH...HHHHHHHH-S-........................................EEEEESS-T CS TRM 11 evtvskkeeVFYNPvqefnRDlsvlvirqlellkkklgek........................................ikvldaLsa 59 + + + +VFYNPvqefnRDls+++ l++++++l+++ +++l+aL+a sp|Q9VK89|TRM1_DROME 24 NAEIVSGGNVFYNPVQEFNRDLSIAA---LNVYRQRLTKErsekalkkqrkkvkeqedekttpvpedppvyeagtryedgLRILEALAA 109 344667799*****************...6666666665566679999***************************************** PP TSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHHHHC...--SS-EEEEEE--SS--GGGHHHHHHHEEE CS TRM 60 sGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneveeivvineddAnvlla...ekkkrfdviDlDPfGspapfldsavqsvke 145 +GlR+iR+a+E++gv++++aND+s +Av++i+ N++ n+vee+++++++dA +l++ + +krfd++DlDP+G p fld a+q++ + sp|Q9VK89|TRM1_DROME 110 TGLRSIRYAQEIAGVRQIVANDLSRQAVASINTNIRHNKVEELIEPSHSDAMTLMYlstQPEKRFDAVDLDPYGCPNRFLDGAMQCLVD 198 ********************************************************998889*************************** PP EEEEEEEE--HHHHTTSSHHHHHHHHSSB---STTHHHHHHHHHHHHHHHHHHTTTEEEEEEEEEEETTEEEEEEEEEE-HHHHHHHHT CS TRM 146 kglLlvTaTDtavLaGaspekclrkYgavalkselcHEvalRillkavareaakYekavepllslsadfYvrvfvkvktgaakvkkvle 234 +glLlvTaTD+avLaG++pe+c+ kYg+v+l++++cHE+alRill+ ++++a++Y k++eplls+sadfY+r+fv+v g+a+ k +++ sp|Q9VK89|TRM1_DROME 199 GGLLLVTATDMAVLAGNAPEACYVKYGSVPLRMKCCHEMALRILLHCIESHANRYGKYIEPLLSISADFYIRIFVRVYVGQAQCKLSMS 287 ***************************************************************************************** PP TEEEEEEETTT..................--EEEE--GGG--..............SB-TTT-SBEEEEEEEE-S-S--HHHHHHHHHH CS TRM 235 klglvlhCsgC..................lseekvkglakfk..............aecehcgaklhlaGPlwlgplkdkefveevlee 291 k+ +++C+gC +cehcg+++hl+GP+w +p+++ efv+ +l + sp|Q9VK89|TRM1_DROME 288 KQSWIYQCTGCetftlqplgitkpnptagN-----------PqqlkfgiptgpavnSQCEHCGHRHHLGGPIWSAPIHNPEFVQDLLTA 365 ***********9888777776555444331...........0233455566667779******************************** PP HCT..-TTS..-HHHHHHHHHHHHHHS--SS-EEHHHHHHHHT-SS---HHHHHHH....TT-EE-TTSTTEEE-SS-HHHHHHHHH CS TRM 292 aekekeeel..kkrvkklLklikeElldaplyydlsqlasvlklsvpplkkvvsalkelgfevsrtHakpkaikTnapleavvevlr 376 +++++ ++l ++r+ ++L++++eEl d+plyy+ ++l+ vlkl+ p+ k++sa++++g++vs +Ha+++++kTnap +++++lr sp|Q9VK89|TRM1_DROME 366 VQETTLQSLgtQRRIVGVLSMVQEELQDVPLYYTPDKLCCVLKLEIVPMLKFRSAILHAGYRVSYSHASKNSLKTNAPPAVLWDILR 452 ***98666667**************************************************************************97 PP >> tr|Q9VEA4|Q9VEA4_DROME Trimethylguanosine synthase 1 OS=Drosophila melanogaster OX=7227 GN=Tgs1 PE=4 SV=3 # 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.075 50 156 .. 292 401 .. 225 410 .. 0.82 Alignments for each domain: == domain 1 score: 11.8 bits; conditional E-value: 1.1e-05 -.......EEEEESS-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHHHHC--SS-EEEEEE--. CS TRM 50 k.......ikvldaLsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneveeivvineddAnvllaekkkrfdviDlDP. 128 + ++da+ + G ai fa +v+a Di + +++ k+N+ + v++ +++ ++d + a++k r +v+ l P tr|Q9VEA4|Q9VEA4_DROME 292 QtarrlacDVIVDAFCGCGGNAIQFANT---CGRVIAVDIDAEKLAMAKHNAGIYGVAHKIEFIHADFLQFAASTKLRPNVVFLSPp 375 355555543589**************87...8999************************99999989999999999999***99853 PP SS--GGGHHHHHHHEEEEEEEEEEE--H CS TRM 129 fGspapfldsavqsvkekglLlvTaTDt 156 +G p ++ +a +++ glL v a + tr|Q9VEA4|Q9VEA4_DROME 376 WGGP-DYQKQATFDIET-GLLPVGASQL 401 7776.688888888877.8888877665 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (376 nodes) Target sequences: 13786 (7394964 residues searched) Passed MSV filter: 422 (0.0306108); expected 275.7 (0.02) Passed bias filter: 304 (0.0220514); expected 275.7 (0.02) Passed Vit filter: 21 (0.00152328); expected 13.8 (0.001) Passed Fwd filter: 2 (0.000145075); expected 0.1 (1e-05) Initial search space (Z): 13786 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.18u 0.01s 00:00:00.19 Elapsed: 00:00:00.06 # Mc/sec: 46341.77 // [ok]