# hmmscan :: search sequence(s) against a profile 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 sequence file: AA/23726_FusoPortal_Gene_1311.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1311 [L=254] Description: # 1425508 # 1426269 # -1 # ID=1_1311;partial=00;start_type=ATG;rbs_motif=AGGAGG;rbs_spacer=5-10bp;gc_cont=0.324 Scores for complete sequence (score includes all domains): --- full sequence --- --- best 1 domain --- -#dom- E-value score bias E-value score bias exp N Model Description ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2e-100 335.0 12.2 2.3e-100 334.8 12.2 1.0 1 DUF3100 Protein of unknown function (DUF3100) ------ inclusion threshold ------ 1.2 9.2 7.7 2.2 8.4 3.9 2.8 3 Phage_holin_3_3 LydA holin phage, holin superfamily III Domain annotation for each model (and alignments): >> DUF3100 Protein of unknown function (DUF3100) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 334.8 12.2 2.7e-104 2.3e-100 1 234 [. 15 246 .. 15 247 .. 0.99 Alignments for each domain: == domain 1 score: 334.8 bits; conditional E-value: 2.7e-104 DUF3100 1 viaeligtlkikigkgkivllpllyalvlglllslkklkvlkkkevkaaskllsvalllliaklgltvgpslkklleagpall 83 vi+e+ig++k+++gkg i+l+p+lya+v+g+l+ k+lk++++k+++ a++l+sv+llll+ak+g+t+gps++kl++a+pal+ 23726_FusoPortal_Gene_1311 15 VISEYIGIVKFNLGKGVIALFPMLYAMVFGVLT--KFLKITNEKDMEDAGSLVSVTLLLLMAKYGTTIGPSFPKLVSASPALI 95 799***************************997..56788******************************************* PP DUF3100 84 lqelGnlgtillalPvAlllglkreaiGatfsiaRepnlaliaekygldspegrGvlavYivGtvfGtlfisllasllaslgi 166 lqe+Gn+gt+ll++P+Al+lglkre+iGat+siaRepn+a+ia+kygldspeg+Gvl+vYivGtvfGt+fi+llas+la++++ 23726_FusoPortal_Gene_1311 96 LQEFGNIGTVLLGVPIALFLGLKRETIGATHSIAREPNIAIIADKYGLDSPEGEGVLGVYIVGTVFGTIFIGLLASFLAAYTP 178 *********************************************************************************** PP DUF3100 167 fhplalamasGvGsgsmmaaasgalaaalpemaeeilalaaasnlltsvtgvyvlifialPlaeklYk 234 +hp++lamasGvGs+smm+a++gal++++p+ma++i+a++a snll++++gvy++i+ialPl+e+lYk 23726_FusoPortal_Gene_1311 179 LHPYSLAMASGVGSASMMTASVGALSTLYPDMADTIAAFGASSNLLSGLDGVYMSIWIALPLTEYLYK 246 *******************************************************************8 PP >> Phage_holin_3_3 LydA holin phage, holin superfamily III # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -3.2 0.0 1 8.6e+03 38 46 .. 40 48 .. 30 59 .. 0.59 2 ? 1.1 0.0 0.047 3.9e+02 27 43 .. 101 117 .. 97 138 .. 0.79 3 ? 8.4 3.9 0.00026 2.2 32 67 .. 157 192 .. 148 202 .. 0.88 Alignments for each domain: == domain 1 score: -3.2 bits; conditional E-value: 1 Phage_holin_3_3 38 SgfvGllvf 46 +++ G+l+ 23726_FusoPortal_Gene_1311 40 AMVFGVLTK 48 444444443 PP == domain 2 score: 1.1 bits; conditional E-value: 0.047 Phage_holin_3_3 27 sikellleivvSgfvGl 43 +i +ll++++++f Gl 23726_FusoPortal_Gene_1311 101 NIGTVLLGVPIALFLGL 117 567799**********8 PP == domain 3 score: 8.4 bits; conditional E-value: 0.00026 Phage_holin_3_3 32 lleivvSgfvGllvfllceelglsqlltlaiaglsG 67 ++ ++ +f+Gll+++l++++ l ++ ++++g++ 23726_FusoPortal_Gene_1311 157 VGTVFGTIFIGLLASFLAAYTPLHPYSLAMASGVGS 192 56677889*************************976 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (254 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 579 (0.0346458); expected 334.2 (0.02) Passed bias filter: 316 (0.0189086); expected 334.2 (0.02) Passed Vit filter: 31 (0.00185495); expected 16.7 (0.001) Passed Fwd filter: 2 (0.000119674); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.25u 0.23s 00:00:00.48 Elapsed: 00:00:00.32 # Mc/sec: 2307.46 // [ok]