breseq  version 0.33.1  revision 8505477f25b3
mutation predictions | marginal predictions | summary statistics | genome diff | command line log

Marginal read alignment evidence (highest frequency 20 of 220 shown, sorted by frequency from high to low)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606720,4660GC55.0% ‑3.9 / 18.4 20A23P (GCC→CCC) phrBdeoxyribodipyrimidine photolyase, FAD‑binding
*REL606720,4640CG51.8% 5.1 / 41.8 27A22G (GCC→GGC) phrBdeoxyribodipyrimidine photolyase, FAD‑binding
*REL6061,118,2090GT44.7% ‑1.3 / 42.9 38intergenic (‑444/‑311)csgD/csgBDNA‑binding transcriptional activator in two‑component regulatory system/curlin nucleator protein, minor subunit in curli complex
*REL6064,392,5560CT44.2% 126.3 / 116.5 154L43L (CTA→TTA) aidBisovaleryl CoA dehydrogenase
*REL6062,735,6080GC42.3% 20.4 / 13.8 26A409P (GCC→CCC) ascFfused cellobiose/arbutin/salicin‑specific PTS enzymes: IIB component/IC component
*REL6062,376,3480TG41.7% 238.7 / 45.6 156T63P (ACC→CCC) folCbifunctional folylpolyglutamate synthase/ dihydrofolate synthase
*REL6062,515,9090AC40.5% 122.3 / 22.4 74V260G (GTG→GGG) ypfIpredicted hydrolase
*REL6064,180,2490TG40.2% 180.4 / 35.5 107G24G (GGT→GGGyjaHhypothetical protein
*REL6061,485,4500AG40.0% 296.5 / 167.8 210G193G (GGA→GGGydcVpredicted spermidine/putrescine transporter subunit
*REL6062,320,0380AC38.6% 344.3 / 41.6 163V5G (GTC→GGC) pmrDpolymyxin resistance protein B
*REL6061,854,4010AC38.5% 152.3 / 40.0 104K88N (AAA→AACyeaOhypothetical protein
*REL6064330TG37.1% 110.7 / 28.7 70V33G (GTG→GGG) thrAbifunctional aspartokinase I/homeserine dehydrogenase I
*REL6062,437,4060GT37.1% 66.7 / 16.2 35D435E (GAC→GAAypdDfused predicted PTS enzymes: Hpr component/enzyme I component/enzyme IIA component
*REL606414,0560TG36.2% 253.1 / 54.7 138T345P (ACC→CCC) yajRpredicted transporter
*REL6064,032,4500TG36.1% 184.2 / 51.6 108T363P (ACC→CCC) glnGfused DNA‑binding response regulator in two‑component regulatory system with GlnL: response regulator/sigma54 interaction protein
*REL6064,474,9710TG35.5% 373.4 / 29.4 166N261T (AAC→ACC) idnTL‑idonate and D‑gluconate transporter
*REL6063,774,9140CG35.4% 303.0 / 124.9 164A20P (GCC→CCC) yicJpredicted transporter
*REL606239,1630TG34.3% 387.5 / 111.5 198V58G (GTG→GGG) dnaQDNA polymerase III subunit epsilon
*REL606406,9580TG34.0% 335.9 / 56.6 156T548P (ACC→CCC) dxs1‑deoxy‑D‑xylulose‑5‑phosphate synthase
*REL606254,0430TG34.0% 84.9 / 17.5 53pseudogene (588/1710 nt)fhiAb0229; flagellar system protein

Marginal new junction evidence (sorted from low to high skew)
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2802968 =225 (1.020)3 (0.020) 3/132 20.3 1.6% coding (961/1302 nt) rumA 23S rRNA (uracil‑5‑)‑methyltransferase
?REL606 = 4211851 171 (0.890)coding (314/1350 nt) lysC aspartate kinase III
* ? REL606 = 203970791 (0.410)4 (0.020) 4/144 20.3 3.7% coding (422/438 nt) vioB VioB, involved in dTDP‑N‑acetylviosamine synthesis
?REL606 = 2040225 121 (0.570)intergenic (‑97/+46) vioB/vioA VioB, involved in dTDP‑N‑acetylviosamine synthesis/VioA, involved in dTDP‑N‑acetylviosamine synthesis
* ? REL606 = 2680757221 (1.010)3 (0.010) 3/138 20.8 1.8% coding (283/897 nt) ECB_02511 ParB family protein
?REL606 = 2681225 116 (0.570)coding (751/897 nt) ECB_02511 ParB family protein
* ? REL606 = 1556432247 (1.120)3 (0.010) 3/138 20.8 1.3% coding (44/1683 nt) ydeN hypothetical protein
?REL606 = 1556925 221 (1.100)coding (714/762 nt) ydeO predicted DNA‑binding transcriptional acfivator
* ? REL606 = 2040186126 (0.570)3 (0.010) 3/140 20.9 2.8% intergenic (‑58/+85) vioB/vioA VioB, involved in dTDP‑N‑acetylviosamine synthesis/VioA, involved in dTDP‑N‑acetylviosamine synthesis
?REL606 = 2040466 94 (0.460)coding (918/1113 nt) vioA VioA, involved in dTDP‑N‑acetylviosamine synthesis
* ? REL606 = 3737261180 (0.820)3 (0.010) 3/142 21.0 2.1% coding (566/1254 nt) waaL Lipid A‑core, surface polymer ligase
?REL606 = 3737695 112 (0.540)coding (132/1254 nt) waaL Lipid A‑core, surface polymer ligase
* ? REL606 2583007 =344 (1.570)73 (0.350) 3/144 21.2 18.2% intergenic (‑177/+275) iscR/yfhQ DNA‑binding transcriptional repressor/predicted methyltransferase
?REL606 = 2583279 327 (1.550)intergenic (‑449/+3) iscR/yfhQ DNA‑binding transcriptional repressor/predicted methyltransferase
* ? REL606 = 1530048271 (1.230)3 (0.010) 3/144 21.2 1.0% coding (171/216 nt) bdm biofilm‑dependent modulation protein
?REL606 3790384 = 311 (1.480)coding (2144/3498 nt) ECB_03531 hypothetical protein
* ? REL606 = 2589475199 (0.910)3 (0.010) 3/148 21.4 1.5% coding (121/891 nt) hcaR DNA‑binding transcriptional activator of 3‑phenylpropionic acid catabolism
?REL606 4277048 = 190 (0.880)coding (1141/2148 nt) fdhF formate dehydrogenase‑H, selenopolypeptide subunit
* ? REL606 = 2948453198 (0.900)3 (0.010) 3/150 21.6 1.3% coding (1853/2145 nt) sbm methylmalonyl‑CoA mutase
?REL606 3953252 = 262 (1.190)coding (270/2547 nt) cyaA adenylate cyclase