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_0009132,001,9110GA37.1% 31.5 / 16.6 35intergenic (‑122/+1)fliA/fliXRNA polymerase sigma factor FliA/small RNA FliX
*NC_0009133,875,5250GT32.4% 54.6 / 12.3 37V30F (GTT→TTT) yidXputative lipoprotein YidX
*NC_0009133,875,5320AT32.4% 36.5 / 21.6 38H32L (CAT→CTT) ‡yidXputative lipoprotein YidX
*NC_000913388,3360AC32.1% 79.2 / 10.1 54N181T (AAC→ACC) tauDtaurine dioxygenase
*NC_0009131,984,1920GC31.6% 77.3 / 21.2 57A293G (GCG→GGG) araGarabinose ABC transporter ATP binding subunit
*NC_0009133,875,5210TC31.4% 44.4 / 14.3 38A28A (GCT→GCCyidXputative lipoprotein YidX
*NC_0009132,828,5060TG30.8% 84.5 / 10.9 69V257G (GTG→GGG) srlDsorbitol‑6‑phosphate 2‑dehydrogenase
*NC_0009134,256,0610CT30.2% 55.9 / 12.4 44E136K (GAA→AAA) plsBglycerol‑3‑phosphate 1‑O‑acyltransferase
*NC_0009133,755,3630AC29.5% 121.7 / 18.4 78G383G (GGT→GGGaldBaldehyde dehydrogenase B
*NC_0009134,256,0690GC28.9% 77.5 / 10.7 47P133R (CCG→CGG) plsBglycerol‑3‑phosphate 1‑O‑acyltransferase
*NC_000913811,7210TG28.4% 145.8 / 11.6 96V67G (GTG→GGG) bioCmalonyl‑acyl carrier protein methyltransferase
*NC_0009131,972,0900AC28.2% 135.9 / 12.1 79V201G (GTG→GGG) tarmethyl‑accepting chemotaxis protein Tar
*NC_0009133,875,5310CT27.8% 57.4 / 10.9 37H32Y (CAT→TAT) ‡yidXputative lipoprotein YidX
*NC_0009134,296,2680TC27.0% 41.7 / 16.3 38intergenic (+474/+168)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
*NC_0009134,561,2150GA26.5% 88.4 / 12.7 49E96K (GAA→AAA) kptARNA 2'‑phosphotransferase
*NC_0009134,502,7610GC25.9% 124.9 / 14.0 60A224G (GCG→GGG) insGIS4 family putative transposase
*NC_0009131,415,3950GA25.4% 122.8 / 12.2 64R664R (CGC→CGTrecEexonuclease VIII, ds DNA exonuclease, 5' ‑‑> 3' specific
*NC_0009132,370,2660TG25.0% 123.9 / 11.3 72G83G (GGT→GGGarnDputative 4‑deoxy‑4‑formamido‑L‑arabinose‑phosphoundecaprenol deformylase ArnD
*NC_0009132,917,0880GC22.9% 171.7 / 11.8 71A678P (GCT→CCT) barAsensor histidine kinase BarA
*NC_0009134,166,8490GA22.6% 27.6 / 10.3 31noncoding (191/1542 nt)rrsB16S ribosomal RNA

Marginal new junction evidence (lowest skew 10 of 24 shown)
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4507459 =NA (NA)17 (0.210) 16/486 3.1 23.9% noncoding (1221/1221 nt) IS30 repeat region
?NC_000913 = 4640419 54 (0.670)coding (18/141 nt) yjjY uncharacterized protein YjjY
* ? NC_000913 = 280111NA (NA)17 (0.210) 16/486 3.1 23.9% noncoding (181/181 nt) IS30 repeat region
?NC_000913 = 4640419 54 (0.670)coding (18/141 nt) yjjY uncharacterized protein YjjY
* ? NC_000913 = 271761NA (NA)17 (0.210) 16/484 3.1 24.0% noncoding (1221/1221 nt) IS30 repeat region
?NC_000913 = 4640419 54 (0.670)coding (18/141 nt) yjjY uncharacterized protein YjjY
* ? NC_000913 = 4508679NA (NA)15 (0.190) 15/488 3.3 21.4% noncoding (1/1221 nt) IS30 repeat region
?NC_000913 4640418 = 55 (0.680)coding (17/141 nt) yjjY uncharacterized protein YjjY
* ? NC_000913 668961 =96 (1.180)13 (0.160) 12/490 3.8 11.9% coding (76/318 nt) rsfS ribosomal silencing factor RsfS
?NC_000913 = 689045 NA (NA)noncoding (1/1195 nt) IS5 repeat region
* ? NC_000913 2804741 =57 (0.700)10 (0.120) 10/490 4.3 14.9% intergenic (+280/‑74) nrdF/proV ribonucleoside‑diphosphate reductase 2 subunit beta/glycine betaine ABC transporter ATP binding subunit ProV
?NC_000913 = 3131340 NA (NA)noncoding (1195/1195 nt) IS5 repeat region
* ? NC_000913 607765 =31 (0.380)8 (0.100) 8/488 4.8 20.5% noncoding (28/60 nt) sokE small RNA SokE
?NC_000913 = 3722075 NA (NA)noncoding (1443/1443 nt) IS150 repeat region
* ? NC_000913 1978503 =NA (NA)7 (0.090) 6/488 5.3 10.1% noncoding (768/768 nt) IS1 repeat region
?NC_000913 4307900 = 62 (0.770)coding (579/696 nt) alsE D‑allulose‑6‑phosphate 3‑epimerase
* ? NC_000913 = 275149NA (NA)5 (0.060) 5/490 5.7 8.5% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 1091173 = 54 (0.670)coding (1117/2424 nt) pgaA partially deacetylated poly‑beta‑1,6‑N‑acetyl‑D‑glucosamine export outer membrane porin
* ? NC_000913 273955 =NA (NA)5 (0.060) 5/490 5.7 9.4% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 = 295085 48 (0.590)coding (58/225 nt) ykgV protein YkgV