# 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_1181.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1181 [L=253] Description: # 1281331 # 1282089 # 1 # ID=1_1181;partial=00;start_type=ATG;rbs_motif=AGGAGG;rbs_spacer=5-10bp;gc_cont=0.274 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.2e-44 150.8 0.0 1.5e-43 148.7 0.0 1.8 1 SIR2 Sir2 family ------ inclusion threshold ------ 0.12 12.0 0.0 0.74 9.5 0.0 2.0 2 TPP_enzyme_M Thiamine pyrophosphate enzyme, central domain Domain annotation for each model (and alignments): >> SIR2 Sir2 family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 148.7 0.0 1.8e-47 1.5e-43 1 177 [] 27 194 .. 27 194 .. 0.93 Alignments for each domain: == domain 1 score: 148.7 bits; conditional E-value: 1.8e-47 -GGGGGGGT---SSSTTSHHHHSCHHHHSSHHHHHHHHHHHHH..HHHHHHHHHHCC.TSS--HHHHHHHHHHHTTCEEEEEE CS SIR2 1 GAGiStesGipdfRskeGlyaklakeelaspealfskkllvee..vfyniarellae.eaepnaaHkalkklekkgkllrliT 80 GAG St+sG++dfR+k+Gly+ l k++ + pe++ s + + ++ +f++ +++ l+ +pn+ H al +lek g l+++iT 23726_FusoPortal_Gene_1181 27 GAGTSTDSGVKDFRGKDGLYKTLYKDK-YRPEEVLSSDFFYSHrdIFMEYVEKELNIkGLKPNKGHMALVELEKIGILKAVIT 108 9***********************999.88999999998877667****99998887699*********************** PP S-SSSHHHHTT-SSEEETTEEEEEEEETTTSBEEECCGGHCHHHCCTT.B-TTTCCBEEEEE--CTSB-S.HHHHHHHHHHHH CS SIR2 81 QNiDgLekkaGskkvvelHGslakakcvkcekkldkeelqeeieaekvpecekcgkllkpdvvfFGealp.dklkkliekvee 162 QNiD+L++ +G+k+v+elHGsl++ c +c+k+ d+ + ++cg++++pdv+++Ge+l+ + ++++i ++e+ 23726_FusoPortal_Gene_1181 109 QNIDDLHQVSGNKNVLELHGSLKRWYCLSCGKTADR------------NFSCECGGVVRPDVTLYGENLNqSVVNEAIYQLEQ 179 *******************************96666............23578*****************9999********* PP -SEEEEES--STCCT CS SIR2 163 adlllviGtSlkVap 177 ad+l+v GtSl+V+p 23726_FusoPortal_Gene_1181 180 ADTLIVAGTSLTVYP 194 *************98 PP >> TPP_enzyme_M Thiamine pyrophosphate enzyme, central domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -0.1 0.0 0.08 6.7e+02 3 29 .. 11 37 .. 9 40 .. 0.81 2 ? 9.5 0.0 8.8e-05 0.74 50 110 .. 152 211 .. 147 224 .. 0.72 Alignments for each domain: == domain 1 score: -0.1 bits; conditional E-value: 0.08 HHHHHHHH-SSEEEEE-CHHCHTTGHH CS TPP_enzyme_M 3 raaellkkakrpvilvGggvsrsgask 29 +++++lk+ k v + G+g+s++ k 23726_FusoPortal_Gene_1181 11 ELVKILKNTKYLVFFGGAGTSTDSGVK 37 689999************999876555 PP == domain 2 score: 9.5 bits; conditional E-value: 8.8e-05 TSS-TTSTTEEEESGGGSCHHHHHHHHH-SEEEEES-SSSTCCCTTTTTTTTTT-EEEEEE CS TPP_enzyme_M 50 gaipedhplylGmlgllgspaanealeeaDlvlllGarlsdirttgsspefapeakiihid 110 g++ d ly+ l ++a le+aD++++ G+ l+ + + ++ ++ + +++i i+ 23726_FusoPortal_Gene_1181 152 GVVRPDVTLYGENLNQSVVNEAIYQLEQADTLIVAGTSLTVYPAA-YYLRYFRGKNLIIIN 211 66677777887788888888889999************9886554.444343344565554 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (253 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 628 (0.0375778); expected 334.2 (0.02) Passed bias filter: 492 (0.0294399); 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.28u 0.27s 00:00:00.55 Elapsed: 00:00:00.36 # Mc/sec: 2043.00 // [ok]