breseq  version 0.35.4  revision f352f80f4bc9
mutation predictions | marginal predictions | summary statistics | genome diff | command line log

Marginal read alignment evidence (highest frequency 20 of 26 shown, sorted by frequency from high to low)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,280,7040AC62.1% 3.5 / 27.3 29T25P (ACC→CCC) radDputative DNA repair helicase RadD
*NC_0009134,296,0600CT50.0% 2.2 / 119.1 80intergenic (+266/+376)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
*NC_0009131,374,0380TC42.9% 7.3 / 21.2 21F273L (TTC→CTC) ‡ycjOputative ABC transporter membrane subunit YcjO
*NC_0009131,374,0340G.42.1% 28.1 / 32.2 19coding (813/882 nt)ycjOputative ABC transporter membrane subunit YcjO
*NC_0009134,414,9060TG40.9% 13.0 / 23.7 22V211G (GTG→GGG) aidBputative acyl‑CoA dehydrogenase AidB
*NC_000913395,2880AT40.6% 26.1 / 44.7 32F334I (TTC→ATC) ampHpeptidoglycan DD‑carboxypeptidase/peptidoglycan DD‑endopeptidase
*NC_0009131,374,0360TG40.0% 11.2 / 20.6 20L272R (CTC→CGC) ycjOputative ABC transporter membrane subunit YcjO
*NC_0009132,349,6210AC40.0% 15.9 / 29.7 25intergenic (‑149/‑66)inaA/yfaHputative lipopolysaccharide kinase InaA/putative uncharacterized protein YfaH
*NC_000913395,2910CG39.4% 24.3 / 44.9 33G333R (GGT→CGT) ampHpeptidoglycan DD‑carboxypeptidase/peptidoglycan DD‑endopeptidase
*NC_000913395,2940CG36.4% 27.0 / 41.6 33G332R (GGC→CGC) ampHpeptidoglycan DD‑carboxypeptidase/peptidoglycan DD‑endopeptidase
*NC_0009134,021,8240TG34.8% 26.4 / 16.1 23G533G (GGT→GGGubiBubiquinone biosynthesis protein UbiB
*NC_0009131,374,0390TG33.3% 36.5 / 22.6 27F273C (TTC→TGC) ‡ycjOputative ABC transporter membrane subunit YcjO
*NC_000913794,9760TC32.4% 48.6 / 12.2 34intergenic (+54/‑113)acrZ/modAmultidrug efflux pump accessory protein AcrZ/molybdate ABC transporter periplasmic binding protein
*NC_0009131,374,0400CT32.1% 24.9 / 22.7 28F273F (TTC→TTT) ‡ycjOputative ABC transporter membrane subunit YcjO
*NC_0009132,814,8960TC31.0% 34.2 / 10.6 29E515G (GAA→GGA) gshAglutamate‑‑cysteine ligase
*NC_0009133,812,3860CA28.6% 25.1 / 17.6 28L58F (TTG→TTTyicRRadC‑like JAB domain‑containing protein YicR
*NC_0009132,794,7870CT27.8% 39.3 / 19.3 36P179S (CCG→TCG) gabP4‑aminobutyrate:H(+) symporter
*NC_0009133,697,7400TG27.0% 72.9 / 16.5 37V428G (GTG→GGG) bcsEc‑di‑GMP‑binding protein BcsE
*NC_0009133,563,4810TG25.0% 46.2 / 12.0 24V490G (GTG→GGG) glpDaerobic glycerol 3‑phosphate dehydrogenase
*NC_0009131,222,2060CA22.2% 69.0 / 10.9 27pseudogene (2606/2648 nt)ycgHputative transporter component, N‑terminal fragment

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 = 120780552 (1.090)10 (0.280) 8/96 1.2 18.5% coding (305/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209602 49 (1.370)pseudogene (54/537 nt) stfE e14 prophage; putative side tail fiber protein fragment