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

Marginal read alignment evidence (highest frequency 20 of 25 shown, sorted by frequency from high to low)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3500CT46.2% 18.9 / 14.0 26S64F (TCT→TTT) bglJDNA‑binding transcriptional regulator BglJ
*NC_0009133,665,7020CA41.9% 28.1 / 14.3 31A37S (GCG→TCG) gadXDNA‑binding transcriptional dual regulator GadX
*NC_0009131,022,9150TG41.9% 46.6 / 10.1 43A323A (GCA→GCCyccSputative transporter YccS
*NC_0009134,296,2680TC40.0% 21.6 / 22.8 25intergenic (+474/+168)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
*NC_0009134,604,3440CT38.5% 25.4 / 11.2 26S62F (TCT→TTT) bglJDNA‑binding transcriptional regulator BglJ
*NC_0009134,056,1150AC38.3% 73.1 / 11.6 48G75G (GGT→GGGglnLprotein histidine kinase NtrB
*NC_0009132,434,8820TG38.0% 83.6 / 11.2 50Y252S (TAC→TCC) truAtRNA pseudouridine(38‑40) synthase
*NC_0009132,326,4500TG37.5% 89.7 / 13.8 56G114G (GGT→GGGatoBacetyl‑CoA acetyltransferase
*NC_000913941,7130TG34.9% 121.4 / 20.2 64V252G (GTG→GGG) dmsAdimethyl sulfoxide reductase subunit A
*NC_0009133,665,6980TC33.3% 44.0 / 10.3 31D38G (GAT→GGT) gadXDNA‑binding transcriptional dual regulator GadX
*NC_0009134,171,5570CG33.3% 60.4 / 46.8 45intergenic (+13/‑80)rrlB/rrfB23S ribosomal RNA/5S ribosomal RNA
*NC_0009133,177,7480CA32.5% 56.1 / 11.5 41E55* (GAG→TAG) nudFADP‑sugar pyrophosphatase
*NC_000913811,7210TG32.0% 125.8 / 12.0 80V67G (GTG→GGG) bioCmalonyl‑acyl carrier protein methyltransferase
*NC_0009133,556,2070TC31.7% 76.2 / 12.0 41S215G (AGT→GGT) rtcARNA 3'‑terminal phosphate cyclase
*NC_0009133,556,2000GA29.3% 66.7 / 12.7 41S217F (TCC→TTC) rtcARNA 3'‑terminal phosphate cyclase
*NC_0009134,482,6480AC28.6% 132.1 / 11.1 70V397G (GTT→GGT) valSvaline‑‑tRNA ligase
*NC_0009132,185,1990TG28.2% 172.0 / 10.6 78T35P (ACC→CCC) thiMhydroxyethylthiazole kinase
*NC_0009134,223,6440CT25.5% 108.4 / 11.3 47intergenic (‑16/‑184)iclR/metHDNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase
*NC_000913941,7270GT25.4% 149.1 / 10.1 67E257* (GAA→TAA) dmsAdimethyl sulfoxide reductase subunit A
*NC_0009132,185,2380AG25.0% 195.9 / 12.4 80S22P (TCC→CCC) thiMhydroxyethylthiazole kinase

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 1411899 =57 (0.970)4 (0.080) 4/442 4.2 9.9% coding (50/936 nt) ttcA tRNA cytosine(32) 2‑sulfurtransferase TtcA
?NC_000913 = 1434958 22 (0.420)intergenic (‑552/+227) ynaM/uspF protein YnaM/universal stress protein F