# 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_679.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_679 [L=150] Description: # 772177 # 772626 # 1 # ID=1_679;partial=00;start_type=ATG;rbs_motif=AGGAG;rbs_spacer=5-10bp;gc_cont=0.309 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.2e-60 202.0 0.6 3.6e-60 201.9 0.6 1.0 1 LacAB_rpiB Ribose/Galactose Isomerase ------ inclusion threshold ------ 0.075 13.0 0.7 0.43 10.5 0.2 2.1 2 Gly_rich_SFCGS Glycine-rich SFCGS 0.085 12.7 0.0 0.26 11.1 0.0 1.7 2 DUF2000 Protein of unknown function (DUF2000) 0.11 13.1 1.9 0.72 10.5 0.5 2.3 2 PTS_IIB PTS system, Lactose/Cellobiose specific IIB s Domain annotation for each model (and alignments): >> LacAB_rpiB Ribose/Galactose Isomerase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 201.9 0.6 8.6e-64 3.6e-60 1 139 [] 3 141 .. 3 141 .. 0.98 Alignments for each domain: == domain 1 score: 201.9 bits; conditional E-value: 8.6e-64 EEEEE-GCGHHHHHHHHHHH.HHTTTEEEEECESSSS---HHHHHHHHHHHHHTTSSSEEEEEESSSHHHHHHHHTSTT--EEE CS LacAB_rpiB 1 iaigsDhagfelkekiieeL.kkkghevidvgteseesvdYpdfakkvaeavasgeadlgiliCgtGiGmsiaankvkgiraal 83 ia+g+Dh+g+elkeki+++L kk+g+evid+gt+s+esvdYp++++ va+ v ++e+d+gil+CgtGiG+siaank+kgiraa 23726_FusoPortal_Gene_679 3 IALGADHGGYELKEKIKQHLaKKEGIEVIDFGTNSTESVDYPKYGHLVAKSVVDKEVDFGILVCGTGIGISIAANKIKGIRAAN 86 89******************6556***********99*********************************************** PP -SSHHHHHHHHHTT--SEEEEETTTS-HHHHHHHHHHHHHS---.HTHHHHHHHHH CS LacAB_rpiB 84 csdeesAklarehnnanvlalGarvigeelakeivkaFlstefeeggrharrvdki 139 c+++++Akl+rehnnan+lalGar++g la +iv++Fls +fe ggrh++rvd+i 23726_FusoPortal_Gene_679 87 CTNTTMAKLTREHNNANILALGARIVGDVLALDIVDEFLSVSFE-GGRHQKRVDQI 141 *******************************************7.*********97 PP >> Gly_rich_SFCGS Glycine-rich SFCGS # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 0.6 0.0 0.12 5.2e+02 84 102 .. 22 40 .. 6 45 .. 0.76 2 ? 10.5 0.2 0.0001 0.43 43 97 .. 55 109 .. 45 127 .. 0.80 Alignments for each domain: == domain 1 score: 0.6 bits; conditional E-value: 0.12 Gly_rich_SFCGS 84 taikegkkvlgfgfldkee 102 a keg +v+ fg +e 23726_FusoPortal_Gene_679 22 LAKKEGIEVIDFGTNSTES 40 4779999999999877775 PP == domain 2 score: 10.5 bits; conditional E-value: 0.0001 Gly_rich_SFCGS 43 kkeeadlgisfcgsggagaltaqtkygykakyglrsveegvtaikegkkvlgfgf 97 ++e d+gi cg+g ++ a++ g +a + + +t ++ +++l +g 23726_FusoPortal_Gene_679 55 VDKEVDFGILVCGTGIGISIAANKIKGIRAANCTNTTMAKLTREHNNANILALGA 109 56789*********99999999999999999999999899998888888888874 PP >> DUF2000 Protein of unknown function (DUF2000) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.1 0.0 6.2e-05 0.26 67 118 .. 17 68 .. 10 83 .. 0.87 2 ? -0.8 0.0 0.31 1.3e+03 10 34 .. 78 102 .. 72 112 .. 0.71 Alignments for each domain: == domain 1 score: 11.1 bits; conditional E-value: 6.2e-05 HHHHHHHHT.T---EEEEGGGGG---HHHHHHHHTT--TTT--EEEEEEEEE CS DUF2000 67 elrekaleregltvvdFtdaaqstksyeeyraalaatpeedleyvGialvGp 118 ++ ++++++eg++v+dF + +++ +y +y + +a++ ++ +Gi ++G 23726_FusoPortal_Gene_679 17 KIKQHLAKKEGIEVIDFGTNSTESVDYPKYGHLVAKSVVDKEVDFGILVCGT 68 4556666779****************************************97 PP == domain 2 score: -0.8 bits; conditional E-value: 0.31 TS-HHHHHHHHHHHHHHHHHH-GGG CS DUF2000 10 dlpaglalNaaavlalglgaaapel 34 +++ +a+N + +++ l+++++++ 23726_FusoPortal_Gene_679 78 KIKGIRAANCTNTTMAKLTREHNNA 102 5556677787777777777777664 PP >> PTS_IIB PTS system, Lactose/Cellobiose specific IIB subunit # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 2.1 0.0 0.073 3e+02 17 40 .. 15 39 .. 11 60 .. 0.66 2 ? 10.5 0.5 0.00017 0.72 3 26 .. 63 85 .. 61 145 .. 0.81 Alignments for each domain: == domain 1 score: 2.1 bits; conditional E-value: 0.073 HHHHHHH.HHHCTEEEEETTEEGGG CS PTS_IIB 17 aekleka.akelgieaegeagvele 40 +ek+++ ak++gie++ + +e 23726_FusoPortal_Gene_679 15 KEKIKQHlAKKEGIEVIDFGTNSTE 39 6777765155555999876665444 PP == domain 2 score: 10.5 bits; conditional E-value: 0.00017 EEEESSSHHHHHHHHHHHHHHHHH CS PTS_IIB 3 llaCgaGmatSlmaaeklekaake 26 +l+Cg+G++ S+ a +k++ + + 23726_FusoPortal_Gene_679 63 ILVCGTGIGISIAA-NKIKGIRAA 85 79**********99.999966555 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (150 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 484 (0.0289612); expected 334.2 (0.02) Passed bias filter: 375 (0.022439); expected 334.2 (0.02) Passed Vit filter: 25 (0.00149593); expected 16.7 (0.001) Passed Fwd filter: 4 (0.000239349); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 4 [number of targets reported over threshold] # CPU time: 0.23u 0.24s 00:00:00.47 Elapsed: 00:00:00.32 # Mc/sec: 1362.67 // [ok]