breseq  version 0.29.0  revision 8f9c342918e4
mutation predictions | marginal predictions | summary statistics | genome diff | command line log

Marginal read alignment evidence (highest frequency 20 of 550 shown)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091315,4100TG51.9% 53.9 / 83.8 53intergenic (+112/‑35)dnaJ/insL1chaperone Hsp40, DnaK co‑chaperone/IS186 transposase
*NC_000913354,2880TG48.5% 98.8 / 16.4 68V528G (GTG→GGG) prpEpropionate‑‑CoA ligase
*NC_0009133,556,2070TC47.8% 74.2 / 26.5 67S215G (AGT→GGT) rtcARNA 3'‑terminal phosphate cyclase
*NC_0009134,381,7280AC46.7% 77.8 / 14.2 60G197G (GGT→GGGfrdAanaerobic fumarate reductase catalytic and NAD/flavoprotein subunit
*NC_000913354,2820TG46.2% 99.8 / 11.6 65V526G (GTT→GGT) prpEpropionate‑‑CoA ligase
*NC_000913925,9000AG45.6% 105.7 / 17.4 69noncoding (72/88 nt)serWtRNA‑Ser
*NC_000913718,6590TA44.6% 106.1 / 36.3 74D601V (GAT→GTT) speFornithine decarboxylase isozyme, inducible
*NC_0009133,556,2000GA43.9% 74.5 / 33.0 66S217F (TCC→TTC) rtcARNA 3'‑terminal phosphate cyclase
*NC_0009134,332,7910AC43.9% 100.8 / 13.6 57V161G (GTG→GGG) basSsensory histidine kinase in two‑component regulatory system with BasR
*NC_0009133,665,7020CA43.6% 147.6 / 56.4 101A37S (GCG→TCG) gadXacid resistance regulon transcriptional activator; autoactivator
*NC_00091315,3910AC43.5% 78.5 / 28.4 62intergenic (+93/‑54)dnaJ/insL1chaperone Hsp40, DnaK co‑chaperone/IS186 transposase
*NC_0009131,153,5730TG43.5% 99.1 / 24.4 62Y92D (TAC→GAC) pabC4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex
*NC_0009134,465,7870GA43.2% 94.2 / 35.8 74R132C (CGC→TGC) treBtrehalose‑specific PTS enzyme: IIB and IIC component
*NC_000913357,8900CA43.2% 47.3 / 25.1 44L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
*NC_0009131,604,2370CT42.5% 49.7 / 20.8 47S64L (TCG→TTG) lsrDautoinducer 2 import system permease protein
*NC_000913436,5160CA41.5% 90.2 / 18.1 53Q295K (CAA→AAA) thiLthiamine monophosphate kinase
*NC_0009131,701,3100AC41.4% 111.4 / 11.2 71T118T (ACA→ACCmalYPLP‑dependent beta‑cystathionase and maltose regulon regulator
*NC_00091315,3950GA41.4% 66.1 / 26.2 59intergenic (+97/‑50)dnaJ/insL1chaperone Hsp40, DnaK co‑chaperone/IS186 transposase
*NC_0009132,349,1450AC41.3% 109.6 / 14.0 63F110V (TTC→GTC) inaAacid‑inducible Kdo/WaaP family putative kinase
*NC_0009132,578,1430AC41.0% 135.1 / 28.5 78L79V (TTA→GTA) maeBmalic enzyme: putative oxidoreductase/phosphotransacetylase

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 =86 (1.140)4 (0.060) 4/246 12.1 6.3% coding (50/936 nt) ttcA tRNA s(2)C32 thioltransferase, iron sulfur cluster protein
?NC_000913 = 1434958 47 (0.750)pseudogene (51/51 nt) ttcC pseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related