# 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_881.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_881 [L=268] Description: # 982055 # 982858 # 1 # ID=1_881;partial=00;start_type=ATG;rbs_motif=AGGA;rbs_spacer=5-10bp;gc_cont=0.257 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.1e-52 177.5 0.4 8.3e-51 172.8 0.2 1.9 2 NAD_kinase ATP-NAD kinase 0.0094 15.5 0.0 0.017 14.6 0.0 1.5 1 DAGK_cat Diacylglycerol kinase catalytic domain Domain annotation for each model (and alignments): >> NAD_kinase ATP-NAD kinase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 4.6 0.0 0.0018 15 1 41 [. 3 41 .. 3 44 .. 0.88 2 ! 172.8 0.2 1e-54 8.3e-51 75 290 .. 39 245 .. 31 246 .. 0.95 Alignments for each domain: == domain 1 score: 4.6 bits; conditional E-value: 0.0018 .EEEEESSCGGTCSHHHHHHHHHHHCT-..EEEEECCGGCC CS NAD_kinase 1 kvgivvnpikkeaieraeevarwLldalsrvlevtveeklg 41 k+ i+ n++k+ ai++ +e++++L+ + + +ev++++++ 23726_FusoPortal_Gene_881 3 KLSIIYNEDKEDAIKIYKELLKYLKSK--KEFEVLDDKNIS 41 5789**********************9..689999888765 PP == domain 2 score: 172.8 bits; conditional E-value: 1e-54 TTTSSSEEEEEE-HHHHHHHHHHHHHCT-EEEEEESSS--SS-SEEGGGHHHHHHHHHHT-SEEEEEEEEEEEEEETTE..... CS NAD_kinase 75 adagvdliivvGGDGTaldvarllqkaeiPvlgvnaGvkgfltefepeaaaellekllegeirvekrevldavimrskgeedia 158 ++++ i+v+GGDGT l+ ++ ++++++ ++++n G++g+lte++ + ++++ e++l+g+i++e+r ++++i ++k 23726_FusoPortal_Gene_881 39 NISQAEYIVVIGGDGTLLRGFKKIKDKKVKIIAINSGTLGYLTEIRKDGYKKIFENILKGKINIEERYFFTVKIGKKKY----- 117 5678999****************************************************************99987777..... PP .EEEEEEESSEEEEEESSSSS-EE.EEEEETTEEEEEEEESEEEEE-TGGGGTHHHHTT--EE-TT---EEEEEESESSS---- CS NAD_kinase 159 diviensvlnevviirgpgstlie.lelyvdtelleadqldglivstPtGstayslsaGgplisPgvlailivpicPhslssrp 241 ++lnev + + +++ i e+yvd+++l+++ dg+i++tPtGstayslsaGgp+++P+++ +li+pi Ph l rp 23726_FusoPortal_Gene_881 118 ------NALNEVFLTKDNIKRNIVsSEIYVDDKFLGKFKGDGVIIATPTGSTAYSLSAGGPIVTPELKLFLITPIAPHNLNTRP 195 ......68**********9998888*********************************************************** PP EEE-TT-EEEEEE-.SS--EEEEETTTEEEEE-TT-EEEEEEESSEEEEE CS NAD_kinase 242 iVvssssklrirva.srseallaldGseelelspgdevrvtasdkkfkli 290 i++s + k+ +++a + ++ +dG+ +++++dev+++ s+++ k++ 23726_FusoPortal_Gene_881 196 IILSGDVKIVLTLAgPSEFGIVNVDGHTHNKINLEDEVEISYSKESLKIV 245 ************99988999************************999887 PP >> DAGK_cat Diacylglycerol kinase catalytic domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 14.6 0.0 2.1e-06 0.017 50 92 .. 37 76 .. 2 99 .. 0.69 Alignments for each domain: == domain 1 score: 14.6 bits; conditional E-value: 2.1e-06 HHHHCC-SEEEEEESHHHHHHHHHHHCTHHCTTSEEEEEE-SS CS DAGK_cat 50 eakeegvdlivvaGGDGtvsevlnglleelatkpplgiiPlGt 92 + + ++++ ivv+GGDGt+ ++ + + +k+++ i Gt 23726_FusoPortal_Gene_881 37 DKNISQAEYIVVIGGDGTLLRGFKKIKD---KKVKIIAINSGT 76 2223488***********9888888888...666666666666 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (268 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 983 (0.05882); expected 334.2 (0.02) Passed bias filter: 588 (0.0351843); expected 334.2 (0.02) Passed Vit filter: 34 (0.00203447); 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.29u 0.25s 00:00:00.54 Elapsed: 00:00:00.34 # Mc/sec: 2291.43 // [ok]