# 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_1060.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1060 [L=310] Description: # 1168200 # 1169129 # 1 # ID=1_1060;partial=00;start_type=ATG;rbs_motif=GGA/GAG/AGG;rbs_spacer=5-10bp;gc_cont=0.223 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1e-17 64.3 0.6 7.4e-16 58.2 0.6 2.1 1 Epimerase NAD dependent epimerase/dehydratase family 1.3e-08 34.2 0.0 8.4e-06 25.0 0.0 2.3 2 RmlD_sub_bind RmlD substrate binding domain 2.2e-06 27.2 1.3 0.00045 19.6 0.1 2.3 2 GDP_Man_Dehyd GDP-mannose 4,6 dehydratase 2.1e-05 24.4 0.0 3e-05 23.9 0.0 1.2 1 NAD_binding_10 NAD(P)H-binding 9.1e-05 21.5 4.2 0.0024 16.8 4.2 2.1 1 3Beta_HSD 3-beta hydroxysteroid dehydrogenase/isomerase ------ inclusion threshold ------ 0.02 14.3 0.0 0.033 13.5 0.0 1.4 1 adh_short short chain dehydrogenase 0.091 12.2 0.7 1.6 8.2 0.0 2.0 2 NmrA NmrA-like family Domain annotation for each model (and alignments): >> Epimerase NAD dependent epimerase/dehydratase family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 58.2 0.6 3.1e-19 7.4e-16 1 240 [. 4 214 .. 4 215 .. 0.89 Alignments for each domain: == domain 1 score: 58.2 bits; conditional E-value: 3.1e-19 EEEETTTSHHHHHHHHHHHHCTSEEEEEES-STTTTCHHCHHCEEEE--TT-HHHHHHHHHHHT-SEEEE--S---HHHHHHS CS Epimerase 1 ilVtGatGfiGsalvkallekgyevigldrlssasntarledlelveldltdrsalerlfkeeqpdavihlAavggvaassed 83 il++G++ f+G+ ++k++lek+y++ +l+r ++++e + + ++d + ++e+++k+ ++d+++ + +a+ e+ 23726_FusoPortal_Gene_1060 4 ILIMGGNQFVGKEIAKNFLEKDYTIYVLNRG----TRKNIEGVFFLKVDRDNLIEMENILKDIEVDIIVDV------SAYTEE 76 89*****************************....66678888**************************99......232222 PP HHHHHHHHTHHHHHHHHHHHHHTCTEEEEEEEGGGGTSBSSSSBSTCSTTS----SSHHHHHHHHHHHHHHHHHHHHT-EEEE CS Epimerase 84 padfieknlqgtlnlleaarkagvkrllflsSsevYGkdaeqpkeettltgplkpnspYaaaKlageklveayaaqyglnaii 166 ++++l+ + k+g k+++ +sS++vY + + +p +e ++tg+ Ya K+++ek + + + y+++++i 23726_FusoPortal_Gene_1060 77 -----------QVDILHKVMKNGFKQYILISSASVYNNIECTPVNEGCQTGENLIWGDYAKNKYLAEKKTIENSNLYNFKYTI 148 ...........677999999*****************99999******999999999************************** PP EEESEEESTTS-TTSTTSHHHHHHHHHHCTSSSEEEETTST-EEEEEEHHHHHHHHHHHHHHHTT-STTEEEEE CS Epimerase 167 lrlfnlyGprngegevsrvipalirkilegkqevtllgdgtalRdflyveDvaeaillvlehvdagdkgeiyni 240 +r f +yG + + +++ ++ ++i+ + ++ + + ++ +f+yveD+a ai +e+ d+ ++i+ni 23726_FusoPortal_Gene_1060 149 FRAFYIYGIG-NNLDRE---NYFFSRIKYNL-PIFIP-SKNNIIQFGYVEDLALAIESSIENSDFY--NQIFNI 214 **********.333333...35678999999.88886.699*********************7774..569998 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.0 0.0 3.5e-09 8.4e-06 2 63 .. 3 75 .. 1 83 [. 0.96 2 ! 6.8 0.0 0.0012 3 125 251 .. 124 252 .. 121 254 .. 0.68 Alignments for each domain: == domain 1 score: 25.0 bits; conditional E-value: 3.5e-09 EEEEESTTSHHHHHHHHHCTTTSEEEEEE-TT...........T--TTSHHHHHHHHHHH--SEEEE------ CS RmlD_sub_bind 2 kilvtGanGqlGseLakqlaeeglevlaltre...........eldltdaeavaellkeakpdvvvnaAAyta 63 kil++G+n +G+e+ak++ e+++ ++ l+r ++d ++ +++++lk+++ d++v++ Ayt+ 23726_FusoPortal_Gene_1060 3 KILIMGGNQFVGKEIAKNFLEKDYTIYVLNRGtrkniegvfflKVDRDNLIEMENILKDIEVDIIVDVSAYTE 75 69***************************999999************************************96 PP == domain 2 score: 6.8 bits; conditional E-value: 0.0012 HHHHHHHHHHHHHH....H-SSEEEEEE-SEESSSSS-H..HHHHHHHHHH-SEEEEE.SS-EE--EEHHHHHHHHHHHHHHC CS RmlD_sub_bind 125 YGrtKlaGEqavla....agkkalivRtswvygekgknf..vktmlrlakereelkvv.aDqlgsPTsardlAevllqlleka 200 Y++ K++ E+ + ++ k++i+R ++yg g+n+ + ++++ k + ++ + + + + +++dlA ++ + +e++ 23726_FusoPortal_Gene_1060 124 YAKNKYLAEKKTIEnsnlYNFKYTIFRAFYIYGI-GNNLdrENYFFSRIKYNLPIFIPsKNNIIQFGYVEDLALAIESSIENS 205 6666666666555533334457777777777776.444411345666666666666641455666778888888888777777 PP HH.GGG-EEEE---BS-EEHHHHHHHHHHHHCHTTSEEEEE-CCCS--SS- CS RmlD_sub_bind 201 aeekelagvyHlanegetswlefArailkeagadseevkpistaeyptpAk 251 +++++ +++++ ++ efA+ i ++ a+++ +k ++t+e+ +A+ 23726_FusoPortal_Gene_1060 206 DF---YNQIFNISGDEYVTMSEFAE-ICGKVMAKKAVIKYVNTEENKIKAR 252 64...799****************7.5566666688999999999998886 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 ! 6.7 0.1 0.0016 3.8 1 69 [. 5 57 .. 5 80 .. 0.68 2 ! 19.6 0.1 1.9e-07 0.00045 122 332 .] 90 287 .. 77 287 .. 0.76 Alignments for each domain: == domain 1 score: 6.7 bits; conditional E-value: 0.0016 EEETTTSHHHHHHHHHHHHCT-CCEEEEEEESS-TTHHHHHCTTT...TTGCCEEEEE--TT-HHHHHH CS GDP_Man_Dehyd 1 LiTGiaGfiGsylaelLlekgyevhgivrrsssfntgriehlykdpeaeleenlklvkgDltDsssler 69 Li G+ f+G +a+ +lek y++++++r + + n ie +++ k+D + ++e+ 23726_FusoPortal_Gene_1060 5 LIMGGNQFVGKEIAKNFLEKDYTIYVLNRGTRK-N---IEG------------VFFLKVDRDNLIEMEN 57 899*********************998776544.3...332............2344444444444444 PP == domain 2 score: 19.6 bits; conditional E-value: 1.9e-07 SEEEEEEEGGGG-SBSSSSB-TTS-----SH...HHHHHHHHHHHHHHHHHHH---EEEEEE-EEESTTS-TTSCCCHHHHHH CS GDP_Man_Dehyd 122 vrfyqaStsEvYGkvqevpqtEttpfyPrsP...YaaaKlyaewivvnyreayglfavngilfNhesprrgetfvtrkitrai 201 ++++ +S++ vY + + +p +E ++ + Ya K ae +++ + y++ ++ + f ++ + + + 23726_FusoPortal_Gene_1060 90 KQYILISSASVYNNIECTPVNEGCQTGENLIwgdYAKNKYLAEKKTIENSNLYNFKYTIFRAFYIYGIGNNL--DRE--NYFF 168 468899****************999888653334*********************************99533..333..3478 PP HHHHTTS-C-EEETTTT-EEEEEEHHHHHHHHHHHHHCSS..---EEE-TT-EEECHHCHHCCHHCCCHHHHHH......... CS GDP_Man_Dehyd 202 akialgkqevlvlGnldakRdwghvkDyveamwlmlqqdk..pedyviatgetysvrefvelaflelgkeiewegkgeeevgr 282 ++i+ ++ + + + +++ +g+v+D + a+ +++ ++ ++i+ e +++ ef+e+ + ++k+ + + 23726_FusoPortal_Gene_1060 169 SRIKYNL--PIFIPSKNNIIQFGYVEDLALAIESSIENSDfyNQIFNISGDEYVTMSEFAEICGKVMAKKAVIKY-------- 241 8999988..445689999******************997522455777777888888888887777666554442........ PP ...HTGGCCCEECCG--SS--SBB-B-THHHHHHH----SS-HHHHHHHH CS GDP_Man_Dehyd 283 akesgkvvvevdpkyfrptevdlllgdatkakeeLgWkpkvsfeelvrem 332 ++++ +++ + p +l+gd++k++ + g++ ++s+ ++++++ 23726_FusoPortal_Gene_1060 242 ---VNTEENKIKARDWFPFREVNLFGDISKLE-NTGFRNTYSLIQGLEKT 287 ...222334445555555556689*******9.89********9999875 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 ! 23.9 0.0 1.3e-08 3e-05 1 143 [. 8 151 .. 8 160 .. 0.77 Alignments for each domain: == domain 1 score: 23.9 bits; conditional E-value: 1.3e-08 TTTSHHHHHHHHHHHHC.T-EEEEEETTHHHCHHHHCTTTCEEEE--TT.CCHHHHHHTT-SE.EEES---TCCHHHHHHHHH CS NAD_binding_10 1 GAtGkvGrllvkellargghevtalvRnpeklaeleaeeglevvegdvtdvedlaealkgada.Vvsalggggkdeegakniv 82 G++ vG++++k l++ ++++ +l R +++ eg+ +++d + ++++ lk++++ +++ +++ e + i+ 23726_FusoPortal_Gene_1060 8 GGNQFVGKEIAKNFLEK-DYTIYVLNRG--TRKNI---EGVFFLKVDRDNLIEMENILKDIEVdIIVDVSAYT---EEQVDIL 81 88889************.5776555555..55555...59*********************987888877644...3456666 PP HHHHHCT--EEEEEEETTTTSG.....CGCHHHHHHHHCCHHHHHHHHHHHHH.HHHS...SSEEEEEEE CS NAD_binding_10 83 daakaagvkRfvlvssigvgde.....vpgefgewnkellkaylaakraaeel.lras...gldyTivrp 143 + ++g k+++l+ss+ v+++ v++ ++ + +++y++ k+ ae+ +++s +++yTi r 23726_FusoPortal_Gene_1060 82 HKVMKNGFKQYILISSASVYNNiectpVNEGCQTGENLIWGDYAKNKYLAEKKtIENSnlyNFKYTIFRA 151 666788**************95323225555666777789********9999548888444567999996 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 ! 16.8 4.2 1e-06 0.0024 2 227 .. 6 198 .. 5 224 .. 0.71 Alignments for each domain: == domain 1 score: 16.8 bits; conditional E-value: 1e-06 3Beta_HSD 2 vtGggGFlGrrivkeLlrleelkevrvlDlresselleelskvnvktkkvlegDvrdkekleealegvdvvihtasivdvagk 84 ++Gg F+G++i k++l+ + ++ vl +++++ + +l++D + ++e+ l++++v i + dv 23726_FusoPortal_Gene_1060 6 IMGGNQFVGKEIAKNFLEKD--YTIYVLNRGTRKNI-------EGV--FFLKVDRDNLIEMENILKDIEVDIIV----DVSAY 73 78888888888888888887..47777655444444.......333..37777777888888888888877766....77664 PP 3Beta_HSD 85 leedeilevnvkGTqnvleacvkagvrvlvYtSsaevvgpnsrgdeilngdEetpyesthkdayaesKalaEklvlkangrtl 167 +ee+ +l +k+g ++ + Ssa+v ++n++ ++ +g + e+ ya+ K laEk ++ + l 23726_FusoPortal_Gene_1060 74 TEEQ----------VDILHKVMKNGFKQYILISSASV-YNNIECTPVNEGCQ--TGENLIWGDYAKNKYLAEKKTIENSN--L 141 4444..........35799999************997.67799999996643..345566777*********98776555..6 PP 3Beta_HSD 168 knggrllTcaLRpagiyGegdkllvpllk.kllksglalfkigekealsdrvYvgNvawah 227 n +T R+ +iyG g++l + + +k +l + i +k+++ ++ Yv+ +a a+ 23726_FusoPortal_Gene_1060 142 YNF--KYT-IFRAFYIYGIGNNLDRENYFfSRIKYNLP-IFIPSKNNIIQFGYVEDLALAI 198 665..345.69**********98877665144455544.45667778899********997 PP >> adh_short short chain dehydrogenase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.5 0.0 1.4e-05 0.033 1 71 [. 2 62 .. 2 106 .. 0.69 Alignments for each domain: == domain 1 score: 13.5 bits; conditional E-value: 1.4e-05 SEEEEETTTSHHHHHHHHHHHHTT-EEEEEESCHHHHHHHHHHHHHTTSEEEEEE-TTTSHHHHHHHHHHH CS adh_short 1 kvvlvtGassGIGraiakelakeGakvvvvdrseekleevakelkeegakalaiqvDvtdreevealveka 71 k +l+ G+ + +G++iak++++++++++v++r +k +++ +++++vD + e+e++++++ 23726_FusoPortal_Gene_1060 2 KKILIMGGNQFVGKEIAKNFLEKDYTIYVLNRGTRK------NIEG----VFFLKVDRDNLIEMENILKDI 62 689**************************9998554......3444....566666666666666666665 PP >> NmrA NmrA-like family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 8.2 0.0 0.00065 1.6 1 99 [. 4 94 .. 4 111 .. 0.91 2 ? 2.6 0.1 0.033 79 166 231 .. 174 241 .. 116 243 .. 0.67 Alignments for each domain: == domain 1 score: 8.2 bits; conditional E-value: 0.00065 EEESSTTSHHHHHHHHHHHHHT-EEEEEES-S-SHHHHHHHHSTTEEEES-TT-HHHHHHHHTT-SEEEE---GGS-HHHHHHHH CS NmrA 1 IlviGaTGrqGksvvrallkaghkvralvRdtkselakslkeaGvelvkgdlddkeslvealkgvdvvfsvtgfyskeiedgkkl 85 Il++G+ Gk +++ +l++++++ +l R t+++ +Gv ++k d d+ ++++lk+ +v ++v s+ +e+++++ 23726_FusoPortal_Gene_1060 4 ILIMGGNQFVGKEIAKNFLEKDYTIYVLNRGTRKNI------EGVFFLKVDRDNLIEMENILKDIEVDIIVDV--SAYTEEQVDI 80 8999999999*********************99875......8***************************999..99******** PP HHHHHHHS-SEEEE CS NmrA 86 vdAakeaGvkhfiy 99 + + + G k++i 23726_FusoPortal_Gene_1060 81 LHKVMKNGFKQYIL 94 ********999985 PP == domain 2 score: 2.6 bits; conditional E-value: 0.033 -EEEEESS-TTS-EE-HHHHHHHHHHHHHHHHHHHTTEEEEE-SEEE-HHHHHHHHHHHHTS-..EEE CS NmrA 166 kvtllgpgnpkavlldaeeDiGafvikilkdprklkgkriklagntlsgneivelfekkigkt..vkY 231 + + p++ ++ ++ +eD + ++ + +++ + ++ +++++ e +e+ k++ k+ +kY 23726_FusoPortal_Gene_1060 174 NLPIFIPSKNNIIQFGYVEDLALAIESSIENSDFYNQIFNISGDEYVTMSEFAEICGKVMAKKavIKY 241 44556666666666688889998888888887777444444457899************999855566 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (310 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 728 (0.0435615); expected 334.2 (0.02) Passed bias filter: 360 (0.0215414); expected 334.2 (0.02) Passed Vit filter: 40 (0.00239349); expected 16.7 (0.001) Passed Fwd filter: 7 (0.000418861); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 7 [number of targets reported over threshold] # CPU time: 0.33u 0.27s 00:00:00.60 Elapsed: 00:00:00.36 # Mc/sec: 2503.28 // [ok]