# 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_626.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_626 [L=191] Description: # 704642 # 705214 # 1 # ID=1_626;partial=00;start_type=ATG;rbs_motif=GGA/GAG/AGG;rbs_spacer=5-10bp;gc_cont=0.225 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.5e-54 183.1 1.1 3.9e-54 182.9 1.1 1.0 1 CoaE Dephospho-CoA kinase 2.8e-08 34.3 1.5 4.2e-08 33.7 1.3 1.4 1 AAA_18 AAA domain 0.00041 20.4 2.8 0.004 17.2 0.2 2.2 2 Cytidylate_kin2 Cytidylate kinase-like family 0.0052 16.9 0.2 0.0084 16.3 0.0 1.6 2 AAA_28 AAA domain ------ inclusion threshold ------ 0.023 14.7 2.1 0.039 14.0 0.2 2.1 2 SKI Shikimate kinase 0.055 13.5 0.7 0.055 13.5 0.7 1.8 2 FokI_N Restriction endonuclease FokI, recognition d 0.1 12.7 0.1 0.18 11.8 0.1 1.8 1 AAA_33 AAA domain 0.22 11.9 6.6 0.19 12.0 2.6 2.1 2 AAA_17 AAA domain 0.36 10.4 5.8 1.3 8.6 5.7 2.1 1 PRK Phosphoribulokinase / Uridine kinase family 1.4 9.8 6.1 0.67 10.8 1.4 2.4 2 DUF4284 Immunity protein 22 1.7 8.6 8.4 1.4 8.9 0.3 3.9 4 SapB_1 Saposin-like type B, region 1 Domain annotation for each model (and alignments): >> CoaE Dephospho-CoA kinase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 182.9 1.1 2.6e-57 3.9e-54 1 178 [. 1 177 [. 1 178 [. 0.97 Alignments for each domain: == domain 1 score: 182.9 bits; conditional E-value: 2.6e-57 EEEEEEESTTSSHHHHHHHHHHTT-EEEEHHHHHHHTTCCTHHHHHHHHHHHGGGCB.TTSSE-HHHHHHHHHTSHHHHHHHHH CS CoaE 1 livgltGgiasGKstvanlfadlgvpivDaDlvArqvvepgsealkkiveeFgeeildedGeLdReaLrkkvFadpeekkaLna 84 +i+gltGgiasGKstv++++a+ g+++ DaD++A+ + e++s + ++i ++Fg++ild++ + d+++L++ vF+++e+ k Lna 23726_FusoPortal_Gene_626 1 MIIGLTGGIASGKSTVSKYLAEKGFKVYDADKIAKDISEKKS-VQEEIISTFGNKILDKNRNVDKKKLKEIVFENKEKLKQLNA 83 69***********************************99876.568899*********************************** PP HHHHHHHHHHHHHHHH--SSEEEEE-TTTTTTTGGGGSSEEEEEE--HHHHHHHHHHCCTSTHHHHHHHHHHS--HHHHHHH-S CS CoaE 85 llhPlirreikkelakleskyvlldvPLLvEkqlekliervlvvdvelelqveRlvardgvseeevkaiiaaQlseeerlalad 168 ++hP + + k+ k +s+++++dvPLL+E++++k ++++lv+ + e q++R+v+rd+++++ +++ii++Qls+eer++ ad 23726_FusoPortal_Gene_626 84 IIHPKVIHFYKELKGKNTSEIIIFDVPLLFESGIDKFCDKILVIISDYEIQLNRIVERDKIDRDLAEKIIKSQLSNEERIKKAD 167 *******999999999999***************************************************************** PP EEEE-.SSHH CS CoaE 169 vvldNdssla 178 vv++N+s+l+ 23726_FusoPortal_Gene_626 168 VVIENNSNLE 177 ******9987 PP >> AAA_18 AAA domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 33.7 1.3 2.8e-11 4.2e-08 1 122 [. 3 150 .. 3 164 .. 0.69 Alignments for each domain: == domain 1 score: 33.7 bits; conditional E-value: 2.8e-11 AAA_18 1 IlltGtpGvGKttlareLskrlgf..........gvnlrdvakenglllevdpe.....aldeskrldedkldklldaleenss 69 I+ltG +++GK+t++++L+++ gf + ++ ++e+ + ld+++ +d++kl +++ + +e + 23726_FusoPortal_Gene_626 3 IGLTGGIASGKSTVSKYLAEK-GFkvydadkiakDISEKKSVQEEII------StfgnkILDKNRNVDKKKLKEIVFENKEKLK 79 89*******************.5655887666652222223333222......1333333555555555555555444444444 PP AAA_18 70 .................aleegenlvvdghlaelfserakdikkvvLtasp.eelleRlekRgyeeakike 122 + + + +++ +++ ++s k+ k+++ s+ e+ l R+ +R++ +++ +e 23726_FusoPortal_Gene_626 80 qlnaiihpkvihfykelKGKNTSEIIIFDVPLLFESGIDKFCDKILVIISDyEIQLNRIVERDKIDRDLAE 150 55559999******77655666667**********99999888855555555***********88877664 PP >> Cytidylate_kin2 Cytidylate kinase-like family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 17.2 0.2 2.6e-06 0.004 1 66 [. 2 73 .. 2 106 .. 0.75 2 ? 3.8 0.2 0.036 54 120 149 .. 127 156 .. 123 173 .. 0.86 Alignments for each domain: == domain 1 score: 17.2 bits; conditional E-value: 2.6e-06 Cytidylate_kin2 1 vITIsrqyGsgGreigekLAekLgipfydr.eiieevaeelgiseeev.....ekyde.ksrlksllskllaa 66 +I +++ +sg++++ + LAek g+++yd +i + ++e++ ++ee++ + +d+ +++ k+ l++++++ 23726_FusoPortal_Gene_626 2 IIGLTGGIASGKSTVSKYLAEK-GFKVYDAdKIAKDISEKKSVQEEIIstfgnKILDKnRNVDKKKLKEIVFE 73 7889999**************9.7******889999***********95444444444444444444444444 PP == domain 2 score: 3.8 bits; conditional E-value: 0.036 Cytidylate_kin2 120 ltapledrvkrvaereglseeearekiket 149 ++++ e ++r+ er +++++ a++ ik++ 23726_FusoPortal_Gene_626 127 IISDYEIQLNRIVERDKIDRDLAEKIIKSQ 156 788999999*********999999999875 PP >> AAA_28 AAA domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 16.3 0.0 5.6e-06 0.0084 3 43 .. 4 46 .. 2 108 .. 0.84 2 ? -3.2 0.0 5.3 8e+03 54 76 .. 147 169 .. 137 172 .. 0.65 Alignments for each domain: == domain 1 score: 16.3 bits; conditional E-value: 5.6e-06 AAA_28 3 vitGgestGKTtlleaLaargypvv..pEaareileeklakge 43 +tGg +GK+t+++ La++g++v+ +a +i e+k+ ++e 23726_FusoPortal_Gene_626 4 GLTGGIASGKSTVSKYLAEKGFKVYdaDKIAKDISEKKSVQEE 46 58*********************99878888888888876665 PP == domain 2 score: -3.2 bits; conditional E-value: 5.3 AAA_28 54 afaervleeqleeeaaaaankvl 76 ++ae++++ ql++e++ ++ +v+ 23726_FusoPortal_Gene_626 147 DLAEKIIKSQLSNEERIKKADVV 169 45666666666666666655555 PP >> SKI Shikimate kinase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 14.0 0.2 2.6e-05 0.039 2 37 .. 10 45 .. 9 85 .. 0.74 2 ? 0.2 0.1 0.44 6.7e+02 100 146 .. 136 160 .. 122 187 .. 0.47 Alignments for each domain: == domain 1 score: 14.0 bits; conditional E-value: 2.6e-05 TSSHHHHHHHHHHHHTSEEEEHHHHHHHHHT.S-HHH CS SKI 2 GaGKstvGkllAkklglefiDlDeliekrtg.ksiae 37 ++GKstv k lA+k g++++D+D++ ++ + ks++e 23726_FusoPortal_Gene_626 10 ASGKSTVSKYLAEK-GFKVYDADKIAKDISEkKSVQE 45 59**********98.9********9977666455543 PP == domain 2 score: 0.2 bits; conditional E-value: 0.44 HHHHHTTSSGGGCSCHHHHHHHHHHHHHHHHHCCSSCEEEECTSSHH CS SKI 100 erlkadkeRPllqekskeeelleelaereplyeeaadvvvdeseesk 146 +r+ + r+++ +++a ++++++ +++ 23726_FusoPortal_Gene_626 136 NRIVE----------------------RDKIDRDLAEKIIKSQLSNE 160 33333......................33333333333333333332 PP >> FokI_N Restriction endonuclease FokI, recognition domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.5 0.7 3.6e-05 0.055 16 72 .. 45 102 .. 33 105 .. 0.86 2 ? -1.1 0.1 1.2 1.9e+03 31 55 .. 147 157 .. 120 185 .. 0.51 Alignments for each domain: == domain 1 score: 13.5 bits; conditional E-value: 3.6e-05 FokI_N 16 kkvvaifvegs.kvhrelveeklpklvsekdGkkelleelearelsiayaelkGkgtp 72 +++++ f ++ +++r++ ++kl+++v e++ k ++l+++ + ++ y+elkGk t+ 23726_FusoPortal_Gene_626 45 EEIISTFGNKIlDKNRNVDKKKLKEIVFENKEKLKQLNAIIHPKVIHFYKELKGKNTS 102 4799999887527899999**************************99********986 PP == domain 2 score: -1.1 bits; conditional E-value: 1.2 FokI_N 31 elveeklpklvsekdGkkelleele 55 l+e+ ++ 23726_FusoPortal_Gene_626 147 DLAEKIIKSQ--------------L 157 2222222222..............2 PP >> AAA_33 AAA domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.8 0.1 0.00012 0.18 5 42 .. 6 42 .. 1 163 [. 0.82 Alignments for each domain: == domain 1 score: 11.8 bits; conditional E-value: 0.00012 AAA_33 5 vGlpGsGKstlaerllkklgavllssdtlrrrlfgtdq 42 G +sGKst+++ l+ + g+ + + d++ + + +++ 23726_FusoPortal_Gene_626 6 TGGIASGKSTVSKYLA-EKGFKVYDADKIAKDISEKKS 42 56678********988.889999999998888776433 PP >> AAA_17 AAA domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.0 2.6 0.00013 0.19 1 87 [. 6 94 .. 6 147 .. 0.65 2 ? 1.7 0.5 0.2 3e+02 46 94 .. 133 181 .. 127 188 .. 0.67 Alignments for each domain: == domain 1 score: 12.0 bits; conditional E-value: 0.00013 EE-TTSSHHHHHHHHHHHHTEEEEEHHHHHHHHHHTTTS..HSSHHHHCTSHHHHHHHHHHHHHHHH....HHHHHSSEEEEEE-S CS AAA_17 1 lGgPGsGKgtqceklvekygfthlsaGdllreeakrggl..akardemrklpleekeelqkeaaeki....aeeakkkgflvdghp 80 +Gg sGK+t+++ l+ek gf+++ d + + ++++++ ++ + + + l++++ ++k+ ++i +e+ k+ ++++ hp 23726_FusoPortal_Gene_626 6 TGGIASGKSTVSKYLAEK-GFKVYD-ADKIAKDISEKKSvqEEIISTFGNKILDKNRNVDKKKLKEIvfenKEKLKQLNAII--HP 87 69999*******666665.******.68888888888763344456666555556666665555555333144447777888..77 PP EEEETTE CS AAA_17 81 rikeqag 87 + + 23726_FusoPortal_Gene_626 88 KVIHFYK 94 6544333 PP == domain 2 score: 1.7 bits; conditional E-value: 0.2 HCTSHHHHHHHHHHHHHHHHHHHHHSSEEEEEE-SEEEETTEEEESSHH CS AAA_17 46 mrklpleekeelqkeaaekiaeeakkkgflvdghprikeqagfldglpk 94 ++ + ++e++++++ +aeki + + +++ +++ ++ e++ +l++l + 23726_FusoPortal_Gene_626 133 IQLNRIVERDKIDRDLAEKIIKSQLSNEERIKKADVVIENNSNLEDLFE 181 5566778889999999999998886655555555555555555665544 PP >> PRK Phosphoribulokinase / Uridine kinase family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 8.6 5.7 0.00086 1.3 1 158 [. 2 149 .. 2 180 .. 0.56 Alignments for each domain: == domain 1 score: 8.6 bits; conditional E-value: 0.00086 PRK 1 vIgvagsSGaGkttvarkleki.fgrekvkaaeiegdsflrldraekdleleerkrakekevsflspeandfdllyevl....k 79 +Ig++g ++Gk+tv++ l + f+ v + d+ ek+ ee ++ +++ ++ ++ d ++l e + 23726_FusoPortal_Gene_626 2 IIGLTGGIASGKSTVSKYLAEKgFK---VY-----DADKIAKDISEKKSVQEEIISTFGNKIL-DKNRNVDKKKLKEIVfenkE 76 8***************998654155...33.....2233344455555555555555544442.33333333333332211114 PP PRK 80 elkegkkvdkpiYsheteeaap..dkeeliegadvlviEGlhalydervaelldlkiyvdisidlelarklqRDiaeRGrs 158 +lk+ +++ p+ h ++e++ e +i++ +l+ G+ ++ d + + +++l+r ++RD+ +R 23726_FusoPortal_Gene_626 77 KLKQLNAIIHPKVIHFYKELKGknTSEIIIFDVPLLFESGI--------DKFCDKILVIISDYEIQLNRIVERDKIDRDLA 149 56667777777777777777662223444555555555555........55555444444555566666666666555444 PP >> DUF4284 Immunity protein 22 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.8 1.4 0.00044 0.67 26 92 .. 12 85 .. 9 102 .. 0.72 2 ? 1.0 0.2 0.48 7.3e+02 53 92 .. 129 170 .. 114 188 .. 0.63 Alignments for each domain: == domain 1 score: 10.8 bits; conditional E-value: 0.00044 DUF4284 26 gkslssqFckdlgiefyDeDfleikfleeesedieelleevsyd.....eei.ieklkekl..keeklnkaNavi 92 gks s+++ + g++ yD D + + e++s++ ee+++ + + ++ ++klke + ++ekl++ Na+i 23726_FusoPortal_Gene_626 12 GKSTVSKYLAEKGFKVYDADKIAKDISEKKSVQ-EEIISTFGNKildknRNVdKKKLKEIVfeNKEKLKQLNAII 85 88999*********************9998888.88887765431111112224445444322466777888887 PP == domain 2 score: 1.0 bits; conditional E-value: 0.48 DUF4284 53 eeesedieelleevsydeeiieklkek..lkeeklnkaNavi 92 + +++++ ++e +d ++ ek++++ +ee+++ka vi 23726_FusoPortal_Gene_626 129 SDYEIQLNRIVERDKIDRDLAEKIIKSqlSNEERIKKADVVI 170 455666677777777777776666543113455566666555 PP >> SapB_1 Saposin-like type B, region 1 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 1.7 0.2 0.17 2.6e+02 15 28 .. 36 49 .. 31 53 .. 0.75 2 ? -0.7 0.0 0.91 1.4e+03 29 37 .. 116 124 .. 114 124 .. 0.79 3 ? 2.0 0.1 0.14 2.1e+02 13 29 .. 149 166 .. 147 166 .. 0.90 4 ? 8.9 0.3 0.00091 1.4 10 36 .. 163 189 .. 163 190 .. 0.93 Alignments for each domain: == domain 1 score: 1.7 bits; conditional E-value: 0.17 HHH-SS--HHHHHH CS SapB_1 15 klLkdNkTeeeIie 28 + ++++ +eeIi+ 23726_FusoPortal_Gene_626 36 DISEKKSVQEEIIS 49 4556889*****97 PP == domain 2 score: -0.7 bits; conditional E-value: 0.91 HHHHHHHHS CS SapB_1 29 aLekvCskL 37 ++k C+k+ 23726_FusoPortal_Gene_626 116 GIDKFCDKI 124 589***985 PP == domain 3 score: 2.0 bits; conditional E-value: 0.14 HHHHH-SS.--HHHHHHH CS SapB_1 13 veklLkdN.kTeeeIiea 29 +ek +k++ ++ee I++a 23726_FusoPortal_Gene_626 149 AEKIIKSQlSNEERIKKA 166 8999***********986 PP == domain 4 score: 8.9 bits; conditional E-value: 0.00091 HHHHHHHH-SS--HHHHHHHHHHHHHH CS SapB_1 10 vkeveklLkdNkTeeeIieaLekvCsk 36 +k+++ +++N++ e++ e +e C++ 23726_FusoPortal_Gene_626 163 IKKADVVIENNSNLEDLFEKVERFCEM 189 5899999******************97 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (191 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 1535 (0.0918502); expected 334.2 (0.02) Passed bias filter: 516 (0.030876); expected 334.2 (0.02) Passed Vit filter: 66 (0.00394926); expected 16.7 (0.001) Passed Fwd filter: 11 (0.00065821); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 11 [number of targets reported over threshold] # CPU time: 0.26u 0.24s 00:00:00.50 Elapsed: 00:00:00.31 # Mc/sec: 1791.11 // [ok]