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

Marginal read alignment evidence (highest frequency 20 of 36 shown, sorted by frequency from high to low)
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913339,9860GA73.2% 62.8 / 25.2 41R73Q (CGA→CAA) yahHputative uncharacterized protein YahH
*NC_000913339,9790GA72.5% 59.0 / 21.5 40V71M (GTG→ATG) yahHputative uncharacterized protein YahH
*NC_000913190,8430CA40.6% 128.0 / 11.1 96intergenic (+244/‑14)rpsB/tsf30S ribosomal subunit protein S2/protein chain elongation factor EF‑Ts
*NC_0009134,136,4670AC35.1% 88.3 / 13.3 60E120D (GAA→GACyijEcystine exporter
*NC_0009131,427,9920GT34.2% 165.2 / 25.3 117P245P (CCC→CCAinsH5IS5 family transposase and trans‑activator
*NC_0009132,477,2650AC33.7% 184.9 / 32.6 104L122* (TTA→TGA) dsdCDNA‑binding transcriptional dual regulator DsdC
*NC_000913225,5790CA32.4% 118.0 / 15.0 105intergenic (+4/‑180)alaV/rrlHtRNA‑Ala/23S ribosomal RNA
*NC_0009132,159,4210AG32.2% 96.6 / 10.8 63I346V (ATT→GTT) mdtCmultidrug efflux pump RND permease subunit MdtC
*NC_000913847,9550GA32.1% 160.6 / 52.7 132A17V (GCG→GTG) glnHL‑glutamine ABC transporter periplasmic binding protein
*NC_0009132,159,4130CT31.7% 88.4 / 20.9 61A343V (GCG→GTG) mdtCmultidrug efflux pump RND permease subunit MdtC
*NC_0009133,054,8860GT31.2% 191.8 / 29.5 112Q187K (CAA→AAA) ygfBPF03695 family protein YgfB
*NC_0009133,054,8940GT31.2% 180.4 / 31.7 113P184Q (CCA→CAA) ygfBPF03695 family protein YgfB
*NC_0009131,922,8150CA29.8% 173.6 / 16.0 114G52V (GGA→GTA) pphAphosphoprotein phosphatase 1
*NC_0009132,767,4450CA29.5% 178.4 / 12.1 96intergenic (+90/‑265)rnlB/yfjPantitoxin RnlB/putative GTP‑binding protein YfjP
*NC_0009133,401,3450GT28.6% 217.9 / 10.1 112intergenic (‑258/+47)mreB/csrDdynamic cytoskeletal protein MreB/regulator of CsrB and CsrC decay
*NC_0009131,334,9660AC28.5% 269.2 / 28.1 158intergenic (+137/‑182)cysB/ymiADNA‑binding transcriptional dual regulator CysB/uncharacterized protein YmiA
*NC_0009133,741,8710GT27.8% 184.4 / 17.5 116S221Y (TCC→TAC) plaRDNA‑binding transcriptional repressor PlaR
*NC_0009132,490,2400GT27.6% 294.3 / 14.0 152intergenic (‑54/+16)yfdV/oxcputative transport protein YfdV/oxalyl‑CoA decarboxylase
*NC_0009132,477,2820GA27.6% 221.0 / 11.7 107C116C (TGC→TGTdsdCDNA‑binding transcriptional dual regulator DsdC
*NC_0009132,365,8400TG27.5% 274.5 / 24.6 178intergenic (‑220/‑70)ais/arnBputative lipopolysaccharide core heptose(II)‑phosphate phosphatase/UDP‑4‑amino‑4‑deoxy‑L‑arabinose aminotransferase

Marginal new junction evidence (lowest skew 10 of 63 shown)
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4353803 =92 (0.870)33 (0.310) 27/598 3.3 26.4% coding (915/1518 nt) lysU lysine‑‑tRNA ligase/Ap4A synthetase/Ap3A synthetase
?NC_000913 4498181 = NA (NA)noncoding (1/1331 nt) IS2 repeat region
* ? NC_000913 3583428 =NA (NA)38 (0.360) 27/598 3.3 37.6% noncoding (1/768 nt) IS1 repeat region
?NC_000913 3889014 = 63 (0.600)coding (285/1416 nt) tnaA tryptophanase
* ? NC_000913 3365556 =NA (NA)29 (0.280) 25/590 3.5 17.8% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 3890222 = 136 (1.290)intergenic (+77/‑14) tnaA/tnaB tryptophanase/tryptophan:H(+) symporter TnaB
* ? NC_000913 = 275149NA (NA)26 (0.250) 23/598 4.0 22.0% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 538086 = 92 (0.870)coding (454/1455 nt) ybbW putative allantoin transporter
* ? NC_000913 2690293 =91 (0.860)21 (0.200) 18/598 5.1 18.8% coding (806/1428 nt) glrK histidine kinase GlrK
?NC_000913 2996361 = NA (NA)noncoding (1331/1331 nt) IS2 repeat region
* ? NC_000913 = 2070271NA (NA)19 (0.180) 17/598 5.3 17.4% noncoding (1/1331 nt) IS2 repeat region
?NC_000913 = 2690297 90 (0.850)coding (802/1428 nt) glrK histidine kinase GlrK
* ? NC_000913 3720632 =NA (NA)17 (0.160) 15/598 5.8 NA intergenic (‑32/‑48) hokA/insJ small toxic polypeptide/insertion element IS150 protein InsA
?NC_000913 = 3722077 NA (NA)intergenic (+28/+251) insK/glyS IS150 family conserved protein InsB/glycine‑‑tRNA ligase subunit beta
* ? NC_000913 = 4606166106 (1.000)9 (0.120) 8/420 6.1 11.4% intergenic (‑1/+34) leuV/leuP tRNA‑Leu/tRNA‑Leu
?NC_000913 4606403 = 66 (0.890)intergenic (‑2/+266) leuQ/rsmC tRNA‑Leu/16S rRNA m(2)G1207 methyltransferase
* ? NC_000913 = 3648146113 (1.070)10 (0.120) 10/482 6.1 9.8% intergenic (+2/‑148) agrA/agrB small RNA AgrA/small regulatory RNA AgrB
?NC_000913 3648378 = 94 (1.100)intergenic (+3/‑150) agrB/arsR small regulatory RNA AgrB/DNA‑binding transcriptional repressor ArsR
* ? NC_000913 = 291401NA (NA)16 (0.150) 14/598 6.1 20.2% noncoding (1/768 nt) IS1 repeat region
?NC_000913 3889014 = 63 (0.600)coding (285/1416 nt) tnaA tryptophanase