# 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_1778.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1778 [L=522] Description: # 1900194 # 1901759 # 1 # ID=1_1778;partial=00;start_type=ATG;rbs_motif=AGGAG;rbs_spacer=5-10bp;gc_cont=0.287 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 0.00051 19.5 1.5 0.2 11.0 0.1 3.1 3 Peptidase_S9 Prolyl oligopeptidase family 0.007 16.1 0.0 0.012 15.4 0.0 1.4 1 Abhydrolase_3 alpha/beta hydrolase fold ------ inclusion threshold ------ 0.048 12.5 0.3 0.16 10.7 0.0 2.1 2 COesterase Carboxylesterase family 0.058 13.1 0.0 0.27 10.9 0.0 2.0 2 Sigma70_r2 Sigma-70 region 2 Domain annotation for each model (and alignments): >> Peptidase_S9 Prolyl oligopeptidase family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 11.0 0.1 4.7e-05 0.2 8 80 .. 128 199 .. 122 206 .. 0.75 2 ? 0.4 0.1 0.082 3.4e+02 83 139 .. 220 270 .. 216 281 .. 0.72 3 ! 5.0 0.0 0.0034 14 131 191 .. 446 506 .. 424 520 .. 0.79 Alignments for each domain: == domain 1 score: 11.0 bits; conditional E-value: 4.7e-05 HHHHCCT-EEEEE--TTSSTTHCHHHHTTSGGGGHHHHHHHHHHHHHHHHTT...S--GGGEEEEEETHHHHHHHH CS Peptidase_S9 8 qlladrGyvvvvanirGsggyGeafhdagkrklgqnefdDfiaaaealakeg...yvdserlaitGgSyGGyltga 80 + l ++Gyvv+++ rG + ++ + gk + + D+ aa+++l ++ d+++++ G+S+GG l++ 23726_FusoPortal_Gene_1778 128 YAL-SKGYVVAAPGARGRTLTDDKGNYIGK---APAAIVDLKAAVRYLYLNDevmPGDANKIISNGTSAGGALSAL 199 444.99999999999999966555444433...456788999***9996555325569*************99875 PP == domain 2 score: 0.4 bits; conditional E-value: 0.082 HH-GGG-SEEEEESE---CCCHHCSTTSTHHHCCCC--TTTCHHHHHHHHHH-CGGG CS Peptidase_S9 83 nqrpdlFkaavaevpvvdvlrydsdtdllteeymewgnpeeneegykelspyspldn 139 + r+d+F a+ a +p++++ + ds + e+m +g + +++ ++ + +++ ++ 23726_FusoPortal_Gene_1778 220 ETRDDIF-AVSAYCPITNLENADSA---Y--EWMYNGVNSYSRMEFTRNTSAQEYND 270 4477888.56799**********99...6..34555555567777777777777665 PP == domain 3 score: 5.0 bits; conditional E-value: 0.0034 HHHH-CGGG-----EEEEEEETT-SSS-THHHHHHHHHHH--S-SEEEEEETTS-SSTTSS CS Peptidase_S9 131 lspyspldnvkaypplLlihGlaDdrvppaealklvaaLqkkgknvlllifpdegHgvgkp 191 ++p+ + d + ++ + hG D+ + +a ++ l+ +L++ gk v++ ++gHg+ + 23726_FusoPortal_Gene_1778 446 MNPMYYIDSNTSTKYWRIRHGAVDKDTSLAIPAILALKLKNSGKVVDFASPWGQGHGGDYD 506 5555555555555555678*************************************88744 PP >> Abhydrolase_3 alpha/beta hydrolase fold # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 15.4 0.0 2.9e-06 0.012 44 91 .. 156 203 .. 155 289 .. 0.78 Alignments for each domain: == domain 1 score: 15.4 bits; conditional E-value: 2.9e-06 -TTHHHHHHHHHHHHHHHTHHHHTETTCEEEEEEETHHHHHHHHHHHH CS Abhydrolase_3 44 pfPaavedalaalrwlaeqawelgadpsrivlaGdSaGGnlaaavalr 91 ++Paa+ d aa+r+l +++ + d+++i+ G+SaGG l+a + 23726_FusoPortal_Gene_1778 156 KAPAAIVDLKAAVRYLYLNDEVMPGDANKIISNGTSAGGALSALLGAS 203 68*****************99999*****************9987654 PP >> COesterase Carboxylesterase family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.7 0.0 3.9e-05 0.16 162 211 .. 160 209 .. 54 314 .. 0.81 2 ? -1.1 0.0 0.15 6.1e+02 311 403 .. 393 485 .. 377 491 .. 0.53 Alignments for each domain: == domain 1 score: 10.7 bits; conditional E-value: 3.9e-05 HHHHHHHHHHHHHHHGGGGTEEEEEEEEEEETHHHHHHHHHHHSGGGTTT CS COesterase 162 GllDqvlALrWVkenIasFGGdpenvtlfGesAGAasvsllllsplsrgl 211 +++D+ +A+r++ n + Gd +++ G+sAG+a+ +ll +s s++ 23726_FusoPortal_Gene_1778 160 AIVDLKAAVRYLYLNDEVMPGDANKIISNGTSAGGALSALLGASGNSQDY 209 789**************************************999998863 PP == domain 2 score: -1.1 bits; conditional E-value: 0.15 EEEEETTTTHHHHCCCSCTT...SSTTCHHHHHHHHHHHHHHHCHHCHHHHHHHHHHHHTTTTTSSTCHHHHHHHHHHHHHHH CS COesterase 311 liGvnkdegllflaylldee...kavkkkeeklsesllkdelyeekekkeelseaileeytdeadrddaeksrkalvelltdv 390 l+G k ++++f ++ +d + ++ ++ k+++++ k + +++ ++++ +++ +y +++ + + r + ++ +d+ 23726_FusoPortal_Gene_1778 393 LFGDKKSDSKHFTKFSFDINnkaAIDYFRNGKFNDKNNKISVPKMA--DKNIIKMMNPMYY-IDSNTSTKYWRIRHGAVDKDT 472 5555666666666655554411122222222333333222333222..5666666666665.344445555555555666777 PP HTHHHHHHHHHHH CS COesterase 391 lfkvpvlrlaklh 403 +++p++ ++kl 23726_FusoPortal_Gene_1778 473 SLAIPAILALKLK 485 8888877777665 PP >> Sigma70_r2 Sigma-70 region 2 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.9 0.0 6.4e-05 0.27 44 67 .. 319 342 .. 313 345 .. 0.83 2 ? -1.5 0.0 0.49 2.1e+03 44 67 .. 401 422 .. 397 424 .. 0.78 Alignments for each domain: == domain 1 score: 10.9 bits; conditional E-value: 6.4e-05 TTCHHHHHHHHHHHHHHHHHHHHT CS Sigma70_r2 44 ergvsfsawllriarnaaidrlRk 67 + + sf+++l i+rn a ++lR 23726_FusoPortal_Gene_1778 319 SGNGSFKTYLSIIIRNSANRALRE 342 4444******************97 PP == domain 2 score: -1.5 bits; conditional E-value: 0.49 TTCHHHHHHHHHHHHHHHHHHHHT CS Sigma70_r2 44 ergvsfsawllriarnaaidrlRk 67 +++ f+ + + i aaid++R 23726_FusoPortal_Gene_1778 401 SKH--FTKFSFDINNKAAIDYFRN 422 444..888888889999*****96 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (522 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 629 (0.0376376); expected 334.2 (0.02) Passed bias filter: 371 (0.0221996); expected 334.2 (0.02) Passed Vit filter: 29 (0.00173528); 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.29u 0.23s 00:00:00.52 Elapsed: 00:00:00.31 # Mc/sec: 4895.07 // [ok]