# 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_174.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_174 [L=509] Description: # 184905 # 186431 # -1 # ID=1_174;partial=00;start_type=ATG;rbs_motif=AGGAG;rbs_spacer=5-10bp;gc_cont=0.248 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.3e-93 312.8 0.0 1.9e-93 312.3 0.0 1.2 1 NeuB NeuB family 2.4e-11 44.0 0.0 9.5e-11 42.1 0.0 2.0 2 SAF SAF domain 2.7e-07 31.3 0.0 8.7e-07 29.6 0.0 1.8 2 Acetyltransf_3 Acetyltransferase (GNAT) domain 0.00014 22.0 0.1 0.00042 20.4 0.1 1.9 1 Acetyltransf_1 Acetyltransferase (GNAT) family ------ inclusion threshold ------ 0.043 14.0 0.1 0.67 10.1 0.0 2.3 2 DUF4352 Domain of unknown function (DUF4352) Domain annotation for each model (and alignments): >> NeuB NeuB family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 312.3 0.0 5.5e-97 1.9e-93 2 239 .. 192 430 .. 191 432 .. 0.99 Alignments for each domain: == domain 1 score: 312.3 bits; conditional E-value: 5.5e-97 HHHHHHHT-SEEEEEE--HHHH--GGGGG--........TT-S-HHHHHHHH---HHHHHHHHHHHHHTT-EEEEEE-SHHHHH CS NeuB 2 ieaAaeagadavKfQtfkaeelvskeakkaeyqkkek..ekeesqlellkklelseeehkelkeyceekgieflstpfdlesvd 83 i++A+e+gadavK+Qt+++++l+ + ++++ +q+k+ ++++ +++l+k++ +++e+hke+++y++e +i ++stpfd+++v+ 23726_FusoPortal_Gene_174 192 IKVAKECGADAVKIQTYTPDTLTLN-CNNEYFQIKDGtiWDGKILYDLYKEAYTPWEWHKEIFDYAKELNICLFSTPFDKTAVN 274 99***********************.99999**999899********************************************* PP HHHHHT-S-EEE-GGGTT-HHHHHHHHTT-S-EEEE-TT--HHHHHHHHHHH.HHHT--EEEEE----SS--GGG--THHHHHH CS NeuB 84 lleslgveayKiaSgeltnlpllekiaktgkpvilStGmatleeieeavevlreagaeklvllhCtseYPapledvnLkaietl 167 lleslg ++yKiaS+e++++pl+e++a+++kpvi+StG+at eei++a+e++++ g+ +++ll+Cts+YPa+led+nL +ie+l 23726_FusoPortal_Gene_174 275 LLESLGNPIYKIASFEINDIPLIEYAASKKKPVIISTGVATEEEIKDAIEACKKVGNYDITLLQCTSQYPAKLEDANLVMIEDL 358 ************************************************************************************ PP HHH-TSEEEEE--SSSSHHHHHHHHHT--EEEEEB-S-TT--STTGGG-B-HHHHHHHHHHHHHHHHHSS-- CS NeuB 168 keafgvpvGySDHtlgiavaiaAvalGAkviEkHftldrnlegpDhkaSlepeelkelveairnvekalgeg 239 +++f+v+ G+SDHt+g va++Ava+GAkv+EkHf+l+r++ gpD+++S+ p+e+ke+v++irnvek++g+ 23726_FusoPortal_Gene_174 359 AKRFSVKSGLSDHTMGHLVATTAVAMGAKVVEKHFILNRSIGGPDSSFSMLPSEFKEMVDNIRNVEKMIGKI 430 *********************************************************************986 PP >> SAF SAF domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -2.1 0.0 1.8 6.1e+03 4 27 .. 376 399 .. 376 403 .. 0.86 2 ! 42.1 0.0 2.8e-14 9.5e-11 2 59 .. 446 499 .. 445 503 .. 0.93 Alignments for each domain: == domain 1 score: -2.1 bits; conditional E-value: 1.8 EEESSTB-TT-BEBTTTEEEECEC CS SAF 4 vvAardikaGeviteddltlkrpi 27 +vA+++++ G ++ e++ l r+i 23726_FusoPortal_Gene_174 376 LVATTAVAMGAKVVEKHFILNRSI 399 699999999999999998888775 PP == domain 2 score: 42.1 bits; conditional E-value: 2.8e-14 EEEEESSTB-TT-BEBTTTEEEECECGGCCTSCCCGGHHHHTTEEESS-B-TT-B-BG CS SAF 2 nVvvAardikaGeviteddltlkrpiashkigllpiaygeviGkratrdikaGepvha 59 +++ +++dik+Gevite +++ +rp +g+ p+ y e+iGk++ rd++ G+p+ 23726_FusoPortal_Gene_174 446 RSLFVSKDIKKGEVITENNIKSVRP----SNGISPKFYYEIIGKKVNRDLEFGTPLLF 499 68999********************....5*************************975 PP >> Acetyltransf_3 Acetyltransferase (GNAT) domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 29.6 0.0 2.6e-10 8.7e-07 18 138 .] 21 139 .. 9 139 .. 0.82 2 ? -2.0 0.0 1.4 4.8e+03 75 111 .. 356 392 .. 289 442 .. 0.62 Alignments for each domain: == domain 1 score: 29.6 bits; conditional E-value: 2.6e-10 Acetyltransf_3 18 lrsdpevmry..lpepltleeaeewlerawaeeseeegygwaivlkdtgv.iGsvglrdleeedeeaelgymlaps.fwgkGya 97 +r+ +ev++y ++ ++ee +wle+ +++ +++++ ++i ++ v iG++++ +++ee e+gy+l++ gkG++ 23726_FusoPortal_Gene_174 21 WRNKKEVRKYmfNNKIISEEEHLKWLES-LKR-ITNKEVFYIIY--NN-VkIGIASIDKISEEI--YEFGYYLNDLvPKGKGIG 97 999999***9855559999********6.443.34455557773..33.45*********9997..8*******88689***** PP Acetyltransf_3 98 teavralleyafeel.glerlvavidkeNeaSirvleklGfe 138 + + + ++f++ +l+++++++ +N aS ++l lG++ 23726_FusoPortal_Gene_174 98 KKVFYYFYLFLFHNYkKLNEIRCEVLNNNLASLNLLFFLGMK 139 ***********955549*******************999986 PP == domain 2 score: -2.0 bits; conditional E-value: 1.4 Acetyltransf_3 75 eeedeeaelgymlapsfwgkGyateavralleyafee 111 e+ +++ + l ++ g+ +at av++ +++ ++ 23726_FusoPortal_Gene_174 356 EDLAKRFSVKSGLSDHTMGHLVATTAVAMGAKVVEKH 392 3333345555555555556666666666555555444 PP >> Acetyltransf_1 Acetyltransferase (GNAT) family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 20.4 0.1 1.3e-07 0.00042 13 112 .. 31 133 .. 20 138 .. 0.82 Alignments for each domain: == domain 1 score: 20.4 bits; conditional E-value: 1.3e-07 HHHHHHHHHHHHHHH....HHSTTSEEEEEEETTEEEEEEEEEEEETTSTTEEEEEEEEESGGG.TTSSHHHHHHHHHHHHHHH CS Acetyltransf_1 13 fppnedaseeefees....fsleaegllvaeedgeliGfisllkipdekggvayieslaVspey.rgkGlGsaLlealeerare 91 ++n+ seee+ ++ ++++++ + + +++ iG +s+ ki + ++y +++ + gkG+G++ + ++++ 23726_FusoPortal_Gene_174 31 MFNNKIISEEEHLKWleslKRITNKEVFYIIYNNVKIGIASIDKISE----EIYEFGYYLNDLVpKGKGIGKKVFYYFYLFLFH 110 455666666774444444445566699999*****************....79999******999******************* PP T..TESEEEEEEETTCHHHHHHH CS Acetyltransf_1 92 e..gikkivlaveedndaairlY 112 + +++i ++v ++n a+ +l+ 23726_FusoPortal_Gene_174 111 NykKLNEIRCEVLNNNLASLNLL 133 *9555**************9986 PP >> DUF4352 Domain of unknown function (DUF4352) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 1.0 0.0 0.13 4.5e+02 69 107 .. 206 242 .. 168 253 .. 0.84 2 ? 10.1 0.0 0.0002 0.67 82 123 .. 411 452 .. 400 456 .. 0.79 Alignments for each domain: == domain 1 score: 1.0 bits; conditional E-value: 0.13 EE-BEEESCCCCTTB--EEEE-TT-EEEEEEEEEESTT- CS DUF4352 69 keydadlsasskaddlflgtlnpGkkvtGklvFdvpkgk 107 ++y++d+ + +++++ +++ G+ +Gk+ +d k++ 23726_FusoPortal_Gene_174 206 QTYTPDTLTLNCNNE--YFQIKDGTIWDGKILYDLYKEA 242 689999999999955..6689999999999999987765 PP == domain 2 score: 10.1 bits; conditional E-value: 0.0002 TB..--EEEE-TT-EEEEEEEEEESTT-CGEEEEE.HBB-TTTC CS DUF4352 82 dd..lflgtlnpGkkvtGklvFdvpkgkkkleLeyqtslsfsgk 123 ++ ++ +++++ +k+ Gk+++++++ k+ +L+++ sl f +k 23726_FusoPortal_Gene_174 411 SEfkEMVDNIRNVEKMIGKISYEISEKKE-ASLKFKRSL-FVSK 452 4444788999*****************99.*****9888.5554 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (509 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 973 (0.0582216); expected 334.2 (0.02) Passed bias filter: 488 (0.0292006); expected 334.2 (0.02) Passed Vit filter: 39 (0.00233365); 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.30u 0.24s 00:00:00.54 Elapsed: 00:00:00.34 # Mc/sec: 4352.00 // [ok]