# 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/Pan_troglodytes.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1.2e-55 188.5 0.4 2.2e-55 187.7 0.4 1.4 1 tr|H2PZ90|H2PZ90_PANTR TYW3 isoform 2 OS=Pan troglodytes OX= Domain annotation for each sequence (and alignments): >> tr|H2PZ90|H2PZ90_PANTR TYW3 isoform 2 OS=Pan troglodytes OX=9598 GN=TYW3 PE=2 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 187.7 0.4 9.6e-60 2.2e-55 3 213 .. 11 193 .. 9 194 .. 0.94 Alignments for each domain: == domain 1 score: 187.7 bits; conditional E-value: 9.6e-60 HHHHHHHH....HHHHHCCHC-GGGHHHHHHHHC-SSEEEEEEES-EEEEE............EESSTT-STT-EEEEEESS---HH CS TYW3 3 KakiLeeLsvpdtdkspkGsvDekirelideinsledlvTtSSCsGRvsvflegrksaeaeekaaksggkkgggkwlfvsHdpvkee 89 Ka++L++ d+s+kGsvDe+++el++ +n +++++TtSSC+GR+ ++ ++ + + +k+++ wl+v+H+ + ++ tr|H2PZ90|H2PZ90_PANTR 11 KAQCLSK-----ADLSRKGSVDEDVVELVQFLNMRDQFFTTSSCAGRILLL---------DQGINGFEVQKQNCCWLLVTHKLCVKD 83 6777776.....6**************************************.........56778888999999************* PP HHHHHHCS--..............SSEEEEEEE--EEEEEESSHHHHHHHHHHHHHCT-TT-EEEEECTTEEEEEEE---EEEEEEE CS TYW3 90 ellaelegleeeeesaksseseeeerlivfkfePmILhvlcrslehAqkllsaAlsaGfreSGivslkkkkvivairsglklevpig 176 +++ +l++ + + ++kfeP++Lhv+cr+l++Aq l+s+A+++Gfr+SGi+ k++k ++a+rs+++levp++ tr|H2PZ90|H2PZ90_PANTR 84 DVIVALKKAN---------------GDATLKFEPFVLHVQCRQLQDAQILHSMAIDSGFRNSGITVGKRGKTMLAVRSTHGLEVPLS 155 *******994...............6789********************************************************** PP SST.S---HHHHHHHHHHHHHHHHHHHHHHHHHHHHC- CS TYW3 177 rkk.llvseeylkllvklanerfeenkkrierfeealk 213 +k l+v+eey+++l+++an+++eenkkrierf + l+ tr|H2PZ90|H2PZ90_PANTR 156 HKGkLMVTEEYIDFLLNVANQKMEENKKRIERFCNCLQ 193 5555*****************************99886 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (214 nodes) Target sequences: 23008 (11942255 residues searched) Passed MSV filter: 892 (0.0387691); expected 460.2 (0.02) Passed bias filter: 547 (0.0237743); expected 460.2 (0.02) Passed Vit filter: 48 (0.00208623); expected 23.0 (0.001) Passed Fwd filter: 3 (0.000130389); expected 0.2 (1e-05) Initial search space (Z): 23008 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.21u 0.01s 00:00:00.22 Elapsed: 00:00:00.08 # Mc/sec: 31945.53 // [ok]