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 | 2,154 | C→T | 47.6% | D606D (GAC→GAT) | thrA → | fused aspartokinase I and homoserine dehydrogenase I |
RA | 70,434 | C→T | 40.5% | N16N (AAC→AAT) | araC → | DNA‑binding transcriptional dual regulator |
RA | 70,581 | C→T | 55.9% | V65V (GTC→GTT) | araC → | DNA‑binding transcriptional dual regulator |
RA | 70,740 | T→G | 56.7% | G118G (GGT→GGG) | araC → | DNA‑binding transcriptional dual regulator |
RA | 71,079 | C→T | 55.8% | R231R (CGC→CGT) | araC → | DNA‑binding transcriptional dual regulator |
RA | 71,082 | T→C | 56.9% | I232I (ATT→ATC) | araC → | DNA‑binding transcriptional dual regulator |
RA | 71,085 | T→C | 56.9% | S233S (AGT→AGC) | araC → | DNA‑binding transcriptional dual regulator |
RA | 71,109 | T→C | 50.0% | T241T (ACT→ACC) | araC → | DNA‑binding transcriptional dual regulator |
RA | 71,175 | A→G | 48.9% | R263R (CGA→CGG) | araC → | DNA‑binding transcriptional dual regulator |
RA | 71,214 | T→C | 26.3% | 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 | 53.7% | D12A (GAT→GCT) | dnaQ → | DNA polymerase III epsilon subunit |
RA | 236,107 | A→C | 52.7% | 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 | 810,113 | G→A | 26.1% | M116I (ATG→ATA) | bioB → | biotin synthase |
RA | 907,713 | G→A | 35.7% | T331I (ACC→ATC) | ybjT ← | conserved hypothetical protein with NAD(P)‑binding Rossmann‑fold domain |
RA | 1,093,686 | T→C | 100% | V130A (GTA→GCA) | ycdT → | predicted diguanylate cyclase |
RA | 1,151,925 | A→C | 30.0% | P207P (CCA→CCC) | fabD → | malonyl‑CoA‑[acyl‑carrier‑protein] transacylase |
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 |
RA | 2,058,034 | C→T | 37.2% | G279G (GGC→GGT) | amn → | AMP nucleosidase |
MC JC | 2,176,983 | Δ1,199 bp | 100% | [gatA]–[gatA] | [gatA], insH, [gatA] | |
RA | 2,286,501 | C→A | 48.2% | A261S (GCG→TCG) | yejK ← | nucleotide associated protein |
MC JC | 2,320,795 | Δ1,336 bp | 100% | insD–[rcsC] | insD, insC, [rcsC] | |
RA | 2,809,581 | G→A | 55.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 |
RA | 3,053,471 | A→G | 21.9% | F9S (TTT→TCT) | pepP ← | proline aminopeptidase P II |
RA | 3,152,234 | G→A | 39.5% | D6N (GAT→AAT) | yghB → | conserved inner membrane protein |
MC JC | 3,261,134 | Δ1,199 bp | 100% | insH–[tdcD] | insH, [tdcD] | |
RA | 3,315,139 | T→A | 29.9% | D244V (GAT→GTT) | infB ← | fused protein chain initiation factor 2, IF2 |
RA | 3,688,001 | A→G | 41.5% | S440G (AGT→GGT) | yifB → | predicted bifunctional protein, enzyme and transcriptional regulator |
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 |
RA | 4,195,173 | A→G | 25.0% | *118* (TAA→TGA) | rplR → | 50S ribosomal subunit protein L18 |
RA | 4,246,990 | C→T | 26.7% | R40H (CGT→CAT) | malG ← | maltose transporter subunit |
MC JC | 4,310,365 | Δ1,199 bp | 100% | insH–[alsK] | insH, [alsK] | |
RA | 4,340,671 | G→A | 33.4% | T347T (ACC→ACT) | adiC ← | arginine:agmatin |
RA | 4,371,273 | Δ2 bp | 100% | intergenic (‑6/+299) | dcuA ← / ← aspA | C4‑dicarboxylate antiporter/aspartate ammonia‑lyase |
RA | 4,416,062 | A→G | 25.8% | N27D (AAT→GAT) | yjfK → | conserved hypothetical protein |
RA | 4,434,349 | T→A | 21.4% | intergenic (+112/‑197) | fklB → / → cycA | FKBP‑type peptidyl‑prolyl cis‑trans isomerase/D‑alanine/D‑serine/glycine transporter |
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 | 171656 | 171695 | 40 | 4 [3] | [3] 4 | fhuB | fused subunits of iron‑hydroxamate transporter and membrane components of ABC superfamily transporter |
* | * | ÷ | W3110S.gb | 283087 | 283146 | 60 | 4 [3] | [3] 4 | yagF | predicted dehydratase |
* | * | ÷ | W3110S.gb | 379235 | 379257 | 23 | 8 [2] | [3] 4 | frmR/yaiO | regulator protein that represses frmRAB operon/hypothetical protein |
* | * | ÷ | W3110S.gb | 573861–575008 | 575999 | 992–2139 | 4 [3] | [0] 68 | insH–[nmpC] | insH,[nmpC] |
* | * | ÷ | W3110S.gb | 721519 | 721529 | 11 | 4 [3] | [2] 4 | kdpE | DNA‑binding response regulator in two‑component regulatory system with KdpD |
* | * | ÷ | W3110S.gb | 1352905 | 1352934–1352917 | 13–30 | 4 [3] | [3] 4 | fabI/ycjD | enoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/conserved hypothetical protein |
* | * | ÷ | W3110S.gb | 1421309 | 1421330 | 22 | 4 [3] | [2] 4 | ydaG/racR | hypothetical protein/predicted DNA‑binding transcriptional regulator |
* | * | ÷ | W3110S.gb | 1928924 | 1928961 | 38 | 4 [3] | [2] 4 | ptrB | protease II |
* | * | ÷ | W3110S.gb | 2876819 | 2876872 | 54 | 4 [3] | [3] 4 | iap/ygbF | aminopeptidase in alkaline phosphatase isozyme conversion/hypothetical protein |
* | * | ÷ | W3110S.gb | 4213633–4215889 | 4216709–4215891 | 3–3077 | 4 [3] | [3] 4 | rrlE–[rrfE] | rrlE,[rrfE] |
* | * | ÷ | W3110S.gb | 4401141 | 4401188 | 48 | 5 [1] | [3] 4 | 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 | = 7063 | 16 (0.590) | 12 (0.480) +GATT |
6/90 | NT | 54.3% | coding (897/1431 nt) | yaaJ | predicted transporter |
? | W3110S.gb | = 2564144 | 6 (0.220) | coding (1397/1404 nt) | eutA | reactivating factor for ethanolamine ammonia lyase | |||||
* | ? | W3110S.gb | 70061 = | 32 (1.180) | 17 (0.810) +CGCGTACAACT |
14/76 | NT | 35.0% | intergenic (‑13/‑326) | araB/araC | L‑ribulokinase/DNA‑binding transcriptional dual regulator |
? | W3110S.gb | 236044 = | 50 (1.840) | intergenic (‑42/‑23) | rnhA/dnaQ | ribonuclease HI, degrades RNA of DNA‑RNA hybrids/DNA polymerase III epsilon subunit | |||||
* | ? | W3110S.gb | = 71265 | 45 (1.650) | 22 (1.140) +AGGGTCCCTAAGTC |
16/70 | NT | 54.0% | coding (879/879 nt) | araC | DNA‑binding transcriptional dual regulator |
? | W3110S.gb | 3464879 = | 8 (0.290) | intergenic (‑255/+46) | murB/rrfB | UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding/5S ribosomal RNA | |||||
* | ? | W3110S.gb | = 236798 | 21 (0.770) | 12 (0.450) | 9/96 | NT | 36.1% | coding (732/732 nt) | dnaQ | DNA polymerase III epsilon subunit |
? | W3110S.gb | 2855720 = | 22 (0.830) | intergenic (‑258/‑29) | ygbA/mutS | hypothetical protein/methyl‑directed mismatch repair protein | |||||
* | ? | W3110S.gb | = 274373 | NA (NA) | 30 (1.100) | 19/98 | 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 = | 25 (0.920) | 25 (0.960) +TC |
22/94 | NT | 67.6% | 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 | = 315707 | NA (NA) | 50 (1.840) | 33/98 | 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) | 27 (0.990) | 18/98 | 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 = | 24 (0.880) | 14 (0.560) | 11/90 | NT | 38.9% | intergenic (‑156/‑29) | ybfF/seqA | conserved hypothetical protein/regulatory protein for replication initiation |
? | W3110S.gb | = 4125340 | 22 (0.880) | coding (836/837 nt) | dam | DNA adenine methylase | |||||
* | ? | W3110S.gb | = 713954 | 24 (0.880) | 30 (1.100) | 21/98 | NT | 48.4% | coding (546/546 nt) | seqA | regulatory protein for replication initiation |
? | W3110S.gb | 2809398 = | 40 (1.470) | intergenic (+63/‑28) | ygaH/mprA | predicted inner membrane protein/DNA‑binding transcriptional regulator | |||||
* | ? | W3110S.gb | 1209367 = | 13 (0.480) | 8 (0.440) | 7/66 | NT | 50.5% | coding (290/630 nt) | ycfK | hypothetical protein |
? | W3110S.gb | 1211196 = | 7 (0.390) | pseudogene (37/537 nt) | stfE | ECK1143:JW5172:b1157; e14 prophage region; predicted side tail fiber protein fragment | |||||
* | ? | W3110S.gb | = 1209382 | 15 (0.550) | 8 (0.440) | 5/66 | NT | 48.5% | coding (305/630 nt) | ycfK | hypothetical protein |
? | W3110S.gb | = 1211179 | 7 (0.390) | pseudogene (54/537 nt) | stfE | ECK1143:JW5172:b1157; e14 prophage region; predicted side tail fiber protein fragment | |||||
* | ? | W3110S.gb | 2018571 = | 3 (0.110) | 20 (0.780) | 16/92 | NT | 87.7% | 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 | 0 (0.000) | 26 (1.000) | 18/94 | NT | 100% | 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) | 30 (1.130) | 23/96 | 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) | 42 (1.580) | 28/96 | 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) | 21 (0.770) | 17/98 | 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 | = 2200411 | 46 (1.690) | 5 (0.230) +AGCCACTTCC |
3/78 | NT | 13.2% | pseudogene (604/825 nt) pseudogene (604/3796 nt) |
molR molR |
DNA‑binding transcriptional regulator, N‑ter fragment; ECK2108:JW2102:b2115 DNA‑binding transcriptional regulator; ECK2108:JW2102+JW5915+JW5916:b4499 |
? | W3110S.gb | 2615024 = | 37 (1.360) | coding (275/1464 nt) | yfgC | predicted peptidase | |||||
* | ? | W3110S.gb | = 2412225 | NA (NA) | 38 (1.400) | 28/98 | NT | 86.4% | intergenic (‑55/‑782) | insA/yfbQ | IS1 element protein/predicted aminotransferase |
? | W3110S.gb | = 3748807 | 6 (0.220) | coding (176/984 nt) | tnaB | tryptophan transporter of low affinity | |||||
* | ? | W3110S.gb | 3748799 = | 4 (0.150) | 33 (1.210) | 24/98 | NT | 89.2% | 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' |