# 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_128.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_128 [L=164] Description: # 145151 # 145642 # -1 # ID=1_128;partial=00;start_type=TTG;rbs_motif=AGGAGG;rbs_spacer=5-10bp;gc_cont=0.335 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.1e-28 97.0 0.2 5.5e-28 96.2 0.2 1.4 1 Ribosomal_S5_C Ribosomal protein S5, C-terminal domain 4.2e-27 93.8 1.5 6.9e-27 93.1 1.5 1.4 1 Ribosomal_S5 Ribosomal protein S5, N-terminal domain ------ inclusion threshold ------ 0.041 13.9 0.7 0.041 13.9 0.7 2.4 2 Antimicrobial_8 Uperin family 0.091 12.9 1.7 0.48 10.6 0.3 2.5 2 Dodecin Dodecin Domain annotation for each model (and alignments): >> Ribosomal_S5_C Ribosomal protein S5, C-terminal domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 96.2 0.2 1.3e-31 5.5e-28 2 72 .] 84 154 .. 83 154 .. 0.99 Alignments for each domain: == domain 1 score: 96.2 bits; conditional E-value: 1.3e-31 EEETTEEEEEEE--TT--EES-HHHHHHHHHHT--SEE---BS-SSHHHHHHHHHHHHHHTTTT-.-GGGG CS Ribosomal_S5_C 2 gkfgavkvlLkpaprGtGiiaggavravlelaGikDvlakslgsrnpinvvkAtfeALkklrspeevaalR 72 gk+ga+++ + pa +GtG+iag a r++lel+G++D+l+k++gsrn++nv++At+eALk lr++e++aalR 23726_FusoPortal_Gene_128 84 GKWGATTLWMAPAYEGTGVIAGSASREILELVGVHDILTKIKGSRNKHNVARATVEALKLLRTAEQIAALR 154 89*******************************************************************98 PP >> Ribosomal_S5 Ribosomal protein S5, N-terminal domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 93.1 1.5 1.7e-30 6.9e-27 2 64 .. 9 71 .. 8 72 .. 0.98 Alignments for each domain: == domain 1 score: 93.1 bits; conditional E-value: 1.7e-30 EEEEEEEEEEEEEETTEEEEEEEEEEEEE-SSSEEEEEEEEESSHHHHHHHHHHHHHCSEEE- CS Ribosomal_S5 2 kekvlkikrvqkqtkgGrrtrfsalvvvGdgnGhvGlGvgkakevatAirkAiekAkknlvpV 64 +ek+lki+rv+k+tkgGr ++fs+l++vGdg+G++GlG+gka+ v++AirkAi +Akkn+v++ 23726_FusoPortal_Gene_128 9 QEKLLKISRVSKTTKGGRTISFSVLAAVGDGEGKIGLGLGKANGVPDAIRKAIAAAKKNIVKI 71 89***********************************************************98 PP >> Antimicrobial_8 Uperin family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.9 0.7 9.7e-06 0.041 1 17 [] 52 69 .. 52 69 .. 0.77 2 ? -0.7 0.1 0.43 1.8e+03 1 9 [. 116 124 .. 116 125 .. 0.87 Alignments for each domain: == domain 1 score: 13.9 bits; conditional E-value: 9.7e-06 Antimicrobial_8 1 GVGDliRKiVs.aIKNVV 17 GV D+iRK+++ a KN+V 23726_FusoPortal_Gene_128 52 GVPDAIRKAIAaAKKNIV 69 99******7650567998 PP == domain 2 score: -0.7 bits; conditional E-value: 0.43 Antimicrobial_8 1 GVGDliRKi 9 GV D+ Ki 23726_FusoPortal_Gene_128 116 GVHDILTKI 124 899999998 PP >> Dodecin Dodecin # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.6 0.3 0.00011 0.48 13 34 .. 51 72 .. 49 74 .. 0.89 2 ? 0.3 0.0 0.19 8.1e+02 25 41 .. 145 161 .. 127 162 .. 0.74 Alignments for each domain: == domain 1 score: 10.6 bits; conditional E-value: 0.00011 S-HHHHHHHHHHHHHHHSSTEE CS Dodecin 13 kSwedAiqnAiaeAsktlrnir 34 +++ dAi++Aia A+k + +i+ 23726_FusoPortal_Gene_128 51 NGVPDAIRKAIAAAKKNIVKIS 72 6788************998885 PP == domain 2 score: 0.3 bits; conditional E-value: 0.19 HHHHHSSTEEEEEEEEE CS Dodecin 25 eAsktlrnirwveVket 41 + ++ + +r+ eVk++ 23726_FusoPortal_Gene_128 145 RTAEQIAALRGLEVKDI 161 56778888888898875 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (164 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 590 (0.035304); expected 334.2 (0.02) Passed bias filter: 458 (0.0274055); expected 334.2 (0.02) Passed Vit filter: 25 (0.00149593); expected 16.7 (0.001) Passed Fwd filter: 4 (0.000239349); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 4 [number of targets reported over threshold] # CPU time: 0.26u 0.27s 00:00:00.53 Elapsed: 00:00:00.34 # Mc/sec: 1402.22 // [ok]