# 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_1505.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1505 [L=257] Description: # 1604759 # 1605529 # 1 # ID=1_1505;partial=00;start_type=ATG;rbs_motif=AGGAG(G)/GGAGG;rbs_spacer=13-15bp;gc_cont=0.222 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 3e-35 121.3 0.4 4.5e-35 120.8 0.4 1.3 1 SpoU_methylase SpoU rRNA Methylase family ------ inclusion threshold ------ 0.073 13.3 0.0 0.36 11.0 0.0 2.2 2 SpoU_sub_bind RNA 2'-O ribose methyltransferase substrate b Domain annotation for each model (and alignments): >> SpoU_methylase SpoU rRNA Methylase family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 120.8 0.4 5.4e-39 4.5e-35 3 142 .] 108 247 .. 106 247 .. 0.93 Alignments for each domain: == domain 1 score: 120.8 bits; conditional E-value: 5.4e-39 EEEEES-SSHHHHHHHHHHHHHTT-SEEEEESGSSCTCGSHHHHHHCTTHHHHHHEEEESSHHHHHHHHCCCTEEEEEESTT- CS SpoU_methylase 3 ivvlveiedpgNlGaiaRtakafgaselvlvkprsadllnpeakraaagalevlelkvvdsleealkelkeegfalvatslkg 85 v+l++++dpgNlG+i+R ++a+++++++l+k ++d +n++++ra++g+ ++l ++ + +e+++ lke+++ +++t l++ 23726_FusoPortal_Gene_1505 108 LVILDDVSDPGNLGTIIRICDATNFKDIILTK-GTVDAYNEKVIRATMGSILNVNL-YYLEKSEIINFLKENNYSIISTYLDK 188 69****************************88.999*******************6.56677889**************9955 PP SEGC..CSHHHTTEEEEEEBTTTBS-HHHHCTSSEEEE----TTSSS--HHHHHHHHHH CS SpoU_methylase 86 ksls..eidlkekvalvfGnEdkGlsnelleladrlvkIpmsgeveSLNvavAvailly 142 + + +i+lkek a++fGnE++G+s+++++ d+ Ip+ ++eSLNvavA++i+ly 23726_FusoPortal_Gene_1505 189 TAILynKIELKEKNAVIFGNEGNGISDDFINITDYKTIIPILSNTESLNVAVATGIILY 247 54444699*************************************************99 PP >> SpoU_sub_bind RNA 2'-O ribose methyltransferase substrate binding # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.0 0.0 4.3e-05 0.36 45 70 .. 69 94 .. 48 99 .. 0.83 2 ? -2.1 0.0 0.56 4.7e+03 33 47 .. 168 182 .. 149 184 .. 0.78 Alignments for each domain: == domain 1 score: 11.0 bits; conditional E-value: 4.3e-05 SpoU_sub_bind 45 kvelvdkekLdklsgtanHQGvvaev 70 k + v+++++++ls+++n QG++ 23726_FusoPortal_Gene_1505 69 KEIVVTRKIFEELSSQENSQGIIILY 94 5578******************9865 PP == domain 2 score: -2.1 bits; conditional E-value: 0.56 SpoU_sub_bind 33 eellklakkrgikve 47 e+++ +k+++ +++ 23726_FusoPortal_Gene_1505 168 SEIINFLKENNYSII 182 489999999988775 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (257 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 1009 (0.0603758); expected 334.2 (0.02) Passed bias filter: 408 (0.0244136); expected 334.2 (0.02) Passed Vit filter: 23 (0.00137626); 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.26u 0.23s 00:00:00.49 Elapsed: 00:00:00.32 # Mc/sec: 2334.71 // [ok]