# 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_41.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_41 [L=379] Description: # 52750 # 53886 # -1 # ID=1_41;partial=00;start_type=ATG;rbs_motif=AGGAGG;rbs_spacer=5-10bp;gc_cont=0.326 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.3e-46 157.1 0.6 4.7e-46 156.6 0.6 1.2 1 Acyl-CoA_dh_1 Acyl-CoA dehydrogenase, C-terminal domain 9.6e-27 92.9 0.2 1.7e-26 92.2 0.2 1.4 1 Acyl-CoA_dh_M Acyl-CoA dehydrogenase, middle domain 2.6e-22 79.5 0.0 5.6e-22 78.4 0.0 1.6 1 Acyl-CoA_dh_N Acyl-CoA dehydrogenase, N-terminal domain 3.1e-19 69.5 0.0 4.7e-19 68.9 0.0 1.3 1 Acyl-CoA_dh_2 Acyl-CoA dehydrogenase, C-terminal domain ------ inclusion threshold ------ 0.042 13.7 0.1 0.13 12.2 0.0 1.7 2 Mor Mor transcription activator family 0.12 12.1 0.0 0.2 11.4 0.0 1.3 1 YiiD_C Putative thioesterase (yiiD_Cterm) Domain annotation for each model (and alignments): >> Acyl-CoA_dh_1 Acyl-CoA dehydrogenase, C-terminal domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 156.6 0.6 1.7e-49 4.7e-46 1 149 [. 227 375 .. 227 376 .. 0.99 Alignments for each domain: == domain 1 score: 156.6 bits; conditional E-value: 1.7e-49 TBHHHHHTTTT--S-----..-SHHHHHHHHHHHHHHHTSB-CC-BGGGSHHHHHHHHHHHHHHHHHHHHHHHHHH--------H CS Acyl-CoA_dh_1 1 gkGfrvametlnheRlavaaaalglarraleealeyarrrkafgkplidhqlvreklaemaaeleaarllvyraaealdagkpdr 85 g+Gf++am+tl++ R +++ a+g+a+r ++ea++y ++r +fgkp+i++q++++k+a+m++++e ar +v +a ++d g p + 23726_FusoPortal_Gene_41 227 GEGFAIAMKTLDQARSWIGCIAVGIAQRGIQEAIAYGKERIQFGKPIIKNQALQFKIADMEIKTETARQMVAHALTKMDLGLPYG 311 68*********************************************************************************** PP HHHHHHHHHHHHHHHHHHHHHHHHHGGGTTSGGGHHHHHHHHTHHHHHSSS-HHHHHHHHHHHH CS Acyl-CoA_dh_1 86 aeasmaKlfaselaqevadlamqllGgaGyleeypverlyrdarvlrigeGtseiqrniiarrl 149 +e+++aK +a ++a+ev ++a+q++Gg+Gy++eypve+l+rd+++++i+eGt+eiqr +ia ++ 23726_FusoPortal_Gene_41 312 KESAIAKCYAGDIAMEVSSEAIQMFGGYGYSREYPVEKLIRDSKIFQIFEGTNEIQRIVIANNV 375 ************************************************************9876 PP >> Acyl-CoA_dh_M Acyl-CoA dehydrogenase, middle domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 92.2 0.2 6.1e-30 1.7e-26 1 97 [] 120 215 .. 120 215 .. 0.93 Alignments for each domain: == domain 1 score: 92.2 bits; conditional E-value: 6.1e-30 EEE---SS-SS-GGG---EEEEETT.EEEEEEEEEEEETTTT-SEEEEEEETS-------EEEEEEE..GGGEEE-----BSS-T CS Acyl-CoA_dh_M 1 afalTEpqaGsdvasvkttaaerdggeyvlnGeKiwitnaaiatlalvlartekpdrarglslflvprdapgvsvrrieeklGlr 85 af+lTEp+aGsd+ + +t +a +dg+eyvlnG+K++itn+a a+++ + a t+k+++ +g+s+f+v+++++g+s + e+k+G+r 23726_FusoPortal_Gene_41 120 AFCLTEPGAGSDASAGRT-TAVKDGDEYVLNGRKCFITNGAVASFYCITAITDKEKGLKGISMFFVEKGTKGLSTGNHEDKMGIR 203 79****************.66667777************************6555589*************************** PP TS-EEEEEEE.. CS Acyl-CoA_dh_M 86 glptgelvfedv 97 + +t+++v+ed+ 23726_FusoPortal_Gene_41 204 TSNTCDVVLEDC 215 **********97 PP >> Acyl-CoA_dh_N Acyl-CoA dehydrogenase, N-terminal domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 78.4 0.0 2e-25 5.6e-22 1 110 [. 6 114 .. 6 117 .. 0.95 Alignments for each domain: == domain 1 score: 78.4 bits; conditional E-value: 2e-25 -HHHHHHHHHHHHHHHHHTTTTHHHHHHHTSGHHHHHHHHHHTTCTSTTS-GGGTS-BS-HHHHHHHHHHHHHHTHHHHHHHHHH CS Acyl-CoA_dh_N 1 teeqealrdtvrefaeeelaphaaewdeegefprelwkklaelGllgitvpeeygGagldyleaaivaeelaradaglalllgvh 85 +ee ++l + v++f+++e++++++e+d++ge p+e+++k+ e+G+ ++ vpe+ygG gl+ ++ a ++ee+a adag+a+++++ 23726_FusoPortal_Gene_41 6 SEEAQDLLKDVKKFCDNEVREQCKEYDKSGEWPKEIYDKAIEQGYQALEVPEKYGGPGLKRVDIAALIEEMAIADAGFATTISA- 89 699*****************************************************************************9998. PP HHTTHHHHHHHS-HHHHHHHHHHHH CS Acyl-CoA_dh_N 86 sslvapailrfGteeqkrryLpkla 110 s+l++++++ G+eeqk+r ++ 23726_FusoPortal_Gene_41 90 SGLAMKPVFISGSEEQKQRMCDLVL 114 78999*************9876665 PP >> Acyl-CoA_dh_2 Acyl-CoA dehydrogenase, C-terminal domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 68.9 0.0 1.7e-22 4.7e-19 2 133 .. 243 364 .. 242 365 .. 0.98 Alignments for each domain: == domain 1 score: 68.9 bits; conditional E-value: 1.7e-22 HHHHHHHHHHHHHHHHHHHHHHHHHCTTT-------HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHTT----HHHHHHHHHH CS Acyl-CoA_dh_2 2 giaavalGaAraalaefvelarervrassgekladdpatqlalaeaaarldAAralleraaaraealaaagkevtvedraearra 86 +i+ +a+G+A++ ++e++++ +er++ +g+++ ++a+q+++a+++ + ++Ar+++++a ++ + g ++ +e+++a 23726_FusoPortal_Gene_41 243 WIGCIAVGIAQRGIQEAIAYGKERIQ--FGKPIIKNQALQFKIADMEIKTETARQMVAHALTK----MDLGLPY----GKESAIA 317 78999*********************..***********************************....99****9....******* PP HHHHHHHHHHHHHHHHHHCGGGGGBTT-HHHHHHHHHHHHCTSGGG- CS Acyl-CoA_dh_2 87 aalavelaiaavdalfeaaGgsalrkesplqRlwRDahaaaaHaaln 133 + +a ++a ++ +++++++Gg+++ e+p++ l RD ++ ++++++n 23726_FusoPortal_Gene_41 318 KCYAGDIAMEVSSEAIQMFGGYGYSREYPVEKLIRDSKIFQIFEGTN 364 ******************************************99988 PP >> Mor Mor transcription activator family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -3.8 0.0 4.3 1.2e+04 69 85 .. 18 34 .. 10 42 .. 0.61 2 ? 12.2 0.0 4.6e-05 0.13 22 74 .. 315 364 .. 309 369 .. 0.90 Alignments for each domain: == domain 1 score: -3.8 bits; conditional E-value: 4.3 H--SS-HHHHHHHHT-- CS Mor 69 efkGnnvaeLarkYgls 85 +f n+v+e ++Y s 23726_FusoPortal_Gene_41 18 KFCDNEVREQCKEYDKS 34 44455556666666555 PP == domain 2 score: 12.2 bits; conditional E-value: 4.6e-05 S.S-....TTSHHHHHHHHHHHH-SS------SHHHHHHHHHHHHHHH--SS- CS Mor 22 geidealaeqigeevadalaahwGGqqiyiPrgvklelelrdreifnefkGnn 74 + i + a i++ev+++ + +GG + r++ +e +rd +if+ f+G n 23726_FusoPortal_Gene_41 315 A-IAKCYAGDIAMEVSSEAIQMFGGYGYS--REYPVEKLIRDSKIFQIFEGTN 364 4.7788999****************9986..9*******************98 PP >> YiiD_C Putative thioesterase (yiiD_Cterm) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.4 0.0 7.3e-05 0.2 18 138 .. 200 322 .. 197 326 .. 0.79 Alignments for each domain: == domain 1 score: 11.4 bits; conditional E-value: 7.3e-05 YiiD_C 18 mgikveeydgkele.lsaplaaNiNdkgtmFaGSlytlatltgW.glvwLqlkeaglkgdivladakikYlkPvtedltaqvslp 100 mgi+ ++ + le +++p +a + ++g Fa ++ tl +W g + + + ++g+++ i + +i++ kP+ ++ q +++ 23726_FusoPortal_Gene_41 200 MGIRTSNTCDVVLEdCRIPASALVGKEGEGFAIAMKTLDQARSWiGCIAVGIAQRGIQEAIAYGKERIQFGKPIIKNQALQFKIA 284 666666666655543566677777778888998999999999998999************************************* PP YiiD_C 101 eeealeedleslarkkkariklevellsgeklaaeftG 138 + e +e +++ + +++l ++ +a + G 23726_FusoPortal_Gene_41 285 DMEIKTETARQMVAHALTKMDLGLPYGKESAIAKCYAG 322 88888888888888888777777766666666666666 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (379 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 481 (0.0287817); expected 334.2 (0.02) Passed bias filter: 426 (0.0254907); expected 334.2 (0.02) Passed Vit filter: 30 (0.00179512); expected 16.7 (0.001) Passed Fwd filter: 6 (0.000359023); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 6 [number of targets reported over threshold] # CPU time: 0.27u 0.24s 00:00:00.51 Elapsed: 00:00:00.32 # Mc/sec: 3443.02 // [ok]