# 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_1475.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1475 [L=108] Description: # 1576140 # 1576463 # 1 # ID=1_1475;partial=00;start_type=ATG;rbs_motif=AGGA;rbs_spacer=5-10bp;gc_cont=0.238 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 8.4e-06 26.2 0.7 1.5e-05 25.4 0.7 1.7 1 FadA Adhesion protein FadA ------ inclusion threshold ------ 0.16 12.2 0.9 0.19 12.0 0.9 1.2 1 Ly49 Ly49-like protein, N-terminal region 0.19 12.1 4.3 0.2 12.0 3.2 1.6 2 Inhibitor_I29 Cathepsin propeptide inhibitor domain (I29) Domain annotation for each model (and alignments): >> FadA Adhesion protein FadA # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 25.4 0.7 2.7e-09 1.5e-05 1 41 [. 1 46 [. 1 107 [. 0.90 Alignments for each domain: == domain 1 score: 25.4 bits; conditional E-value: 2.7e-09 FadA 1 MKKillvllgvlvvssisfaaekt.....asleseldaleaqleeL 41 MKKi+l+l++vl++ + ++++ asl +++ al++++ ++ 23726_FusoPortal_Gene_1475 1 MKKIILMLVSVLIINACTSTKNAPfneveASLNQKYGALSNEYYKM 46 ****************99999998889999***********99765 PP >> Ly49 Ly49-like protein, N-terminal region # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.0 0.9 3.4e-05 0.19 17 98 .. 4 91 .. 1 103 [. 0.71 Alignments for each domain: == domain 1 score: 12.0 bits; conditional E-value: 3.4e-05 XXXXXXXXXXXXXXXXXXXXXX.XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX..HH.........HHHHHHHHHHHHHHHHH CS Ly49 17 llLvtvavLvtkifqysqekheleetlnlsqnysimqndinlkeelLrnkslecsilkd.........lleslkreqkrlysktk 92 ++L+ v+vL+++ +++ e +l+q+y ++ n+ + ++L+n +e ++ + ++ lk+++k+ s++ 23726_FusoPortal_Gene_1475 4 IILMLVSVLIINACTSTKNAPFNEVEASLNQKYGALSNEYY---KMLENPIVEKDRRSIlnkfesfrtEVRDLKKNRKNPSSNET 85 57999*******999988887777777**********9988...68999999887655333455555555556665555555555 PP TTTSXS CS Ly49 93 tvlksl 98 vl+s+ 23726_FusoPortal_Gene_1475 86 RVLNSF 91 555555 PP >> Inhibitor_I29 Cathepsin propeptide inhibitor domain (I29) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -1.6 0.2 0.63 3.5e+03 43 49 .. 25 34 .. 13 45 .. 0.56 2 ? 12.0 3.2 3.5e-05 0.2 13 56 .. 54 97 .. 33 99 .. 0.86 Alignments for each domain: == domain 1 score: -1.6 bits; conditional E-value: 0.63 Inhibitor_I29 43 f...tlglNk 49 f + +lN+ 23726_FusoPortal_Gene_1475 25 FnevEASLNQ 34 3333666774 PP == domain 2 score: 12.0 bits; conditional E-value: 3.5e-05 Inhibitor_I29 13 rhsseinerfenFkeNlrqiqehNanknvtftlglNkfsDlshe 56 ++ + i+ +fe+F++ +r ++ ++n+++ t lN f D s++ 23726_FusoPortal_Gene_1475 54 KDRRSILNKFESFRTEVRDLKKNRKNPSSNETRVLNSFIDKSST 97 4778899*****************9999999*********9986 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (108 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 829 (0.0496051); expected 334.2 (0.02) Passed bias filter: 437 (0.0261489); expected 334.2 (0.02) Passed Vit filter: 28 (0.00167544); expected 16.7 (0.001) Passed Fwd filter: 3 (0.000179512); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 3 [number of targets reported over threshold] # CPU time: 0.25u 0.26s 00:00:00.51 Elapsed: 00:00:00.34 # Mc/sec: 923.41 // [ok]