# 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_1603.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1603 [L=392] Description: # 1722932 # 1724107 # -1 # ID=1_1603;partial=00;start_type=ATG;rbs_motif=AGGAG;rbs_spacer=5-10bp;gc_cont=0.318 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1.6e-46 158.4 5.6 4.9e-44 150.3 1.4 2.1 2 ETF Electron transfer flavoprotein domain 4.8e-35 119.3 1.0 1.7e-34 117.6 0.7 2.0 2 ETF_alpha Electron transfer flavoprotein FAD-binding d ------ inclusion threshold ------ 0.042 14.5 0.1 0.11 13.1 0.1 1.7 1 DXP_reductoisom 1-deoxy-D-xylulose 5-phosphate reductoisomer Domain annotation for each model (and alignments): >> ETF Electron transfer flavoprotein domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 10.4 0.3 7.4e-05 0.41 2 48 .. 10 55 .. 9 70 .. 0.77 2 ! 150.3 1.4 8.8e-48 4.9e-44 11 180 .. 70 243 .. 63 245 .. 0.92 Alignments for each domain: == domain 1 score: 10.4 bits; conditional E-value: 7.4e-05 EEEE-EE.TTEE-HHHHHHHHHHHHHHHHHT--EEEEEEES....-- CS ETF 2 lvlaehgvdgklnpvdlealeaArrlkekaggevtavvvGppeaeea 48 lv+ae++ dg+l++v le+l++A +l+ + + ++++ +G++ a+ a 23726_FusoPortal_Gene_1603 10 LVFAEQR-DGILQNVGLELLGKATELAYEINKQIALKDAGDELADFA 55 89****9.******************998778888777777666555 PP == domain 2 score: 150.3 bits; conditional E-value: 8.8e-48 TEE-HHHHHHHHHHHHHHHH..HT--EEEEEEES....--..THHHHGGGT-SEEEEEE-GGGSS--HHHHHHHHHHHHHHH- CS ETF 11 gklnpvdlealeaArrlkek..aggevtavvvGppeaeea..lrealaamGadkvilvddpalagydaeatakalaalikked 89 ++l + d+e+ ++ +k++ + ++vta++vG+ + l++ l + Gadkv++vd p+l+ yd+ea+a++l a i++++ 23726_FusoPortal_Gene_1603 70 ATLEEDDQELKNKVADIKKNnpDAAKVTALLVGH----NIknLAQDLIKAGADKVLVVDKPELEVYDTEAYAQVLTATINDQK 148 5667778888888888888886559*********....555699*************************************** PP -SEEEEESSHHHHHHHHHHHHHTT---EEEESEEEEE.EE...EEEEEECTTTTEEEEE--S.SSSEEEEE-TTSS----S-T CS ETF 90 pdlvlagataigkqvaprlAalLgvplitdvteleveggk...ltvtRpieggnglaevvea.dlpavltvrpgafepaaima 168 p++vl+gat+ g+++apr++++ ++l++d+t+le ++k l +tRp++ggn +a++v++ ++p+++tvrpg++++++++ 23726_FusoPortal_Gene_1603 149 PEIVLFGATTLGRDLAPRVSSRIATGLTADCTKLELLKDKerqLGMTRPAFGGNLMATIVSPdHRPQMATVRPGVMKKLPKSD 231 *********************************9996566799**************************************** PP T---EEEE--HH CS ETF 169 akkkeveeveva 180 +++ev+e++v+ 23726_FusoPortal_Gene_1603 232 DRTGEVVEFPVT 243 ********9986 PP >> ETF_alpha Electron transfer flavoprotein FAD-binding domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -1.5 0.0 0.44 2.5e+03 30 43 .. 186 199 .. 162 206 .. 0.67 2 ! 117.6 0.7 3e-38 1.7e-34 1 84 [] 264 347 .. 264 347 .. 0.99 Alignments for each domain: == domain 1 score: -1.5 bits; conditional E-value: 0.44 ETF_alpha 30 ealgaavgasRavv 43 + + ++g++R++ 23726_FusoPortal_Gene_1603 186 KDKERQLGMTRPAF 199 33355677777765 PP == domain 2 score: 117.6 bits; conditional E-value: 3e-38 ETF_alpha 1 vdLeeakvvvsgGrGlkseenlelleelaealgaavgasRavvdagwipeerqvGqtGktvaPklyiavGisGaiqhlaGmke 83 vd++eak++vsgGrG++s++n+elle+la+++g+ v+ sRa vdag +p++rqvGqtGktv+P++y+a+GisGaiqh+aGm+e 23726_FusoPortal_Gene_1603 264 VDISEAKILVSGGRGVGSKQNFELLEDLATEIGGIVSSSRAQVDAGNMPHDRQVGQTGKTVRPEVYFACGISGAIQHVAGMEE 346 699******************************************************************************99 PP ETF_alpha 84 s 84 s 23726_FusoPortal_Gene_1603 347 S 347 7 PP >> DXP_reductoisom 1-deoxy-D-xylulose 5-phosphate reductoisomerase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.1 0.1 2e-05 0.11 16 66 .. 85 135 .. 80 161 .. 0.82 Alignments for each domain: == domain 1 score: 13.1 bits; conditional E-value: 2e-05 HHHHHCTTTEEEEEEEESSTHHHHHHHHHHHT-SEEEESSHHHHHHHHHHH CS DXP_reductoisom 16 dvieenpdkfevvalaagsnvellaeqakefkpklvavadeeaaeelkeal 66 d+ ++npd+ +v+al +g+n+++la+ + + ++ v v d+ + e ++ 23726_FusoPortal_Gene_1603 85 DIKKNNPDAAKVTALLVGHNIKNLAQDLIKAGADKVLVVDKPELEVYDTEA 135 678999*************************99999999998877666554 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (392 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 644 (0.0385352); expected 334.2 (0.02) Passed bias filter: 458 (0.0274055); expected 334.2 (0.02) Passed Vit filter: 39 (0.00233365); 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.32u 0.24s 00:00:00.56 Elapsed: 00:00:00.34 # Mc/sec: 3351.64 // [ok]