# 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/Drosophila_melanogaster.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 7.9e-44 148.3 0.0 2e-43 147.0 0.0 1.7 1 tr|A1Z759|A1Z759_DROME Uncharacterized protein, isoform A OS 1.5e-09 37.7 0.0 3e-09 36.8 0.0 1.5 1 tr|A1Z758|A1Z758_DROME TA01342p OS=Drosophila melanogaster O Domain annotation for each sequence (and alignments): >> tr|A1Z759|A1Z759_DROME Uncharacterized protein, isoform A OS=Drosophila melanogaster OX=7227 GN=CG12822 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 147.0 0.0 3e-47 2e-43 2 119 .] 116 233 .. 115 233 .. 0.95 Alignments for each domain: == domain 1 score: 147.0 bits; conditional E-value: 3e-47 S-SSGGG--.-EEEEEE-HHH....HGGGTTGGG-SEEEEEEE-TTS...-.SS--E-TT-.SS-E.EGGGS--S-SSS-EEEEEEE CS TrmO 2 iPrqpglaeeaeaeielepey..nseealegleefshlwllfvfhenaekklkakvrpprlggnervGvfatRsphRpnpiglsvvk 86 +Prq+ + +++++ i+l++ + n+e++legle+fshlwl+++fh+n+++ kakv pprlgg ervGvf+tRsphRp+pigls+v+ tr|A1Z759|A1Z759_DROME 116 VPRQSIVGSRLRGIIQLNDGVftNPEHSLEGLEDFSHLWLIYHFHRNNSHP-KAKVAPPRLGG-ERVGVFSTRSPHRPCPIGLSLVE 200 8*****************76545899*******************996666.9********99.999******************** PP EEEEETTEEEEES----TT-EEEEEEE--HHHH CS TrmO 87 lekvegnvlevsgvDlldgtpvlDiKPYvpyaD 119 +ek+e+ ++ + g+D++dgtpvlDiKPY+p++D tr|A1Z759|A1Z759_DROME 201 IEKIENATISFFGTDMVDGTPVLDIKPYIPHYD 233 *******************************99 PP >> tr|A1Z758|A1Z758_DROME TA01342p OS=Drosophila melanogaster OX=7227 GN=143303_at PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 36.8 0.0 4.3e-13 3e-09 2 52 .. 116 167 .. 115 181 .. 0.84 Alignments for each domain: == domain 1 score: 36.8 bits; conditional E-value: 4.3e-13 S-SSGGG--.-EEEEEE-HHH....HGGGTTGGG-SEEEEEEE-TTS...-.S CS TrmO 2 iPrqpglaeeaeaeielepey..nseealegleefshlwllfvfhenaekklk 52 +Prq+ + +++++ i+l++ + n+e++legle+fshlwl+++fh+n+++ k tr|A1Z758|A1Z758_DROME 116 VPRQSIVGSRLRGIIQLNDGVftNPEHSLEGLEDFSHLWLIYHFHRNNSHP-K 167 8*****************76545899*******************884443.3 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (119 nodes) Target sequences: 13786 (7394964 residues searched) Passed MSV filter: 230 (0.0166836); expected 275.7 (0.02) Passed bias filter: 199 (0.0144349); expected 275.7 (0.02) Passed Vit filter: 22 (0.00159582); 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.08u 0.00s 00:00:00.08 Elapsed: 00:00:00.04 # Mc/sec: 22000.02 // [ok]