# 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/Amphimedon_queenslandica.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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1.5e-22 81.4 4.0 2.4e-22 80.7 4.0 1.3 1 tr|A0A1X7VIY9|A0A1X7VIY9_AMPQE Uncharacterized protein OS=Amphi ------ inclusion threshold ------ 0.039 15.9 3.0 0.039 15.9 3.0 5.6 7 tr|A0A1X7V6P6|A0A1X7V6P6_AMPQE Uncharacterized protein OS=Amphi Domain annotation for each sequence (and alignments): >> tr|A0A1X7VIY9|A0A1X7VIY9_AMPQE Uncharacterized protein OS=Amphimedon queenslandica OX=400682 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 80.7 4.0 1.1e-26 2.4e-22 1 82 [] 200 277 .. 200 277 .. 0.91 Alignments for each domain: == domain 1 score: 80.7 bits; conditional E-value: 1.1e-26 WBS_methylT 1 LfaGkseqqqlpkglgedgeeeseeaeqvkyekrrrkrkkrkkkkkkavksreWilkKKerrrkkGkeVkrdSKYTGRk 79 Lf+G +++ +lpk+lg+++ ++ e++q+++ ++r+k k r++k+k +vksr+Wi++KKer++++GkeV++dSKYTGRk tr|A0A1X7VIY9|A0A1X7VIY9_AMPQE 200 LFTG-GAPVALPKALGTEDLKS--EQNQISFVEARKKLK-RRGKDKVPVKSRDWIIAKKERKQRQGKEVRPDSKYTGRK 274 7999.89************765..46*******999988.577777799****************************** PP WBS_methylT 80 Rka 82 Rk+ tr|A0A1X7VIY9|A0A1X7VIY9_AMPQE 275 RKP 277 *95 PP >> tr|A0A1X7V6P6|A0A1X7V6P6_AMPQE Uncharacterized protein OS=Amphimedon queenslandica OX=400682 PE=4 SV=1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 15.9 3.0 1.8e-06 0.039 19 81 .. 181 248 .. 164 249 .. 0.75 2 ? -0.4 0.2 0.22 4.8e+03 24 45 .. 282 311 .. 265 328 .. 0.54 3 ? 0.7 6.3 0.1 2.2e+03 8 59 .. 414 466 .. 408 472 .. 0.76 4 ? -4.0 0.2 2 4.3e+04 29 40 .. 539 550 .. 527 573 .. 0.57 5 ? -2.8 0.7 1.3 2.8e+04 7 20 .. 787 800 .. 775 824 .. 0.53 6 ? -2.0 0.4 0.71 1.5e+04 10 44 .. 904 938 .. 898 968 .. 0.70 7 ? -1.8 0.4 0.59 1.3e+04 24 41 .. 1013 1030 .. 995 1054 .. 0.52 Alignments for each domain: == domain 1 score: 15.9 bits; conditional E-value: 1.8e-06 WBS_methylT 19 geeeseeaeqvkyekrrrkrkkr.....kkkkkkavksreWilkKKerrrkkGkeVkrdSKYTGRkRk 81 +++++ee+e+v++++++r ++ +++++++++ ++++ +K + +++k keVk +S + GRkR+ tr|A0A1X7V6P6|A0A1X7V6P6_AMPQE 181 ESDKEEEEEKVQIKEQERVIETIiidddNDDDDDDEEVKDYLKQKDDFEQRKMKEVKYTSPFRGRKRP 248 23344478999999999966665554334444455559*****************************7 PP == domain 2 score: -0.4 bits; conditional E-value: 0.22 WBS_methylT 24 e......eaeqvkyekrrr..krkkrkkkk 45 + +a+qv++ +++ kkrk + tr|A0A1X7V6P6|A0A1X7V6P6_AMPQE 282 KettlkpSAKQVHVHVQENkeGTKKRKTVR 311 333444544455444444400344333333 PP == domain 3 score: 0.7 bits; conditional E-value: 0.1 WBS_methylT 8 qqqlpkglgedgeeese.eaeqvkyekrrrkrkkrkkkkkkavksreWilkKK 59 +++ ++++ ++ +e++ +e+v++++r+++++k+ k +++ k+re + +K tr|A0A1X7V6P6|A0A1X7V6P6_AMPQE 414 REEIEEAKEIKELSERQkVKERVDIKRREEMKQKEDMKLREEIKRREEVNEKE 466 4566666666655544447899********99999999999999999887775 PP == domain 4 score: -4.0 bits; conditional E-value: 2 WBS_methylT 29 vkyekrrrkrkk 40 ++ kr+++rk+ tr|A0A1X7V6P6|A0A1X7V6P6_AMPQE 539 IDEGKRTKRRKS 550 333333333332 PP == domain 5 score: -2.8 bits; conditional E-value: 1.3 WBS_methylT 7 eqqqlpkglgedge 20 ++++ g+g +g+ tr|A0A1X7V6P6|A0A1X7V6P6_AMPQE 787 KKKSSESGKGSSGT 800 23333344444443 PP == domain 6 score: -2.0 bits; conditional E-value: 0.71 WBS_methylT 10 qlpkglgedgeeeseeaeqvkyekrrrkrkkrkkk 44 +l g+++++++++ +++ +++kr+ + +++kk+ tr|A0A1X7V6P6|A0A1X7V6P6_AMPQE 904 KLATGRKKQASATEPKKKLLDISKRKPPTSTEKKQ 938 55666666665555577888888888877776655 PP == domain 7 score: -1.8 bits; conditional E-value: 0.59 WBS_methylT 24 eeaeqvkyekrrrkrkkr 41 ++++ v + +++++++++ tr|A0A1X7V6P6|A0A1X7V6P6_AMPQE 1013 NNNSIVGIITKKEENDAK 1030 234444444444443333 PP Internal pipeline statistics summary: ------------------------------------- Query model(s): 1 (82 nodes) Target sequences: 43435 (14716620 residues searched) Passed MSV filter: 3362 (0.077403); expected 868.7 (0.02) Passed bias filter: 1125 (0.0259008); expected 868.7 (0.02) Passed Vit filter: 102 (0.00234834); expected 43.4 (0.001) Passed Fwd filter: 5 (0.000115115); expected 0.4 (1e-05) Initial search space (Z): 43435 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.24u 0.01s 00:00:00.25 Elapsed: 00:00:00.12 # Mc/sec: 10056.36 // [ok]