# 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]