breseq  version 0.33.1  revision 8505477f25b3
mutation predictions | marginal predictions | summary statistics | genome diff | command line log

Predicted mutations
evidence position mutation freq annotation gene description
RA 70,434 C→T 61.9% N16N (AAC→AAT araC → DNA‑binding transcriptional dual regulator
RA 70,581 C→T 55.6% V65V (GTC→GTT araC → DNA‑binding transcriptional dual regulator
RA 70,740 T→G 66.4% G118G (GGT→GGG araC → DNA‑binding transcriptional dual regulator
RA 71,079 C→T 57.3% R231R (CGC→CGT araC → DNA‑binding transcriptional dual regulator
RA 71,082 T→C 56.5% I232I (ATT→ATC araC → DNA‑binding transcriptional dual regulator
RA 71,085 T→C 58.1% S233S (AGT→AGC araC → DNA‑binding transcriptional dual regulator
RA 71,109 T→C 55.9% T241T (ACT→ACC araC → DNA‑binding transcriptional dual regulator
RA 71,175 A→G 51.0% R263R (CGA→CGG araC → DNA‑binding transcriptional dual regulator
RA 71,214 T→C 63.4% F276F (TTT→TTC araC → DNA‑binding transcriptional dual regulator
RA 155,463 A→T 100% *247K (TAA→AAA)  ecpD ← predicted periplasmic pilin chaperone
RA 236,101 A→C 57.0% D12A (GAT→GCT)  dnaQ → DNA polymerase III epsilon subunit
RA 236,107 A→C 57.1% E14A (GAA→GCA)  dnaQ → DNA polymerase III epsilon subunit
RA 547,694 A→G 100% pseudogene (114/261 nt)
pseudogene (114/1259 nt)
ylbE →
ylbE →
ECK0512:JW0507:b4507; hypothetical protein, N‑ter fragment
ECK0512:JW0508+JW0507:b4572; hypothetical protein
RA 547,831:1 +G 100% pseudogene (251/261 nt)
pseudogene (251/1259 nt)
ylbE →
ylbE →
ECK0512:JW0507:b4507; hypothetical protein, N‑ter fragment
ECK0512:JW0508+JW0507:b4572; hypothetical protein
RA 556,858 A→T 100% L36Q (CTG→CAG)  folD ← bifunctional 5,10‑methylene‑tetrahydrofolate dehydrogenase and 5,10‑methylene‑tetrahydrofolate cyclohydrolase
MC JC 654,214 Δ1,199 bp 100% [dcuC][dcuC] [dcuC], insH, [dcuC]
RA 1,093,686 T→C 100% V130A (GTA→GCA)  ycdT → predicted diguanylate cyclase
MC JC 1,301,041 Δ1,336 bp 100% insCinsD insC, insD
RA 1,303,631 C→T 100% Q246Y (CAG→TAT)  oppA → oligopeptide transporter subunit
RA 1,303,633 G→T 100% Q246Y (CAG→TAT oppA → oligopeptide transporter subunit
MC JC 2,176,983 Δ1,199 bp 100% [gatA][gatA] [gatA], insH, [gatA]
MC JC 2,320,795 Δ1,336 bp 100% insD[rcsC] insD, insC, [rcsC]
RA 2,809,581 G→A 66.7% G52G (GGG→GGA mprA → DNA‑binding transcriptional regulator
RA 2,866,110 T→A 100% intergenic (‑20/+160) rpoS ← / ← nlpD RNA polymerase, sigma S (sigma 38) factor/predicted outer membrane lipoprotein
RA 3,249,944 C→A 8.8% P89T (CCG→ACG)  yqjG → predicted S‑transferase
MC JC 3,742,938 Δ1,336 bp 100% insDinsC insD, insC
MC JC 3,746,913 Δ1,199 bp 100% [tnaB][tnaB] [tnaB], insH, [tnaB]
RA 3,824,280 T→C 100% pseudogene (413/687 nt) rph → ECK3633:JW3618:b3643; ribonuclease PH
RA 3,824,326 T→C 100% pseudogene (459/687 nt) rph → ECK3633:JW3618:b3643; ribonuclease PH
RA 3,824,362 A→T 100% pseudogene (495/687 nt) rph → ECK3633:JW3618:b3643; ribonuclease PH
RA 3,824,368 G→A 100% pseudogene (501/687 nt) rph → ECK3633:JW3618:b3643; ribonuclease PH
RA 3,824,374 C→T 100% pseudogene (507/687 nt) rph → ECK3633:JW3618:b3643; ribonuclease PH
RA 3,824,380 2 bp→AA 100% pseudogene (513‑514/687 nt) rph → ECK3633:JW3618:b3643; ribonuclease PH
RA 3,824,410 T→A 100% pseudogene (543/687 nt) rph → ECK3633:JW3618:b3643; ribonuclease PH
RA 3,824,514 T→C 100% pseudogene (647/687 nt) rph → ECK3633:JW3618:b3643; ribonuclease PH
RA 3,824,524 T→G 100% pseudogene (657/687 nt) rph → ECK3633:JW3618:b3643; ribonuclease PH
RA 3,824,533:1 +G 100% pseudogene (666/687 nt) rph → ECK3633:JW3618:b3643; ribonuclease PH
MC JC 4,310,365 Δ1,199 bp 100% insH[alsK] insH, [alsK]
RA 4,371,273 Δ2 bp 100% intergenic (‑6/+299) dcuA ← / ← aspA C4‑dicarboxylate antiporter/aspartate ammonia‑lyase
RA 4,563,589 T→A 100% K207M (AAG→ATG)  iadA ← isoaspartyl dipeptidase

Unassigned missing coverage evidence
   seq id start end size ←reads reads→ gene description
* * ÷ W3110S.gb 511879 511894 16 16 [8] [7] 11 ybaT predicted transporter
* * ÷ W3110S.gb 573864–575008 575999 992–2136 10 [9] [0] 145 insH–[nmpC] insH,[nmpC]
* * ÷ W3110S.gb 721516 721537 22 10 [8] [8] 11 kdpE DNA‑binding response regulator in two‑component regulatory system with KdpD
* * ÷ W3110S.gb 2296383 2296401 19 10 [6] [7] 11 ccmF heme lyase, CcmF subunit
* * ÷ W3110S.gb 4213632–4216719 4216719 1–3088 10 [9] [9] 12 rrlE–[rrfE] rrlE,[rrfE]
* * ÷ W3110S.gb 4329957 4329971 15 12 [8] [9] 10 phnC/phnB phosphonate/organophosphate ester transporter subunit/conserved hypothetical protein

Unassigned new junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? W3110S.gb 70061 =48 (0.890)74 (1.760)
+CGCGTACAACT
37/78 NT 59.9% intergenic (‑13/‑326) araB/araC L‑ribulokinase/DNA‑binding transcriptional dual regulator
?W3110S.gb 236044 = 80 (1.480)intergenic (‑42/‑23) rnhA/dnaQ ribonuclease HI, degrades RNA of DNA‑RNA hybrids/DNA polymerase III epsilon subunit
* ? W3110S.gb = 7126570 (1.290)41 (1.060)
+AGGGTCCCTAAGTC
23/72 NT 57.2% coding (879/879 nt) araC DNA‑binding transcriptional dual regulator
?W3110S.gb 3464879 = 16 (0.300)intergenic (‑255/+46) murB/rrfB UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding/5S ribosomal RNA
* ? W3110S.gb = 23679843 (0.790)43 (1.770)
+27 bp
27/46 NT 67.3% coding (732/732 nt) dnaQ DNA polymerase III epsilon subunit
?W3110S.gb 4124504 = 50 (0.920)intergenic (+106/‑1) damX/dam hypothetical protein/DNA adenine methylase
* ? W3110S.gb = 274373NA (NA)75 (1.380) 30/100 NT 100% intergenic (‑32/‑176) insH/mmuP IS5 transposase and trans‑activator/predicted S‑methylmethionine transporter
?W3110S.gb 576000 = 0 (0.000)pseudogene (109/1128 nt) nmpC DLP12 prophage region; ECK0544:JW5078:b0553; truncated outer membrane porin
* ? W3110S.gb 314453 =49 (0.900)43 (0.830)
+TC
28/96 NT 64.7% coding (873/888 nt) eaeH attaching and effacing protein, pathogenesis factor
?W3110S.gb 3932930 = 0 (0.000)coding (219/1608 nt) dppA dipeptide transporter
* ? W3110S.gb = 315707NA (NA)67 (1.240) 40/100 NT 100% coding (654/687 nt) ykgA predicted DNA‑binding transcriptional regulator
?W3110S.gb 576911 = 0 (0.000)coding (76/498 nt) ybcS predicted lysozyme
* ? W3110S.gb 566000 =NA (NA)56 (1.030) 35/100 NT 100% intergenic (+90/‑56) ybcD/insE DLP12 prophage region; ECK0532:JW0527:b4508; predicted replication protein fragment/IS3 element protein InsE
?W3110S.gb = 576913 0 (0.000)coding (78/498 nt) ybcS predicted lysozyme
* ? W3110S.gb 713380 =51 (0.940)73 (1.470) 40/92 NT 61.4% intergenic (‑156/‑29) ybfF/seqA conserved hypothetical protein/regulatory protein for replication initiation
?W3110S.gb = 4125340 45 (0.900)coding (836/837 nt) dam DNA adenine methylase
* ? W3110S.gb = 71395474 (1.370)89 (1.640) 47/100 NT 55.1% coding (546/546 nt) seqA regulatory protein for replication initiation
?W3110S.gb 2809398 = 71 (1.310)intergenic (+63/‑28) ygaH/mprA predicted inner membrane protein/DNA‑binding transcriptional regulator
* ? W3110S.gb = 80909932 (0.590)4 (0.100)
+GAATGTGAATAA
4/76 NT 13.1% coding (581/1290 nt) bioA 7,8‑diaminopelargonic acid synthase, PLP‑dependent
?W3110S.gb 2239936 = 38 (0.700)intergenic (+102/+141) yeiA/mglC predicted oxidoreductase/methyl‑galactoside transporter subunit
* ? W3110S.gb 1209367 =24 (0.440)20 (0.550) 17/68 NT 57.8% coding (290/630 nt) ycfK hypothetical protein
?W3110S.gb 1211196 = 13 (0.360)pseudogene (37/537 nt) stfE ECK1143:JW5172:b1157; e14 prophage region; predicted side tail fiber protein fragment
* ? W3110S.gb = 120938221 (0.390)9 (0.250) 7/68 NT 39.9% coding (305/630 nt) ycfK hypothetical protein
?W3110S.gb = 1211179 13 (0.360)pseudogene (54/537 nt) stfE ECK1143:JW5172:b1157; e14 prophage region; predicted side tail fiber protein fragment
* ? W3110S.gb 1439508 =33 (0.610)4 (0.080) 3/92 NT 8.7% coding (2991/3525 nt) ydbK fused predicted Fe‑S subunit of pyruvate‑flavodoxin oxidoreductase
?W3110S.gb 1439539 = 54 (1.080)coding (2960/3525 nt) ydbK fused predicted Fe‑S subunit of pyruvate‑flavodoxin oxidoreductase
* ? W3110S.gb 2018571 =6 (0.110)33 (0.650) 23/94 NT 85.4% coding (567/687 nt) fliH flagellar biosynthesis protein
?W3110S.gb 2519723 = NA (NA)intergenic (+33/‑54) nupC/insL nucleoside (except guanosine) transporter/predicted transposase
* ? W3110S.gb = 20185834 (0.070)41 (0.790) 31/96 NT 91.4% coding (579/687 nt) fliH flagellar biosynthesis protein
?W3110S.gb = 2521059 NA (NA)intergenic (+170/+30) insL/yfeA predicted transposase/predicted diguanylate cyclase
* ? W3110S.gb 2172308 =NA (NA)69 (1.300) 43/98 NT 100% pseudogene (1563/2041 nt) gatR DNA‑binding transcriptional regulator; ECK2083:JW5340+JW2074:b4498
?W3110S.gb = 3920281 0 (0.000)coding (88/213 nt) cspA major cold shock protein
* ? W3110S.gb = 2173564NA (NA)63 (1.190) 32/98 NT 100% pseudogene (307/2041 nt)
pseudogene (307/339 nt)
gatR
gatR
DNA‑binding transcriptional regulator; ECK2083:JW5340+JW2074:b4498
DNA‑binding transcriptional regulator, C‑ter fragment; ECK2083:JW2074:b2090
?W3110S.gb 3920278 = 0 (0.000)coding (91/213 nt) cspA major cold shock protein
* ? W3110S.gb = 2173565NA (NA)63 (1.160) 39/100 NT 100% pseudogene (306/2041 nt)
pseudogene (306/339 nt)
gatR
gatR
DNA‑binding transcriptional regulator; ECK2083:JW5340+JW2074:b4498
DNA‑binding transcriptional regulator, C‑ter fragment; ECK2083:JW2074:b2090
?W3110S.gb = 3932933 0 (0.000)coding (222/1608 nt) dppA dipeptide transporter
* ? W3110S.gb = 2412225NA (NA)97 (1.790) 49/100 NT 87.4% intergenic (‑55/‑782) insA/yfbQ IS1 element protein/predicted aminotransferase
?W3110S.gb = 3748807 14 (0.260)coding (176/984 nt) tnaB tryptophan transporter of low affinity
* ? W3110S.gb 2564137 =11 (0.200)13 (0.280)
+ATGCGGT
7/86 NT 32.2% coding (1404/1404 nt) eutA reactivating factor for ethanolamine ammonia lyase
?W3110S.gb = 3997968 53 (0.980)coding (684/753 nt) yhiQ predicted SAM‑dependent methyltransferase
* ? W3110S.gb = 32611330 (0.000)49 (0.980) 33/92 NT 100% intergenic (‑59/‑69) tdcE/insH pyruvate formate‑lyase 4/2‑ketobutyrate formate‑lyase/IS5 element protein
?W3110S.gb 3262333 = 0 (0.000)coding (1183/1296 nt) tdcD propionate kinase/acetate kinase C, anaerobic
* ? W3110S.gb 3748799 =3 (0.060)70 (1.290) 36/100 NT 95.9% coding (184/984 nt) tnaB tryptophan transporter of low affinity
?W3110S.gb 4056222 = NA (NA)intergenic (‑564/+15) yrhB/insB hypothetical protein/IS1 transposase InsAB'
* ? W3110S.gb = 37505920 (0.000)80 (1.610) 43/92 NT 100% intergenic (‑104/‑69) tnaA/insH tryptophanase/L‑cysteine desulfhydrase, PLP‑dependent/IS5 element protein
?W3110S.gb 3751792 = 0 (0.000)intergenic (+151/+116) insH/tnaC IS5 element protein/tryptophanase leader peptide
* ? W3110S.gb = 387404980 (1.480)4 (0.080) 3/94 NT 5.4% coding (32/843 nt) yibA lyase containing HEAT‑repeat
?W3110S.gb = 4016519 65 (1.280)pseudogene (12/246 nt) yrhC ECK3468:JW5678:b4552; hypothetical protein fragment
* ? W3110S.gb 3921545 =73 (1.350)3 (0.060) 3/88 NT 5.4% coding (173/711 nt) yiaF conserved hypothetical protein
?W3110S.gb 3921580 = 41 (0.860)coding (208/711 nt) yiaF conserved hypothetical protein
* ? W3110S.gb 4547364 =35 (0.650)13 (0.290) 9/82 NT 27.5% intergenic (+49/‑433) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin
?W3110S.gb 4547678 = 40 (0.900)intergenic (+363/‑119) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin
* ? W3110S.gb = 454737232 (0.590)8 (0.180) 7/82 NT 19.5% intergenic (+57/‑425) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin
?W3110S.gb = 4547668 40 (0.900)intergenic (+353/‑129) fimE/fimA tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin