# 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/Caenorhabditis_elegans.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1.4e-154 515.0 0.0 1.6e-154 514.8 0.0 1.0 1 sp|Q23270|TRM1_CAEEL tRNA (guanine(26)-N(2))-dimethyltransfe 1.1e-06 28.3 0.0 2.1e-06 27.3 0.0 1.4 1 sp|A8WHT1|TRM5_CAEEL tRNA (guanine(37)-N1)-methyltransferase Domain annotation for each sequence (and alignments): >> sp|Q23270|TRM1_CAEEL tRNA (guanine(26)-N(2))-dimethyltransferase OS=Caenorhabditis elegans OX=6239 GN=trm-1 PE=1 SV= # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 514.8 0.0 1.6e-158 1.6e-154 1 376 [] 9 442 .. 9 442 .. 0.95 Alignments for each domain: == domain 1 score: 514.8 bits; conditional E-value: 1.6e-158 TTEEEEE----.SS........--TT--SS--GGGHHHHHHHHHH...HHHHHHHH-S-........................EEEEES CS TRM 1 Geakikvedke.vt........vskkeeVFYNPvqefnRDlsvlvirqlellkkklgek........................ikvlda 56 G+++ik+edke vt v + +VFYNPvqefnRDl+v+v+rq++++++k++e+ i++lda sp|Q23270|TRM1_CAEEL 9 GDSAIKSEDKEeVTviqegqakVGFHGPVFYNPVQEFNRDLTVTVLRQFSADHQKWAEEqkqlkteeeppkkknklainedgkIRILDA 97 7788888888755588888889***********************************99****************************** PP S-TTSHHHHHHHHH-S-EEEEEEE-S-HHHHHHHHHHHHHTT--GGEEEE-S-HHHHHC...--SS-EEEEEE--SS--GGGHHHHHHH CS TRM 57 LsasGlRaiRfalEvpgvkevlaNDisekAvelikeNvklneveeivvineddAnvlla...ekkkrfdviDlDPfGspapfldsavqs 142 LsasGlRa+Rf++Evp+v ++aND+s++Av++i+eNvkln ve+iv+++ +dA ++++ ++krf+++DlDP+G++++fldsavq+ sp|Q23270|TRM1_CAEEL 98 LSASGLRALRFSKEVPNVGFIMANDFSDNAVASIQENVKLNGVEDIVEAHFGDAVMTMMehrGIDKRFHAVDLDPYGTASTFLDSAVQC 186 **********************************************************9777789************************ PP EEEEEEEEEEE--HHHHTTSSHHHHHHHHSSB---STTHHHHHHHHHHHHHHHHHHTTTEEEEEEEEEEETTEEEEEEEEEE-HHHHHH CS TRM 143 vkekglLlvTaTDtavLaGaspekclrkYgavalkselcHEvalRillkavareaakYekavepllslsadfYvrvfvkvktgaakvkk 231 v+++g+L+vT+TD+avL+G++pe+c++kY+av++++++cHEv+lRill+a++++a++Y++++epl+s+s+dfYvrvfv+v+tga+++k+ sp|Q23270|TRM1_CAEEL 187 VADRGILMVTCTDMAVLCGNTPEACYNKYDAVTTRMKCCHEVGLRILLRAIDSAANRYTRYIEPLVSISVDFYVRVFVRVHTGAFQAKQ 275 ***************************************************************************************** PP HHTTEEEEEEETTT.............--EEEE--GGG--.......SB-TTT-SBEEEEEEEE-S-S--HHHHHHHHHHHCT...-TT CS TRM 232 vleklglvlhCsgC.............lseekvkglakfk.......aecehcgaklhlaGPlwlgplkdkefveevleeaekek.eee 299 + +k+g+vl+CsgC ++++++++ +++ ++c+hc lh++GP++l+p+++k fv+++le+++++ +e+ sp|Q23270|TRM1_CAEEL 276 SGTKVGTVLVCSGChsmeplvllkrgeGNQQSKYSIPTVRhsisgpgNRCIHCLLPLHQIGPIYLAPIHSKPFVTSLLERLKSTPeAER 364 **************9999999977666555677777777799*****************************************984444 PP S..-HHHHHHHHHHHHHHS--SS-EEHHHHHHHHT-SS---HHHHHHH....TT-EE-TTSTTEEE-SS-HHHHHHHHH CS TRM 300 l..kkrvkklLklikeElldaplyydlsqlasvlklsvpplkkvvsalkelgfevsrtHakpkaikTnapleavvevlr 376 l ++r++++L++++eE ld++lyy+++q+a+v+k+svp++++v+sa++++gf+vs++H++p+aikTnap++ +++++r sp|Q23270|TRM1_CAEEL 365 LgtHGRLQGVLTMVNEE-LDDVLYYEHNQMANVVKVSVPKSQSVRSAILNAGFKVSGSHCNPRAIKTNAPMHLLWDIYR 442 488**************.9*********************************************************986 PP >> sp|A8WHT1|TRM5_CAEEL tRNA (guanine(37)-N1)-methyltransferase OS=Caenorhabditis elegans OX=6239 GN=C53A5.17 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 27.3 0.0 2.1e-10 2.1e-06 74 117 .. 267 311 .. 243 319 .. 0.89 Alignments for each domain: == domain 1 score: 27.3 bits; conditional E-value: 2.1e-10 EEEEEEE-S-HHHHHHHHHHHHHTT--G.GEEEE-S-HHHHHC-- CS TRM 74 vkevlaNDisekAvelikeNvklnevee.ivvineddAnvllaek 117 k+v+aND++++ v+ +k Nv ln+++e +++i+++dA+++++e+ sp|A8WHT1|TRM5_CAEEL 267 PKRVVANDLNPESVKWLKVNVGLNKIKEeRIEIHNMDAKMFIKEN 311 579***********************955************9754 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (376 nodes) Target sequences: 20008 (8229220 residues searched) Passed MSV filter: 706 (0.0352859); expected 400.2 (0.02) Passed bias filter: 406 (0.0202919); expected 400.2 (0.02) Passed Vit filter: 22 (0.00109956); expected 20.0 (0.001) Passed Fwd filter: 2 (9.996e-05); expected 0.2 (1e-05) Initial search space (Z): 20008 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.21u 0.01s 00:00:00.22 Elapsed: 00:00:00.08 # Mc/sec: 38677.33 // [ok]