# 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_1830.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1830 [L=151] Description: # 1959514 # 1959966 # -1 # ID=1_1830;partial=00;start_type=ATG;rbs_motif=AGGAGG;rbs_spacer=5-10bp;gc_cont=0.274 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.6e-17 63.4 0.3 3.1e-17 62.5 0.3 1.4 1 SepF Cell division protein SepF ------ inclusion threshold ------ 0.075 13.0 0.0 0.2 11.6 0.0 1.9 1 HHH_2 Helix-hairpin-helix motif 1.2 9.1 5.8 1.8 8.5 0.0 2.6 2 DUF2335 Predicted membrane protein (DUF2335) Domain annotation for each model (and alignments): >> SepF Cell division protein SepF # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 62.5 0.3 5.5e-21 3.1e-17 3 72 .] 57 126 .. 56 126 .. 0.98 Alignments for each domain: == domain 1 score: 62.5 bits; conditional E-value: 5.5e-21 EEE-SSGGGHHHHHHHHHTT-EEEEE-CCCTCHHHHHHHHHHHHHHHHCT-EEEEECTTEEEEE-TT-EE CS SepF 3 vveprsfedaqeiadelkegkiVilnledldkeeakrivdflsgvvyaldGdiqkvgkkvflitPknvdi 72 ++ p++fed+++ia ++++ k++ +nle++ +++a+ri+dfl+g++ ++++++++k+v+ i+P+n+++ 23726_FusoPortal_Gene_1830 57 FIDPKQFEDCKKIATYIEKEKMITINLENIGPNVAQRIMDFLAGAMEIKNASFAQIAKHVYTIVPENMKV 126 689****************************************************************997 PP >> HHH_2 Helix-hairpin-helix motif # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.6 0.0 3.7e-05 0.2 38 56 .. 81 99 .. 61 107 .. 0.80 Alignments for each domain: == domain 1 score: 11.6 bits; conditional E-value: 3.7e-05 HHH_2 38 levddiGeeiAqsivefFe 56 +++++iG+ +Aq i +f++ 23726_FusoPortal_Gene_1830 81 INLENIGPNVAQRIMDFLA 99 3568*************97 PP >> DUF2335 Predicted membrane protein (DUF2335) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 0.4 3.2 0.11 5.9e+02 21 41 .. 22 42 .. 19 49 .. 0.75 2 ? 8.5 0.0 0.00033 1.8 9 23 .. 83 97 .. 81 105 .. 0.88 Alignments for each domain: == domain 1 score: 0.4 bits; conditional E-value: 0.11 DUF2335 21 iamaekEqehRheleekllda 41 i+++e+Eq+++ ++ ++a 23726_FusoPortal_Gene_1830 22 IEEMEQEQTSKALSKRQKMEA 42 789999999877666666655 PP == domain 2 score: 8.5 bits; conditional E-value: 0.00033 DUF2335 9 YesilPGaAeRIiam 23 e+i P+ A+RI++ 23726_FusoPortal_Gene_1830 83 LENIGPNVAQRIMDF 97 799**********85 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (151 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 1526 (0.0913116); expected 334.2 (0.02) Passed bias filter: 612 (0.0366204); expected 334.2 (0.02) Passed Vit filter: 36 (0.00215414); expected 16.7 (0.001) Passed Fwd filter: 3 (0.000179512); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 3 [number of targets reported over threshold] # CPU time: 0.24u 0.24s 00:00:00.48 Elapsed: 00:00:00.32 # Mc/sec: 1371.76 // [ok]