# 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_1657.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1657 [L=397] Description: # 1779634 # 1780824 # -1 # ID=1_1657;partial=00;start_type=ATG;rbs_motif=AGGAG;rbs_spacer=5-10bp;gc_cont=0.234 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 4.6e-119 398.1 9.4 5.8e-119 397.8 9.4 1.0 1 HIGH_NTase1 HIGH Nucleotidyl Transferase 9.5e-06 25.8 0.0 4.4e-05 23.6 0.0 2.1 2 CTP_transf_like Cytidylyltransferase-like 0.0084 15.7 0.2 0.37 10.3 0.1 2.5 2 Citrate_ly_lig Citrate lyase ligase C-terminal domain Domain annotation for each model (and alignments): >> HIGH_NTase1 HIGH Nucleotidyl Transferase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 397.8 9.4 1e-122 5.8e-119 2 389 .. 4 383 .. 3 388 .. 0.94 Alignments for each domain: == domain 1 score: 397.8 bits; conditional E-value: 1e-122 HIGH_NTase1 2 kvvGiIaEYNPfHnGHlyhlekakeltgadvviavmSGnFvqRGepAildKwtRaemALenGvDlViElPvvyavqsAelFAk 84 +vvG+I+EYNPfHnGHl+h++++ +l+++++ iavmSG++vqRGep++++K +++++AL++G+D+ViElP++y++qsAe+FAk 23726_FusoPortal_Gene_1657 4 NVVGLIVEYNPFHNGHLHHIQEIDRLFEDNIKIAVMSGDYVQRGEPSLINKLEKTKIALSQGIDIVIELPTFYSTQSAEIFAK 86 68********************************************************************************* PP HIGH_NTase1 85 gaVslLnalgvvdslvFGsEsgdiealkkiaeilaeepeeykealkeelkkglsypkarekalkeykkeeelllsspNniLgi 167 g+V+lLn+l ++++vFGsEs+d+++lk+ia+i ++e++ +l+e+l +g+syp+a++k+l + +N+iL++ 23726_FusoPortal_Gene_1657 87 GSVNLLNKLS-CSHIVFGSESNDLDKLKRIATIS--LTKEFELSLREFLAEGFSYPTAFSKVLFD-------EKLGSNDILAL 159 **********.********************765..5689999*********************7.......45689****** PP HIGH_NTase1 168 eYlkalkklnskikpvtikRegagyhdeelseeeeiaSAtaiRkal.keek....kleelkklvPestlellkeelkkkalls 245 eYl+a+k++nski+ ++ikRe++gy d+e ++++SAt+iRk+l ++++ kl+++k+lvPe ++++l+e+ +++ + 23726_FusoPortal_Gene_1657 160 EYLRAIKTINSKIEAYCIKREKTGYYDDEK---DNFSSATYIRKILlDSNEkkenKLNKIKNLVPEFSYKILEEN--FGVFSC 237 ****************************54...58***********6444467766777**************99..567999 PP HIGH_NTase1 246 eedlskllkykllteskeelkkildvsegLenrikkalkkaesleelleklktKrytytriqRillhillnikkedlkklkkk 328 ++d+++l+ky+++++ + elk+i+d++ gLenr++k ++ ++e++++++ tKr+t +r+qRillh l+n++++ ++k+k++ 23726_FusoPortal_Gene_1657 238 LSDFYDLIKYNIIKNCS-ELKNIQDLEIGLENRLYKYSLENLKFEDFFNEVLTKRITISRLQRILLHSLFNLTENITEKVKNE 319 ****************7.5************************************************************9987 PP HIGH_NTase1 329 kpylrvLGfnkkgrklLkeikkks...slplitklskakkklkeeakklleldiratdlyslls 389 p +++LGf++kg+k+L+ +kk++ + +++t+ ++ k+ l+ee+ kl++ + ++++y++ s 23726_FusoPortal_Gene_1657 320 VPFVKILGFSTKGQKYLNYLKKSKnynERKILTSNRNLKEILNEEEAKLFNFNELCSQIYRIKS 383 68**************99999986444789******************************7754 PP >> CTP_transf_like Cytidylyltransferase-like # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 23.6 0.0 7.9e-09 4.4e-05 6 142 .. 12 200 .. 11 201 .. 0.76 2 ? -2.7 0.0 1 5.7e+03 100 133 .. 313 346 .. 290 363 .. 0.68 Alignments for each domain: == domain 1 score: 23.6 bits; conditional E-value: 7.9e-09 -TT--HHHHHHHHHHHHHSSESEEEEEESSTHCCST-SCSSTHHHHHHHHHHHTTCTTEEEEEEESEHH..........HHHH CS CTP_transf_like 6 FDpiHlGHlrlleqakklleedlivgvpsdeppkkkkkplasaeeRlemlelakevdevisvdpweltr..........evlk 78 ++p+H+GHl+ ++++ +l+e + ++v+s + +++ +++l++ e+++ + l + d vi+ +++t+ ++l+ 23726_FusoPortal_Gene_1657 12 YNPFHNGHLHHIQEIDRLFEDNIKIAVMSGDYVQRGEPSLINKLEKTK-IALSQGIDIVIELPTFYSTQsaeifakgsvNLLN 93 9*******************998899********87566665555544.3444899*************************** PP HTTTTEEEEEEEHHHHHHHHCHHHHHHHHHHHHHEEEEEE........................................... CS CTP_transf_like 79 klkpdklylvigadvledfsewkeleeilgnvklvvvsrp........................................... 118 kl + +v+g ++ d+++ k++ +i +++ r+ 23726_FusoPortal_Gene_1657 94 KLSCS--HIVFGSES-NDLDKLKRIATISLTKEFELSLREflaegfsyptafskvlfdeklgsndilaleylraiktinskie 173 *****..99999998.5666666666666666666666666666666668888888888888888888888888888877777 PP ...HTTCEEEEEEECCTTTCSHHHHHH CS CTP_transf_like 119 ...ekqanvllkaekrtngiSstkire 142 k+++ + ++++ n+ S+t+ir+ 23726_FusoPortal_Gene_1657 174 aycIKREKTGYYDDEKDNFSSATYIRK 200 776333344455666777888888887 PP == domain 2 score: -2.7 bits; conditional E-value: 1 HHHHHHHHHHHHHEEEEEEHTTCEEEEEEECCTT CS CTP_transf_like 100 wkeleeilgnvklvvvsrpekqanvllkaekrtn 133 +++++ + vk+ ++s + +++ +lk++k++n 23726_FusoPortal_Gene_1657 313 TEKVKNEVPFVKILGFSTKGQKYLNYLKKSKNYN 346 4466666666666666666555555555555554 PP >> Citrate_ly_lig Citrate lyase ligase C-terminal domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 2.9 0.0 0.013 72 4 22 .. 8 26 .. 5 75 .. 0.81 2 ! 10.3 0.1 6.6e-05 0.37 136 181 .. 170 227 .. 154 228 .. 0.74 Alignments for each domain: == domain 1 score: 2.9 bits; conditional E-value: 0.013 Citrate_ly_lig 4 ivmnanPftlGhryLveka 22 +++ nPf +Gh + +++ 23726_FusoPortal_Gene_1657 8 LIVEYNPFHNGHLHHIQEI 26 78999*******9988765 PP == domain 2 score: 10.3 bits; conditional E-value: 6.6e-05 Citrate_ly_lig 136 aeievveierk......ekeeelisaskvrkll......keknleklaelvPkttley 181 ++ie i+r ++++++ sa+ +rk+l ke++l+k+++lvP+ +++ 23726_FusoPortal_Gene_1657 170 SKIEAYCIKREktgyydDEKDNFSSATYIRKILldsnekKENKLNKIKNLVPEFSYKI 227 347777777765555556789999*******985555556778999*******99986 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (397 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 892 (0.0533748); expected 334.2 (0.02) Passed bias filter: 372 (0.0222595); expected 334.2 (0.02) Passed Vit filter: 29 (0.00173528); 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.28u 0.24s 00:00:00.52 Elapsed: 00:00:00.33 # Mc/sec: 3497.25 // [ok]