# 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_715.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_715 [L=203] Description: # 808026 # 808634 # -1 # ID=1_715;partial=00;start_type=ATG;rbs_motif=AGGA;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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 5.4e-54 182.6 0.1 6.3e-54 182.4 0.1 1.0 1 Dak2 DAK2 domain Domain annotation for each model (and alignments): >> Dak2 DAK2 domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 182.4 0.1 3.8e-58 6.3e-54 1 173 [. 26 199 .. 26 200 .. 0.98 Alignments for each domain: == domain 1 score: 182.4 bits; conditional E-value: 3.8e-58 TTS-S-HHHHHHHHHHHHHHHHHC---SSHHHHHHHHHHHHHHH--STHHHHHHHHHHHHHHHC.....TT---HHHHHHHHHH CS Dak2 1 pvGDGDtGtnmargaeavlkaleeleeaslaevlkalakallsgagGtsGvlysqflralakalkekeeldaedlaaalkaale 84 ++GDGD+G n+arg+++++++l++ ++ +++++++++++ll+++gG+sG++y+++++++++ lk+k+e++ + l + l++++e 23726_FusoPortal_Gene_715 26 EIGDGDHGVNLARGFTEIKNQLNNFKDLPVSDIFTKMGMILLTKVGGASGAIYGTAFMSAGTYLKGKTEFNNQILLETLNSMIE 109 59********************************************************************************** PP HHHHHH---TTSSSTHHHHHHHHHHHHT............T---HHHHHHHHHHCCS---SSGGGGGGGGGTTS--HHHHHHHH CS Dak2 85 aiqavgkarpgdkTmlDvlrpavealekakd.gedlkealekaakaaekgaeaTkellaklgrAsylgerslgvvDpGAvglal 167 +iq++gka +kTmlD++ p+ + lek+ + g++lk++ +++++ a++ +eaTk++ a++grAsylgers+g++DpGA++++l 23726_FusoPortal_Gene_715 110 GIQKRGKAVLNEKTMLDTIMPTYNFLEKSFNeGKNLKDIKNEVIEVAKNSMEATKDIIATKGRASYLGERSVGHIDPGAMSSYL 193 ****************************8888**************************************************** PP HHHHHH CS Dak2 168 ileala 173 ++++++ 23726_FusoPortal_Gene_715 194 MIKTVC 199 **9986 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (203 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 820 (0.0490665); expected 334.2 (0.02) Passed bias filter: 528 (0.0315941); expected 334.2 (0.02) Passed Vit filter: 33 (0.00197463); expected 16.7 (0.001) Passed Fwd filter: 1 (5.98372e-05); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 1 [number of targets reported over threshold] # CPU time: 0.23u 0.23s 00:00:00.46 Elapsed: 00:00:00.30 # Mc/sec: 1967.09 // [ok]