# 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_1356.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1356 [L=259] Description: # 1463816 # 1464592 # 1 # ID=1_1356;partial=00;start_type=ATG;rbs_motif=AGGAGG;rbs_spacer=5-10bp;gc_cont=0.310 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 8.6e-70 234.8 0.3 1.1e-69 234.5 0.3 1.0 1 ECH_1 Enoyl-CoA hydratase/isomerase 1.5e-37 129.8 0.3 6.2e-35 121.2 0.2 2.0 2 ECH_2 Enoyl-CoA hydratase/isomerase Domain annotation for each model (and alignments): >> ECH_1 Enoyl-CoA hydratase/isomerase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 234.5 0.3 1.3e-73 1.1e-69 3 245 .. 9 251 .. 7 257 .. 0.98 Alignments for each domain: == domain 1 score: 234.5 bits; conditional E-value: 1.3e-73 ETTEEEEEE--CGCTTEE-HHHHHHHHHHHHHHHHSTTECEEEEEE.BTSEEEEEB-HHHHHHHHCHHHHHHHCCHHHHHHHH CS ECH_1 3 eagvaviklnrPeavNalsaelltEliqalekleddksvravvLtG.gdkaFaaGadlkemaeeekakqallkeevlelwtkl 84 e+ + ++++rPea+Nal+ ++l+El + ++ ++ + + r+v LtG g k F aGad+ em + + +++ + +e ++k+ 23726_FusoPortal_Gene_1356 9 ENFIGIVTIERPEALNALNSQVLDELSATFAGINLE-TTRVVLLTGsGSKSFVAGADIAEMSTLNSDEGTKFGYKGNEIFRKI 90 567889*****************************9.9**************************99999************** PP HHSSSEEEEEE-SEEETHHHHHHHHSSEEEEETTTEEE-GGGGGTS--CTTCHHHHHHHHHHHHHHHHHHHT-EEEHHHHHHC CS ECH_1 85 ytvpkpviaavnGaaigGGlelalacDiviaadnasfglnelklGlipgaggtqilarkvGkskAeeilllgkkisaqealkv 167 t p+pvia +nG+a+gGG+ela+ cD++i+++na fg++e+ lG++pg+ggtq+lar +G kA+e++++++ i+a+eal++ 23726_FusoPortal_Gene_1356 91 ETFPLPVIAVINGFALGGGCELAMSCDFRICSENAIFGQPEVGLGITPGFGGTQRLARLIGLGKAKEMIYTANTIKADEALNI 173 *********************************************************************************** PP TSSSEEESGGGHHHHHHHHHHHHHHSCHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHCTSHHHHHHHHHHHTCSH CS ECH_1 168 GlvdkvvpeeqlleealelaeklaekspaalrqlkealnakledalaeqleaeaklfestfskddikeGlqaflEkrk 245 Glv++v p+e+l+eea++la k+a+++p a+r+ k+a+n+ ++ +++ +e klf s f+++d k G++afl k k 23726_FusoPortal_Gene_1356 174 GLVNHVYPQETLMEEAMKLAGKIAKNAPFAVRACKKAINEGIDTDMDRAIIIEEKLFGSCFATEDQKVGMKAFLDKVK 251 ***************************************************************************987 PP >> ECH_2 Enoyl-CoA hydratase/isomerase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 121.2 0.2 7.4e-39 6.2e-35 1 180 [. 12 188 .. 12 195 .. 0.94 2 ! 8.3 0.0 0.00016 1.3 249 311 .. 189 249 .. 185 258 .. 0.89 Alignments for each domain: == domain 1 score: 121.2 bits; conditional E-value: 7.4e-39 .EEEEE--GGGTT-B-HHHHHHHHHHHHHHCC-TT-SEEEEEESSSSECE--B--HHHHHHHHCS-THHHHHHHHHHHHHHHH CS ECH_2 1 vgvitLnrpkalNaLslemvralkealkewekddevklvvlkgagekaFCAGgDvralyesakegkkeeeeffreeYklnyli 83 +g++t++rp+alNaL+ +++++l++++ + + +++++v+l+g+g+k+F AG D++++ + ++++ ++ f + ++ i 23726_FusoPortal_Gene_1356 12 IGIVTIERPEALNALNSQVLDELSATFAGI-NLETTRVVLLTGSGSKSFVAGADIAEMSTLNSDEGTK---FGYKGNEIFRKI 90 589*************************98.67789*************************9997755...************ PP HC-SS-EEEEE-BEEETHHHHHHTTSSEEEE-TT-EEE--GGGGT----TTHHHHHHCSST.THHHHHHHH---B-HHHHHHT CS ECH_2 84 atypkpyvalmdGivmGgGvGlsvhgsfrvvtertvfAmPEtaIGlfPDvGgsyflsrlpg.elGlylalTGarlkgadalaa 165 +t+p p++a+++G+++GgG+ l++++ fr+ +e+++f PE+++G++P Gg+ l+rl g ++ + T++++k+ +al++ 23726_FusoPortal_Gene_1356 91 ETFPLPVIAVINGFALGGGCELAMSCDFRICSENAIFGQPEVGLGITPGFGGTQRLARLIGlGKAKEMIYTANTIKADEALNI 173 ***********************************************************887889999*************** PP TS--EE--GGGHHHH CS ECH_2 166 glAthyvpsekleal 180 gl +h+ p+e l + 23726_FusoPortal_Gene_1356 174 GLVNHVYPQETLMEE 188 *********998665 PP == domain 2 score: 8.3 bits; conditional E-value: 0.00016 HHHHHHHHHHS-HHHHHHHHHHHHHHTTT--HHHHHHHHHHHHHHHHHHSSHHHHHHHHTTS- CS ECH_2 249 akktletlkkrsptslkvtleqlrrgrklsslaevlrmElrlaqrllaekDfveGvrallvdk 311 a k + +++k++p++++++ +++++g ++ +++ ++ E +l +++a++D G +a l dk 23726_FusoPortal_Gene_1356 189 AMKLAGKIAKNAPFAVRACKKAINEGIDT-DMDRAIIIEEKLFGSCFATEDQKVGMKAFL-DK 249 6778899**********************.5************************99976.55 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (259 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 404 (0.0241742); expected 334.2 (0.02) Passed bias filter: 364 (0.0217808); expected 334.2 (0.02) Passed Vit filter: 23 (0.00137626); expected 16.7 (0.001) Passed Fwd filter: 2 (0.000119674); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 2 [number of targets reported over threshold] # CPU time: 0.28u 0.26s 00:00:00.54 Elapsed: 00:00:00.34 # Mc/sec: 2214.47 // [ok]