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 | 121,985 | Δ1 bp | 6.4% | intergenic (‑434/‑107) | aroP ← / → pdhR | aromatic amino acid transporter/pyruvate dehydrogenase complex repressor; autorepressor |
MC JC | 257,908 | Δ776 bp | 100% | [crl] | [crl] | |
JC | 417,368 | Δ5 bp | 100% | coding (227‑231/690 nt) | phoB → | response regulator in two‑component regulatory system with PhoR |
RA | 701,378:1 | +T | 100% | coding (217/1221 nt) | nagC ← | N‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor |
RA | 704,008 | C→A | 100% | P22Q (CCG→CAG) | nagE → | N‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components |
RA | 720,627 | A→T | 7.8% | L15F (TTA→TTT) | ybfK → | uncharacterized protein |
RA | 1,130,125 | T→A | 5.7% | S2C (AGT→TGT) | flgM ← | anti‑sigma factor for FliA (sigma 28) |
RA | 1,157,896:1 | +A | 100% | coding (28/1434 nt) | ptsG → | fused glucose‑specific PTS enzymes: IIB component/IIC component |
RA | 1,269,450 | C→T | 7.1% | intergenic (+61/+253) | rdlA → / ← ldrB | sRNA antisense regulator affects LdrA translation; proposed addiction module in LDR‑A repeat, with toxic peptide LdrA/toxic polypeptide, small |
RA | 1,269,451 | G→T | 6.4% | intergenic (+62/+252) | rdlA → / ← ldrB | sRNA antisense regulator affects LdrA translation; proposed addiction module in LDR‑A repeat, with toxic peptide LdrA/toxic polypeptide, small |
RA | 1,466,258 | G→C | 6.5% | pseudogene (867/2513 nt) | ydbA → | pseudogene, autotransporter homolog; interrupted by IS2 and IS30 |
RA | 1,637,431 | A→G | 14.7% | K129E (AAA→GAA) | ydfO → | Qin prophage; uncharacterized protein |
RA | 1,708,844 | C→G | 7.2% | A642A (GCC→GCG) | rsxC → | SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex |
RA | 1,726,169 | T→G | 100% | F49L (TTT→TTG) | nemR → | transcriptional repressor for the nemRA‑gloA operon, quinone‑, glyoxal‑, and HOCl‑activated |
JC JC | 1,755,755 | IS5 (+) +4 bp | 100% | coding (58‑61/1413 nt) | pykF → | pyruvate kinase I |
MC JC | 1,978,503 | Δ776 bp | 100% | insB1–insA | insB1, insA | |
JC JC | 1,979,486 | IS5 (+) +4 bp | 100% | intergenic (‑271/‑264) | insA ← / → uspC | IS1 repressor TnpA/universal stress protein |
RA | 2,132,787 | A→C | 100% | I204S (ATC→AGC) | wcaA ← | putative glycosyl transferase |
RA | 2,173,361 | Δ2 bp | 100% | pseudogene (1‑2/442 nt) | gatC ← | pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC |
RA | 2,304,537 | A→G | 8.0% | intergenic (+144/+571) | eco → / ← mqo | ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain |
RA | 2,304,549 | C→T | 8.9% | intergenic (+156/+559) | eco → / ← mqo | ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain |
RA | 2,567,332 | Δ1 bp | 5.6% | coding (776/1227 nt) | eutH ← | ethanolamine transporter |
RA | 3,560,455:1 | +G | 100% | intergenic (‑2/+1) | glpR ← / ← glpR | pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon |
MC JC | 3,815,859 | Δ82 bp | 100% | [rph]–[rph] | [rph], [rph] | |
MC JC | 3,909,312 | Δ774 bp | 100% | [pstA]–[pstC] | [pstA], [pstC] | |
MC JC | 4,001,645 | Δ5 bp | 100% | coding (220‑224/951 nt) | corA → | magnesium/nickel/cobalt transporter |
RA | 4,171,648 | A→C | 100% | noncoding (12/120 nt) | rrfB → | 5S ribosomal RNA of rrnB operon |
RA | 4,184,543 | C→A | 100% | P1100Q (CCG→CAG) | rpoB → | RNA polymerase, beta subunit |
RA | 4,296,060 | C→T | 9.5% | intergenic (+266/+376) | gltP → / ← yjcO | glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein |
RA | 4,296,190 | A→G | 26.5% | intergenic (+396/+246) | gltP → / ← yjcO | glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein |
RA | 4,296,191 | A→C | 28.3% | intergenic (+397/+245) | gltP → / ← yjcO | glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein |
RA | 4,296,268 | T→C | 21.5% | intergenic (+474/+168) | gltP → / ← yjcO | glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein |
RA | 4,296,286 | C→T | 13.6% | intergenic (+492/+150) | gltP → / ← yjcO | glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein |
RA | 4,296,304 | A→C | 6.5% | intergenic (+510/+132) | gltP → / ← yjcO | glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein |
RA | 4,296,380:1 | +CG | 100% | intergenic (+586/+56) | gltP → / ← yjcO | glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein |
RA | 4,326,088 | T→G | 6.6% | intergenic (‑347/+311) | yjdN ← / ← yjdM | metalloprotein superfamily protein/zinc‑ribbon family protein |
RA | 4,326,147 | T→C | 9.8% | intergenic (‑406/+252) | yjdN ← / ← yjdM | metalloprotein superfamily protein/zinc‑ribbon family protein |
RA | 4,326,148 | G→A | 6.5% | intergenic (‑407/+251) | yjdN ← / ← yjdM | metalloprotein superfamily protein/zinc‑ribbon family protein |
Unassigned missing coverage evidence | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | start | end | size | ←reads | reads→ | gene | description | |||
* | * | ÷ | NC_000913 | 1293040 | 1303173 | 10134 | 87 [0] | [8] 102 | tdk–[oppB] | tdk,insZ,insZ,insZ,adhE,ychE,oppA,[oppB] |
* | * | ÷ | NC_000913 | 3423885–3424233 | 3424401–3424239 | 7–517 | 33 [32] | [30] 33 | rrlD | 23S ribosomal RNA of rrnD operon |
* | * | ÷ | NC_000913 | 4110668 | 4111588 | 921 | 33 [32] | [32] 46 | [cdh]–tpiA | [cdh],tpiA |
Unassigned new junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | NC_000913 | 257908 = | NA (NA) | 79 (0.830) | 59/514 | NT | 100% | noncoding (768/768 nt) | IS1 | repeat region |
? | NC_000913 | = 1293039 | 0 (0.000) | intergenic (‑117/‑488) | hns/tdk | global DNA‑binding transcriptional dual regulator H‑NS/thymidine kinase/deoxyuridine kinase | |||||
* | ? | NC_000913 | = 258675 | NA (NA) | 92 (0.970) | 66/512 | NT | 92.9% | noncoding (1/768 nt) | IS1 | repeat region |
? | NC_000913 | 1303174 = | 7 (0.070) | coding (276/921 nt) | oppB | oligopeptide ABC transporter permease | |||||
* | ? | NC_000913 | 1207790 = | 43 (0.450) | 41 (0.460) | 37/482 | NT | 50.5% | coding (290/630 nt) | stfP | e14 prophage; uncharacterized protein |
? | NC_000913 | 1209619 = | 40 (0.450) | pseudogene (1/501 nt) | stfE | pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related | |||||
* | ? | NC_000913 | = 1207805 | 40 (0.420) | 59 (0.660) | 44/482 | NT | 60.4% | coding (305/630 nt) | stfP | e14 prophage; uncharacterized protein |
? | NC_000913 | = 1209602 | 40 (0.450) | pseudogene (18/501 nt) | stfE | pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related | |||||
* | ? | NC_000913 | 1411899 = | 96 (1.000) | 6 (0.070) | 6/462 | NT | 8.2% | coding (50/936 nt) | ttcA | tRNA s(2)C32 thioltransferase, iron sulfur cluster protein |
? | NC_000913 | = 1434958 | 48 (0.560) | pseudogene (51/51 nt) | ttcC | pseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related | |||||
* | ? | NC_000913 | 4316637 = | 48 (0.500) | 4 (0.040) | 4/490 | NT | 7.1% | coding (582/1137 nt) | phnM | ribophosphonate triphosphate hydrolase |
? | NC_000913 | 4316823 = | 59 (0.650) | coding (396/1137 nt) | phnM | ribophosphonate triphosphate hydrolase |