# 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_902.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_902 [L=240] Description: # 1004551 # 1005270 # 1 # ID=1_902;partial=00;start_type=ATG;rbs_motif=AGGAGG;rbs_spacer=5-10bp;gc_cont=0.212 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- ------ inclusion threshold ------ 0.031 14.3 0.0 0.053 13.6 0.0 1.4 1 AbiEi_4 Transcriptional regulator, AbiEi antitoxin 0.053 13.4 1.9 0.19 11.6 0.1 2.3 2 Peptidase_M99_m beta-barrel domain of carboxypeptidase M99 Domain annotation for each model (and alignments): >> AbiEi_4 Transcriptional regulator, AbiEi antitoxin # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.6 0.0 6.3e-06 0.053 13 49 .] 19 59 .. 15 59 .. 0.92 Alignments for each domain: == domain 1 score: 13.6 bits; conditional E-value: 6.3e-06 AbiEi_4 13 lfttreleaaGispsa....lrrLvrkGwLvrlrrGvYali 49 lf+ +++ + +is++a l r + G++ r+ +G+Y l+ 23726_FusoPortal_Gene_902 19 LFSCDDFFEKNISRTAvlksLERIEKTGKIKRVEKGFYYLP 59 788899999*****************************986 PP >> Peptidase_M99_m beta-barrel domain of carboxypeptidase M99 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.6 0.1 2.3e-05 0.19 22 65 .. 52 96 .. 47 101 .. 0.77 2 ? 0.6 0.1 0.062 5.2e+02 24 63 .. 188 230 .. 167 235 .. 0.76 Alignments for each domain: == domain 1 score: 11.6 bits; conditional E-value: 2.3e-05 EEEEEEEETT--GGGS-EEESSSSEEEEEETTE.EEEEETTEEEE CS Peptidase_M99_m 22 klnyfPlPknskissiPlesrayilGllkkdqe.vllkyGnklmt 65 + ++ lP +s +iP+ s+ yi +l kd+e + Gn l 23726_FusoPortal_Gene_902 52 EKGFYYLPIKSIFGEIPISSNEYIKKFLYKDKEkIGYISGNNLYN 96 556889*****99***********998877765144556776665 PP == domain 2 score: 0.6 bits; conditional E-value: 0.062 EEEEEETT--GGGS-EEESSSSEEEEEE...TTEEEEEETTEE CS Peptidase_M99_m 24 nyfPlPknskissiPlesrayilGllkk...dqevllkyGnkl 63 +y+Pl + + si e r + +lk+ +++ ++k G kl 23726_FusoPortal_Gene_902 188 EYYPLFVSAILGSILEEKRFKKVQILKEkiaKKKTFFKLGIKL 230 5888888888888888888888888886222577888888887 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (240 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 735 (0.0439804); expected 334.2 (0.02) Passed bias filter: 306 (0.0183102); expected 334.2 (0.02) Passed Vit filter: 20 (0.00119674); 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.24u 0.24s 00:00:00.48 Elapsed: 00:00:00.32 # Mc/sec: 2180.27 // [ok]