# 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_699.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_699 [L=218] Description: # 792400 # 793053 # 1 # ID=1_699;partial=00;start_type=ATG;rbs_motif=AGGAGG;rbs_spacer=5-10bp;gc_cont=0.320 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.8e-70 236.6 3.5 2e-70 236.5 3.5 1.0 1 CoA_trans Coenzyme A transferase 2.1e-06 27.7 1.2 4.5e-05 23.4 0.0 2.4 2 AcetylCoA_hydro Acetyl-CoA hydrolase/transferase N-terminal 0.006 16.2 0.2 0.0097 15.5 0.2 1.5 1 FAA_hydrolase Fumarylacetoacetate (FAA) hydrolase family Domain annotation for each model (and alignments): >> CoA_trans Coenzyme A transferase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 236.5 3.5 3.6e-74 2e-70 3 216 .. 7 215 .. 5 216 .. 0.99 Alignments for each domain: == domain 1 score: 236.5 bits; conditional E-value: 3.6e-74 HHHHHHCTT--TTEEEEE-SBTTBTHHHHHHHHHHHHTTTSEEEEESTTTE..ECCHHHHHTTTEEEEEESCTTBBCEHHHHHH CS CoA_trans 3 taaeavakeikdGmtvavGGFgligiPealiaalarsgvkdltviiseagv..lglgplllngsvkkviasyggetvnprfgkq 84 +++ea++ ++kdGmtv++GGF ++g+Pe++i+al+++gvkdlt++ +++g+ g+g+l++n++vkkvias++g+ np++g++ 23726_FusoPortal_Gene_699 7 SMEEAIS-HVKDGMTVHIGGFIACGTPESIITALIEKGVKDLTIVANDTGLidKGIGRLVVNNQVKKVIASHIGT--NPETGRR 87 89*****.*******************************************9***********************..******* PP HHTTSSEEEEHHHHHHHHHHHHHHHTEETTEEEESTTTTSHGGTTESEEEEETTEEEEEEESSSEEEEEEECCEEETTCSEESS CS CoA_trans 85 flsgelefellpqgtlaerlraGgaglpelgfltntgigtlvapklkkvpgfgGaeyilepalkadvalikaskadeegnlvfe 168 + sge+e+el+pqgtlaer ra+g gl+ g+lt+tg+gt+v+++ k++ +++G++y+le+++kadvali +k de gn+ e 23726_FusoPortal_Gene_699 88 MQSGEMEVELVPQGTLAERVRAAGYGLG--GILTPTGLGTIVQEG-KQIINVDGKDYLLEKPIKADVALIFGTKVDELGNVICE 168 ****************************..***************.9************************************* PP GGGHTTHSHHHHHHEEEEEEEEEEEECTTSSECCCESBTGGGHSEEEE CS CoA_trans 169 ksagnfnkpaiaeaakvtileveeivekgkldpekvliPgvlvdavve 216 k+++nfn p +a+aa v+i+e +eiv+ g+l+pe+ +i ++++d++v+ 23726_FusoPortal_Gene_699 169 KTTKNFN-PLMATAADVVIVEALEIVPAGSLSPEHLDISRIFIDYIVK 215 *******.***************************************8 PP >> AcetylCoA_hydro Acetyl-CoA hydrolase/transferase N-terminal domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 23.4 0.0 8.2e-09 4.5e-05 7 46 .. 3 42 .. 1 54 [. 0.90 2 ? 3.0 0.2 0.015 84 125 187 .. 149 214 .. 125 217 .. 0.68 Alignments for each domain: == domain 1 score: 23.4 bits; conditional E-value: 8.2e-09 AcetylCoA_hydro 7 kklkepeeaaelvkngshielsgftavgtakaipkylakr 46 +k+ + eea++ vk+g+ ++++gf+a gt+ i ++l+ + 23726_FusoPortal_Gene_699 3 QKIVSMEEAISHVKDGMTVHIGGFIACGTPESIITALIEK 42 788899*******************************965 PP == domain 2 score: 3.0 bits; conditional E-value: 0.015 AcetylCoA_hydro 125 iDvalieataldedGylnlG.vgvrakKviiEvaelvvivvkknpPrvngydelialdkv..dyiv 187 Dvali t +de G + + +++ ++ a +v+++ + +P +e ++++ dyiv 23726_FusoPortal_Gene_699 149 ADVALIFGTKVDELGNVICEkTTKNFNPLMATAADVVIVEALEIVPAGSLSPEHLDISRIfiDYIV 214 689999999999999666552444555555555666666666777766666677777776677776 PP >> FAA_hydrolase Fumarylacetoacetate (FAA) hydrolase family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 15.5 0.2 1.7e-06 0.0097 155 217 .. 3 107 .. 1 128 [. 0.68 Alignments for each domain: == domain 1 score: 15.5 bits; conditional E-value: 1.7e-06 GTTTSHHHHHHHHHHHH-B-TT-EEEEEES..................................................S--E CS FAA_hydrolase 155 dmifspaeliahlsagltLepGdliltGTpsGv.........ga..................................eeeslk 195 + i+s++e i+h+ g+t++ G i+ GTp++ ++ ++e+ + 23726_FusoPortal_Gene_699 3 QKIVSMEEAISHVKDGMTVHIGGFIACGTPESIitaliekgvKDltivandtglidkgigrlvvnnqvkkviashigtNPETGR 86 6699**************************887655555444004444444444444444555555555555555555666666 PP E--TT-EEEEEETTTEEEEEEE CS FAA_hydrolase 196 flkdgDeveveieglGelrntV 217 ++ g+ +eve+ G+l+ +V 23726_FusoPortal_Gene_699 87 RMQSGE-MEVELVPQGTLAERV 107 666664.555664455555555 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (218 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 327 (0.0195668); expected 334.2 (0.02) Passed bias filter: 272 (0.0162757); expected 334.2 (0.02) Passed Vit filter: 22 (0.00131642); 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.23u 0.23s 00:00:00.46 Elapsed: 00:00:00.31 # Mc/sec: 2044.30 // [ok]