# 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/Strongylocentrotus_purpuratus.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 4.4e-18 66.5 6.9 7.1e-18 65.8 6.9 1.3 1 tr|W4Z511|W4Z511_STRPU Uncharacterized protein OS=Strongyloc ------ inclusion threshold ------ 1.6 10.2 3.8 0.75 11.2 0.2 2.3 2 tr|W4XD29|W4XD29_STRPU Uncharacterized protein OS=Strongyloc Domain annotation for each sequence (and alignments): >> tr|W4Z511|W4Z511_STRPU Uncharacterized protein OS=Strongylocentrotus purpuratus OX=7668 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 65.8 6.9 5e-22 7.1e-18 1 81 [. 203 279 .. 203 280 .. 0.86 Alignments for each domain: == domain 1 score: 65.8 bits; conditional E-value: 5e-22 WBS_methylT 1 LfaGkseqqqlpkglgedgeeeseeaeqvkyekrrrkrkkrkkkkkkavksreWilkKKerrrkkGk.eVkrdSKYTGRkRk 81 LfaG ++ +lpk+lg+ ++ +q ++++ r + ++ k+k + +ksreWi +KKerrr++Gk eV++dSKYTGRkR tr|W4Z511|W4Z511_STRPU 203 LFAG-VSNPTLPKALGTGNQRGVA--TQASFTDSRMRYRNMKGK--SVKKSREWIKEKKERRRRQGKgEVRPDSKYTGRKRS 279 8999.89999*******9976655..888888888877755555..4568*******************************6 PP >> tr|W4XD29|W4XD29_STRPU Uncharacterized protein OS=Strongylocentrotus purpuratus OX=7668 PE=3 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.2 0.2 5.2e-05 0.75 37 71 .. 303 337 .. 275 339 .. 0.62 2 ? -2.4 0.3 0.92 1.3e+04 18 41 .. 353 378 .. 345 384 .. 0.56 Alignments for each domain: == domain 1 score: 11.2 bits; conditional E-value: 5.2e-05 WBS_methylT 37 krkkrkkkkkkavk.....sreWilkKKerrrkkGkeVkr 71 + + ++ + sr+W +++Ke+ +++G ++k+ tr|W4XD29|W4XD29_STRPU 303 RGQGDGGD-----TegkpdSRDWSFAQKELLKHQGIDIKE 337 22222222.....123356****************99986 PP == domain 2 score: -2.4 bits; conditional E-value: 0.92 WBS_methylT 18 dgeeese..eaeqvkyekrrrkrkkr 41 ++ee ++ e+++ +ye++ + +k+ tr|W4XD29|W4XD29_STRPU 353 EKEEVDQllEQQRLDYESKLQALQKQ 378 22233335677777777776655554 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (82 nodes) Target sequences: 28594 (14266044 residues searched) Passed MSV filter: 3043 (0.106421); expected 571.9 (0.02) Passed bias filter: 836 (0.0292369); expected 571.9 (0.02) Passed Vit filter: 131 (0.00458138); expected 28.6 (0.001) Passed Fwd filter: 11 (0.000384696); expected 0.3 (1e-05) Initial search space (Z): 28594 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.28u 0.01s 00:00:00.29 Elapsed: 00:00:00.10 # Mc/sec: 11698.16 // [ok]