# 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/Methyltransf_15.hmm.txt # target sequence database: proteomes/Drosophila_melanogaster.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: Methyltransf_15 [M=165] Accession: PF09445.9 Description: RNA cap guanine-N2 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.7e-44 151.1 0.0 3.6e-44 150.0 0.0 1.6 1 tr|Q9VEA4|Q9VEA4_DROME Trimethylguanosine synthase 1 OS=Dros 2.3e-05 23.8 0.0 4.5e-05 22.8 0.0 1.5 1 tr|Q9VVV7|Q9VVV7_DROME Uncharacterized protein OS=Drosophila ------ inclusion threshold ------ 0.023 14.0 0.0 0.048 13.0 0.0 1.5 1 sp|Q8IRE4|TRM5_DROME tRNA (guanine(37)-N1)-methyltransfera Domain annotation for each sequence (and alignments): >> 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 ! 150.0 0.0 7.9e-48 3.6e-44 2 164 .. 300 455 .. 299 456 .. 0.95 Alignments for each domain: == domain 1 score: 150.0 bits; conditional E-value: 7.9e-48 Methyltransf_15 2 evvlDlfgGaGGnviqfanvfekviaidinlsrlklaqhNaevYgvsdkvtlvlgdvfellaklkleeikvdlvflSpPwGGpsYkk 88 v++D+f+G+GGn+iqfan+ +via+di+ ++l++a+hNa +Ygv++k++++++d+ + a++kl+ +vflSpPwGGp+Y+k tr|Q9VEA4|Q9VEA4_DROME 300 DVIVDAFCGCGGNAIQFANTCGRVIAVDIDAEKLAMAKHNAGIYGVAHKIEFIHADFLQFAASTKLRPN---VVFLSPPWGGPDYQK 383 689**********************************************************99998877...*************** PP Methyltransf_15 89 kekfdlekkleplklekllkesvklskniilfLPrnsdlkqlaaltaevlenaskvkvvklkenglvkallaylge 164 + +fd+e+ l p+ ++l++ s +l + + +fLPrn+++kq+ al + ++++v + + ++ al+ay+g+ tr|Q9VEA4|Q9VEA4_DROME 384 QATFDIETGLLPVGASQLMQLSRSLASDVAFFLPRNANMKQVVAL----SGVGQQCEVEHNYLDTRMVALTAYYGK 455 ******************************************998....577788999999999999999999986 PP >> tr|Q9VVV7|Q9VVV7_DROME Uncharacterized protein OS=Drosophila melanogaster OX=7227 GN=cg3808 PE=1 SV=3 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 22.8 0.0 9.8e-09 4.5e-05 3 65 .. 418 479 .. 416 526 .. 0.84 Alignments for each domain: == domain 1 score: 22.8 bits; conditional E-value: 9.8e-09 Methyltransf_15 3 vvlDlfgGaGGnviqfanvfekviaidinlsrlklaqhNaevYgvsdkvtlvlgdvfellakl 65 ++lD+++G+G + +fa+ +kv++++i ++ +k a+ Nae+ g+++ ++ +g+ + ++++ tr|Q9VVV7|Q9VVV7_DROME 418 TMLDICCGTGTITLAFAKHCKKVMGVEIVPDAIKDAEFNAEANGIKN-AKFFTGNADDFIKSM 479 79*****************************************9964.456666666655544 PP >> sp|Q8IRE4|TRM5_DROME tRNA (guanine(37)-N1)-methyltransferase OS=Drosophila melanogaster OX=7227 GN=CG32281 PE=2 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.0 0.0 1e-05 0.048 3 60 .. 252 309 .. 250 341 .. 0.86 Alignments for each domain: == domain 1 score: 13.0 bits; conditional E-value: 1e-05 Methyltransf_15 3 vvlDlfgGaGGnviqfanvfekviaidinlsrlklaqhNaevYgvsdkvtlvlgdvfe 60 v+ D+f+G G +i a+ ++v a d n++ + qhNa+ ++++ d + sp|Q8IRE4|TRM5_DROME 252 VLYDVFAGVGPFSIPAAKKRCHVLANDLNPESFRWLQHNAKRNKCLPNIKMSNKDGRQ 309 567**************************************99999999887777665 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (165 nodes) Target sequences: 13786 (7394964 residues searched) Passed MSV filter: 375 (0.0272015); expected 275.7 (0.02) Passed bias filter: 300 (0.0217612); expected 275.7 (0.02) Passed Vit filter: 32 (0.0023212); expected 13.8 (0.001) Passed Fwd filter: 3 (0.000217612); expected 0.1 (1e-05) Initial search space (Z): 13786 [actual number of targets] Domain search space (domZ): 3 [number of targets reported over threshold] # CPU time: 0.10u 0.00s 00:00:00.10 Elapsed: 00:00:00.04 # Mc/sec: 30504.23 // [ok]