breseq  version 0.32.1  revision aa8f2595a244
mutation predictions | marginal predictions | summary statistics | genome diff | command line log

Marginal read alignment evidence (highest frequency 20 of 24 shown, sorted by frequency from high to low)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,286,9650AG60.9% 28.7 / 63.8 69intergenic (+439/+101)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
*NC_0009132,857,1440CT60.0% 2.8 / 49.3 32L18F (CTC→TTC) ‡mutSmethyl‑directed mismatch repair protein
*NC_0009132,857,1450TA59.4% 11.5 / 42.9 32L18H (CTC→CAC) ‡mutSmethyl‑directed mismatch repair protein
*NC_0009132,857,1460CG59.4% 16.7 / 48.6 32L18L (CTC→CTG) ‡mutSmethyl‑directed mismatch repair protein
*NC_0009132,857,1470AT59.4% 56.0 / 42.4 32coding (55/2562 nt)mutSmethyl‑directed mismatch repair protein
*NC_0009132,857,1480GC58.1% 53.4 / 44.6 32coding (56/2562 nt)mutSmethyl‑directed mismatch repair protein
*NC_0009132,857,1420AG56.2% 5.3 / 33.2 32Y17C (TAT→TGT) mutSmethyl‑directed mismatch repair protein
*NC_0009132,857,1430TA56.2% 2.4 / 47.1 32Y17* (TAT→TAAmutSmethyl‑directed mismatch repair protein
*NC_0009134,171,5570CA55.6% ‑0.7 / 47.0 27intergenic (+13/‑80)rrlB/rrfB23S ribosomal RNA of rrnB operon/5S ribosomal RNA of rrnB operon
*NC_0009134,296,0600CT49.2% 17.2 / 84.1 65intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
*NC_0009132,857,1411.G43.8% ‑5.3 / 63.5 32coding (49/2562 nt)mutSmethyl‑directed mismatch repair protein
*NC_000913566,2770CT40.0% 75.4 / 91.2 90pseudogene (100/279 nt)exoDpseudogene, DLP12 prophage, phage‑type exonuclease family,Phage or Prophage Related
*NC_000913566,2450GA39.7% 61.2 / 71.5 73pseudogene (132/279 nt)exoDpseudogene, DLP12 prophage, phage‑type exonuclease family,Phage or Prophage Related
*NC_000913566,1730CG39.4% 71.7 / 102.2 71pseudogene (204/279 nt)exoDpseudogene, DLP12 prophage, phage‑type exonuclease family,Phage or Prophage Related
*NC_0009132,857,1431.A37.5% 10.5 / 55.9 32coding (51/2562 nt)mutSmethyl‑directed mismatch repair protein
*NC_000913566,3230CT36.7% 79.8 / 70.8 80pseudogene (54/279 nt)exoDpseudogene, DLP12 prophage, phage‑type exonuclease family,Phage or Prophage Related
*NC_000913566,2050TC34.4% 50.4 / 57.6 61pseudogene (172/279 nt)exoDpseudogene, DLP12 prophage, phage‑type exonuclease family,Phage or Prophage Related
*NC_000913566,3260TC32.5% 84.3 / 69.3 83pseudogene (51/279 nt)exoDpseudogene, DLP12 prophage, phage‑type exonuclease family,Phage or Prophage Related
*NC_000913566,3320TG31.7% 120.4 / 96.2 82pseudogene (45/279 nt)exoDpseudogene, DLP12 prophage, phage‑type exonuclease family,Phage or Prophage Related
*NC_0009133,728,1102.C29.2% 32.9 / 25.9 24coding (1262/1455 nt)xylBxylulokinase

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 = 120780532 (0.680)7 (0.180) 6/178 3.0 18.2% coding (305/630 nt) stfP e14 prophage, uncharacterized protein
?NC_000913 = 1209602 36 (0.910)pseudogene (18/501 nt) stfE pseudogene, e14 prophage, side tail fiber protein fragment family,Phage or Prophage Related
* ? NC_000913 = 132314752 (1.110)13 (0.280) 5/210 3.7 33.3% intergenic (‑65/‑73) trpL/yciV trp operon leader peptide/PHP domain protein
?NC_000913 = 3742633 0 (0.000)intergenic (‑101/‑100) yiaJ/yiaK transcriptional repressor for the yiaKLMNO‑lyxK‑sgbHUE operon/2,3‑diketo‑L‑gulonate reductase, NADH‑dependent
* ? NC_000913 2765212 =76 (1.620)3 (0.070) 3/194 4.3 4.1% intergenic (‑59/+301) yfjL/yfjM CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein/CP4‑57 prophage, uncharacterized protein
?NC_000913 2765237 = 70 (1.610)intergenic (‑84/+276) yfjL/yfjM CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein/CP4‑57 prophage, uncharacterized protein
* ? NC_000913 = 20665384 (0.090)3 (0.070) 3/198 4.4 44.3% noncoding (816/1195 nt) IS5 repeat region
?NC_000913 = 2102132 NA (NA)noncoding (812/1195 nt) IS5 repeat region
* ? NC_000913 = 274334NA (NA)3 (0.070) 3/198 4.4 NA noncoding (816/1195 nt) IS5 repeat region
?NC_000913 = 274338 NA (NA)noncoding (812/1195 nt) IS5 repeat region
* ? NC_000913 = 159861250 (1.060)5 (0.110) 3/210 4.5 9.2% intergenic (‑625/+5) yneO/lsrK pseudogene, AidA homolog/autoinducer‑2 (AI‑2) kinase
?NC_000913 1598618 = 49 (1.040)coding (1592/1593 nt) lsrK autoinducer‑2 (AI‑2) kinase