# 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/TYW3.hmm.txt # target sequence database: proteomes/Macaca_mulatta.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: TYW3 [M=214] Accession: PF02676.13 Description: Methyltransferase TYW3 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 3.5e-43 147.7 0.2 9.1e-23 80.9 0.0 2.3 2 tr|F7GQ32|F7GQ32_MACMU Uncharacterized protein OS=Macaca mul Domain annotation for each sequence (and alignments): >> tr|F7GQ32|F7GQ32_MACMU Uncharacterized protein OS=Macaca mulatta OX=9544 GN=TYW3 PE=4 SV=2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 65.8 0.0 1.8e-22 3.9e-18 3 90 .. 11 84 .. 9 86 .. 0.87 2 ! 80.9 0.0 4.3e-27 9.1e-23 138 213 .. 84 160 .. 83 161 .. 0.95 Alignments for each domain: == domain 1 score: 65.8 bits; conditional E-value: 1.8e-22 HHHHHHHH....HHHHHCCHC-GGGHHHHHHHHC-SSEEEEEEES-EEEEE............EESSTT-STT-EEEEEESS---HHH CS TYW3 3 KakiLeeLsvpdtdkspkGsvDekirelideinsledlvTtSSCsGRvsvflegrksaeaeekaaksggkkgggkwlfvsHdpvkeee 90 Ka++L++ d+s+kGsvD++++el++ +n++++++TtSSC+GR+ ++ + + + +k+++ wl+v+H+p+++++ tr|F7GQ32|F7GQ32_MACMU 11 KAQCLSK-----ADLSRKGSVDKDVVELVQFLNTRDQFFTTSSCAGRILLL---------DGGINGFEVQKQNCGWLLVTHKPCEKDD 84 6777776.....6**************************************.........344567777788899*********8876 PP == domain 2 score: 80.9 bits; conditional E-value: 4.3e-27 HHHHHHHHCT-TT-EEEEECTTEEEEEEE---EEEEEEESST.S---HHHHHHHHHHHHHHHHHHHHHHHHHHHHC- CS TYW3 138 kllsaAlsaGfreSGivslkkkkvivairsglklevpigrkk.llvseeylkllvklanerfeenkkrierfeealk 213 +++s+A+++Gfr+SGi+ k++k+++a+rs+++levp+++k l+v+eey+++l+++an+++eenkkrierf++ l+ tr|F7GQ32|F7GQ32_MACMU 84 DVHSVAIDSGFRNSGITVGKRGKIMLAVRSTHGLEVPLSHKGkLMVTEEYIDFLLNVANQKMEENKKRIERFYNCLH 160 5899***********************************5555******************************9886 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (214 nodes) Target sequences: 21211 (10888574 residues searched) Passed MSV filter: 832 (0.0392249); expected 424.2 (0.02) Passed bias filter: 492 (0.0231955); expected 424.2 (0.02) Passed Vit filter: 39 (0.00183867); expected 21.2 (0.001) Passed Fwd filter: 3 (0.000141436); expected 0.2 (1e-05) Initial search space (Z): 21211 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.19u 0.01s 00:00:00.20 Elapsed: 00:00:00.08 # Mc/sec: 29126.94 // [ok]