# 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_1224.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1224 [L=164] Description: # 1334647 # 1335138 # 1 # ID=1_1224;partial=00;start_type=ATG;rbs_motif=GGxGG;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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 9.8e-13 46.7 5.9 0.002 18.0 0.6 3.7 3 MORN_2 MORN repeat variant 0.0066 16.5 0.2 0.074 13.2 0.0 2.6 2 PPR PPR repeat ------ inclusion threshold ------ 0.019 14.6 0.1 0.032 13.8 0.1 1.4 1 PPR_1 PPR repeat Domain annotation for each model (and alignments): >> MORN_2 MORN repeat variant # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 13.8 0.1 8.4e-06 0.047 1 17 [. 70 86 .. 70 89 .. 0.88 2 ! 15.8 0.1 1.9e-06 0.01 2 18 .. 95 111 .. 94 114 .. 0.84 3 ! 18.0 0.6 3.7e-07 0.002 3 21 .. 120 138 .. 118 139 .. 0.94 Alignments for each domain: == domain 1 score: 13.8 bits; conditional E-value: 8.4e-06 MORN_2 1 leGeekeYyenGklkee 17 ++Ge+++Yy nG+l+ + 23726_FusoPortal_Gene_1224 70 KDGEFRRYYLNGNLLMR 86 69************865 PP == domain 2 score: 15.8 bits; conditional E-value: 1.9e-06 MORN_2 2 eGeekeYyenGklkeel 18 +G+++ Yy+nGk++ + 23726_FusoPortal_Gene_1224 95 DGIWEDYYPNGKVLMSG 111 8***********98665 PP == domain 3 score: 18.0 bits; conditional E-value: 3.7e-07 MORN_2 3 GeekeYyenGklkeeltYk 21 G +k+Y e+G+l+ e++Y+ 23726_FusoPortal_Gene_1224 120 GNWKYYNESGQLLGEVPYD 138 99****************7 PP >> PPR PPR repeat # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 13.2 0.0 1.3e-05 0.074 4 18 .. 82 96 .. 82 99 .. 0.89 2 ? 0.2 0.0 0.19 1.1e+03 5 14 .. 107 116 .. 106 118 .. 0.82 Alignments for each domain: == domain 1 score: 13.2 bits; conditional E-value: 1.3e-05 PPR 4 nslIsgycknGklee 18 n l++gyck+Gk++ 23726_FusoPortal_Gene_1224 82 NLLMRGYCKAGKVDG 96 89**********986 PP == domain 2 score: 0.2 bits; conditional E-value: 0.19 PPR 5 slIsgycknG 14 +l+sgy k + 23726_FusoPortal_Gene_1224 107 VLMSGYMKGN 116 69****9965 PP >> PPR_1 PPR repeat # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.8 0.1 5.8e-06 0.032 11 24 .. 82 95 .. 82 96 .. 0.93 Alignments for each domain: == domain 1 score: 13.8 bits; conditional E-value: 5.8e-06 HHHHHHHHHTS-HH CS PPR_1 11 NtLIdGlCraGrvd 24 N L+ G+C+aG+vd 23726_FusoPortal_Gene_1224 82 NLLMRGYCKAGKVD 95 7899*********9 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (164 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 849 (0.0508018); expected 334.2 (0.02) Passed bias filter: 641 (0.0383557); expected 334.2 (0.02) Passed Vit filter: 46 (0.00275251); 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.23s 00:00:00.47 Elapsed: 00:00:00.31 # Mc/sec: 1537.91 // [ok]