mutation predictions | marginal predictions | summary statistics | genome diff | command line log
|
breseq version 0.39.0
mutation predictions | marginal predictions | summary statistics | genome diff | command line log |
| Marginal read alignment evidence (highest frequency 20 of 28 shown, sorted by frequency from high to low) | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,878,947 | 0 | C | T | 40.0% | 25.4 / 12.8 | 30 | E276K (GAA→AAA) | cas1 | multifunctional nuclease Cas1 |
| * | NC_000913 | 3,548,925 | 0 | A | C | 32.1% | 73.2 / 11.1 | 57 | G382G (GGT→GGG) | malQ | 4‑alpha‑glucanotransferase |
| * | NC_000913 | 2,667,126 | 0 | A | C | 31.8% | 97.9 / 10.8 | 66 | Y241D (TAT→GAT) | hcaT | putative 3‑phenylpropionate transporter |
| * | NC_000913 | 1,061,930 | 0 | G | A | 31.0% | 88.3 / 19.1 | 58 | L44L (CTG→CTA) | torD | trimethylamine‑N‑oxide reductase‑specific chaperone |
| * | NC_000913 | 613,283 | 0 | A | G | 27.2% | 125.8 / 18.3 | 84 | I183V (ATT→GTT) | fes | ferric enterobactin esterase |
| * | NC_000913 | 1,715,620 | 0 | A | C | 26.9% | 77.6 / 10.9 | 52 | G90G (GGT→GGG) | pdxY | pyridoxal kinase 2 |
| * | NC_000913 | 300,120 | 0 | C | G | 26.8% | 94.6 / 12.7 | 57 | R272P (CGC→CCC) | paoC | aldehyde dehydrogenase: molybdenum cofactor‑binding subunit |
| * | NC_000913 | 2,950,173 | 0 | A | G | 24.5% | 101.0 / 11.0 | 53 | E311G (GAG→GGG) | argA | N‑acetylglutamate synthase |
| * | NC_000913 | 3,393,679 | 0 | A | T | 24.0% | 96.7 / 10.2 | 50 | S860R (AGT→AGA) | yhdP | phospholipid transport factor YhdP |
| * | NC_000913 | 566,344 | 0 | A | G | 23.1% | 112.7 / 31.8 | 66 | intergenic (‑366/‑340) | intD/renD | putative integrase/protein RenD |
| * | NC_000913 | 566,275 | 0 | A | G | 21.7% | 131.5 / 33.1 | 70 | intergenic (‑297/‑409) | intD/renD | putative integrase/protein RenD |
| * | NC_000913 | 566,283 | 0 | G | C | 21.5% | 152.5 / 44.7 | 67 | intergenic (‑305/‑401) | intD/renD | putative integrase/protein RenD |
| * | NC_000913 | 566,320 | 0 | A | C | 21.5% | 131.0 / 28.1 | 66 | intergenic (‑342/‑364) | intD/renD | putative integrase/protein RenD |
| * | NC_000913 | 640,880 | 0 | T | G | 21.5% | 271.7 / 11.2 | 122 | G376G (GGT→GGG) | ahpF | alkyl hydroperoxide reductase, AhpF component |
| * | NC_000913 | 3,914,523 | 0 | C | T | 21.4% | 138.5 / 11.8 | 71 | V226V (GTG→GTA) | glmU | fused N‑acetylglucosamine‑1‑phosphate uridyltransferase and glucosamine‑1‑phosphate acetyltransferase |
| * | NC_000913 | 566,284 | 0 | C | G | 21.2% | 158.6 / 45.8 | 67 | intergenic (‑306/‑400) | intD/renD | putative integrase/protein RenD |
| * | NC_000913 | 566,350 | 0 | C | G | 21.2% | 155.9 / 39.2 | 67 | intergenic (‑372/‑334) | intD/renD | putative integrase/protein RenD |
| * | NC_000913 | 566,190 | 0 | C | A | 21.0% | 174.0 / 32.9 | 81 | intergenic (‑212/‑494) | intD/renD | putative integrase/protein RenD |
| * | NC_000913 | 566,332 | 0 | T | G | 20.9% | 123.9 / 44.2 | 68 | intergenic (‑354/‑352) | intD/renD | putative integrase/protein RenD |
| * | NC_000913 | 566,347 | 0 | G | A | 20.9% | 134.3 / 28.9 | 67 | intergenic (‑369/‑337) | intD/renD | putative integrase/protein RenD |
| Marginal new junction evidence (sorted from low to high skew) | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
| * | ? | NC_000913 | = 3451295 | 68 (1.190) | 3 (0.050) | 3/470 | 4.5 | 3.8% | coding (70/822 nt) | rplB | 50S ribosomal subunit protein L2 |
| ? | NC_000913 | 3451798 = | 88 (1.610) | coding (489/606 nt) | rplD | 50S ribosomal subunit protein L4 | |||||
| * | ? | NC_000913 | = 2290113 | NA (NA) | 3 (0.050) | 3/490 | 4.6 | 10.7% | noncoding (1/1195 nt) | IS5 | repeat region |
| ? | NC_000913 | 2904758 = | 25 (0.440) | coding (661/672 nt) | queE | putative 7‑carboxy‑7‑deazaguanine synthase QueE | |||||