# 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_209.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_209 [L=266] Description: # 222207 # 223004 # 1 # ID=1_209;partial=00;start_type=ATG;rbs_motif=GGAGG;rbs_spacer=5-10bp;gc_cont=0.318 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 3.6e-28 97.2 0.1 7.6e-28 96.2 0.1 1.6 1 RepA_N Replication initiator protein A (RepA) N-terminus 0.0056 16.5 0.0 0.019 14.8 0.0 1.9 1 HTH_36 Helix-turn-helix domain ------ inclusion threshold ------ 0.052 13.3 0.1 0.085 12.6 0.1 1.2 1 Telethonin Telethonin protein Domain annotation for each model (and alignments): >> RepA_N Replication initiator protein A (RepA) N-terminus # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 96.2 0.1 1.4e-31 7.6e-28 2 76 .] 15 89 .. 14 89 .. 0.99 Alignments for each domain: == domain 1 score: 96.2 bits; conditional E-value: 1.4e-31 RepA_N 2 sfyrlpkvLledeeykklsaeakvlYalLkdRlelSiknnwiDeegrvYliftieelaellncskqkvikilkeL 76 +f+++p++L++ ee+k lsaea++lY++L++R+++S+knnwiD+egrv+++f +ee++++ n+sk ++ik+l+eL 23726_FusoPortal_Gene_209 15 NFLKVPEILVDGEEFKGLSAEAIILYSMLLKRTGMSFKNNWIDKEGRVFIYFAVEEIMKRRNISKPTAIKTLDEL 89 7************************************************************************98 PP >> HTH_36 Helix-turn-helix domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 14.8 0.0 3.4e-06 0.019 1 50 [. 30 90 .. 30 94 .. 0.85 Alignments for each domain: == domain 1 score: 14.8 bits; conditional E-value: 3.4e-06 HTH_36 1 nrLspsaklllgslasfa.........nknggcfp..sqetlakrlglsknTvsraikeLe 50 ++Ls +a++l+ +l++ +k+g +f +e ++kr ++sk T+++ ++eL 23726_FusoPortal_Gene_209 30 KGLSAEAIILYSMLLKRTgmsfknnwiDKEGRVFIyfAVEEIMKRRNISKPTAIKTLDELD 90 579**************977777777777777666556*********************95 PP >> Telethonin Telethonin protein # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.6 0.1 1.5e-05 0.085 7 56 .. 121 170 .. 118 192 .. 0.77 Alignments for each domain: == domain 1 score: 12.6 bits; conditional E-value: 1.5e-05 Telethonin 7 svreeneakresytaeWkdlslktrpeercslkeedesrretlkrqgqar 56 +v+e++ +k +++t+e kd+ l ++ +e +++k++d+ e k + r 23726_FusoPortal_Gene_209 121 QVKENDFQKLKNFTSEVKDVDLRSKENELQEVKNVDRNYIENNKSKYSKR 170 79**********************************99888765544443 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (266 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 1360 (0.0813787); expected 334.2 (0.02) Passed bias filter: 525 (0.0314146); expected 334.2 (0.02) Passed Vit filter: 49 (0.00293202); 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.28u 0.24s 00:00:00.52 Elapsed: 00:00:00.33 # Mc/sec: 2343.24 // [ok]