# 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_1287.faa # target HMM database: Pfam-A.hmm # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Query: 23726_FusoPortal_Gene_1287 [L=442] Description: # 1404471 # 1405796 # -1 # ID=1_1287;partial=00;start_type=GTG;rbs_motif=GGAGG;rbs_spacer=5-10bp;gc_cont=0.205 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.7e-17 63.5 18.7 4.9e-16 59.3 18.7 2.1 1 AAA_21 AAA domain, putative AbiEii toxin, Type IV TA syste 5.8e-12 45.8 53.5 8.7e-12 45.3 53.5 1.3 1 AAA_15 AAA ATPase domain ------ inclusion threshold ------ 0.014 14.7 16.5 1.3 8.3 16.5 2.6 1 SMC_N RecF/RecN/SMC N terminal domain 0.04 14.1 0.0 12 6.0 0.0 2.4 2 AAA_22 AAA domain 0.081 12.6 0.0 0.22 11.2 0.0 1.7 1 AAA_29 P-loop containing region of AAA domain 0.13 12.5 0.0 0.3 11.4 0.0 1.6 2 AAA_16 AAA ATPase domain 0.17 11.4 1.1 0.46 10.1 0.1 2.3 2 Cecropin Cecropin family 0.71 10.3 5.8 0.15 12.5 0.4 2.6 3 ABC_tran ABC transporter Domain annotation for each model (and alignments): >> AAA_21 AAA domain, putative AbiEii toxin, Type IV TA system # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 59.3 18.7 2.3e-19 4.9e-16 3 303 .] 32 377 .. 30 377 .. 0.79 Alignments for each domain: == domain 1 score: 59.3 bits; conditional E-value: 2.3e-19 AAA_21 3 vliGpNgsGKSnLLdaldllarl..............................drdvegltdelsggggirkipslldlydpe 55 ++G Ng GKSn+ d ++ l+r+ + + g+ ++++ +++ ++i+ ++ +++ + 23726_FusoPortal_Gene_1287 32 SIYGENGIGKSNFVDSFYTLKRIistrtinekvriltekqkdlqsndfdkalyFFGQLGSIIKNGFFSDSTDIITESKTINSK 114 59***************************************977666666666666666666666666666666666666666 PP AAA_21 56 ipfefeieeflsqgvryryglrlersqveeklssvetllesgllldkrsleeedvfnpeargirggvsvfslkelsvsepndl 138 +++e+ f ++++ y+++ ++ ++ + +++++l ++++ + e + + ++ + +++++ ke++ +++ + 23726_FusoPortal_Gene_1287 115 DNMVIEVG-FKIKSKSGIYRIETDDTDI-------ISEKLDFTLNKNKVNFFEITKKEKNLNESVFIDNEYKKEILSIIEKYW 189 66666666.6666666666666666666.......788888888888888888888888888888888888888888888888 PP AAA_21 139 qfgkaltditsrawsllllldglrkerigllnifadyfgfakekrllqrlvrslklveielsdlgegeeliasvderlrer.. 219 + ++l +i ++++ +++ + + ++ if +f ++ l + + +k++e e + g +++ ++ ++e+ 23726_FusoPortal_Gene_1287 190 GKHSLLSLIAYEIEDKKENYVK----KRIFNGIFEVINFFLSLSILSRNKMEVFKDIEKEKLFYGTL-SINEEKKLTNIENvi 267 8888888888888875555555....77777777777777777777777777777777777776533.333444444444489 PP AAA_21 220 .............................gvvflengdgqelparelSdGtkrllalltallralkkgsllliDEpEngLHPk 273 + + +++ ++ ++ +S Gtk +l++l++l++a+k ++++iDE++ng+H 23726_FusoPortal_Gene_1287 268 ntffialfsdikqayykkkidndkinyilCFKKNIYNKLIDIEYNRESTGTKNILKILPYLISAAKG-KTVIIDEIDNGIHDL 349 999**************************************************************98.*************97 PP AAA_21 274 llrrllellkelsrngaQviltTHSpllld 303 l+ ++le l+e ++Q+i+tTH++lll+ 23726_FusoPortal_Gene_1287 350 LMLKILENLSED--LKGQLIITTHNTLLLE 377 777777555444..559**********985 PP >> AAA_15 AAA ATPase domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ! 45.3 53.5 4.2e-15 8.7e-12 5 382 .. 4 373 .. 1 376 [. 0.72 Alignments for each domain: == domain 1 score: 45.3 bits; conditional E-value: 4.2e-15 AAA_15 5 kikikNFRslkdteidldee......ltvliGkNnsGKSsilealdlflnnkekknenkvikiekedlekfeniseskeikId 81 ik+kN++sl + e+dl ++ l + G+N++GKS++++ ++ ++ + +i+ + l ++ ++ s++++ + 23726_FusoPortal_Gene_1287 4 YIKLKNYKSLIELEVDLTKKentpkkLISIYGENGIGKSNFVDSFYTLKRII----STRTINEKVRILTEKQKDLQSNDFDKA 82 69************************************************88....344555555556666666666666666 PP AAA_15 82 vefliefeyledkkkilklkskkkskeyestlkskgkkkadlllkkkekdlkklikdkklklkkneysdekisvekkkekkel 164 ++f +i++ ++ s+ + +sk ++++d+++++ +++k + + ++++ +++++ek + +k+k + 23726_FusoPortal_Gene_1287 83 LYF-----FGQLGSIIKNGFF---SDSTDIITESKTINSKDNMVIEVGFKIKSKSGIYRIETDDTDIISEKLDFTLNKNKVNF 157 666.....1111333333333...3333345567788999999999999999999999999999******************* PP AAA_15 165 ekkkeeklekekaeikidkedlkklikkleeslnkeelknvlksaiiklllisadRnvsd.....edsekkssllgkllekil 242 ++ +++++++++ + id+e k++++ +e++ k++l +++ +i++ + + ++ + + + + s+l++ + + + 23726_FusoPortal_Gene_1287 158 FEITKKEKNLNESVF-IDNEYKKEILSIIEKYWGKHSLLSLIAYEIEDKKENYVKKRIFNgifevINFFLSLSILSRNKMEVF 239 *99999998888888.9****************99994444444444444443333333312222112223344455555555 PP AAA_15 243 eealqekl.....eeeleeklkkieeklkeidkelaeelk......kilekkidkikkfgykgkkneeieik..........d 304 ++ +ekl ++++e+kl +ie+ ++++ l ++ k ki ++ki+ 23726_FusoPortal_Gene_1287 240 KDIEKEKLfygtlSINEEKKLTNIENVINTFFIALFSDIKqayykkKIDNDKIN-----------------YilcfkkniynK 305 555455557777777777777777777777777776666665433333333333.................133334444446 PP AAA_15 305 geelplnynGsGvqrlillillileaekkdeneekinilaiEEPEahlHPqlqkvlikllkelannnktQviitTHSp 382 ++++n+ ++G+++++ ++ +++a k+ +++i+E+++ H + k l +l+++ k+Q+iitTH+ 23726_FusoPortal_Gene_1287 306 LIDIEYNRESTGTKNILKILPYLISAAKG-------KTVIIDEIDNGIHD---LLMLKILENLSEDLKGQLIITTHNT 373 679**************************.......*************5...679********************97 PP >> SMC_N RecF/RecN/SMC N terminal domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 8.3 16.5 0.00062 1.3 4 205 .. 4 379 .. 1 396 [. 0.64 Alignments for each domain: == domain 1 score: 8.3 bits; conditional E-value: 0.00062 EEEEESCSSEECEEEE....E-TTSSEEEEEESTTSSHHHHHHHHHHHHTT.TSSS--HHHHCCCT--.SSSSSTTCCEEEEC CS xxxxxxxxxxxxxxxx....xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF SMC_N 4 rleienFksyagktvi....lpfsksftaivGpNGSGKsnlldAilFvLGersakklrseklsdlihkskekasvkkaeveit 82 ++++n+ks+ + +v+ + k++ +i+G+NG GKsn++d+ ++L ++ +ek++ l k+k+ +s ++ + + 23726_FusoPortal_Gene_1287 4 YIKLKNYKSLIELEVDltkkENTPKKLISIYGENGIGKSNFVDSF-YTLKRIISTRTINEKVRILTEKQKDLQSNDFDKALYF 85 5889999999865554111156778899****************9.9998665555555666666665444444444444444 PP E...................EECCEE-B----BEEEEEEEETTT-.EEEEE..TTEEE-HHHHHHHHHCTT--TT........ CS x...................xxxxxxxxxxxxxxxxxxxxxxxxx.xxxxx..xxxxxxxxxxxxxxxxxxxxxx........ RF SMC_N 83 f...................dnedkklkiekeevsisrrvkrkge.seyki..ngkevtkkevsellesagiske........ 135 f ++++ + +k++++i++ k k + y+i +++++ ++++ l++ +++ 23726_FusoPortal_Gene_1287 86 FgqlgsiikngffsdstdiiTESK--TINSKDNMVIEVGFKIKSKsGIYRIetDDTDIISEKLDFTLNKNKVNFFeitkkekn 166 466665555554443332211111..222223333333222222202333322333333344444444444333367777888 PP ................................................................................... CS ................................................................................... RF SMC_N 136 ................................................................................... 135 23726_FusoPortal_Gene_1287 167 lnesvfidneykkeilsiiekywgkhsllsliayeiedkkenyvkkrifngifevinfflslsilsrnkmevfkdiekeklfy 249 88888888888888888888888888888888888888888888888888888888888888888888888888888888888 PP ...............................................................CS-HHHHHHHHHHHHHHHHT CS ...............................................................xxxxxxxxxxxxxxxxxxxx RF SMC_N 136 ...............................................................lLSgGektLvalaLlfAiqk 155 S G k ++ + i + 23726_FusoPortal_Gene_1287 250 gtlsineekkltnienvintffialfsdikqayykkkidndkinyilcfkkniynklidieynRESTGTKNILKILPY-LISA 331 888888888888888888888888888888888888888888888888888888888888887778888877665322.2333 PP CS--SEEEECSTTTTS-HHHHHHHHHHHHHHCCTSEEEEE-SSHHHHTT- CS xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF SMC_N 156 vkpaplyllDeidaaLDeknvkkvaellkeksknaQfivislreemleka 205 +k + +++Deid + + + k+ e l e+ + Q+i++++++ +le+ 23726_FusoPortal_Gene_1287 332 AKGK-TVIIDEIDNGIHDLLMLKILENLSEDL-KGQLIITTHNTLLLEEE 379 3333.357788888887777777777777766.67888888887777665 PP >> AAA_22 AAA domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 6.0 0.0 0.0059 12 5 25 .. 28 48 .. 23 84 .. 0.84 2 ? 5.7 0.0 0.0074 15 77 117 .. 320 360 .. 269 374 .. 0.69 Alignments for each domain: == domain 1 score: 6.0 bits; conditional E-value: 0.0059 AAA_22 5 agilvvtGesGtGKTtllrrl 25 ++++ ++Ge G+GK+ + + 23726_FusoPortal_Gene_1287 28 KKLISIYGENGIGKSNFVDSF 48 7899***********777665 PP == domain 2 score: 5.7 bits; conditional E-value: 0.0074 AAA_22 77 ellealeqalraravielliiDEaqh.lrs.eealeeLrslln 117 +l+ l ++a++ + +iiDE+++ + + +l+ L +l++ 23726_FusoPortal_Gene_1287 320 NILKILPYLISAAKG-KTVIIDEIDNgI-HdLLMLKILENLSE 360 444555555556666.899*******88.55444444444444 PP >> AAA_29 P-loop containing region of AAA domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.2 0.0 0.00011 0.22 24 45 .. 30 51 .. 4 53 .. 0.85 Alignments for each domain: == domain 1 score: 11.2 bits; conditional E-value: 0.00011 AAA_29 24 vtlltGpsGsGKSTllDAiqal 45 + + G++G GKS ++D++ +l 23726_FusoPortal_Gene_1287 30 LISIYGENGIGKSNFVDSFYTL 51 56799************99886 PP >> AAA_16 AAA ATPase domain # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 11.4 0.0 0.00014 0.3 22 49 .. 25 52 .. 16 112 .. 0.72 2 ? -4.0 0.0 7.3 1.5e+04 36 45 .. 334 343 .. 331 348 .. 0.80 Alignments for each domain: == domain 1 score: 11.4 bits; conditional E-value: 0.00014 AAA_16 22 grp.elvllsGesGvGKsslvrellrale 49 ++p l+ ++Ge G+GKs v+++++ l+ 23726_FusoPortal_Gene_1287 25 NTPkKLISIYGENGIGKSNFVDSFYT-LK 52 4556899**************99995.33 PP == domain 2 score: -4.0 bits; conditional E-value: 7.3 AAA_16 36 GKsslvrell 45 GK+ +++e+ 23726_FusoPortal_Gene_1287 334 GKTVIIDEID 343 9999999986 PP >> Cecropin Cecropin family # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 10.1 0.1 0.00022 0.46 7 17 .. 207 218 .. 204 222 .. 0.85 2 ? -2.8 0.1 2.4 5e+03 3 8 .. 239 244 .. 238 244 .. 0.89 Alignments for each domain: == domain 1 score: 10.1 bits; conditional E-value: 0.00022 Cecropin 7 En.aGKrireAA 17 En ++Kri++++ 23726_FusoPortal_Gene_1287 207 ENyVKKRIFNGI 218 9**********4 PP == domain 2 score: -2.8 bits; conditional E-value: 2.4 Cecropin 3 gKKLEn 8 +K +E+ 23726_FusoPortal_Gene_1287 239 FKDIEK 244 9****8 PP >> ABC_tran ABC transporter # score bias c-Evalue i-Evalue hmmfrom hmm to alifrom ali to envfrom env to acc --- ------ ----- --------- --------- ------- ------- ------- ------- ------- ------- ---- 1 ? 12.5 0.4 7e-05 0.15 10 44 .. 27 66 .. 22 186 .. 0.78 2 ? -0.2 0.1 0.6 1.3e+03 97 126 .. 302 331 .. 188 344 .. 0.44 3 ? -2.1 0.1 2.3 4.8e+03 28 70 .. 352 394 .. 350 409 .. 0.73 Alignments for each domain: == domain 1 score: 12.5 bits; conditional E-value: 7e-05 TTSEEEEEESTTSSHHHHHHHHTTSS--.....SEEEEEE CS xxxxxxxxxxxxxxxxxxxxxxxxxxxx.....xxxxxxx RF ABC_tran 10 egekvaivGenGaGKStLlkllagllkp.....teGeill 44 + +++ i GenG GKS ++ ++ l + ++++++ 23726_FusoPortal_Gene_1287 27 PKKLISIYGENGIGKSNFVDSFYTLKRIistrtINEKVRI 66 568899**********999988888776333333333333 PP == domain 2 score: -0.2 bits; conditional E-value: 0.6 TTCBGESTTSGGS-HHHHHHHHHHHHHHHT CS xxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF ABC_tran 97 lkdtvvksspssLSgGqkqrvalarallkk 126 + ++ ++ + + S G k +++ l+++ 23726_FusoPortal_Gene_1287 302 IYNKLIDIEYNRESTGTKNILKILPYLISA 331 334444445566666666666666666554 PP == domain 3 score: -2.1 bits; conditional E-value: 2.3 HHHHTTSS--SEEEEEETTEEGGGB-HHHCCHCEEEE-SS--S CS xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF ABC_tran 28 lkllagllkpteGeilldgkdlkeqeleslrkeigvlpqepql 70 lk+l +l+++ +G++ + +++ e+e+ + +i +++ +++ 23726_FusoPortal_Gene_1287 352 LKILENLSEDLKGQLIITTHNTLLLEEEFIKDSIYIFKVDENA 394 7889999999999998887777776777777777666555544 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (442 residues searched) Target model(s): 16712 (2907032 nodes) Passed MSV filter: 884 (0.0528961); expected 334.2 (0.02) Passed bias filter: 294 (0.0175921); expected 334.2 (0.02) Passed Vit filter: 22 (0.00131642); expected 16.7 (0.001) Passed Fwd filter: 8 (0.000478698); expected 0.2 (1e-05) Initial search space (Z): 16712 [actual number of targets] Domain search space (domZ): 8 [number of targets reported over threshold] # CPU time: 0.31u 0.24s 00:00:00.55 Elapsed: 00:00:00.33 # Mc/sec: 3893.66 // [ok]