# 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_781.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_781 [L=166] Description: # 880709 # 881206 # 1 # ID=1_781;partial=00;start_type=ATG;rbs_motif=AGGA;rbs_spacer=5-10bp;gc_cont=0.217 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.9e-09 37.7 0.0 2.3e-09 37.4 0.0 1.2 1 Acetyltransf_1 Acetyltransferase (GNAT) family 2.2e-06 27.5 0.0 2.9e-06 27.2 0.0 1.2 1 Acetyltransf_10 Acetyltransferase (GNAT) domain 4.8e-06 26.8 0.5 5.8e-06 26.5 0.5 1.1 1 Acetyltransf_4 Acetyltransferase (GNAT) domain 0.0022 18.3 0.3 0.0072 16.6 0.1 1.8 2 Acetyltransf_7 Acetyltransferase (GNAT) domain ------ inclusion threshold ------ 0.061 13.3 0.0 0.12 12.4 0.0 1.6 1 Ribosomal_L14 Ribosomal protein L14p/L23e Domain annotation for each model (and alignments): >> Acetyltransf_1 Acetyltransferase (GNAT) family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 37.4 0.0 6.9e-13 2.3e-09 35 117 .] 60 141 .. 23 141 .. 0.86 Alignments for each domain: == domain 1 score: 37.4 bits; conditional E-value: 6.9e-13 EEEEEETTEEEEEEEEEEEETTSTTEEEEEEEEESGGGTTSSHHHHHHHHHHHHHHHTTESEEEEEEETTCHHHHHHHHHTTE CS Acetyltransf_1 35 llvaeedgeliGfisllkipdekggvayieslaVspeyrgkGlGsaLlealeerareegikkivlaveedndaairlYeklGF 117 + ++++++e++G + +l +++ +++++ i ++ ++ + +++GlG++ +le+ ++++ + ki + v +dn+ +++k F 23726_FusoPortal_Gene_781 60 YRLISFNDEILGVVDYLTDYP-EKNTILIGFFIIKNDKQKQGLGTKIFGYLEKSFKNKKFLKIRIGVLVDNEIGLSFWKKQNF 141 334679***************.8******************************99999**999***************88887 PP >> Acetyltransf_10 Acetyltransferase (GNAT) domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 27.2 0.0 8.5e-10 2.9e-06 29 126 .. 56 164 .. 32 165 .. 0.78 Alignments for each domain: == domain 1 score: 27.2 bits; conditional E-value: 8.5e-10 Acetyltransf_10 29 gealllvaeeggkivGviel.rdgs.....hisllfVdpeyqgqGigkaLleaveeeaekdglkl.eltvnasp.yavpfYekl 104 +++ + + +++i+Gv+ + +d i+ ++++++ q+qG+g++++ ++e+ +++++ + ++ v ++ f +k+ 23726_FusoPortal_Gene_781 56 NQKNYRLISFNDEILGVVDYlTDYPekntiLIGFFIIKNDKQKQGLGTKIFGYLEKSFKNKKFLKiRIGVLVDNeIGLSFWKKQ 139 4444444444678*****7767776777889**************************999999999******99999******* PP Acetyltransf_10 105 GFkavgee...qeenGirfvpMeke 126 Fk+++++ +e++ +++ Meke 23726_FusoPortal_Gene_781 140 NFKEIERKflkFEKSEKEVIVMEKE 164 *****77511156666677777665 PP >> Acetyltransf_4 Acetyltransferase (GNAT) domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 26.5 0.5 1.7e-09 5.8e-06 31 133 .. 36 142 .. 9 161 .. 0.82 Alignments for each domain: == domain 1 score: 26.5 bits; conditional E-value: 1.7e-09 Acetyltransf_4 31 vssveefekrlesl.d..sekqylflvfe.sgrlvGyatlrqfkk.akdwkaelsvylardnnkeglgrellnalvevlkemni 109 +s+ +ek +e++ + +++q + + +++++G + ++ k+ + ++ +++d++k+glg ++ +l++ +k+ + 23726_FusoPortal_Gene_781 36 PPSISNVEKDIEAIpNgvQKNQKNYRLISfNDEILGVVDYL-TDYpEKNTILIGFFIIKNDKQKQGLGTKIFGYLEKSFKNKKF 118 388999999999998668999999999999********998.56645777777788889************************* PP Acetyltransf_4 110 enllaeiasnninaivflkalGFe 133 +++ + + ++n+ f+k+ F+ 23726_FusoPortal_Gene_781 119 LKIRIGVLVDNEIGLSFWKKQNFK 142 *******************98775 PP >> Acetyltransf_7 Acetyltransferase (GNAT) domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -0.9 0.0 0.65 2.2e+03 58 74 .. 12 27 .. 3 28 .. 0.64 2 ! 16.6 0.1 2.2e-06 0.0072 11 76 .] 66 143 .. 48 143 .. 0.66 Alignments for each domain: == domain 1 score: -0.9 bits; conditional E-value: 0.65 EEETHHHHHHHHHHTTE CS Acetyltransf_7 58 lelnaekaikfYeklGF 74 el++ +a++ Y k+ + 23726_FusoPortal_Gene_781 12 NELKN-EALNIYLKNDY 27 23433.79999999876 PP == domain 2 score: 16.6 bits; conditional E-value: 2.2e-06 TTEEEEEEEEEEES.....TCTEEEEE-GGGTTSSHHHHHHHH.......TTSE.EEEE.EEEETHHHHHHHHHHTTEEE CS Acetyltransf_7 11 dgkivGfialsled.....aeitrlaVdpeyrgkGiGraLlea.......aekg.iile.qlelnaekaikfYeklGFee 76 +++i+G + + + i+++ ++ + +++G+G+++ + + + +++ + +++e ++f++k F+e 23726_FusoPortal_Gene_781 66 NDEILGVVDYLTDYpekntILIGFFIIKNDKQKQGLGTKIFGYleksfknK--KfLKIRiGVLVDNEIGLSFWKKQNFKE 143 47888877766666776533333388**************98744333321..221122233444446677*****9986 PP >> Ribosomal_L14 Ribosomal protein L14p/L23e # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.4 0.0 3.5e-05 0.12 88 114 .. 95 121 .. 56 124 .. 0.89 Alignments for each domain: == domain 1 score: 12.4 bits; conditional E-value: 3.5e-05 -TTS-BSS-.B-SBEEHHHHHH.-HHH CS Ribosomal_L14 88 nkkgeplGtrifgpvarelrekkflki 114 +k+++ lGt+ifg + ++kkflki 23726_FusoPortal_Gene_781 95 DKQKQGLGTKIFGYLEKSFKNKKFLKI 121 4568899*******************9 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (166 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 798 (0.0477501); expected 334.2 (0.02) Passed bias filter: 363 (0.0217209); expected 334.2 (0.02) Passed Vit filter: 18 (0.00107707); 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.23u 0.22s 00:00:00.45 Elapsed: 00:00:00.31 # Mc/sec: 1556.67 // [ok]