# 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_1491.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1491 [L=210] Description: # 1588330 # 1588959 # 1 # ID=1_1491;partial=00;start_type=ATG;rbs_motif=AGGA;rbs_spacer=5-10bp;gc_cont=0.214 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 4e-25 88.1 0.1 6.2e-25 87.5 0.1 1.3 1 TPK_catalytic Thiamin pyrophosphokinase, catalytic domain 9.8e-09 34.7 0.0 1.6e-08 34.0 0.0 1.4 1 TPK_B1_binding Thiamin pyrophosphokinase, vitamin B1 binding Domain annotation for each model (and alignments): >> TPK_catalytic Thiamin pyrophosphokinase, catalytic domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 87.5 0.1 7.4e-29 6.2e-25 11 111 .. 30 121 .. 21 127 .. 0.95 Alignments for each domain: == domain 1 score: 87.5 bits; conditional E-value: 7.4e-29 TPK_catalytic 11 rvcaDgGanrlyeeekerlellpdvivGDlDSirpevkefyeskgteivkkp.dqdstDlekaialvlekeeedknksivvlg 92 + caDgGan +y +l l+p+ i GDlDSi++evkefy++k+++++k + ++d+tD e+ ++ +++k+ + i+++ 23726_FusoPortal_Gene_1491 30 IYCADGGANFCY-----ELTLIPKEIYGDLDSIKDEVKEFYQEKKVKFIKFKiEKDYTDSELLLNEIQNKY-----DVIYCIA 102 57**********.....****************************************************99.....99***** PP TPK_catalytic 93 alggRfDhtlaninvllkl 111 +lgg +Dh+l+nin+l k+ 23726_FusoPortal_Gene_1491 103 GLGGSIDHELTNINLLAKY 121 ****************986 PP >> TPK_B1_binding Thiamin pyrophosphokinase, vitamin B1 binding domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 34.0 0.0 1.9e-12 1.6e-08 13 67 .. 147 200 .. 125 202 .. 0.87 Alignments for each domain: == domain 1 score: 34.0 bits; conditional E-value: 1.9e-12 .-EEEEEESSSSEEEEEEESBSS-EEEEEESTTCETEEEEEBS.SSEEEEEESBS CS TPK_B1_binding 13 gkyvgliPlggevtglttkGlkwnLedwelefggllstSNeveeeeevtietdsg 67 ++ ++++ ++++v+glt+kG+k+ +e+ +++ g+ +++SN ++ e+++ ++sg 23726_FusoPortal_Gene_1491 147 NTKISFVIFSDQVKGLTLKGFKYSIENLDIKKGEARCISNIIV-ENKANLLIKSG 200 5779***************************************.77777776666 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (210 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 887 (0.0530756); expected 334.2 (0.02) Passed bias filter: 404 (0.0241742); expected 334.2 (0.02) Passed Vit filter: 21 (0.00125658); 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.27u 0.23s 00:00:00.50 Elapsed: 00:00:00.33 # Mc/sec: 1849.93 // [ok]