# 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_1609.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1609 [L=295] Description: # 1728323 # 1729207 # 1 # ID=1_1609;partial=00;start_type=ATG;rbs_motif=GGA/GAG/AGG;rbs_spacer=5-10bp;gc_cont=0.237 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 1e-27 96.8 0.0 1.6e-27 96.2 0.0 1.3 1 cobW CobW/HypB/UreG, nucleotide-binding domain 2.7e-07 30.3 0.0 4.6e-07 29.6 0.0 1.4 1 CobW_C Cobalamin synthesis protein cobW C-terminal 0.006 16.5 0.0 0.011 15.6 0.0 1.6 1 TsaE Threonylcarbamoyl adenosine biosynthesis pro ------ inclusion threshold ------ 0.013 15.9 0.0 0.03 14.7 0.0 1.7 1 AAA_18 AAA domain 0.025 14.2 0.0 0.04 13.5 0.0 1.3 1 Viral_helicase1 Viral (Superfamily 1) RNA helicase 0.034 14.4 0.0 0.042 14.1 0.0 1.4 1 AAA_16 AAA ATPase domain 0.17 11.9 0.0 0.23 11.5 0.0 1.2 1 AAA_33 AAA domain Domain annotation for each model (and alignments): >> cobW CobW/HypB/UreG, nucleotide-binding domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 96.2 0.0 6.6e-31 1.6e-27 3 176 .. 4 175 .. 2 176 .. 0.92 Alignments for each domain: == domain 1 score: 96.2 bits; conditional E-value: 6.6e-31 EEEEESTTCCHHHHHHHHHCT--SS-EEEEESS-SSCHCCCCCCHCCCT.EEEETTTTECCHCTT-HHHHHHHHHHHHTT-SE CS cobW 3 tvltGflGsGKTtLlnhllenraglklavivNefgevsiDaellretgveivelenGciCctlredllkelealeereerldv 85 ++++GflG+GKTt++++l +n +l+ v++Ne+ +++iD+++l+e++ ++ e+++GciCc+++ d+ +++++ ++ ++++ 23726_FusoPortal_Gene_1609 4 LLVSGFLGAGKTTFIKELPKN-INLEFVVLENEYADIGIDRDFLDEKNLNVWEMSEGCICCSMKGDFKSSIKRIYSEI-NPEY 84 789**************9866.899****************************************************9.9*** PP EEEEECTTTTHHHHHTT.SHHHHHCCCECEEEEEEETTTH.CCHCCHHHHHHHHHHTSSEEEEE-HHHHCHC..HHHHHHHHH CS cobW 86 iliEttGlaePapvaft.eeelksdltldavvavvDaknleaeeekieekaedqlaaadllvlnKtDlldeeallekleerlr 167 ++iE tG+ ++ + e ++d+++ + ++++D + +++ ++ ++ + d l+ + ++l+K + ++ ++e++++++ 23726_FusoPortal_Gene_1609 85 LIIEPTGVGMLSSIIENiREINNNDIEILSPLTLIDITSFNEYLKTFNDFFYDNLKNTGKVILTKLESFNSF-DVESIKSEIL 166 *******998888866634445569***************9999************************9987.679****999 PP HHHTTS-EE CS cobW 168 qlnpeakvl 176 + n++++++ 23726_FusoPortal_Gene_1609 167 KTNSNLEII 175 999999886 PP >> CobW_C Cobalamin synthesis protein cobW C-terminal domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 29.6 0.0 1.9e-10 4.6e-07 17 86 .. 232 292 .. 224 294 .. 0.88 Alignments for each domain: == domain 1 score: 29.6 bits; conditional E-value: 1.9e-10 CobW_C 17 lealleelpegilRaKGivwlagkpdalvlqgvggrlrlepagawwaaededrrseLVfIGrdldkeale 86 l++l++ ++i+RaKGiv+++g + ++ v + +++e + +++++++V+IG +ld + l+ 23726_FusoPortal_Gene_1609 232 LNRLVNGDFGKIYRAKGIVKIDGY--WGKFNLVYKNFEME-------PITDAKGTKIVIIGNNLDIKNLK 292 5666667899**************..**************.......334689************98886 PP >> TsaE Threonylcarbamoyl adenosine biosynthesis protein TsaE # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 15.6 0.0 4.6e-06 0.011 21 46 .. 3 28 .. 1 58 [. 0.90 Alignments for each domain: == domain 1 score: 15.6 bits; conditional E-value: 4.6e-06 EEEEEESTTSSHHHHHHHHHHHTT-- CS TsaE 21 vilLeGdLGaGKTtlvrglakalgve 46 ++l +G LGaGKTt++++l k++++e 23726_FusoPortal_Gene_1609 3 ILLVSGFLGAGKTTFIKELPKNINLE 28 57899**************9988776 PP >> AAA_18 AAA domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 14.7 0.0 1.3e-05 0.03 2 38 .. 5 37 .. 4 88 .. 0.79 Alignments for each domain: == domain 1 score: 14.7 bits; conditional E-value: 1.3e-05 AAA_18 2 lltGtpGvGKttlareLskrlgfgvnlrdvakengll 38 l++G G+GKtt+ +eL k +nl +v+ en++ 23726_FusoPortal_Gene_1609 5 LVSGFLGAGKTTFIKELPKN----INLEFVVLENEYA 37 689**************887....5666666665544 PP >> Viral_helicase1 Viral (Superfamily 1) RNA helicase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.5 0.0 1.7e-05 0.04 2 47 .. 5 53 .. 4 66 .. 0.61 Alignments for each domain: == domain 1 score: 13.5 bits; conditional E-value: 1.7e-05 Viral_helicase1 2 vvhGvpGcGKStlirkllrts..liltvnrptaelltegk.pklpnlsv 47 +v+G+ G+GK+t i++l ++ ++ ++ a + ++ +++nl+v 23726_FusoPortal_Gene_1609 5 LVSGFLGAGKTTFIKELPKNInlEFVVLENEYADIGIDRDfLDEKNLNV 53 89***************99945222222233333333333444444444 PP >> AAA_16 AAA ATPase domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 14.1 0.0 1.8e-05 0.042 25 62 .. 2 38 .. 1 157 [. 0.91 Alignments for each domain: == domain 1 score: 14.1 bits; conditional E-value: 1.8e-05 AAA_16 25 elvllsGesGvGKsslvrellralekdggiflrgkcdq 62 +++l+sG G GK++ ++el++ ++ + +++l+ ++ + 23726_FusoPortal_Gene_1609 2 RILLVSGFLGAGKTTFIKELPKNINLE-FVVLENEYAD 38 689*******************99998.7777666555 PP >> AAA_33 AAA domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.5 0.0 9.8e-05 0.23 2 30 .. 4 34 .. 3 80 .. 0.85 Alignments for each domain: == domain 1 score: 11.5 bits; conditional E-value: 9.8e-05 AAA_33 2 vllvGlpGsGKstlaerllkklg..avllss 30 +l+ G G+GK+t++++l k+ +v+l+ 23726_FusoPortal_Gene_1609 4 LLVSGFLGAGKTTFIKELPKNINleFVVLEN 34 6789**************9987755666654 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (295 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 613 (0.0366802); expected 334.2 (0.02) Passed bias filter: 350 (0.020943); expected 334.2 (0.02) Passed Vit filter: 56 (0.00335089); expected 16.7 (0.001) Passed Fwd filter: 7 (0.000418861); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 7 [number of targets reported over threshold] # CPU time: 0.29u 0.24s 00:00:00.53 Elapsed: 00:00:00.34 # Mc/sec: 2522.28 // [ok]