# 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_563.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_563 [L=187] Description: # 622509 # 623069 # 1 # ID=1_563;partial=00;start_type=ATG;rbs_motif=GGA/GAG/AGG;rbs_spacer=11-12bp;gc_cont=0.251 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 3.8e-25 88.3 0.4 4.7e-25 88.0 0.4 1.0 1 GDYXXLXY GDYXXLXY protein ------ inclusion threshold ------ 0.52 10.8 7.5 2.7 8.5 7.5 2.3 1 DUF3139 Protein of unknown function (DUF3139) Domain annotation for each model (and alignments): >> GDYXXLXY GDYXXLXY protein # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 88.0 0.4 5.7e-29 4.7e-25 3 135 .. 19 163 .. 17 177 .. 0.89 Alignments for each domain: == domain 1 score: 88.0 bits; conditional E-value: 5.7e-29 GDYXXLXY 3 lvllalivqaestlatGkevlLetaPvDPrdllrGdYvaLrYdisr...................ke..kaekedeleagttvy 65 ++ ++ ++++e+ ++ ++Le+aPvDPr++l+GdY++L+Ydi + +e + ++ ++l++ +++y 23726_FusoPortal_Gene_563 19 VITGFSAMKEENYKKLDSYFYLELAPVDPRSILQGDYMTLNYDITDkvsdfiynnrtyiydgeneNEveEIRELRKLADAKKAY 102 5778899999****************************************************986224689999********** PP GDYXXLXY 66 vilekdadglaqvarvsrsktawkplpanevaLkgkvrdgrveyglerYfvpegqaeeleqalreealyv 135 ++++ d++++a++ +v+++kt +++ +++++k ++ v+++ ++Y ++eg+++++++a +y 23726_FusoPortal_Gene_563 103 IAVRLDENKVAKFVKVTKEKT----DEKDLLFIAYKTDGFDVNINANSYLFQEGTGNKYQDA-----RYS 163 **********99999999887....457999****99999*********************4.....443 PP >> DUF3139 Protein of unknown function (DUF3139) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 8.5 7.5 0.00033 2.7 2 37 .. 6 56 .. 5 158 .. 0.77 Alignments for each domain: == domain 1 score: 8.5 bits; conditional E-value: 0.00033 DUF3139 2 kkiiiilillillivillvgskpikylgkkidkyl.................e 37 kki +i+ ++il+++++++ k+ +y kk+d+y+ 23726_FusoPortal_Gene_563 6 KKILLIINIIILFVITGFSAMKEENY--KKLDSYFylelapvdprsilqgdyM 56 89999999999999996555544334..******9777777766666655441 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (187 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 1008 (0.0603159); expected 334.2 (0.02) Passed bias filter: 527 (0.0315342); expected 334.2 (0.02) Passed Vit filter: 45 (0.00269268); 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.24u 0.23s 00:00:00.47 Elapsed: 00:00:00.31 # Mc/sec: 1753.60 // [ok]