# 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_694.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_694 [L=264] Description: # 786098 # 786889 # 1 # ID=1_694;partial=00;start_type=ATG;rbs_motif=AGGAGG;rbs_spacer=5-10bp;gc_cont=0.317 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.7e-32 111.6 0.0 3e-32 110.8 0.0 1.4 1 OAM_dimer Dimerisation domain of d-ornithine 4,5-aminomuta 6.8e-19 68.0 0.2 1.2e-18 67.2 0.2 1.3 1 B12-binding B12 binding domain ------ inclusion threshold ------ 0.013 14.6 0.0 0.015 14.3 0.0 1.2 1 DUF1749 Protein of unknown function (DUF1749) 0.041 13.4 0.0 0.11 12.0 0.0 1.7 2 AA_kinase Amino acid kinase family Domain annotation for each model (and alignments): >> OAM_dimer Dimerisation domain of d-ornithine 4,5-aminomutase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 110.8 0.0 7.2e-36 3e-32 1 77 [. 21 97 .. 21 98 .. 0.98 Alignments for each domain: == domain 1 score: 110.8 bits; conditional E-value: 7.2e-36 EE---SSSTT--EEEEEEEESS-SCHHHHHHHHHHHHCTTSEEEEEEEEEEEETTEEEEEEEEE--S-EEGGG---- CS OAM_dimer 1 likPygdtlnDGvvqlsftLPvpaeeeAeeAAlelakkmGledpeVvhkeelheeftffvvyGkvshsvdlselkvp 77 +i+Pygdt+nDG+vq+sftLPv ++e++ eAAl+la+kmG+ +p+V+++e+l++ef+f+vvyG +s+svd++++kv+ 23726_FusoPortal_Gene_694 21 QIRPYGDTMNDGKVQMSFTLPVACNEKGIEAALQLARKMGFVNPAVAFSEALDKEFSFYVVYGATSFSVDYTAIKVQ 97 59*************************************************************************97 PP >> B12-binding B12 binding domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 67.2 0.2 2.8e-22 1.2e-18 3 117 .. 123 246 .. 121 249 .. 0.93 Alignments for each domain: == domain 1 score: 67.2 bits; conditional E-value: 2.8e-22 B12-binding 3 vllatvggdlhplGlavlaaalrdagf.........eVqllganvenedivaairaekpdvvglSalltsqleavkelvel... 74 ++a+ g+d+h++G++++ +++ ag+ + ++lg +v+ne+++++++++k+d++++S+++t++ ++++l++l 23726_FusoPortal_Gene_694 123 MVGASTGTDAHTVGIDAIMNMKGYAGHyglerykgvRAYNLGSQVPNEEFIKKAIELKADALLVSQTVTQKDVHIENLTNLvel 206 5899***********************************************************************999999999 PP B12-binding 75 ..lkrirprvkvvvGGahpsldpeetleak.lgiddvvvgegedal 117 +++ r+++++++GG ++++ ++ ak lg +d+ +g g++a 23726_FusoPortal_Gene_694 207 leAEGLRDKIILIAGG--ARITNDL---AKeLG-YDAGFGPGKYAD 246 ****************..5677777...555**.9*******9885 PP >> DUF1749 Protein of unknown function (DUF1749) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 14.3 0.0 3.6e-06 0.015 56 120 .. 69 131 .. 59 217 .. 0.77 Alignments for each domain: == domain 1 score: 14.3 bits; conditional E-value: 3.6e-06 DUF1749 56 adalde.agwslvq..lqlsSsyggwGtgsLdrDteeleklveYlkktksGeekkekvvlmGHStGsq 120 ++ald+ ++ +v +s y++ + +L+ Dt +++++ +Y++++ G + vv++G StG+ 23726_FusoPortal_Gene_694 69 SEALDKeFSFYVVYgaTSFSVDYTAIKVQALEIDTMDMHECEKYIEEN-FG----REVVMVGASTGTD 131 556644134444431134556799999*******************99.56....89*********86 PP >> AA_kinase Amino acid kinase family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -1.9 0.0 0.46 1.9e+03 33 57 .. 74 96 .. 45 143 .. 0.61 2 ? 12.0 0.0 2.6e-05 0.11 11 73 .. 189 250 .. 185 261 .. 0.77 Alignments for each domain: == domain 1 score: -1.9 bits; conditional E-value: 0.46 HTTEEEEEE--HHHHHHHHHHTCCC CS AA_kinase 33 sgnklvvVhGggsftdgllalksgl 57 ++ vV G+ sf + a+ ++ 23726_FusoPortal_Gene_694 74 KEFSFYVVYGATSFSVDYTAI--KV 96 444455666665555555555..33 PP == domain 2 score: 12.0 bits; conditional E-value: 2.6e-05 CCCCTCC..HHHHHHHHHHHHHHHHTTEEEEEE--HHHHHHHHHHTCCCCCHHHC....E..HHH CS AA_kinase 11 sltdkee..aslrrlaeqiaalkesgnklvvVhGggsftdgllalksglssgelaaglrstleea 73 +t+k+ ++l +l+e ++a k++++ Gg +t+ l+++ g+++g+ g++ +++++ 23726_FusoPortal_Gene_694 189 TVTQKDVhiENLTNLVELLEAE-GLRDKIILIAGGARITNDLAKE-LGYDAGFGP-GKYADDVAT 250 5677777777888888888888.888*******************.****99998.555555444 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (264 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 628 (0.0375778); expected 334.2 (0.02) Passed bias filter: 513 (0.0306965); expected 334.2 (0.02) Passed Vit filter: 46 (0.00275251); expected 16.7 (0.001) Passed Fwd filter: 4 (0.000239349); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 4 [number of targets reported over threshold] # CPU time: 0.25u 0.23s 00:00:00.48 Elapsed: 00:00:00.31 # Mc/sec: 2475.67 // [ok]