# 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_262.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_262 [L=207] Description: # 279707 # 280327 # -1 # ID=1_262;partial=00;start_type=ATG;rbs_motif=AGGAG;rbs_spacer=5-10bp;gc_cont=0.279 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.6e-67 226.3 0.2 1.8e-67 226.1 0.2 1.0 1 TMP-TENI Thiamine monophosphate synthase ------ inclusion threshold ------ 0.016 14.5 0.1 0.019 14.3 0.1 1.1 1 NMO Nitronate monooxygenase 0.024 13.7 0.1 0.028 13.5 0.1 1.3 1 NanE Putative N-acetylmannosamine-6-phosphate epimer 0.043 13.2 0.4 1.1 8.6 0.0 2.2 2 His_biosynth Histidine biosynthesis protein 0.078 12.6 0.1 0.15 11.8 0.1 1.6 2 QRPTase_C Quinolinate phosphoribosyl transferase, C-termi Domain annotation for each model (and alignments): >> TMP-TENI Thiamine monophosphate synthase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 226.1 0.2 5.4e-71 1.8e-67 1 181 [] 8 188 .. 8 188 .. 0.99 Alignments for each domain: == domain 1 score: 226.1 bits; conditional E-value: 5.4e-71 TMP-TENI 1 lylvtdskldkneellevveealkggvtlvqlRekeldarellelaeklaalcrkygvpliindrvdlalavgadgvHlgqddl 84 +ylvtd+k+++++++ +++ee +kggv++vqlRek++++++++++a k++++c++y+v +iindr+d+ +av+adgvHlgq+d+ 23726_FusoPortal_Gene_262 8 IYLVTDEKACNGKDFYKCIEESIKGGVKIVQLREKNISTKDFYKKALKVKEICKNYEVLFIINDRLDITQAVEADGVHLGQSDM 91 6*********************************************************************************** PP TMP-TENI 85 pvaearellgpdliiGlSvhsleeaeeaealgadYigvGpifpTdtKkdaeplgleglkelaealeiPvvaiGGitlenvaevl 168 p+++are+l+++ +iG+++ ++eeae+a+ lgadYig G+if T+tK++a+ l++e lk+++++++iPv+aiGGi+ +nv ++ 23726_FusoPortal_Gene_262 92 PIEKAREILKDKFLIGATARNIEEAEKAQLLGADYIGSGAIFGTSTKDNAKRLEMEDLKKIVNSVKIPVFAIGGININNVWMLK 175 ************************************************************************************ PP TMP-TENI 169 eagadgvavvsai 181 + g++gv+ vs+i 23726_FusoPortal_Gene_262 176 NIGLQGVCSVSGI 188 ***********98 PP >> NMO Nitronate monooxygenase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 14.3 0.1 5.6e-06 0.019 81 203 .. 48 167 .. 32 200 .. 0.68 Alignments for each domain: == domain 1 score: 14.3 bits; conditional E-value: 5.6e-06 NMO 81 ekykkel..eeaalelnieegvrfyrkvevvldekvevvsfgfglpteevvarlkeagvvliptissvkearkvearGaDaliv 162 + ykk l +e ++++++ ++ ++++ ++++ v +g + + e + + ++ t+ ++ea k++ +GaD i 23726_FusoPortal_Gene_262 48 DFYKKALkvKEICKNYEVLFIIN--DRLDITQAVEADGVHLGQSDMPIEKAREILKDKFLIGATARNIEEAEKAQLLGADY-IG 128 44444431034455555555555..7777777777777777666665555555555566666788899************7.56 PP NMO 163 qGpeaGGhqGtpevgdkglfaLvpqvldavdiPviAaGGii 203 G G + + +++ +++ ++++++ v+iPv A GGi+ 23726_FusoPortal_Gene_262 129 SGAIFGTSTK--DNAKRLEMEDLKKIVNSVKIPVFAIGGIN 167 6666666655..3333678888999**************97 PP >> NanE Putative N-acetylmannosamine-6-phosphate epimerase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.5 0.1 8.4e-06 0.028 82 155 .. 94 168 .. 28 201 .. 0.78 Alignments for each domain: == domain 1 score: 13.5 bits; conditional E-value: 8.4e-06 HHHHHHHHHH--EEEEE-SSHHHHHHHHHTT-SEEE-.TTTTSSTTS-.SS--HHHHHHHHHT.T--EEEES---S CS NanE 82 eelikkikekyqllmaDistleeglvakklGvdivGt.tlsGytkeeekaepDfalvkalsea.glkviaeGrint 155 e+ ++i + l+ a ++ee+ a+ lG+d++G+ ++ G ++ + + +++ +k+++++ ++ v a G in 23726_FusoPortal_Gene_262 94 EKA-REILKDKFLIGATARNIEEAEKAQLLGADYIGSgAIFGTSTKDNAKRLEMEDLKKIVNSvKIPVFAIGGINI 168 443.445555567889999*****************9458898888888888887777766543899999998875 PP >> His_biosynth Histidine biosynthesis protein # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 2.6 0.1 0.022 72 64 104 .. 92 136 .. 72 146 .. 0.60 2 ? 8.6 0.0 0.00032 1.1 180 218 .. 146 183 .. 140 194 .. 0.83 Alignments for each domain: == domain 1 score: 2.6 bits; conditional E-value: 0.022 HHHHHHHH-SSEEEEESS--SHHHHHHHHHCT-S....EEEESHH CS His_biosynth 64 lveeiaeevfipvqvGGGIrsledveelleagad....kviigta 104 +e+ +e ++ ++ +G r++e++e++ +gad i gt+ 23726_FusoPortal_Gene_262 92 PIEKAREILKDKFLIGATARNIEEAEKAQLLGADyigsGAIFGTS 136 366666666677777777777777777777777733333344444 PP == domain 2 score: 8.6 bits; conditional E-value: 0.00032 HHHHHHHHHHHSSEEEEESS--SHHHHHHHHHTTCSEEE CS His_biosynth 180 leltkelaeavkipviasGGvgsledlkelksegvdavi 218 +e +k++ ++vkipv a GG+ +++ + lk+ g ++v 23726_FusoPortal_Gene_262 146 MEDLKKIVNSVKIPVFAIGGI-NINNVWMLKNIGLQGVC 183 5668999*************9.89999999988888875 PP >> QRPTase_C Quinolinate phosphoribosyl transferase, C-terminal domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -3.0 0.0 1.5 4.9e+03 94 112 .. 27 47 .. 21 67 .. 0.52 2 ? 11.8 0.1 4.3e-05 0.15 64 150 .. 87 180 .. 69 194 .. 0.80 Alignments for each domain: == domain 1 score: -3.0 bits; conditional E-value: 1.5 HHHHHHT-SEEEEE..S--HH CS QRPTase_C 94 eevleaeadiimLD..nfsve 112 ee ++ ++ i++L n+s++ 23726_FusoPortal_Gene_262 27 EESIKGGVKIVQLRekNISTK 47 555555666666631133333 PP == domain 2 score: 11.8 bits; conditional E-value: 4.3e-05 SSHHHHHHHHHHHHTTTSEEEEEESSHHHHHHHHHHT-SEEEEES--HHHHHHHHHHHHHHSTT.......SEEEEESS--TTT CS QRPTase_C 64 gsiveavksvreavdkevkievevesleeleevleaeadiimLDnfsveelkkaveelkeknkr.......vlievSGgvtlen 140 g+ ++++re ++++ i +++ ee+e++ +ad i + + +k+ +++l++ + + + +++ Gg++ +n 23726_FusoPortal_Gene_262 87 GQSDMPIEKAREILKDKFLIGATARNIEEAEKAQLLGADYIGSGAIFGTSTKDNAKRLEMEDLKkivnsvkIPVFAIGGININN 170 555566789999998887799999*****************999997777777777776666656667766999********** PP HHHHHCTT-S CS QRPTase_C 141 veeyaksgvd 150 v ++++g++ 23726_FusoPortal_Gene_262 171 VWMLKNIGLQ 180 **99999876 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (207 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 816 (0.0488272); expected 334.2 (0.02) Passed bias filter: 536 (0.0320728); expected 334.2 (0.02) Passed Vit filter: 44 (0.00263284); expected 16.7 (0.001) Passed Fwd filter: 5 (0.000299186); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 5 [number of targets reported over threshold] # CPU time: 0.31u 0.27s 00:00:00.58 Elapsed: 00:00:00.36 # Mc/sec: 1671.54 // [ok]