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 | 77.4% | N16N (AAC→AAT) | araC → | DNA‑binding transcriptional dual regulator |
RA | 70,581 | C→T | 75.0% | V65V (GTC→GTT) | araC → | DNA‑binding transcriptional dual regulator |
RA | 70,740 | T→G | 85.3% | G118G (GGT→GGG) | araC → | DNA‑binding transcriptional dual regulator |
RA | 71,079 | C→T | 81.4% | R231R (CGC→CGT) | araC → | DNA‑binding transcriptional dual regulator |
RA | 71,082 | T→C | 81.0% | I232I (ATT→ATC) | araC → | DNA‑binding transcriptional dual regulator |
RA | 71,085 | T→C | 80.8% | S233S (AGT→AGC) | araC → | DNA‑binding transcriptional dual regulator |
RA | 71,109 | T→C | 78.0% | T241T (ACT→ACC) | araC → | DNA‑binding transcriptional dual regulator |
RA | 71,175 | A→G | 82.4% | R263R (CGA→CGG) | araC → | DNA‑binding transcriptional dual regulator |
RA | 71,214 | T→C | 80.8% | F276F (TTT→TTC) | araC → | DNA‑binding transcriptional dual regulator |
RA | 155,463 | A→T | 100% | *247K (TAA→AAA) | ecpD ← | predicted periplasmic pilin chaperone |
RA | 219,687 | A→T | 14.3% | A103A (GCT→GCA) | rcsF ← | predicted outer membrane protein, signal |
RA | 236,101 | A→C | 74.8% | D12A (GAT→GCT) | dnaQ → | DNA polymerase III epsilon subunit |
RA | 236,107 | A→C | 78.0% | 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 |
RA | 1,184,032 | G→A | 16.7% | Y125Y (TAC→TAT) | potD ← | polyamine transporter subunit |
RA | 1,184,053 | A→C | 16.1% | P118P (CCT→CCG) | potD ← | polyamine transporter subunit |
MC JC | 1,301,041 | Δ1,336 bp | 100% | insC–insD | 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 | 79.8% | 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 |
MC JC | 3,261,134 | Δ1,199 bp | 100% | insH–[tdcD] | insH, [tdcD] | |
MC JC | 3,742,938 | Δ1,336 bp | 100% | insD–insC | insD, insC | |
MC JC | 3,746,913 | Δ1,199 bp | 100% | [tnaB]–[tnaB] | [tnaB], insH, [tnaB] | |
MC JC | 3,750,593 | Δ1,199 bp | 100% | insH | insH | |
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 | 573859–575008 | 575999 | 992–2141 | 5 [3] | [0] 44 | insH–[nmpC] | insH,[nmpC] |
* | * | ÷ | W3110S.gb | 907991 | 908010 | 20 | 4 [3] | [3] 4 | ybjT | conserved hypothetical protein with NAD(P)‑binding Rossmann‑fold domain |
* | * | ÷ | W3110S.gb | 1173190 | 1173208 | 19 | 4 [3] | [3] 4 | mfd | transcription‑repair coupling factor |
* | * | ÷ | W3110S.gb | 1328313 | 1328336 | 24 | 4 [3] | [3] 4 | yciQ | predicted inner membrane protein |
* | * | ÷ | W3110S.gb | 1828084 | 1828124 | 41 | 4 [2] | [3] 5 | astE | succinylglutamate desuccinylase |
* | * | ÷ | W3110S.gb | 2041778 | 2041800 | 23 | 5 [3] | [3] 4 | yedY | predicted reductase |
* | * | ÷ | W3110S.gb | 2569112 | 2569133 | 22 | 5 [3] | [3] 4 | eutE | predicted aldehyde dehydrogenase, ethanolamine utilization protein |
* | * | ÷ | W3110S.gb | 2777796 | 2777828 | 33 | 5 [3] | [3] 4 | ypjA | adhesin‑like autotransporter |
* | * | ÷ | W3110S.gb | 2827379 | 2827422 | 44 | 4 [3] | [3] 7 | gutM | DNA‑binding transcriptional activator |
* | * | ÷ | W3110S.gb | 3044667 | 3044716 | 50 | 4 [3] | [3] 4 | ygfF/gcvP | predicted NAD(P)‑binding oxidoreductase with NAD(P)‑binding Rossmann‑fold domain/glycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex |
* | * | ÷ | W3110S.gb | 3799043 | 3799077 | 35 | 4 [3] | [3] 4 | yicM | predicted transporter |
* | * | ÷ | W3110S.gb | 3859252 | 3859284 | 33 | 4 [3] | [3] 4 | yibK/lldD | predicted rRNA methylase/L‑lactate dehydrogenase, FMN‑linked |
* | * | ÷ | W3110S.gb | 4135552 | 4135608 | 57 | 5 [3] | [2] 4 | yhfS/frlR | conserved hypothetical protein/predicted DNA‑binding transcriptional regulator |
* | * | ÷ | W3110S.gb | 4213633–4215889 | 4216706–4215891 | 3–3074 | 4 [3] | [2] 4 | rrlE–[rrfE] | rrlE,[rrfE] |
* | * | ÷ | W3110S.gb | 4401124 | 4401179 | 56 | 5 [3] | [3] 5 | amiB | N‑acetylmuramoyl‑l‑alanine amidase II |
Unassigned new junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | W3110S.gb | 70061 = | 20 (0.880) | 83 (4.650) +CGCGTACAACT |
50/80 | NT | 80.9% | intergenic (‑13/‑326) | araB/araC | L‑ribulokinase/DNA‑binding transcriptional dual regulator |
? | W3110S.gb | 236044 = | 30 (1.320) | intergenic (‑42/‑23) | rnhA/dnaQ | ribonuclease HI, degrades RNA of DNA‑RNA hybrids/DNA polymerase III epsilon subunit | |||||
* | ? | W3110S.gb | = 71265 | 19 (0.830) | 32 (1.940) +AGGGTCCCTAAGTC |
23/74 | NT | 77.2% | coding (879/879 nt) | araC | DNA‑binding transcriptional dual regulator |
? | W3110S.gb | 3464879 = | 7 (0.310) | intergenic (‑255/+46) | murB/rrfB | UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding/5S ribosomal RNA | |||||
* | ? | W3110S.gb | 97398 = | 38 (1.670) | 3 (0.140) | 3/98 | NT | 9.3% | coding (312/1317 nt) | murD | UDP‑N‑acetylmuramoyl‑L‑alanine:D‑glutamate ligase |
? | W3110S.gb | 97426 = | 22 (1.010) | coding (340/1317 nt) | murD | UDP‑N‑acetylmuramoyl‑L‑alanine:D‑glutamate ligase | |||||
* | ? | W3110S.gb | 225154 = | NA (NA) | 3 (0.150) | 3/92 | NT | NA | noncoding (1384/1542 nt) | rrsH | 16S ribosomal RNA |
? | W3110S.gb | 225176 = | NA (NA) | noncoding (1406/1542 nt) | rrsH | 16S ribosomal RNA | |||||
* | ? | W3110S.gb | = 225915 | NA (NA) | 8 (0.420) | 5/86 | NT | NA | noncoding (157/2904 nt) | rrlH | 23S ribosomal RNA |
? | W3110S.gb | = 225926 | NA (NA) | noncoding (168/2904 nt) | rrlH | 23S ribosomal RNA | |||||
* | ? | W3110S.gb | = 236798 | 23 (1.010) | 43 (4.010) +27 bp |
27/48 | NT | 78.2% | coding (732/732 nt) | dnaQ | DNA polymerase III epsilon subunit |
? | W3110S.gb | 4124504 = | 28 (1.230) | intergenic (+106/‑1) | damX/dam | hypothetical protein/DNA adenine methylase | |||||
* | ? | W3110S.gb | 244046 = | 17 (0.750) | 3 (0.130) | 3/100 | NT | 13.1% | coding (504/579 nt) | lpcA | D‑sedoheptulose 7‑phosphate isomerase |
? | W3110S.gb | 244098 = | 23 (1.030) | coding (556/579 nt) | lpcA | D‑sedoheptulose 7‑phosphate isomerase | |||||
* | ? | W3110S.gb | = 274373 | NA (NA) | 18 (0.790) | 15/102 | 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 = | 13 (0.570) | 14 (0.640) +TC |
13/98 | NT | 67.5% | coding (873/888 nt) | eaeH | attaching and effacing protein, pathogenesis factor |
? | W3110S.gb | 3932930 = | 1 (0.040) | coding (219/1608 nt) | dppA | dipeptide transporter | |||||
* | ? | W3110S.gb | 315232 = | NA (NA) | 3 (0.150) | 3/92 | NT | NA | coding (422/867 nt) | insF | IS3 element protein InsF |
? | W3110S.gb | 315255 = | NA (NA) | coding (445/867 nt) | insF | IS3 element protein InsF | |||||
* | ? | W3110S.gb | = 315707 | NA (NA) | 24 (1.050) | 19/102 | 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 | 468217 = | 21 (0.920) | 4 (0.200) | 3/90 | NT | 18.0% | coding (123/1773 nt) | mdlA | fused predicted multidrug transporter subunits and ATP‑binding components of ABC superfamily |
? | W3110S.gb | 468244 = | 18 (0.900) | coding (150/1773 nt) | mdlA | fused predicted multidrug transporter subunits and ATP‑binding components of ABC superfamily | |||||
* | ? | W3110S.gb | 566000 = | NA (NA) | 33 (1.450) | 24/102 | 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 = | 20 (0.880) | 69 (3.290) | 47/94 | NT | 81.4% | intergenic (‑156/‑29) | ybfF/seqA | conserved hypothetical protein/regulatory protein for replication initiation |
? | W3110S.gb | = 4125340 | 13 (0.620) | coding (836/837 nt) | dam | DNA adenine methylase | |||||
* | ? | W3110S.gb | = 713954 | 34 (1.490) | 86 (3.770) | 51/102 | NT | 75.4% | coding (546/546 nt) | seqA | regulatory protein for replication initiation |
? | W3110S.gb | 2809398 = | 22 (0.970) | intergenic (+63/‑28) | ygaH/mprA | predicted inner membrane protein/DNA‑binding transcriptional regulator | |||||
* | ? | W3110S.gb | 894129 = | 13 (0.570) | 3 (0.150) | 3/88 | NT | 20.5% | intergenic (+274/‑77) | ybjN/potF | predicted oxidoreductase/putrescine transporter subunit |
? | W3110S.gb | 894150 = | 12 (0.610) | intergenic (+295/‑56) | ybjN/potF | predicted oxidoreductase/putrescine transporter subunit | |||||
* | ? | W3110S.gb | 1209367 = | 6 (0.260) | 10 (0.640) | 9/70 | NT | 58.6% | coding (290/630 nt) | ycfK | hypothetical protein |
? | W3110S.gb | 1211196 = | 10 (0.640) | pseudogene (37/537 nt) | stfE | ECK1143:JW5172:b1157; e14 prophage region; predicted side tail fiber protein fragment | |||||
* | ? | W3110S.gb | = 1209382 | 7 (0.310) | 9 (0.580) | 9/70 | NT | 54.9% | coding (305/630 nt) | ycfK | hypothetical protein |
? | W3110S.gb | = 1211179 | 10 (0.640) | pseudogene (54/537 nt) | stfE | ECK1143:JW5172:b1157; e14 prophage region; predicted side tail fiber protein fragment | |||||
* | ? | W3110S.gb | = 1286086 | 20 (0.880) | 3 (0.150) | 3/90 | NT | 13.4% | coding (906/1539 nt) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
? | W3110S.gb | = 1286104 | 21 (1.040) | coding (924/1539 nt) | narH | nitrate reductase 1, beta (Fe‑S) subunit | |||||
* | ? | W3110S.gb | = 1311245 | 17 (0.750) | 3 (0.150) | 3/88 | NT | 18.4% | coding (739/1254 nt) | kch | voltage‑gated potassium channel |
? | W3110S.gb | = 1311272 | 12 (0.610) | coding (712/1254 nt) | kch | voltage‑gated potassium channel | |||||
* | ? | W3110S.gb | 1324223 = | 14 (0.610) | 4 (0.200) | 3/88 | NT | 27.5% | coding (438/1563 nt) | trpE | component I of anthranilate synthase |
? | W3110S.gb | 1324247 = | 9 (0.460) | coding (414/1563 nt) | trpE | component I of anthranilate synthase | |||||
* | ? | W3110S.gb | = 1499413 | 31 (1.360) | 3 (0.150) | 3/92 | NT | 9.5% | coding (814/1626 nt) | mdoD | glucan biosynthesis protein, periplasmic |
? | W3110S.gb | = 1499426 | 29 (1.410) | coding (827/1626 nt) | mdoD | glucan biosynthesis protein, periplasmic | |||||
* | ? | W3110S.gb | 1883832 = | 24 (1.050) | 3 (0.150) | 3/92 | NT | 15.5% | coding (207/1086 nt) | yeaU | predicted dehydrogenase |
? | W3110S.gb | 1883858 = | 11 (0.540) | coding (233/1086 nt) | yeaU | predicted dehydrogenase | |||||
* | ? | W3110S.gb | 2018571 = | 1 (0.040) | 16 (0.750) | 15/96 | NT | 94.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 | = 2018583 | 2 (0.090) | 23 (1.050) | 16/98 | NT | 92.3% | 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) | 25 (1.120) | 20/100 | 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 | = 2173564 | NA (NA) | 17 (0.760) | 12/100 | 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 | = 2173565 | NA (NA) | 29 (1.270) | 23/102 | NT | 96.7% | 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 | 1 (0.040) | coding (222/1608 nt) | dppA | dipeptide transporter | |||||
* | ? | W3110S.gb | = 2304640 | 16 (0.700) | 3 (0.150) | 3/92 | NT | 15.2% | coding (1448/2487 nt) | napA | nitrate reductase, periplasmic, large subunit |
? | W3110S.gb | = 2304649 | 19 (0.920) | coding (1439/2487 nt) | napA | nitrate reductase, periplasmic, large subunit | |||||
* | ? | W3110S.gb | = 2412225 | NA (NA) | 31 (1.360) | 25/102 | NT | 93.9% | intergenic (‑55/‑782) | insA/yfbQ | IS1 element protein/predicted aminotransferase |
? | W3110S.gb | = 3748807 | 2 (0.090) | coding (176/984 nt) | tnaB | tryptophan transporter of low affinity | |||||
* | ? | W3110S.gb | = 3098949 | 30 (1.320) | 3 (0.150) | 3/92 | NT | 10.0% | coding (436/1047 nt) | ansB | periplasmic L‑asparaginase II |
? | W3110S.gb | = 3098965 | 27 (1.310) | coding (420/1047 nt) | ansB | periplasmic L‑asparaginase II | |||||
* | ? | W3110S.gb | 3265071 = | 35 (1.540) | 3 (0.150) | 3/92 | NT | 9.3% | coding (813/990 nt) | tdcB | catabolic threonine dehydratase, PLP‑dependent |
? | W3110S.gb | 3265090 = | 27 (1.310) | coding (794/990 nt) | tdcB | catabolic threonine dehydratase, PLP‑dependent | |||||
* | ? | W3110S.gb | 3748799 = | 1 (0.040) | 24 (1.050) | 21/102 | NT | 96.0% | 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 | 4212589 = | NA (NA) | 3 (0.140) | 3/94 | NT | NA | noncoding (934/1542 nt) | rrsE | 16S ribosomal RNA |
? | W3110S.gb | 4212609 = | NA (NA) | noncoding (954/1542 nt) | rrsE | 16S ribosomal RNA | |||||
* | ? | W3110S.gb | 4292798 = | 25 (1.100) | 3 (0.150) | 3/88 | NT | 13.8% | intergenic (+6/‑39) | nrfA/nrfB | nitrite reductase, formate‑dependent, cytochrome/nitrite reductase, formate‑dependent, penta‑heme cytochrome c |
? | W3110S.gb | 4292820 = | 16 (0.810) | intergenic (+28/‑17) | nrfA/nrfB | nitrite reductase, formate‑dependent, cytochrome/nitrite reductase, formate‑dependent, penta‑heme cytochrome c | |||||
* | ? | W3110S.gb | 4480518 = | 28 (1.230) | 4 (0.190) | 3/96 | NT | 11.7% | coding (400/1815 nt) | yjgL | hypothetical protein |
? | W3110S.gb | 4480541 = | 34 (1.590) | coding (423/1815 nt) | yjgL | hypothetical protein | |||||
* | ? | W3110S.gb | 4547364 = | 14 (0.610) | 4 (0.210) | 4/84 | NT | 21.9% | intergenic (+49/‑433) | fimE/fimA | tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin |
? | W3110S.gb | 4547678 = | 17 (0.910) | intergenic (+363/‑119) | fimE/fimA | tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin |