# 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_1767.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1767 [L=353] Description: # 1888523 # 1889581 # 1 # ID=1_1767;partial=00;start_type=ATG;rbs_motif=GGAG/GAGG;rbs_spacer=5-10bp;gc_cont=0.216 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.1e-16 60.8 2.4 4.1e-13 49.2 1.6 2.3 2 SIS SIS domain ------ inclusion threshold ------ 0.014 15.3 0.0 0.031 14.2 0.0 1.6 1 SIS_2 SIS domain 0.055 13.6 0.2 0.11 12.6 0.2 1.5 1 PK_C Pyruvate kinase, alpha/beta domain Domain annotation for each model (and alignments): >> SIS SIS domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 49.2 1.6 7.4e-17 4.1e-13 6 127 .. 40 161 .. 37 167 .. 0.83 2 ! 10.0 0.0 9.5e-05 0.53 3 67 .. 211 275 .. 209 285 .. 0.91 Alignments for each domain: == domain 1 score: 49.2 bits; conditional E-value: 7.4e-17 SEEEEEESHHHHHHHHHHHHHHHHHHTSEEEEEEHHHHHHTGGGG-STTCEEEEEESSSHHHHHHHHHHHHHHTT.SEEEEEE CS SIS 6 kriyvlGaGtsyaiakelalkleeiggivvevedaselkhgvlalideddlvilisfsgetkdlvaaaelakerg.akiiaIT 87 k+i++l++G+s ++a+ ++ + i +i+ e+ + +++ ++ dlvi+is+sg +++++ a++++k+ + + iaIT 23726_FusoPortal_Gene_1767 40 KNILILATGSSMNAALITKYFISDILNINIEIKEPFN-YYNYEKINENIDLVIAISQSGKSASTISALKYVKKCKnIPSIAIT 121 89******************************99988.56645666777********************8888662666**** PP SSTTSHHHHHSSEEEEEEEECHCSSSHCHHHHHHHHHHHH CS SIS 88 dskgsplareadlalyinaepetlvastksltaqlaalda 127 ++++s++a+e+ +l++ + e tk ++a+ + l + 23726_FusoPortal_Gene_1767 122 SNNMSIIAKESNMILDLGIGIEQVGFVTKGFSATVLNLFL 161 ***********99999999999944445555555544443 PP == domain 2 score: 10.0 bits; conditional E-value: 9.5e-05 CTSSEEEEEESHHHHHHHHHHHHHHHHHHTSEEEEEEHHHHHHTGGGG-STTCEEEEEESSSHHH CS SIS 3 akakriyvlGaGtsyaiakelalkleeiggivvevedaselkhgvlalideddlvilisfsgetk 67 +akr + +G G+ + + ke+ k+ e+ ++ + + ++++hg+ ++++++++ ++ g + 23726_FusoPortal_Gene_1767 211 LNAKRFIGIGYGACFGLVKEFETKFTETIRLPSQGFELEAYMHGPYLEANKEHIIFYFDNKGKLS 275 689999*******************9999999999999999*****************9998766 PP >> SIS_2 SIS domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 14.2 0.0 5.5e-06 0.031 102 138 .] 85 122 .. 75 122 .. 0.84 Alignments for each domain: == domain 1 score: 14.2 bits; conditional E-value: 5.5e-06 -TT-EEEEE-SS--SHHHHHHHHHHHHH.T-EEEEEE- CS SIS_2 102 kpgDvLlviSnsGnnavvveaaeaaker.gmkviAlts 138 ++ D+ + iS+sG++a+ + a++++k++ ++ iA+ts 23726_FusoPortal_Gene_1767 85 ENIDLVIAISQSGKSASTISALKYVKKCkNIPSIAITS 122 4569*********************9962678888887 PP >> PK_C Pyruvate kinase, alpha/beta domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.6 0.2 2.1e-05 0.11 15 98 .. 86 180 .. 81 205 .. 0.70 Alignments for each domain: == domain 1 score: 12.6 bits; conditional E-value: 2.1e-05 T-SEEEEE-SSSHHHHHHH..HT-TS...SEEEEEESHHHH.....HHHHH.HCTTEEHHH.HHHHH...TTSHHHHHHHHHH CS PK_C 15 gakaIvvltesGstArlis..kyrpk...vpIiavtpnelt.....arqla.lyrGvypvlvekaktsestdellreavkvak 86 ++++++++++sG++A is ky +k +p ia+t+n+ ++ + l G+ v + + s ++ +l+ a+ ak 23726_FusoPortal_Gene_1767 86 NIDLVIAISQSGKSASTISalKYVKKcknIPSIAITSNNMSiiakeSNMILdLGIGIEQVGFVTKGFSATVLNLFLLAIILAK 168 7999************98833444447779999***98888786643333346667777766643334555555666688888 PP HTTSS-STSEEE CS PK_C 87 eaglikkgdevv 98 e++li+++++ v 23726_FusoPortal_Gene_1767 169 EKKLISANQKEV 180 888887766544 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (353 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 639 (0.038236); expected 334.2 (0.02) Passed bias filter: 244 (0.0146003); expected 334.2 (0.02) Passed Vit filter: 15 (0.000897559); 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.25u 0.22s 00:00:00.47 Elapsed: 00:00:00.30 # Mc/sec: 3420.61 // [ok]