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

Marginal read alignment evidence (highest frequency 20 of 313 shown, sorted by frequency from high to low)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,650,2660GA66.2% 441.7 / inf 345P10P (CCG→CCAarsCarsenate reductase
*NC_0009133,650,2480TC66.0% 389.5 / inf 341I4I (ATT→ATCarsCarsenate reductase
*NC_0009133,650,1910TC62.1% 254.3 / inf 298A419A (GCT→GCCarsBarsenite/antimonite:H(+) antiporter
*NC_0009133,650,1820GA60.9% 233.4 / inf 300V416V (GTG→GTAarsBarsenite/antimonite:H(+) antiporter
*NC_0009133,650,3690GC59.6% 250.2 / inf 301V45L (GTC→CTC) arsCarsenate reductase
*NC_0009133,650,3830CG59.5% 250.3 / inf 296A49A (GCC→GCGarsCarsenate reductase
*NC_0009133,650,3410GT58.8% 163.8 / inf 301L35L (CTG→CTTarsCarsenate reductase
*NC_0009133,650,3440AG58.8% 179.8 / inf 301E36E (GAA→GAGarsCarsenate reductase
*NC_0009133,650,3540AT58.8% 243.5 / inf 302T40S (ACG→TCG) ‡arsCarsenate reductase
*NC_0009133,650,3340AT58.6% 245.6 / inf 303H33L (CAT→CTT) arsCarsenate reductase
*NC_0009133,650,3560GA58.5% 213.3 / inf 302T40T (ACG→ACA) ‡arsCarsenate reductase
*NC_0009133,650,3650AG58.4% 175.7 / inf 303E43E (GAA→GAGarsCarsenate reductase
*NC_0009133,650,3530AT58.3% 230.3 / inf 300P39P (CCA→CCTarsCarsenate reductase
*NC_0009133,650,3270AG58.2% 167.2 / inf 304I31V (ATT→GTT) arsCarsenate reductase
*NC_0009133,650,3500GA58.1% 196.4 / inf 302P38P (CCG→CCAarsCarsenate reductase
*NC_0009133,650,3260TC58.1% 155.0 / inf 303T30T (ACT→ACCarsCarsenate reductase
*NC_0009133,650,3320CT58.1% 216.7 / inf 303I32I (ATC→ATTarsCarsenate reductase
*NC_0009133,650,3380TC57.9% 126.3 / inf 303Y34Y (TAT→TACarsCarsenate reductase
*NC_0009131,634,3820TC57.8% 6.2 / 58.0 45K168K (AAA→AAGtfaQputative tail fiber assembly protein TfaQ
*NC_0009133,650,4080CT57.6% 171.2 / inf 264L58L (CTG→TTG) arsCarsenate reductase

Marginal new junction evidence (lowest skew 10 of 15 shown)
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 397085054 (0.530)38 (0.370) 26/598 3.2 41.3% coding (718/1131 nt) wecB UDP‑N‑acetylglucosamine 2‑epimerase
?NC_000913 4498181 = NA (NA)noncoding (1/1331 nt) IS2 repeat region
* ? NC_000913 = 3187426NA (NA)29 (0.290) 23/596 3.7 32.6% noncoding (1331/1331 nt) IS2 repeat region
?NC_000913 3970846 = 60 (0.590)coding (714/1131 nt) wecB UDP‑N‑acetylglucosamine 2‑epimerase
* ? NC_000913 = 566774101 (0.990)46 (0.450) 19/594 4.5 30.4% pseudogene (91/93 nt) renD protein RenD
?NC_000913 568036 = 110 (1.090)pseudogene (2/213 nt) renD protein RenD
* ? NC_000913 = 69674166 (0.650)12 (0.160) 11/442 5.0 18.1% intergenic (‑1/+15) metU/glnW tRNA‑Met/tRNA‑Gln
?NC_000913 697135 = 60 (0.800)intergenic (‑2/+378) metT/asnB tRNA‑Met/asparagine synthetase B
* ? NC_000913 = 4606166104 (1.020)8 (0.110) 8/420 5.8 10.8% intergenic (‑1/+34) leuV/leuP tRNA‑Leu/tRNA‑Leu
?NC_000913 4606403 = 59 (0.820)intergenic (‑2/+266) leuQ/rsmC tRNA‑Leu/16S rRNA m(2)G1207 methyltransferase
* ? NC_000913 = 294746789 (0.870)6 (0.080) 6/432 6.9 8.2% intergenic (+4/‑30) metZ/metW tRNA‑Met/tRNA‑Met
?NC_000913 2947688 = 70 (0.950)intergenic (+5/+69) metV/amiC tRNA‑Met/N‑acetylmuramoyl‑L‑alanine amidase C
* ? NC_000913 = 3648146137 (1.340)7 (0.090) 6/482 7.5 6.4% intergenic (+2/‑148) agrA/agrB small RNA AgrA/small regulatory RNA AgrB
?NC_000913 3648378 = 94 (1.140)intergenic (+3/‑150) agrB/arsR small regulatory RNA AgrB/DNA‑binding transcriptional repressor ArsR
* ? NC_000913 = 128732959 (0.580)5 (0.070) 4/424 7.9 11.1% intergenic (‑1/+209) tyrV/tyrT tRNA‑Tyr/tRNA‑Tyr
?NC_000913 1287624 = 38 (0.530)intergenic (‑2/+158) tyrT/purU tRNA‑Tyr/formyltetrahydrofolate deformylase
* ? NC_000913 = 78144650 (0.490)4 (0.050) 4/442 8.1 8.9% intergenic (+2/‑131) lysZ/lysQ tRNA‑Lys/tRNA‑Lys
?NC_000913 781655 = 45 (0.600)intergenic (+3/‑430) lysQ/nadA tRNA‑Lys/quinolinate synthase
* ? NC_000913 = 128690654 (0.530)7 (0.070) 6/558 8.4 11.9% coding (198/276 nt) ychS putative uncharacterized protein YchS
?NC_000913 1287557 = 53 (0.560)noncoding (66/85 nt) tyrT tRNA‑Tyr