breseq  version 0.39.0  
mutation predictions | marginal predictions | summary statistics | genome diff | command line log

Marginal read alignment evidence
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,396,7720TC76.9% 66.3 / 14.7 40F236F (TTT→TTCamiBN‑acetylmuramoyl‑L‑alanine amidase B
*NC_0009134,392,9970GA65.7% 35.6 / 26.8 35V357V (GTC→GTTqueGepoxyqueuosine reductase
*NC_0009134,393,5520GA64.3% 23.4 / 24.8 28F172F (TTC→TTTqueGepoxyqueuosine reductase
*NC_000913734,6820AG52.0% 2.4 / 24.5 29intergenic (+580/‑968)ybfB/ybfCuncharacterized protein YbfB/uncharacterized protein YbfC
*NC_000913266,3580TG40.9% 14.0 / 10.8 23T66P (ACC→CCC) yafYinner membrane lipoprotein YafY
*NC_000913734,6360CT38.5% 11.4 / 18.9 27intergenic (+534/‑1014)ybfB/ybfCuncharacterized protein YbfB/uncharacterized protein YbfC
*NC_0009132,972,6060TC36.7% 20.6 / 26.1 31V337A (GTG→GCG) tasputative NADP(H)‑dependent aldo‑keto reductase Tas
*NC_0009131,111,9680CT30.3% 34.0 / 22.6 33A369V (GCC→GTC) opgHosmoregulated periplasmic glucans biosynthesis protein H
*NC_000913732,7650CT28.2% 46.2 / 20.8 39T1061T (ACC→ACTrhsCrhs element protein RhsC
*NC_0009134,405,5520AG27.6% 41.7 / 14.8 29E289G (GAG→GGG) purAadenylosuccinate synthetase
*NC_0009132,980,7920CG25.8% 57.4 / 21.1 31G464A (GGC→GCC) araEarabinose:H(+) symporter
*NC_0009134,421,8350AG25.0% 75.3 / 15.7 40Y43C (TAC→TGC) ulaCL‑ascorbate specific PTS enzyme IIA component
*NC_0009134,424,8640TC25.0% 56.6 / 20.8 36intergenic (‑73/‑254)yjfY/rpsFDUF1471 domain‑containing protein YjfY/30S ribosomal subunit protein S6
*NC_0009133,717,7710AG24.0% 30.3 / 11.2 25G154G (GGA→GGGghrBglyoxylate reductase
*NC_000913366,5260TG23.6% 1427.7 / 87.4 651T329P (ACC→CCC) lacIDNA‑binding transcriptional repressor LacI
*NC_0009133,474,8220AG21.9% 75.3 / 17.9 41R48R (CGT→CGCtusBsulfurtransferase complex subunit TusB
*NC_0009133,358,1200GA21.1% 75.4 / 15.5 38L1132L (CTG→CTAgltBglutamate synthase subunit GltB
*NC_000913829,7251.G20.0% 67.5 / 20.8 30coding (248/999 nt)ybhGHlyD_D23 family protein YbhG
*NC_0009134,606,0941.G20.0% 42.1 / 15.1 20noncoding (72/87 nt)leuVtRNA‑Leu

Marginal new junction evidence (sorted from low to high skew)
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 445696933 (1.270)5 (0.200) 5/478 2.3 13.3% coding (1185/1374 nt) mpl UDP‑N‑acetylmuramate‑‑L‑alanyl‑gamma‑D‑glutamyl‑ meso‑2,6‑diaminoheptanedioate ligase
?NC_000913 4456977 = 33 (1.310)coding (1193/1374 nt) mpl UDP‑N‑acetylmuramate‑‑L‑alanyl‑gamma‑D‑glutamyl‑ meso‑2,6‑diaminoheptanedioate ligase