# 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_310.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_310 [L=455] Description: # 337412 # 338776 # -1 # ID=1_310;partial=00;start_type=ATG;rbs_motif=GGAG/GAGG;rbs_spacer=5-10bp;gc_cont=0.221 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 ------- ------ ----- ------- ------ ----- ---- -- -------- ----------- 2.3e-27 96.3 10.5 2.3e-27 96.3 10.5 2.3 3 ATPase_2 ATPase domain predominantly from Archaea 2.6e-26 91.9 3.3 2.6e-26 91.9 3.3 3.0 2 DUF234 Archaea bacterial proteins of unknown functi 3.2e-13 49.8 4.3 3.2e-13 49.8 4.3 2.6 2 AAA_14 AAA domain 3.1e-07 30.6 0.0 6.8e-07 29.5 0.0 1.5 1 AAA_22 AAA domain 3.3e-07 30.7 0.0 2.9e-06 27.7 0.0 2.1 1 AAA_16 AAA ATPase domain 0.0001 22.4 0.1 0.00036 20.7 0.1 1.9 1 Roc Ras of Complex, Roc, domain of DAPkinase 0.0002 21.2 0.0 0.00045 20.1 0.0 1.5 1 Mrr_cat Restriction endonuclease 0.0014 18.3 0.1 0.0035 17.0 0.1 1.6 1 AAA_24 AAA domain 0.0023 18.3 0.1 0.0061 16.9 0.1 1.7 1 AAA ATPase family associated with various cellul 0.0049 16.2 0.3 0.0095 15.2 0.3 1.4 1 Guanylate_cyc_2 Guanylylate cyclase 0.0064 16.4 1.1 0.0064 16.4 1.1 2.8 2 HTH_45 Winged helix-turn-helix ------ inclusion threshold ------ 0.011 14.8 0.4 0.033 13.3 0.4 1.8 1 TniB Bacterial TniB protein 0.019 15.1 0.1 0.045 13.9 0.1 1.6 1 UPF0102 Uncharacterised protein family UPF0102 0.12 12.2 1.5 0.16 11.8 0.2 1.9 2 HTH_5 Bacterial regulatory protein, arsR family Domain annotation for each model (and alignments): >> ATPase_2 ATPase domain predominantly from Archaea # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 96.3 10.5 1.9e-30 2.3e-27 1 232 [. 3 204 .. 3 205 .. 0.91 2 ? -3.4 0.1 5.7 6.8e+03 63 63 .. 330 330 .. 285 357 .. 0.53 3 ? 0.3 0.4 0.43 5.2e+02 25 61 .. 389 425 .. 376 453 .. 0.73 Alignments for each domain: == domain 1 score: 96.3 bits; conditional E-value: 1.9e-30 ATPase_2 1 FfdRekelkeleeilerepneiyfiYGpegsGKtaLlkeiienledkdkyviyyialreesiskldefieeeeekkklekalgi 84 F++Rekel+ l++ +++ n+++++YG++++GKt+L+ke+i+ + + +y++a ++ ++ ++++f k+++++++++ 23726_FusoPortal_Gene_310 3 FINREKELETLNKEYKK-DNSFVVLYGRRRVGKTTLIKEFIKGK-----KAFYFFADKQNENLQIERF------KNQISEYFKD 74 99*************99.8999*******************988.....668999999999999****......7777777755 PP ATPase_2 85 avpkaeleeilkeksklvfkaielllekikkkgkkiiiiiDelqklidlglngaellvkeLlnliesltkekhlcvviclsSeg 168 ++ k+ ++k+ + +++ll ki+++ kiii+iDe+q+l +n +++ + +++++++++k+k++++++c+s + 23726_FusoPortal_Gene_310 75 EFLKK-------IEIKDWDTVFDYLLTKISNE--KIIIVIDEFQYLC--MIN--KDFSSIFQRIYDEKLKDKNIMIILCGSLIS 145 55442.......2469999**********999..9************..999..999*************************** PP ATPase_2 169 lfie.rvyneselegrhky.ilvddldketalnflkegfkeleelskedfeeiyryvGGkPkalve 232 ++ + +es+l+gr++ i+++ ++++ +f+++ +s +++ e+y+++GG Pk++ 23726_FusoPortal_Gene_310 146 MMYSeTLAYESPLYGRRTAqIKLQAIKFKYYSKFFRD-------KSTQELIELYSITGGVPKYILS 204 *987699***********99***************99.......4889***************976 PP == domain 2 score: -3.4 bits; conditional E-value: 5.7 ATPase_2 63 s 63 + 23726_FusoPortal_Gene_310 330 T 330 2 PP == domain 3 score: 0.3 bits; conditional E-value: 0.43 ATPase_2 25 iYGpegsGKtaLlkeiienledkdkyviyyialrees 61 ++G + K ++ +i+++l++k+k+vi+ ++ ree+ 23726_FusoPortal_Gene_310 389 VFGECKYSKKQVGLNILNELKEKSKKVIWKNDKREEY 425 6677777777777889999999999999998888765 PP >> DUF234 Archaea bacterial proteins of unknown function # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? -0.1 0.5 1 1.2e+03 16 36 .. 56 76 .. 41 107 .. 0.62 2 ! 91.9 3.3 2.2e-29 2.6e-26 1 92 [. 313 398 .. 313 416 .. 0.94 Alignments for each domain: == domain 1 score: -0.1 bits; conditional E-value: 1 DUF234 16 eseevlekikeelneyvgkvF 36 +++ +e++k++++ey+ +F 23726_FusoPortal_Gene_310 56 NENLQIERFKNQISEYFKDEF 76 233333666666666665555 PP == domain 2 score: 91.9 bits; conditional E-value: 2.2e-29 DUF234 1 wfrfvypnrseielgeseevlekikeelneyvgkvFEelvrelllelnkagklpfrftkvGrwWkkeeeiDivaldeeekkalf 84 wf++vyp++s++e+++ +++ +ki++e++ yv+k +E+l+re ++ ++++f++ kvGrwW+k++eiDival e++k ++f 23726_FusoPortal_Gene_310 313 WFSYVYPYQSYLEIENLTYIKNKIENEFDLYVSKIYEDLARESIW-----ENMKFPLVKVGRWWDKDTEIDIVALGEDNK-IVF 390 *********************************************.....568899***********************6.*** PP DUF234 85 gEvKwkns 92 gE+K++++ 23726_FusoPortal_Gene_310 391 GECKYSKK 398 *****994 PP >> AAA_14 AAA domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 49.8 4.3 2.7e-16 3.2e-13 3 127 .. 22 173 .. 20 177 .. 0.66 2 ? -3.2 0.0 6.3 7.5e+03 61 83 .. 227 248 .. 201 260 .. 0.49 Alignments for each domain: == domain 1 score: 49.8 bits; conditional E-value: 2.7e-16 AAA_14 3 kilvitGaRrvGKTtllkqlikellvkkesilyinldderd....................kkladpe.....kllelileklk 61 +v+ G RrvGKTtl+k++ik + +y+ d++++ kk+ e ++ +++l k++ 23726_FusoPortal_Gene_310 22 SFVVLYGRRRVGKTTLIKEFIKGKK-----AFYFFADKQNEnlqierfknqiseyfkdeflKKI---EikdwdTVFDYLLTKIS 97 679******************8887.....3333333322222333322233333333332222...02223355566666555 PP AAA_14 62 eekktylflDEiqev....pewekalkslvdeg...knldvivtGSnskllskeia...tslagrva.eielyplsf 127 ++k ++++DE q++ ++++++ ++++de+ kn+ +i+ GS+ +++++e+ + l gr + +i+l ++f 23726_FusoPortal_Gene_310 98 -NEKIIIVIDEFQYLcminKDFSSIFQRIYDEKlkdKNIMIILCGSLISMMYSETLayeSPLYGRRTaQIKLQAIKF 173 .59************9***************6666699*********999988764344555555443677766666 PP == domain 2 score: -3.2 bits; conditional E-value: 6.3 AAA_14 61 keekktylflDEiqevpewekal 83 + ++ ++l+ E+ +++ +++l 23726_FusoPortal_Gene_310 227 Y-SEPKFLLQEEVNDLSRYFSIL 248 1.344555555555555555555 PP >> AAA_22 AAA domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 29.5 0.0 5.7e-10 6.8e-07 3 129 .. 19 141 .. 15 144 .. 0.80 Alignments for each domain: == domain 1 score: 29.5 bits; conditional E-value: 5.7e-10 AAA_22 3 ktagilvvtGesGtGKTtllrrlleqldeeakaeeerpevvfvtlpsst.sskdllrellralglelka..erlsseelleale 83 k+ +++v++G +GKTtl++++ + + ++ ++f + ++++ + + +++ + +++e + e ++ + ++++l 23726_FusoPortal_Gene_310 19 KDNSFVVLYGRRRVGKTTLIKEFIKGK---------KAFYFFADKQNENlQIERFKNQISEYFKDEFLKkiEIKDWDTVFDYLL 93 57899*********************9.........777777777777777777777777776666555566679999999999 PP AAA_22 84 qalraravielliiDEaqhlr..seealeeLrsllnl..srkgvqliLvG 129 +++ + +++iDE+q+l +++ + + ++++ ++k++ +iL+G 23726_FusoPortal_Gene_310 94 TKISNEKI--IIVIDEFQYLCmiNKDFSSIFQRIYDEklKDKNIMIILCG 141 99999988..*********9764455555558888889999******999 PP >> AAA_16 AAA ATPase domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 27.7 0.0 2.5e-09 2.9e-06 2 99 .. 3 87 .. 2 129 .. 0.66 Alignments for each domain: == domain 1 score: 27.7 bits; conditional E-value: 2.5e-09 AAA_16 2 lvGReeeleqLlsalervrsgrpelvllsGesGvGKsslvrellralekdggiflrgkcdqlesaipy....salaqalareellr 83 ++ Re+ele+L+ +++ + + v+l+G + vGK++l++e+ + k+ ++f++ k ++ + ++ +++ ++e+l+ 23726_FusoPortal_Gene_310 3 FINREKELETLNKEYKKDN----SFVVLYGRRRVGKTTLIKEFIK--GKKAFYFFADKQNEN---LQIerfkNQISEYF-KDEFLK 78 788************9665....79*******************7..555599988888888...33211112233333.223333 PP AAA_16 84 qlldeaesaeleawre 99 +++ +w++ 23726_FusoPortal_Gene_310 79 KIE-------IKDWDT 87 332.......333333 PP >> Roc Ras of Complex, Roc, domain of DAPkinase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 20.7 0.1 3e-07 0.00036 2 63 .. 24 87 .. 24 108 .. 0.86 Alignments for each domain: == domain 1 score: 20.7 bits; conditional E-value: 3e-07 EEEESSTTSSHHHHHHCHCT-C.....CC-EEETTTTEEEEEEECGGSGGCTTBTTTTCEEEEEEEE CS Roc 2 ivlvGdsgvGKtsllkrfvedk.....fsekekstigvdfltkevtkikglendkkvkkvklniWDt 63 +vl G +vGKt+l+k+f+++k f +k+++ ++++ ++++++ ++ +d+ +kk++++ WDt 23726_FusoPortal_Gene_310 24 VVLYGRRRVGKTTLIKEFIKGKkafyfFADKQNENLQIERFKNQIS---EYFKDEFLKKIEIKDWDT 87 7899****************8744444579999999**99998875...455678899********9 PP >> Mrr_cat Restriction endonuclease # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 20.1 0.0 3.8e-07 0.00045 16 105 .. 356 450 .. 345 453 .. 0.82 Alignments for each domain: == domain 1 score: 20.1 bits; conditional E-value: 3.8e-07 HHHHHSSSEEE---B..-TTSSB--EEEETTEEEEEEEEEEE--TS-B-HHHHHHHHHHCTT......S..SEEEEEESS-B-H CS Mrr_cat 16 llealGykeveltgr..sgDggiDliakkdplgserilvqaKryknkvgvkeVqelvgalkk......rgadkgifvttsdftk 91 ++e++ + + +gr +D iD++a +++ + ++ ++K +k++vg + ++el kk ++ + i+ + s+f++ 23726_FusoPortal_Gene_310 356 IWENMKFP-LVKVGRwwDKDTEIDIVALGEDN--KIVFGECKYSKKQVGLNILNELKEKSKKviwkndKREEYYIIFSKSGFSQ 436 67777777.5666663346778****876553..599*******************99999889988777999999999***** PP HHHHHHHTT--EEE CS Mrr_cat 92 aAreeaeelnkieL 105 + +e+a++ n+i+L 23726_FusoPortal_Gene_310 437 DLIELAKKKNNIIL 450 ********999998 PP >> AAA_24 AAA domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 17.0 0.1 2.9e-06 0.0035 5 74 .. 24 108 .. 21 124 .. 0.72 Alignments for each domain: == domain 1 score: 17.0 bits; conditional E-value: 2.9e-06 AAA_24 5 iliyGpsgiGKTtlakel..gpkplfidtekGlksldadkfpdie............swedfldevvel.aaeladyktlviDt 73 +++yG +GKTtl+ke+ g+k++++ +k +l+ +f++ ++d+ +++l + ++ ++viD 23726_FusoPortal_Gene_310 24 VVLYGRRRVGKTTLIKEFikGKKAFYFFADKQNENLQIERFKNQIseyfkdeflkkiEIKDWDTVFDYLlTKISNEKIIIVIDE 107 79******************9999999999988889988887544566555555555444444444444555555555666665 PP AAA_24 74 i 74 23726_FusoPortal_Gene_310 108 F 108 5 PP >> AAA ATPase family associated with various cellular activities (AAA) # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 16.9 0.1 5.1e-06 0.0061 1 115 [. 24 147 .. 24 156 .. 0.59 Alignments for each domain: == domain 1 score: 16.9 bits; conditional E-value: 5.1e-06 EEEESSTTSSHHHHHHHHHHH..............HTHEEEEEETTTCCSSSTTHHHHHHHHHHHHHHHTS............. CS AAA 1 vllyGppGtGKTllakavake..............lgvefleisgsellskyvgesekkirelfkeakeka............. 57 v+lyG +GKT+l+k +k l ++ +++i+e+fk+ k+ 23726_FusoPortal_Gene_310 24 VVLYGRRRVGKTTLIKEFIKGkkafyffadkqnenL---------------QIERFKNQISEYFKDEFLKKieikdwdtvfdyl 92 689*****************9433333322222222...............344444444444444211112233344444444 PP ......TSEEEEEETGGGHCHHHSSCSTHHHHHHHHHHHHHHHHTTTSSSTEEEEEEESSGGGS CS AAA 58 ......kpsilfiDEidalaksrsgseseeeervvnqLlteldgvkkkeskvivigatnrpdkl 115 ++ i++iDE+ l+ + ++ + + + d+ k ++++++i+ + ++++ 23726_FusoPortal_Gene_310 93 ltkisnEKIIIVIDEFQYLCMIN--------KDFSSIFQRIYDEKLK-DKNIMIILCGSLISMM 147 44445667************999........4555555555555555.5568888777666555 PP >> Guanylate_cyc_2 Guanylylate cyclase # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 15.2 0.3 7.9e-06 0.0095 88 196 .. 63 171 .. 41 177 .. 0.84 Alignments for each domain: == domain 1 score: 15.2 bits; conditional E-value: 7.9e-06 Guanylate_cyc_2 88 kdekRvnklfkkakeagievekrsvsveelvehlasgkvaivLvnaalLkcdlcksnkllselrkllkrkseytGhyvvlcGy. 170 + ++ +++ fk+ ++ie+++ ++ l++ ++++k++iv+ + ++L c + k ++s++++++++k + + ++lcG 23726_FusoPortal_Gene_310 63 RFKNQISEYFKDEFLKKIEIKDWDTVFDYLLTKISNEKIIIVIDEFQYL-CMINK--DFSSIFQRIYDEKLKDKNIMIILCGSl 143 4456889999*******************************99999999.65554..588889*******************62 PP Guanylate_cyc_2 171 ..dakekkllyrnPassdkvcsisveal 196 + +++l y++P + +++ +i+++a+ 23726_FusoPortal_Gene_310 144 isMMYSETLAYESPLYGRRTAQIKLQAI 171 22455788999***99999999988765 PP >> HTH_45 Winged helix-turn-helix # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 16.4 1.1 5.4e-06 0.0064 9 53 .. 245 289 .. 241 311 .. 0.87 2 ? -0.7 0.1 1.1 1.3e+03 49 72 .. 388 411 .. 379 415 .. 0.68 Alignments for each domain: == domain 1 score: 16.4 bits; conditional E-value: 5.4e-06 HTH_45 9 adILeaisegeakkTkimykanLsykllnkYLdmLvekglikksk 53 ++IL+ais g++k+++i ++ l++Y+++L++ ++++k++ 23726_FusoPortal_Gene_310 245 FSILNAISIGHTKMSAISSYLQINAGGLSSYISKLIDLDILEKES 289 79**************************************97654 PP == domain 2 score: -0.7 bits; conditional E-value: 1.1 HTH_45 49 ikkskekYvlTekGlelleklkel 72 i + kY+ + Gl++l++lke 23726_FusoPortal_Gene_310 388 IVFGECKYSKKQVGLNILNELKEK 411 556666777777888888888764 PP >> TniB Bacterial TniB protein # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.3 0.4 2.8e-05 0.033 22 130 .. 9 110 .. 4 141 .. 0.67 Alignments for each domain: == domain 1 score: 13.3 bits; conditional E-value: 2.8e-05 TniB 22 rLeeLleapkrkrmpnllivGdsnnGKtmivekFarlhppekd.edaekvpvlvvqappkpdekrlykaiLealgaplkararv 104 +Le+L + k++ + ++++G+ +GKt+++++F + +++ +d+++ ++ + +r+ ++i e ++ ++ ++ ++ 23726_FusoPortal_Gene_310 9 ELETLNKEYKKDN-SFVVLYGRRRVGKTTLIKEFIKGKKAFYFfADKQN---------ENLQIERFKNQISEYFKDEFLKKIEI 82 6777766666666.8899*****************99987643133332.........33444555555555555555555555 PP TniB 105 skleqqvlellrelkvr..vLiiDEihn 130 ++ + ll +++ + +++iDE++ 23726_FusoPortal_Gene_310 83 KDWDTVFDYLLTKISNEkiIIVIDEFQY 110 5555555555555543322556666665 PP >> UPF0102 Uncharacterised protein family UPF0102 # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 13.9 0.1 3.7e-05 0.045 21 54 .. 370 406 .. 349 426 .. 0.74 Alignments for each domain: == domain 1 score: 13.9 bits; conditional E-value: 3.7e-05 EEETTEEEEEEE...TTEEEEEEEEE-.......---- CS UPF0102 21 wrcragEiDiia...tevlvfvEVKtRssekfg.spae 54 w + EiDi+a ++ +vf E K+ ++++g + ++ 23726_FusoPortal_Gene_310 370 WWDKDTEIDIVAlgeDNKIVFGECKYS-KKQVGlNILN 406 334568*********889********5.4444424444 PP >> HTH_5 Bacterial regulatory protein, arsR family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.8 0.2 0.00013 0.16 7 47 .] 247 287 .. 245 287 .. 0.96 2 ? -2.3 0.1 3.3 3.9e+03 5 18 .. 402 415 .. 402 415 .. 0.83 Alignments for each domain: == domain 1 score: 11.8 bits; conditional E-value: 0.00013 HHHHHTTSSEEHHHHHHHHTS-HHHHHHHHHHHHHTTSEEE CS HTH_5 7 ILklLaegelcvcelaeelglsqstvShHLkkLreaGLVek 47 IL++ + g ++ + ++ +l+++ + +S + +kL +++++ek 23726_FusoPortal_Gene_310 247 ILNAISIGHTKMSAISSYLQINAGGLSSYISKLIDLDILEK 287 9************************************9986 PP == domain 2 score: -2.3 bits; conditional E-value: 3.3 HHHHHHHTTSSEEH CS HTH_5 5 lkILklLaegelcv 18 l IL+ L e++++v 23726_FusoPortal_Gene_310 402 LNILNELKEKSKKV 415 68999999988876 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (455 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 942 (0.0563667); expected 334.2 (0.02) Passed bias filter: 339 (0.0202848); expected 334.2 (0.02) Passed Vit filter: 47 (0.00281235); expected 16.7 (0.001) Passed Fwd filter: 14 (0.000837721); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 14 [number of targets reported over threshold] # CPU time: 0.37u 0.28s 00:00:00.65 Elapsed: 00:00:00.38 # Mc/sec: 3480.79 // [ok]