# 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_712.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_712 [L=258] Description: # 804956 # 805729 # -1 # ID=1_712;partial=00;start_type=ATG;rbs_motif=AGGA;rbs_spacer=5-10bp;gc_cont=0.256 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 6.9e-41 139.8 0.0 8.9e-41 139.4 0.0 1.1 1 adh_short short chain dehydrogenase 4.3e-20 72.2 0.0 5.2e-20 71.9 0.0 1.1 1 adh_short_C2 Enoyl-(Acyl carrier protein) reductase 6.4e-09 35.5 0.0 2e-08 33.9 0.0 1.9 2 Epimerase NAD dependent epimerase/dehydratase family 9.4e-08 32.1 0.0 1.4e-07 31.5 0.0 1.3 1 KR KR domain 3e-07 29.8 0.0 5.1e-07 29.0 0.0 1.5 1 Polysacc_synt_2 Polysaccharide biosynthesis protein 3.4e-06 26.3 0.0 5.6e-06 25.6 0.0 1.4 1 RmlD_sub_bind RmlD substrate binding domain 5.2e-06 26.4 0.0 8.1e-06 25.8 0.0 1.3 1 NAD_binding_10 NAD(P)H-binding 1e-05 25.0 0.0 1.4e-05 24.6 0.0 1.2 1 GDP_Man_Dehyd GDP-mannose 4,6 dehydratase 8.3e-05 22.2 0.0 0.0026 17.2 0.0 2.1 2 NmrA NmrA-like family 0.00019 20.4 0.1 0.00037 19.5 0.1 1.5 1 3Beta_HSD 3-beta hydroxysteroid dehydrogenase/isomeras ------ inclusion threshold ------ 0.081 13.4 0.3 0.48 11.0 0.1 2.2 2 T_Ag_DNA_bind Origin of replication binding protein 0.18 11.6 0.0 0.76 9.6 0.0 1.9 2 Glyco_tran_WecB Glycosyl transferase WecB/TagA/CpsF family Domain annotation for each model (and alignments): >> adh_short short chain dehydrogenase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 139.4 0.0 6.4e-44 8.9e-41 2 190 .. 3 190 .. 2 194 .. 0.97 Alignments for each domain: == domain 1 score: 139.4 bits; conditional E-value: 6.4e-44 EEEEETTTSHHHHHHHHHHHHTT-EEEEEESCHHHHHHHHHHHHHT..TSEEEEEE-TTTSHHHHHHHHHHHHHHHSSESEEEE CS adh_short 2 vvlvtGassGIGraiakelakeGakvvvvdrseekleevakelkee..gakalaiqvDvtdreevealvekavkklgkldvLvn 83 +l+tGassGIG ++ ++la++ k+++++rs +kl+ ++kel+e+ + + ++i++D+td +++e++ve++ +d+++n 23726_FusoPortal_Gene_712 3 KILITGASSGIGEELTRNLANKSKKLFLLARSLDKLNLLKKELEEKfsSLECVCIKYDLTDINNLENIVENC-----DVDLVIN 81 69********************************************99999*********************.....9****** PP ---SS--SSGGGGSHHHHHHHHHHHTHHHHHHHHHHHHHHHHCTCEEEEEEEEGGGTSB-TTBHHHHHHHHHHHHHHHHHHHHH CS adh_short 84 NAGitgsgpfselseeewervievNltgvflltravlpamlkqssGrIvnvsSvagkvpvagvaaYsAsKaavagltrslakel 167 +AG + ++fs+ls++e ++i+vN+ +++ lt+ + +++l++ +G+I n++S+a++++ + +a+Ys++K a+ +++ +l +el 23726_FusoPortal_Gene_712 82 CAGFGKITDFSKLSDKEDLDTINVNFISPLILTKKFSEKFLQKGQGTILNICSTAALYQHPYMAVYSSAKSALLHYSLALDEEL 165 ************************************************************************************ PP GGCT..EEEEEEEE-SBSSHHHHHH CS adh_short 168 aptg..irvnavaPGlvdTdmakel 190 + ++ +rv +v+PG++++++ ++ 23726_FusoPortal_Gene_712 166 SHKNknVRVLSVCPGPTASNFFEKD 190 998888************9987765 PP >> adh_short_C2 Enoyl-(Acyl carrier protein) reductase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 71.9 0.0 3.7e-23 5.2e-20 1 183 [. 8 192 .. 8 218 .. 0.87 Alignments for each domain: == domain 1 score: 71.9 bits; conditional E-value: 3.7e-23 TTTTSHHHHHHHHHHHTT-EEEEEESSGGH.HHHHHHHHHHTTSE..E.EE--TTSHHHHHHHHHHHHHHHSSEEEEEE----- CS adh_short_C2 1 GaasgIGraiAkalaeeGaevvltdlneea.kkrveelaeelgae..v.lecDVtdeeqvealveevvekfgklDilvnnagia 80 Ga sgIG+++ + la++ ++ l ++ ++ + +el e++++ v ++ D td +++e++ve++ +D ++n+ag+ 23726_FusoPortal_Gene_712 8 GASSGIGEELTRNLANKSKKLFLLARSLDKlNLLKKELEEKFSSLecVcIKYDLTDINNLENIVENC-----DVDLVINCAGFG 86 899**********************999996666667777777444457***************997.....59********** PP EGCSS-GGGH-HHHHHHHHHHHTHHHHHHHHHHHH.HCTSC.CEEEEEEEGGGTSBSTTTTHHHHHHHHHHHHHHHHHHHHGG. CS adh_short_C2 81 pklkgpfletsreefdkaldvnlyslfllakaalp.lmkeg.gsivnlssvgaervvpnynaygaaKaalealtrylAvelgp. 161 + ++f ++s++e ++++vn+ s + l+k + + +++g g i+n+ s++a p ++y+ aK al + + +l el + 23726_FusoPortal_Gene_712 87 KI--TDFSKLSDKEDLDTINVNFISPLILTKKFSEkFLQKGqGTILNICSTAALYQHPYMAVYSSAKSALLHYSLALDEELSHk 168 **..******************************99777777**************************************9986 PP .GTEEEEEEEE--BHCTG.GGGST CS adh_short_C2 162 .kgiRvNavspgpikTla.asgik 183 k++Rv +v+pgp++ + ++i+ 23726_FusoPortal_Gene_712 169 nKNVRVLSVCPGPTASNFfEKDIQ 192 689***********9877444454 PP >> Epimerase NAD dependent epimerase/dehydratase family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 33.9 0.0 1.4e-11 2e-08 1 103 [. 4 118 .. 4 177 .. 0.82 2 ? -1.7 0.0 1.1 1.5e+03 213 229 .. 203 219 .. 194 225 .. 0.79 Alignments for each domain: == domain 1 score: 33.9 bits; conditional E-value: 1.4e-11 EEEETTTSHHHHHHHHHHHHCTSEEEEEES-STTTTCHHCHHC........EEEE--TT-HHHHHHHHHHHT-SEEEE--S--- CS Epimerase 1 ilVtGatGfiGsalvkallekgyevigldrlssasntarledl........elveldltdrsalerlfkeeqpdavihlAavgg 76 il+tGa+ iG l+++l++k ++ l+r ++ n ++e ++ dltd ++le+++++ ++d vi++A+ g 23726_FusoPortal_Gene_712 4 ILITGASSGIGEELTRNLANKSKKLFLLARSLDKLNLLKKELEekfsslecVCIKYDLTDINNLENIVENCDVDLVINCAGFGK 87 8***********************999999977777665554489*****999****************99*********9888 PP HHH....HHHSHHHHHHHHTHHHHHHHHHHH CS Epimerase 77 vaa....ssedpadfieknlqgtlnlleaar 103 + + s ++ d+i+ n++ l l + + 23726_FusoPortal_Gene_712 88 ITDfsklSDKEDLDTINVNFISPLILTKKFS 118 7652222556778888888888777766555 PP == domain 2 score: -1.7 bits; conditional E-value: 1.1 EEHHHHHHHHHHHHHHH CS Epimerase 213 lyveDvaeaillvlehv 229 + +eDva+ i+ v+e+ 23726_FusoPortal_Gene_712 203 MSSEDVAKRIIKVIENK 219 56899999999999994 PP >> KR KR domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 31.5 0.0 1e-10 1.4e-07 4 154 .. 5 152 .. 3 168 .. 0.77 Alignments for each domain: == domain 1 score: 31.5 bits; conditional E-value: 1e-10 EEETTTSHHHHHHHHHHHHHHS-EEEEE-SS..--HHHHHHHHHHHHT..T-EEEE----TT-HHHHHHHHHHHHHHSS---EE CS KR 4 litGGlGgLGlelarwlaekgarhlvllsRsaaadedskaliaeleal..gaevevvaaDvsdkdavealleeikaeeppikGv 85 litG+ g+G el+r la+k ++l ll+Rs d+ + l +ele++ e+ +++ D++d +++e+++e+++ + v 23726_FusoPortal_Gene_712 5 LITGASSGIGEELTRNLANKS-KKLFLLARS--LDKLNL-LKKELEEKfsSLECVCIKYDLTDINNLENIVENCD-----VDLV 79 9*****************996.6899*****..344433.667888875577999*******************9.....666* PP EE---------TTT--HHHHHHHHTTTTHHHHHHHHHTTT....S--SEEEEEEEHHHHH--TT-HHHHHHHH CS KR 86 ihaAgvlrdklleemtaedleqvlapKvkGawnLhelled....esldffvlfSSiagllgsagqanYaaaNa 154 i++Ag+ + + ++++++++ ++ + L + +++ + ++ + S a+l+ ++ a Y++a + 23726_FusoPortal_Gene_712 80 INCAGFGKITDFSKLSDKEDLDTINVNFISPLILTKKFSEkflqKGQGTILNICSTAALYQHPYMAVYSSAKS 152 ******9999999987765555555555555555555544333355567777888888888888888888865 PP >> Polysacc_synt_2 Polysaccharide biosynthesis protein # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 29.0 0.0 3.6e-10 5.1e-07 1 88 [. 4 86 .. 4 124 .. 0.88 Alignments for each domain: == domain 1 score: 29.0 bits; conditional E-value: 3.6e-10 EEEETTTSHHHHHHHHHHHHHB-SEEEEEES-HHHHHHHHHH......-TTEEEE....E--TT-HHHHHHHTT..T-SEEEE--- CS Polysacc_synt_2 1 vlvtGagGsiGselvrqilkekpkkiillerdElklyeiekeleeelkkkklkveivpvigdvrdkerlekaleelkvdvvfhaAA 86 +l+tGa+ iG+el+r++++ k kk++ll r+ kl ++kelee++++ l e+v++ d++d ++le+++e+ +vd+v+++A 23726_FusoPortal_Gene_712 4 ILITGASSGIGEELTRNLAN-KSKKLFLLARSLDKLNLLKKELEEKFSS--L--ECVCIKYDLTDINNLENIVENCDVDLVINCAG 84 89*************99765.6799******************888876..5..56688899*********************997 PP -- CS Polysacc_synt_2 87 lK 88 + 23726_FusoPortal_Gene_712 85 FG 86 65 PP >> RmlD_sub_bind RmlD substrate binding domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 25.6 0.0 4e-09 5.6e-06 2 64 .. 3 88 .. 1 118 [. 0.87 Alignments for each domain: == domain 1 score: 25.6 bits; conditional E-value: 4e-09 EEEEESTTSHHHHHHHHHCTTTSEEEEEE-TT.......................T--TTSHHHHHHHHHHH--SEEEE------H CS RmlD_sub_bind 2 kilvtGanGqlGseLakqlaeeglevlaltre.......................eldltdaeavaellkeakpdvvvnaAAytav 64 kil+tGa+ +G eL+++la++ ++ l r+ + dltd +++++++++ + d+v+n+A + ++ 23726_FusoPortal_Gene_712 3 KILITGASSGIGEELTRNLANKSKKLFLLARSldklnllkkeleekfsslecvciKYDLTDINNLENIVENCDVDLVINCAGFGKI 88 79*********************9**9999999************************************************98543 PP >> NAD_binding_10 NAD(P)H-binding # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 25.8 0.0 5.8e-09 8.1e-06 1 183 [. 8 218 .. 8 219 .. 0.77 Alignments for each domain: == domain 1 score: 25.8 bits; conditional E-value: 5.8e-09 TTTSHHHHHHHHHHHHC.T-EEEEEETTHHHCH....HHHCT.T..TCEEEE--TT.CCHHHHHHTT-SE..EEES---TCC.. CS NAD_binding_10 1 GAtGkvGrllvkellargghevtalvRnpekla....eleae.e..glevvegdvtdvedlaealkgada..Vvsalggggk.. 73 GA+ ++G +l++ l++++ ++ +l+R+ +kl+ ele++ + ++ ++ d+td + l++ ++++d+ V+ ++g g+ 23726_FusoPortal_Gene_712 8 GASSGIGEELTRNLANKS-KKLFLLARSLDKLNllkkELEEKfSslECVCIKYDLTDINNLENIVENCDVdlVINCAGFGKItd 90 9****************7.8************844444443333225557899*************999665999997777767 PP .....HHHH................HHHHHHHHHHCT--EEEEEEETTTTSGCGCHHHHHHHHCCHHHHHHHHHHHHHHHH..S CS NAD_binding_10 74 .....deeg................aknivdaakaagvkRfvlvssigvgdevpgefgewnkellkaylaakraaeellra..s 134 d +k+ + ++g ++ + s + +++ ++++ +a l a +e+l + + 23726_FusoPortal_Gene_712 91 fsklsD--KedldtinvnfisplilTKKFSEKFLQKGQGTILNICSTAALY--QHPYMAVYSSAKSALLHYSLALDEELSHknK 170 666552..1446778999999999999999999999988887665555553..4667788888878888888899999865337 PP SSEEEEEEE-CEE-ST.--SEEEEETTSS----EEEHHHHHHHHHHHCT- CS NAD_binding_10 135 gldyTivrpggLtdep.etkyrlvekgepfkggsisradVAdvlvaaldd 183 ++ v+pg + e++ + ek + ++ ++s +dVA+ +++++++ 23726_FusoPortal_Gene_712 171 NVRVLSVCPGPTASNFfEKD--IQEKFGSSQKFMMSSEDVAKRIIKVIEN 218 99999999999988864654..5566666678899***********9998 PP >> GDP_Man_Dehyd GDP-mannose 4,6 dehydratase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 24.6 0.0 9.8e-09 1.4e-05 1 133 [. 5 139 .. 5 145 .. 0.74 Alignments for each domain: == domain 1 score: 24.6 bits; conditional E-value: 9.8e-09 EEETTTSHHHHHHHHHHHHCT-CCEEEEEEESS-TTHHHHHCTTT...TTGCCEEEEE--TT-HHHHHHHHHHHT-SEEEE--S CS GDP_Man_Dehyd 1 LiTGiaGfiGsylaelLlekgyevhgivrrsssfntgriehlykdpeaeleenlklvkgDltDssslerllaevqpdevynLAA 84 LiTG + iG l++ L +k + + r +++n+ + e l++ + ++ + +k DltD ++le++++++++d v+n A+ 23726_FusoPortal_Gene_712 5 LITGASSGIGEELTRNLANKSKKLFLLARSLDKLNLLKKE-LEE--K-FSSLECVCIKYDLTDINNLENIVENCDVDLVINCAG 84 9***********************9999888887754443.333..2.2345899****************************9 PP -....--HHHHTTSHHCHHHHHTHHHHHHHHHHHHC..T..TTSEEEEEEEGGGG CS GDP_Man_Dehyd 85 q....shVdvsfekPeytaevnvlGtlrLLeairal..glekkvrfyqaStsEvY 133 + + s ++ t++vn + l L + + + + + + + St+ +Y 23726_FusoPortal_Gene_712 85 FgkitDFSKLSDKEDLDTINVNFISPLILTKKFSEKflQKGQGTILNICSTAALY 139 7222134456667777888888888877777665541222345666666666666 PP >> NmrA NmrA-like family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 17.2 0.0 1.9e-06 0.0026 1 67 [. 4 76 .. 4 96 .. 0.90 2 ? 2.8 0.0 0.049 69 169 200 .. 191 222 .. 179 231 .. 0.84 Alignments for each domain: == domain 1 score: 17.2 bits; conditional E-value: 1.9e-06 EEESSTTSHHHHHHHHHHHHHT-EEEEEES-S......-SHHHHHHHHSTTEEEES-TT-HHHHHHHHTT-SE CS NmrA 1 IlviGaTGrqGksvvrallkaghkvralvRdt......kselakslkeaGvelvkgdlddkeslvealkgvdv 67 Il++Ga+ +G+ ++r l+++ k + l R + k+el +++++ + +k dl d ++l++++++ dv 23726_FusoPortal_Gene_712 4 ILITGASSGIGEELTRNLANKSKKLFLLARSLdklnllKKELEEKFSSLECVCIKYDLTDINNLENIVENCDV 76 89******************************************************************97554 PP == domain 2 score: 2.8 bits; conditional E-value: 0.049 EEESS-TTS-EE-HHHHHHHHHHHHHHHHHHH CS NmrA 169 llgpgnpkavlldaeeDiGafvikilkdprkl 200 ++ + ++++++++ eD+++ +ik++++++ + 23726_FusoPortal_Gene_712 191 IQEKFGSSQKFMMSSEDVAKRIIKVIENKKRF 222 5666667788999**************98876 PP >> 3Beta_HSD 3-beta hydroxysteroid dehydrogenase/isomerase family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 19.5 0.1 2.7e-07 0.00037 1 96 [. 5 106 .. 5 114 .. 0.72 Alignments for each domain: == domain 1 score: 19.5 bits; conditional E-value: 2.7e-07 3Beta_HSD 1 vvtGggGFlGrrivkeLlrlee..lkevrvlDlre..sselleelskvnvktkkvlegDvrdkekleealegvdv..vihtas. 77 ++tG++ +G+++ ++L + + r+lD+ + ++el e++s+++ +++ D++d ++le+ +e+ dv vi a+ 23726_FusoPortal_Gene_712 5 LITGASSGIGEELTRNLANKSKklFLLARSLDKLNllKKELEEKFSSLECV---CIKYDLTDINNLENIVENCDVdlVINCAGf 85 69****************9887544567999*7651145555566666665...999999999999999997665116666653 PP 3Beta_HSD 78 ..ivdvagkleedeilevnvk 96 i+d + ++ ++ +nv+ 23726_FusoPortal_Gene_712 86 gkITDFSKLSDKEDLDTINVN 106 334566666677778888876 PP >> T_Ag_DNA_bind Origin of replication binding protein # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.0 0.1 0.00035 0.48 15 48 .. 20 53 .. 16 71 .. 0.85 2 ? -0.6 0.0 1.4 2e+03 37 72 .. 188 224 .. 167 230 .. 0.74 Alignments for each domain: == domain 1 score: 11.0 bits; conditional E-value: 0.00035 T_Ag_DNA_bind 15 ifsnktleafliyttkeKlkeLykkllekfkvkf 48 + nk + fl+ +++Kl+ L+k+l+ekf+ 23726_FusoPortal_Gene_712 20 NLANKSKKLFLLARSLDKLNLLKKELEEKFSSLE 53 6789**************************9754 PP == domain 2 score: -0.6 bits; conditional E-value: 1.4 T_Ag_DNA_bind 37 ykkllekfkvkfkllakyeekav.lllitakkhRvsa 72 +k ++ekf + k +++ e+ a ++ ++++k+R s 23726_FusoPortal_Gene_712 188 EKDIQEKFGSSQKFMMSSEDVAKrIIKVIENKKRFSI 224 6778899988777776666555546666677778775 PP >> Glyco_tran_WecB Glycosyl transferase WecB/TagA/CpsF family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 9.6 0.0 0.00054 0.76 34 108 .. 15 87 .. 6 97 .. 0.86 2 ? -0.9 0.0 0.89 1.2e+03 44 74 .. 196 225 .. 159 232 .. 0.57 Alignments for each domain: == domain 1 score: 9.6 bits; conditional E-value: 0.00054 Glyco_tran_WecB 34 eelleraeeeglrvfllGgkpevlekaaerlrekypglriagthdgykeeeeeaiieainaskpdlllvglGapk 108 eel ++++++++++fll + + l+ ++++l+ek+ l+ + ++ ++ +++ + +++ + dl++ + G+ k 23726_FusoPortal_Gene_712 15 EELTRNLANKSKKLFLLARSLDKLNLLKKELEEKFSSLECVCIKYDL--TDINNLENIVENCDVDLVINCAGFGK 87 67788888889*****************************9887775..78888888889999999999988876 PP == domain 2 score: -0.9 bits; conditional E-value: 0.89 Glyco_tran_WecB 44 glrvfllGgkpevlekaaerlrekypglria 74 g + ++ ++++v++++++ +++k ++ i 23726_FusoPortal_Gene_712 196 GSSQKFMMSSEDVAKRIIKVIENK-KRFSII 225 344444555556666666666555.344454 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (258 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 864 (0.0516994); expected 334.2 (0.02) Passed bias filter: 526 (0.0314744); expected 334.2 (0.02) Passed Vit filter: 57 (0.00341072); expected 16.7 (0.001) Passed Fwd filter: 12 (0.000718047); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 12 [number of targets reported over threshold] # CPU time: 0.26u 0.23s 00:00:00.49 Elapsed: 00:00:00.30 # Mc/sec: 2500.05 // [ok]