# 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/WBS_methylT.hmm.txt # target sequence database: proteomes/Trypanosoma_brucei.fasta # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: WBS_methylT [M=82] Accession: PF12589.7 Description: Methyltransferase involved in Williams-Beuren syndrome 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 6.9e-08 32.1 7.9 1.7e-07 30.9 7.9 1.7 1 tr|Q583K3|Q583K3_TRYB2 Uncharacterized protein OS=Trypanosom ------ inclusion threshold ------ 0.22 11.3 3.4 2 8.2 0.1 2.5 2 tr|Q38DU6|Q38DU6_TRYB2 Elongation factor 1-alpha, putative O Domain annotation for each sequence (and alignments): >> tr|Q583K3|Q583K3_TRYB2 Uncharacterized protein OS=Trypanosoma brucei brucei (strain 927/4 GUTat10.1) OX=185431 GN=Tb # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 30.9 7.9 3.9e-11 1.7e-07 33 81 .. 328 376 .. 282 377 .. 0.64 Alignments for each domain: == domain 1 score: 30.9 bits; conditional E-value: 3.9e-11 WBS_methylT 33 krrrkrkkrkkkkkkavksreWilkKKerrrkkGkeVkrdSKYTGRkRk 81 + r+ +r+k+++ + ++eW+l KK +rrk G + +d+KYT R+R+ tr|Q583K3|Q583K3_TRYB2 328 GGRQAGARRRKDNRPVTGTKEWVLLKKAERRKFGLKTSEDTKYTMRRRR 376 222233333444444445******************************8 PP >> tr|Q38DU6|Q38DU6_TRYB2 Elongation factor 1-alpha, putative OS=Trypanosoma brucei brucei (strain 927/4 GUTat10.1) OX= # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 2.6 0.9 0.026 1.1e+02 25 48 .. 164 187 .. 136 198 .. 0.50 2 ? 8.2 0.1 0.00047 2 29 71 .. 226 267 .. 208 270 .. 0.75 Alignments for each domain: == domain 1 score: 2.6 bits; conditional E-value: 0.026 WBS_methylT 25 eaeqvkyekrrrkrkkrkkkkkka 48 + ++ + +++r+++++++k+kk+ tr|Q38DU6|Q38DU6_TRYB2 164 SPSSAQLQQQRQQQQQQHKGKKDM 187 444455555555555554444443 PP == domain 2 score: 8.2 bits; conditional E-value: 0.00047 WBS_methylT 29 vkyekrrrkrkkrkkkkkkavksreWilkKKerrrkkGkeVkr 71 v ++ +r++k+ ++ kk++ + W+l + e++r++G + + tr|Q38DU6|Q38DU6_TRYB2 226 VTQSEVDRNEKNARQMKKESF-KFAWLLDQSEEERRRGVTIDA 267 555556666665555555433.789*************88765 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (82 nodes) Target sequences: 8561 (4332363 residues searched) Passed MSV filter: 753 (0.087957); expected 171.2 (0.02) Passed bias filter: 198 (0.0231281); expected 171.2 (0.02) Passed Vit filter: 18 (0.00210256); expected 8.6 (0.001) Passed Fwd filter: 4 (0.000467235); expected 0.1 (1e-05) Initial search space (Z): 8561 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.08u 0.01s 00:00:00.09 Elapsed: 00:00:00.04 # Mc/sec: 8881.34 // [ok]