# 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/TRM13.hmm.txt # target sequence database: proteomes/Drosophila_melanogaster.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: TRM13 [M=262] Accession: PF05206.13 Description: Methyltransferase TRM13 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 4.4e-91 304.7 1.9 8.6e-91 303.7 1.9 1.5 1 tr|Q8SX24|Q8SX24_DROME RE67941p OS=Drosophila melanogaster O 0.0029 16.8 3.3 1.6 7.8 0.0 3.0 3 sp|Q9VJ34|GSTCD_DROME Glutathione S-transferase C-terminal Domain annotation for each sequence (and alignments): >> tr|Q8SX24|Q8SX24_DROME RE67941p OS=Drosophila melanogaster OX=7227 GN=Dmel\CG18048 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 303.7 1.9 1.2e-94 8.6e-91 2 261 .. 184 428 .. 183 429 .. 0.97 Alignments for each domain: == domain 1 score: 303.7 bits; conditional E-value: 1.2e-94 TRM13 2 qqsslignleelkllekssayvEfGaGkgelsryvnqaleedkkenskyvLidrasnrlkaDrklrkdeselaikRlkidIkdLnls 88 q++sl+g+le+ ++l +++ y+EfGaGkg+l+ +++++l+e+k ++s++vLidr+s r+k+D+kl+++e +++R+++dI dL+ls tr|Q8SX24|Q8SX24_DROME 184 QTTSLLGILEHDHQLMDHTSYIEFGAGKGQLAYFLATVLQEQKLSHSQVVLIDRMSLRHKKDNKLANRE---VVQRIRADIADLKLS 267 899******************************************************************...8************** PP TRM13 89 aleelekkkkvvavsKHLCGaATDLtLrcllnsekaskkakleglliAlCChhvCswkeyvnkeyleelgitkeefqiltkmvSWav 175 al+el+k++++va sKHLCGaATDLtLrc+l +as +++liAlCChh+Csw++yv++++l+e gi +ef iltkmvSWav tr|Q8SX24|Q8SX24_DROME 268 ALPELKKTQRTVAFSKHLCGAATDLTLRCILGDGNAS----SDYVLIALCCHHRCSWRSYVGRKFLQEAGIGPREFVILTKMVSWAV 350 *********************************9986....59******************************************** PP TRM13 176 sgkrkseaeeedekeekekeekkeeeeseeekksklsseereeiGlkakrlidegRllalkekgfeaelvkYvekevslEnvlLla 261 +g+ ++e++k++++++ + +e+++++l+++ere+iG+++kr++d+gRl++l+++g++ael++Yv+++v+lEnv+Lla tr|Q8SX24|Q8SX24_DROME 351 CGTG--------LSRERRKAMESADFQLTETNTQRLTRQEREQIGQQCKRVLDYGRLEHLRSHGYQAELKFYVPRDVTLENVVLLA 428 ****........56788899999999**********************************************************98 PP >> sp|Q9VJ34|GSTCD_DROME Glutathione S-transferase C-terminal domain-containing protein homolog OS=Drosophila melanogas # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -0.4 0.0 0.076 5.2e+02 16 42 .. 384 410 .. 368 433 .. 0.70 2 ! 5.7 2.3 0.0011 7.3 100 122 .. 454 476 .. 449 489 .. 0.84 3 ! 7.8 0.0 0.00023 1.6 218 260 .. 536 578 .. 511 580 .. 0.89 Alignments for each domain: == domain 1 score: -0.4 bits; conditional E-value: 0.076 TRM13 16 lekssayvEfGaGkgelsryvnqalee 42 ++ +v+f++G g+l+ ++ l + sp|Q9VJ34|GSTCD_DROME 384 AQPGDRIVDFCSGTGHLAILLALKLPN 410 3677899**********9988876644 PP == domain 2 score: 5.7 bits; conditional E-value: 0.0011 TRM13 100 vavsKHLCGaATDLtLrcllnse 122 ++ s H CG+ATD+ L+ + + + sp|Q9VJ34|GSTCD_DROME 454 IGASLHACGTATDIVLQQCRRAK 476 57799***********9887665 PP == domain 3 score: 7.8 bits; conditional E-value: 0.00023 TRM13 218 eiGlkakrlidegRllalkekgfeaelvkYvekevslEnvlLl 260 Gl + ++d+ Rll +e g+++ l+ +++ ++ n+lL+ sp|Q9VJ34|GSTCD_DROME 536 LQGLHCMSVVDTDRLLQAEEAGYQVILTRLKPEQCTPKNHLLV 578 4699************************************997 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (262 nodes) Target sequences: 13786 (7394964 residues searched) Passed MSV filter: 1357 (0.0984332); expected 275.7 (0.02) Passed bias filter: 511 (0.0370666); expected 275.7 (0.02) Passed Vit filter: 50 (0.00362687); 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): 2 [number of targets reported over threshold] # CPU time: 0.21u 0.01s 00:00:00.22 Elapsed: 00:00:00.07 # Mc/sec: 27678.29 // [ok]