breseq  version 0.27.1  revision 87c22d663cc3
mutation predictions | marginal predictions | summary statistics | genome diff | command line log

Marginal read alignment evidence (highest frequency 20 of 46 shown)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,151,7590GA52.1% 56.1 / inf 341E49K (GAA→AAA) acpPacyl carrier protein (ACP)
*NC_0009132,610,1070GA49.1% 16.8 / inf 390E541K (GAG→AAG) hyfGhydrogenase 4, subunit
*NC_0009134,186,4250CT48.7% 31.7 / inf 415P359L (CCA→CTA) rpoCRNA polymerase, beta prime subunit
*NC_0009133,473,6120AC46.5% 37.9 / inf 406L157* (TTA→TGA) rpsG30S ribosomal subunit protein S7
*NC_0009131,907,3270GA46.3% 89.9 / inf 341Q37* (CAG→TAG) cspCstress protein, member of the CspA‑family
*NC_0009134,187,6670TA43.0% 212.1 / inf 377F773Y (TTC→TAC) rpoCRNA polymerase, beta prime subunit
*NC_000913363,4000A.42.6% 706.0 / inf 325coding (2906/3075 nt)lacZbeta‑D‑galactosidase
*NC_0009134,213,0510AC42.2% 63.2 / 169.8 116noncoding (12/120 nt)rrfE5S ribosomal RNA of rrnE operon
*NC_0009131,674,7710G.38.2% 489.0 / inf 213coding (859/1035 nt)tqsApheromone AI‑2 transporter
*NC_0009131,674,7700T.38.0% 469.8 / inf 214coding (858/1035 nt)tqsApheromone AI‑2 transporter
*NC_0009131,674,7720G.38.0% 499.3 / inf 213coding (860/1035 nt)tqsApheromone AI‑2 transporter
*NC_0009131,674,7730G.38.0% 500.7 / inf 213coding (861/1035 nt)tqsApheromone AI‑2 transporter
*NC_0009131,674,7680C.37.7% 509.4 / inf 215coding (856/1035 nt)tqsApheromone AI‑2 transporter
*NC_0009131,674,7690G.37.7% 504.9 / inf 215coding (857/1035 nt)tqsApheromone AI‑2 transporter
*NC_0009131,674,7740C.37.2% 504.3 / inf 215coding (862/1035 nt)tqsApheromone AI‑2 transporter
*NC_0009132,947,5910A.35.8% 157.7 / 90.8 67intergenic (+18/‑16)metW/metVtRNA‑Met/tRNA‑Met
*NC_0009131,674,7660G.35.1% 568.5 / inf 233coding (854/1035 nt)tqsApheromone AI‑2 transporter
*NC_0009131,674,7650G.34.9% 566.1 / inf 233coding (853/1035 nt)tqsApheromone AI‑2 transporter
*NC_0009131,674,7670G.34.8% 571.4 / inf 233coding (855/1035 nt)tqsApheromone AI‑2 transporter
*NC_0009131,674,7640G.34.3% 572.4 / inf 233coding (852/1035 nt)tqsApheromone AI‑2 transporter

Marginal new junction evidence (lowest skew 10 of 59 shown)
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2288919 =NA (NA)254 (0.670) 132/300 3.4 52.2% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 = 3393819 233 (0.610)coding (2440/3801 nt) yhdP DUF3971‑AsmA2 domains protein
* ? NC_000913 = 1674768134 (0.350)160 (0.440) 126/288 3.5 55.3% coding (856/1035 nt) tqsA pheromone AI‑2 transporter
?NC_000913 1674781 = 130 (0.360)coding (869/1035 nt) tqsA pheromone AI‑2 transporter
* ? NC_000913 = 1104115154 (0.400)170 (0.450) 122/298 4.3 51.7% coding (165/456 nt) csgB curlin nucleator protein, minor subunit in curli complex
?NC_000913 1104420 = 165 (0.440)intergenic (+14/‑27) csgB/csgA curlin nucleator protein, minor subunit in curli complex/curlin subunit, amyloid curli fibers, cryptic
* ? NC_000913 2288919 =NA (NA)271 (0.720) 114/298 5.2 55.3% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 = 2413641 219 (0.580)coding (172/1203 nt) ackA acetate kinase A and propionate kinase 2
* ? NC_000913 = 2290113NA (NA)237 (0.620) 107/300 6.2 50.9% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 2413638 = 229 (0.600)coding (169/1203 nt) ackA acetate kinase A and propionate kinase 2
* ? NC_000913 = 4542690419 (1.100)81 (0.230) 71/282 11.1 16.7% intergenic (+57/‑425) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
?NC_000913 = 4542986 417 (1.160)intergenic (+353/‑129) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
* ? NC_000913 4542682 =419 (1.100)84 (0.230) 69/282 11.5 17.0% intergenic (+49/‑433) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
?NC_000913 4542996 = 426 (1.190)intergenic (+363/‑119) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
* ? NC_000913 = 1287084148 (0.390)84 (0.250) 45/260 16.1 31.4% intergenic (+558/+3) narI/tpr nitrate reductase 1, gamma (cytochrome b(NR)) subunit/protamine‑like protein
?NC_000913 1287557 = 239 (0.720)noncoding (66/85 nt) tyrT tRNA‑Tyr
* ? NC_000913 780554 =428 (1.120)101 (0.270)
+GAAG
52/292 16.3 26.2% noncoding (1/76 nt) lysT tRNA‑Lys
?NC_000913 = 2521006 156 (0.410)noncoding (76/76 nt) valU tRNA‑Val
* ? NC_000913 781577 =258 (0.680)179 (0.480)
+GAAG
51/292 16.5 47.0% noncoding (1/76 nt) lysQ tRNA‑Lys
?NC_000913 = 2521006 156 (0.410)noncoding (76/76 nt) valU tRNA‑Val