# 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/TrmO.hmm.txt # target sequence database: proteomes/Dictyostelium_discoideum.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: TrmO [M=119] Accession: PF01980.16 Description: tRNA-methyltransferase O 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.1e-39 134.9 0.0 2.2e-39 133.9 0.0 1.5 1 tr|Q54CD3|Q54CD3_DICDI Uncharacterized protein OS=Dictyostel 1.6e-38 131.1 0.0 8.6e-26 90.0 0.1 2.5 2 tr|Q54N07|Q54N07_DICDI Uncharacterized protein OS=Dictyostel Domain annotation for each sequence (and alignments): >> tr|Q54CD3|Q54CD3_DICDI Uncharacterized protein OS=Dictyostelium discoideum OX=44689 GN=DDB0191740 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 133.9 0.0 3.5e-43 2.2e-39 1 115 [. 112 236 .. 112 239 .. 0.93 Alignments for each domain: == domain 1 score: 133.9 bits; conditional E-value: 3.5e-43 TS-SSGGG--.-EEEEEE-HHH..HGGGTTGGG-SEEEEEEE-TTS..............-.SS--E-TT-.SS-E.EGGGS--S-S CS TrmO 1 giPrqpglaeeaeaeielepeynseealegleefshlwllfvfhen...........aekklkakvrpprlggnervGvfatRsphR 76 g+Prq+ ++++ +a ++l++++n++ +legl+++sh+wl+f+fhen ++++kakv pprl+g +vG+f+tR+phR tr|Q54CD3|Q54CD3_DICDI 112 GTPRQGLIVTKGKALLKLKTSVNPTYSLEGLDQYSHCWLIFIFHENtntlrsvrsdrPASGVKAKVAPPRLDG-IKVGLFSTRTPHR 197 79******************************************997777776655544447**********9.899********** PP SS-EEEEEEEEEEEETTEEEEES----TT-EEEEEEE-- CS TrmO 77 pnpiglsvvklekvegnvlevsgvDlldgtpvlDiKPYv 115 pnpiglsv+k++kv g+++++sg+Dl+dgtp+lD+KPY tr|Q54CD3|Q54CD3_DICDI 198 PNPIGLSVAKIDKVVGDTIYLSGIDLIDGTPILDVKPYS 236 **************************************6 PP >> tr|Q54N07|Q54N07_DICDI Uncharacterized protein OS=Dictyostelium discoideum OX=44689 GN=DDB0186581 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 38.6 0.0 1.2e-13 7.4e-10 1 56 [. 41 95 .. 41 114 .. 0.86 2 ! 90.0 0.1 1.3e-29 8.6e-26 51 119 .] 159 225 .. 128 225 .. 0.88 Alignments for each domain: == domain 1 score: 38.6 bits; conditional E-value: 1.2e-13 TS-SSGG.G--.-EEEEEE-HHH..HGGGTTGGG-SEEEEEEE-TTS...-.SS--E CS TrmO 1 giPrqpg.laeeaeaeielepeynseealegleefshlwllfvfhenaekklkakvr 56 g+Prq+ l+++++++i+++ + ++al+gle++sh+w+lf+fh+n ++++k+k++ tr|Q54N07|Q54N07_DICDI 41 GTPRQSYkLVPDSSGKIKIDCKH-ADDALRGLEGYSHVWVLFWFHKN-NEDYKKKIK 95 79****769***********998.89*******************99.444466665 PP == domain 2 score: 90.0 bits; conditional E-value: 1.3e-29 .SS--E-TT-.SS-E.EGGGS--S-SSS-EEEEEEEEEEEETTEEEEES----TT-EEEEEEE--HHHH CS TrmO 51 lkakvrpprlggnervGvfatRsphRpnpiglsvvklekvegnvlevsgvDlldgtpvlDiKPYvpyaD 119 + +v+pprlg +r+Gvf+tRsphR n igls++k+ek+e +++++sgvD+++gtp++DiKPY+p++D tr|Q54N07|Q54N07_DICDI 159 RHVTVKPPRLG--KRIGVFSTRSPHRYNDIGLSLCKIEKIEYDTIYLSGVDIINGTPIIDIKPYIPQYD 225 467999****9..6*****************************************************99 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (119 nodes) Target sequences: 12739 (6856946 residues searched) Passed MSV filter: 324 (0.0254337); expected 254.8 (0.02) Passed bias filter: 224 (0.0175838); expected 254.8 (0.02) Passed Vit filter: 11 (0.00086349); expected 12.7 (0.001) Passed Fwd filter: 2 (0.000156998); expected 0.1 (1e-05) Initial search space (Z): 12739 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.07u 0.00s 00:00:00.07 Elapsed: 00:00:00.04 # Mc/sec: 20399.41 // [ok]