# 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_1881.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1881 [L=336] Description: # 2016907 # 2017914 # -1 # ID=1_1881;partial=00;start_type=ATG;rbs_motif=GGAGG;rbs_spacer=5-10bp;gc_cont=0.276 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 5e-58 195.3 0.0 8.5e-58 194.5 0.0 1.4 1 2-Hacid_dh_C D-isomer specific 2-hydroxyacid dehydrogenase, 3.6e-26 91.3 0.1 5.1e-26 90.8 0.1 1.1 1 2-Hacid_dh D-isomer specific 2-hydroxyacid dehydrogenase, 5.2e-05 23.3 0.3 0.00012 22.1 0.1 1.7 2 AdoHcyase_NAD S-adenosyl-L-homocysteine hydrolase, NAD bindi 8.1e-05 22.7 0.1 0.00017 21.7 0.1 1.5 1 NAD_binding_2 NAD binding domain of 6-phosphogluconate dehyd 0.0018 17.8 0.1 0.0041 16.6 0.0 1.6 2 IlvN Acetohydroxy acid isomeroreductase, NADPH-bind 0.0059 16.4 0.0 1.4 8.7 0.0 2.3 2 ADH_zinc_N Zinc-binding dehydrogenase ------ inclusion threshold ------ 0.013 15.4 0.0 0.046 13.7 0.0 2.0 1 Shikimate_DH Shikimate / quinate 5-dehydrogenase 0.021 14.4 0.0 2.3 7.8 0.0 2.2 2 CHASE2 CHASE2 domain Domain annotation for each model (and alignments): >> 2-Hacid_dh_C D-isomer specific 2-hydroxyacid dehydrogenase, NAD binding domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 194.5 0.0 4.1e-61 8.5e-58 1 178 [] 112 300 .. 112 300 .. 0.95 Alignments for each domain: == domain 1 score: 194.5 bits; conditional E-value: 4.1e-61 HHHHHHHHTTHHHHHHHHHTTCT-TTGCGSB--GTTSEEEEES-SHHHHHHHHHHHHTT-EEEEESSSHHHHHCHHHHTSEEC CS 2-Hacid_dh_C 1 lalllallrrlaeadeevregewksekallgkelsgktvGiiGlGrIGqavakrlkafgmkviaydrskkkeeeeeelgvelv 83 ++l+la++r++++a+ ++reg++ s + l+g +l++kt GiiG+G+IGq ++k+l++f+mkviayd +++++ ++elg+e+v 23726_FusoPortal_Gene_1881 112 VGLILAVNRKINKAYVRTREGNF-SINGLMGIDLYEKTAGIIGTGKIGQILIKILRGFDMKVIAYDLFPNQKV-ADELGFEYV 192 589********************.99********************************************888.99******* PP SHHHHHHH-SEEEE-S--STTTTTSBSHHHHHHSTTTEEEEESS-GGGB-HHHHHHHHHCTSSCEEEES--TTS......... CS 2-Hacid_dh_C 84 sleellaesDivslhlpltpetrhlinaeelakmkkgailiNtaRGglvdeeaLleaLksgkiagaalDvfeee......... 157 sl+el+a+sDi+sl++plt++t+ +in++++ kmk+g+il+Nt+RG l+d ++L+eaLk++ki+++alDv+eee 23726_FusoPortal_Gene_1881 193 SLDELYANSDIISLNCPLTKDTKYMINRRSMLKMKDGVILVNTGRGMLIDSADLVEALKDKKIGAVALDVYEEEenyffedks 275 **************************************************************************999999988 PP .SSGCHH...HHHHSTTEEE-SS-T CS 2-Hacid_dh_C 158 .plpeds...pllelpnviltPHia 178 ++ ed+ +ll++ nv++t+H a 23726_FusoPortal_Gene_1881 276 tQVIEDDilgRLLSFYNVLITSHQA 300 855555433388***********87 PP >> 2-Hacid_dh D-isomer specific 2-hydroxyacid dehydrogenase, catalytic domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 90.8 0.1 2.4e-29 5.1e-26 2 133 .. 7 331 .. 6 332 .. 0.97 Alignments for each domain: == domain 1 score: 90.8 bits; conditional E-value: 2.4e-29 EEEST.CHHHHHHHHH....TTE.EEEEESCSSHHHCHHHHTTESEEEE-TTS-BSHHHHHHH..TT--EEEESSSSCTTB-H CS xxxxxxxxxxxxxxxx....xxx.xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx..xxxxxxxxxxxxxxxxxx RF 2-Hacid_dh 2 lileplreeelellke....egv.evevkdellteellekakdadalivrsntkvtaevleal..pkLkviaragvGvDnvDl 77 +++ +++++++e++k+ e+++ + ltee++ k++d+++ + n+++++e+++ + +++k+ a +++G++nv l 23726_FusoPortal_Gene_1881 7 IFF-DIKDYDKEFFKKygadY-NfEMTFLKVRLTEETANLTKGYDVVCGFANDNINKETIDIMaeNGIKLLAMRCAGFNNVSL 87 455.589*********99642.24888888***************************************************** PP HHHHHTT-EEEE-TTTTHHHHHHHH.......................................................... CS xxxxxxxxxxxxxxxxxxxxxxxxx.......................................................... RF 2-Hacid_dh 78 daakerGilVtnvpgystesvAElt.......................................................... 102 + +e+ ++V++vp+ys++++AE+t 23726_FusoPortal_Gene_1881 88 KDVNEK-FKVVRVPAYSPHAIAEYTvglilavnrkinkayvrtregnfsinglmgidlyektagiigtgkigqilikilrgfd 169 ******.**************************************************************************** PP ................................................................................... CS ................................................................................... RF 2-Hacid_dh 103 ................................................................................... 102 23726_FusoPortal_Gene_1881 170 mkviaydlfpnqkvadelgfeyvsldelyansdiislncpltkdtkyminrrsmlkmkdgvilvntgrgmlidsadlvealkd 252 *********************************************************************************** PP ................................................T-BHHHHHHHHHHHHHHHHHHHTTCCGTTBC CS ................................................xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF 2-Hacid_dh 103 ................................................aaTeeaqeniaeeaaenlvaflkgespanav 133 ++T+ea+ +i ++++n+++f++g+++ n+v 23726_FusoPortal_Gene_1881 253 kkigavaldvyeeeenyffedkstqvieddilgrllsfynvlitshqaYFTKEAVGAITVTTLNNIKDFVEGRPLVNEV 331 ***********************************************999***********************999988 PP >> AdoHcyase_NAD S-adenosyl-L-homocysteine hydrolase, NAD binding domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -2.2 0.0 1.7 3.5e+03 58 75 .. 73 90 .. 38 101 .. 0.54 2 ! 22.1 0.1 5.7e-08 0.00012 23 145 .. 146 276 .. 135 293 .. 0.73 Alignments for each domain: == domain 1 score: -2.2 bits; conditional E-value: 1.7 HHHHHHCTT-EE--HHHC CS AdoHcyase_NAD 58 calqaameGfevvtleev 75 l Gf+ v l++v 23726_FusoPortal_Gene_1881 73 KLLAMRCAGFNNVSLKDV 90 223333344444444444 PP == domain 2 score: 22.1 bits; conditional E-value: 5.7e-08 TSEEEEE--SHHHHHHHHHHHHCT-EEEEE-S-HHHHHHHHCTT-EE--HHHCTTT-SEEE.E...-SSSSSSB-HHHHCCS- CS AdoHcyase_NAD 23 GkvavvaGyGdvGkGcaaslkglGarvivtevdPicalqaameGfevvtleevvkkadivv.t...ttGnkdiitlehlkkmk 101 k+a ++G G +G+ + l+g+ +vi ++ P + a Gfe v l+e+ +++di+ t +k +i++ + kmk 23726_FusoPortal_Gene_1881 146 EKTAGIIGTGKIGQILIKILRGFDMKVIAYDLFPNQKV-ADELGFEYVSLDELYANSDIISlNcplTKDTKYMINRRSMLKMK 227 68999*****************************7766.5556***************985232224566779********** PP TTEEEEE-SSTTTTBT....HHHHHTT..T-EEEEEETTEEEEE-TTS-E CS AdoHcyase_NAD 102 ddaivcniGhfdseid....vaalean..alkkeevkpqvdrvtlkdgkk 145 d +i+ n G id v+al+ + ++ +v + +++ ++d+ + 23726_FusoPortal_Gene_1881 228 DGVILVNTGRGM-LIDsadlVEALKDKkiGAVALDVYEEEENYFFEDKST 276 ********9743.3440111444544422334445556666666666554 PP >> NAD_binding_2 NAD binding domain of 6-phosphogluconate dehydrogenase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 21.7 0.1 8.2e-08 0.00017 3 113 .. 150 259 .. 148 263 .. 0.91 Alignments for each domain: == domain 1 score: 21.7 bits; conditional E-value: 8.2e-08 EEE--SHHHHHHHHHHHHTT-EEEEE-SSHHHHHHHHHTTEEE-SSHHHHHHHCSEEEE-SSSHHHHHHHHHS.HCCC--TT- CS NAD_binding_2 3 gfiGlGvMGsnmalnllkagykvavynrtkekveelvaegakaaesieelvasldvvitmvkagkavdevieg.llealkkgd 84 g+iG+G++G+ +++ l ++kv++y++ ++ +++++e + s+ el a+ d++ + + +k+++ +i++ + ++k g 23726_FusoPortal_Gene_1881 150 GIIGTGKIGQILIKILRGFDMKVIAYDLFPN--QKVADELGFEYVSLDELYANSDIISLNCPLTKDTKYMINRrSMLKMKDGV 230 99*************************9854..577777677888999*************************98889***** PP EEEE-S---HHHHHHHHHHHHHTT-EEEE CS NAD_binding_2 85 ilidgstsspedtrelakelkekglkfld 113 il++++ d+ +l ++lk k++ + 23726_FusoPortal_Gene_1881 231 ILVNTGRGMLIDSADLVEALKDKKIGAVA 259 ***********************998765 PP >> IlvN Acetohydroxy acid isomeroreductase, NADPH-binding domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -2.6 0.0 1.6 3.4e+03 66 114 .. 53 74 .. 38 96 .. 0.52 2 ! 16.6 0.0 2e-06 0.0041 2 92 .. 144 233 .. 143 259 .. 0.79 Alignments for each domain: == domain 1 score: -2.6 bits; conditional E-value: 1.6 E-S-HHHHHHHHHHHTGGG--TT-EEEESSSHHHHTTSS---TTSEEEE CS IlvN 66 iliPDelqkevyeeeiepnlkegkalafahGfnihfglikppkdvdvvl 114 + D ni+ ++i++ ++ + l 23726_FusoPortal_Gene_1881 53 GFAND---------------------------NINKETIDIMAENGIKL 74 22222...........................34444444444444444 PP == domain 2 score: 16.6 bits; conditional E-value: 2e-06 HHTS-EEEES-SHHHHHHHHHHHHTT--EEEEE-TT-HHHHHHHHTT-EEEEHHHHHHT-SEEEE-S-HHHH.HHHHHHHTGG CS IlvN 2 lkgkkiaviGyGsqGhaqalnlrdsgldvvvglregsksvkkAkeegfevltvaeavkkadvvmiliPDelq.kevyeeeiep 83 l +k+ +iG G G+ + lr ++v+ ++ + k A e gfe ++++e+ +++d++ + P + k +++++ + 23726_FusoPortal_Gene_1881 144 LYEKTAGIIGTGKIGQILIKILRGFDMKVI--AYDLFPNQKVADELGFEYVSLDELYANSDIISLNCPLTKDtKYMINRRSML 224 568999*****************7777765..4555556799************************99765514456666677 PP G--TT-EEE CS IlvN 84 nlkegkala 92 ++k+g +l+ 23726_FusoPortal_Gene_1881 225 KMKDGVILV 233 777777654 PP >> ADH_zinc_N Zinc-binding dehydrogenase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 8.7 0.0 0.00067 1.4 15 83 .. 4 71 .. 1 78 [. 0.85 2 ! 5.3 0.0 0.0078 16 2 39 .. 157 193 .. 157 207 .. 0.92 Alignments for each domain: == domain 1 score: 8.7 bits; conditional E-value: 0.00067 SEEEEEESSHHHHHHHHHCTESEEEETTTSC..HHHHHHHHCTTSSEEEEEESSCSCHHHHHHHHHEEEEE CS ADH_zinc_N 15 arviavdsseeklelakelGAdhvinskded..fvkaikeltggkgvdvvfecvGapatfdqalsllapgG 83 ++i d+++ e+ k++GAd+ + + + ++++ ++lt+g +dvv+ ++ + ++ ++++a++G 23726_FusoPortal_Gene_1881 4 TKIIFFDIKDYDKEFFKKYGADYNFEMTFLKvrLTEETANLTKG--YDVVCGFAN-DNINKETIDIMAENG 71 5899999***************9999987777788999999887..****99988.777777777777776 PP == domain 2 score: 5.3 bits; conditional E-value: 0.0078 HHHHHHHHHHHTTSEEEEEESSHHHHHHHHHCTESEEE CS ADH_zinc_N 2 vGlaavqlAkalGarviavdsseeklelakelGAdhvi 39 +G++++++++ + +via d + + + a elG ++v 23726_FusoPortal_Gene_1881 157 IGQILIKILRGFDMKVIAYDLFP-NQKVADELGFEYVS 193 8********************99.99*******99886 PP >> Shikimate_DH Shikimate / quinate 5-dehydrogenase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.7 0.0 2.2e-05 0.046 9 79 .. 143 204 .. 137 244 .. 0.76 Alignments for each domain: == domain 1 score: 13.7 bits; conditional E-value: 2.2e-05 ..TT-EEEEES-SHHHHHHHHHHHCCT-CEEEEEESSHHHHHHHHHHHT.ET-EEEEGCGHHHHHHT-SEE CS Shikimate_DH 9 slkekkvlliGaGemaelvakhLlakgvkkvvvaNRtlerakelaeelkgeeiealkleelkellaeadvv 79 +l ek++ +iG+G++++ ++k L +k v+ +l +++a+el+ e +++l+el a++d++ 23726_FusoPortal_Gene_1881 143 DLYEKTAGIIGTGKIGQILIKILRGFDMK---VIAYDLFPNQKVADELGFE---YVSLDELY---ANSDII 204 68899******************976666...8888999999999999555...99999988...566666 PP >> CHASE2 CHASE2 domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 4.6 0.0 0.01 21 67 177 .. 3 208 .. 1 239 [. 0.44 2 ? 7.8 0.0 0.0011 2.3 57 95 .. 244 292 .. 242 328 .. 0.71 Alignments for each domain: == domain 1 score: 4.6 bits; conditional E-value: 0.01 CHASE2 67 gakaiglDilfsepeeydealaeal.ak........................gpnvv................p.......aa 101 + k+i +Di + yd+++ + a+ v + a 23726_FusoPortal_Gene_1881 3 KTKIIFFDI-----KDYDKEFFKKYgADynfemtflkvrlteetanltkgydV---VcgfandninketidimAengikllAM 77 678999999.....335555555552223333333334444443333333220...133333333333333331333444466 PP CHASE2 102 aglGfvnv.vpDaDg...vvRrpllvese...lipslalalalrylkae............................klrlpv 149 + Gf+nv +D + vvR p + + + l la+ ++ + + 23726_FusoPortal_Gene_1881 78 RCAGFNNVsLKDVNEkfkVVRVPAYSPHAiaeYTVGLILAVNRKIN--KayvrtregnfsinglmgidlyektagiiG----T 154 8889999965677776666777655333366643355555555333..122344456666666655555543222220....1 PP CHASE2 150 desgglllnyrgg.sgt.................fptvsladvle........v 177 + g++l++ +g ++ f++vsl ++++ + 23726_FusoPortal_Gene_1881 155 GKIGQILIKILRGfDMKviaydlfpnqkvadelgFEYVSLDELYAnsdiislnC 208 112446664444432111222222233344444466677766653344444441 PP == domain 2 score: 7.8 bits; conditional E-value: 0.0011 CHASE2 57 arllerLeaagakaiglDilfsepeey...........dealaealakgp 95 a+l+e L++ + a++lD+ ++e+e+y d+ l + l+ ++ 23726_FusoPortal_Gene_1881 244 ADLVEALKDKKIGAVALDV-YEEEENYffedkstqvieDDILGRLLSFYN 292 89*****************.666766788777777777666666664443 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (336 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 634 (0.0379368); expected 334.2 (0.02) Passed bias filter: 470 (0.0281235); expected 334.2 (0.02) Passed Vit filter: 48 (0.00287219); expected 16.7 (0.001) Passed Fwd filter: 8 (0.000478698); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 8 [number of targets reported over threshold] # CPU time: 0.27u 0.23s 00:00:00.50 Elapsed: 00:00:00.30 # Mc/sec: 3255.88 // [ok]