breseq version 0.27.1 revision 87c22d663cc3
mutation predictions | marginal predictions | summary statistics | genome diff | command line log |
Marginal read alignment evidence (highest frequency 20 of 56 shown) | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,151,759 | 0 | G | A | 47.9% | 49.4 / inf | 383 | E49K (GAA→AAA) | acpP | acyl carrier protein (ACP) |
* | NC_000913 | 4,187,667 | 0 | T | A | 47.8% | 73.2 / inf | 471 | F773Y (TTC→TAC) | rpoC | RNA polymerase, beta prime subunit |
* | NC_000913 | 4,186,425 | 0 | C | T | 46.7% | 125.7 / inf | 486 | P359L (CCA→CTA) | rpoC | RNA polymerase, beta prime subunit |
* | NC_000913 | 2,610,107 | 0 | G | A | 45.8% | 136.9 / inf | 417 | E541K (GAG→AAG) | hyfG | hydrogenase 4, subunit |
* | NC_000913 | 1,907,327 | 0 | G | A | 45.6% | 135.0 / inf | 453 | Q37* (CAG→TAG) | cspC | stress protein, member of the CspA‑family |
* | NC_000913 | 4,213,051 | 0 | A | C | 44.4% | 30.1 / 207.3 | 126 | noncoding (12/120 nt) | rrfE | 5S ribosomal RNA of rrnE operon |
* | NC_000913 | 3,473,612 | 0 | A | C | 43.6% | 116.6 / inf | 496 | L157* (TTA→TGA) | rpsG | 30S ribosomal subunit protein S7 |
* | NC_000913 | 2,947,591 | 0 | A | . | 40.8% | 202.8 / 159.2 | 98 | intergenic (+18/‑16) | metW/metV | tRNA‑Met/tRNA‑Met |
* | NC_000913 | 3,269,994 | 0 | G | T | 36.0% | 188.7 / 41.1 | 196 | intergenic (+392/+631) | yhaC/garK | pentapetide repeats‑related protein/glycerate kinase I |
* | NC_000913 | 1,674,771 | 0 | G | . | 34.1% | 570.6 / 313.0 | 228 | coding (859/1035 nt) | tqsA | pheromone AI‑2 transporter |
* | NC_000913 | 1,674,773 | 0 | G | . | 33.9% | 576.6 / 313.0 | 228 | coding (861/1035 nt) | tqsA | pheromone AI‑2 transporter |
* | NC_000913 | 1,674,772 | 0 | G | . | 33.8% | 576.9 / 313.0 | 228 | coding (860/1035 nt) | tqsA | pheromone AI‑2 transporter |
* | NC_000913 | 1,674,774 | 0 | C | . | 33.3% | 584.3 / 312.0 | 231 | coding (862/1035 nt) | tqsA | pheromone AI‑2 transporter |
* | NC_000913 | 1,674,768 | 0 | C | . | 33.2% | 591.6 / 311.8 | 232 | coding (856/1035 nt) | tqsA | pheromone AI‑2 transporter |
* | NC_000913 | 1,674,770 | 0 | T | . | 33.2% | 554.9 / 302.2 | 232 | coding (858/1035 nt) | tqsA | pheromone AI‑2 transporter |
* | NC_000913 | 1,674,769 | 0 | G | . | 33.1% | 591.9 / 312.1 | 233 | coding (857/1035 nt) | tqsA | pheromone AI‑2 transporter |
* | NC_000913 | 1,674,764 | 0 | G | . | 31.7% | 637.9 / 310.3 | 244 | coding (852/1035 nt) | tqsA | pheromone AI‑2 transporter |
* | NC_000913 | 1,674,765 | 0 | G | . | 31.7% | 631.0 / 310.3 | 244 | coding (853/1035 nt) | tqsA | pheromone AI‑2 transporter |
* | NC_000913 | 1,674,766 | 0 | G | . | 31.6% | 634.2 / 310.3 | 244 | coding (854/1035 nt) | tqsA | pheromone AI‑2 transporter |
* | NC_000913 | 1,674,767 | 0 | G | . | 31.6% | 646.9 / 310.3 | 244 | coding (855/1035 nt) | tqsA | pheromone AI‑2 transporter |
Marginal new junction evidence (lowest skew 10 of 50 shown) | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | NC_000913 | = 1104115 | 191 (0.430) | 203 (0.460) | 147/298 | 3.2 | 51.1% | coding (165/456 nt) | csgB | curlin nucleator protein, minor subunit in curli complex |
? | NC_000913 | 1104420 = | 198 (0.450) | intergenic (+14/‑27) | csgB/csgA | curlin nucleator protein, minor subunit in curli complex/curlin subunit, amyloid curli fibers, cryptic | |||||
* | ? | NC_000913 | = 3366750 | NA (NA) | 222 (0.500) | 135/300 | 4.4 | 42.7% | noncoding (1/1195 nt) | IS5 | repeat region |
? | NC_000913 | 3393816 = | 298 (0.680) | coding (2443/3801 nt) | yhdP | DUF3971‑AsmA2 domains protein | |||||
* | ? | NC_000913 | = 2290113 | NA (NA) | 294 (0.670) | 133/300 | 4.6 | 51.5% | noncoding (1/1195 nt) | IS5 | repeat region |
? | NC_000913 | 2413638 = | 277 (0.630) | coding (169/1203 nt) | ackA | acetate kinase A and propionate kinase 2 | |||||
* | ? | NC_000913 | 2288919 = | NA (NA) | 227 (0.520) | 132/298 | 4.6 | 45.6% | noncoding (1195/1195 nt) | IS5 | repeat region |
? | NC_000913 | = 2413641 | 271 (0.620) | coding (172/1203 nt) | ackA | acetate kinase A and propionate kinase 2 | |||||
* | ? | NC_000913 | = 1674768 | 158 (0.360) | 178 (0.420) | 126/288 | 4.8 | 54.5% | coding (856/1035 nt) | tqsA | pheromone AI‑2 transporter |
? | NC_000913 | 1674781 = | 146 (0.350) | coding (869/1035 nt) | tqsA | pheromone AI‑2 transporter | |||||
* | ? | NC_000913 | 4542682 = | 476 (1.080) | 68 (0.160) | 61/282 | 15.8 | 12.9% | intergenic (+49/‑433) | fimE/fimA | tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin) |
? | NC_000913 | 4542996 = | 469 (1.130) | intergenic (+363/‑119) | fimE/fimA | tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin) | |||||
* | ? | NC_000913 | = 4542690 | 469 (1.060) | 66 (0.160) | 60/282 | 16.1 | 12.7% | intergenic (+57/‑425) | fimE/fimA | tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin) |
? | NC_000913 | = 4542986 | 466 (1.130) | intergenic (+353/‑129) | fimE/fimA | tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin) | |||||
* | ? | NC_000913 | 696863 = | 137 (0.310) | 108 (0.270) +GAAGAAACAATC |
48/276 | 19.0 | 30.0% | intergenic (‑33/+2) | glnW/glnU | tRNA‑Gln/tRNA‑Gln |
? | NC_000913 | = 697134 | 411 (0.930) | intergenic (‑1/+379) | metT/asnB | tRNA‑Met/asparagine synthetase B | |||||
* | ? | NC_000913 | = 2067353 | NA (NA) | 214 (0.500) | 49/292 | 19.9 | 42.5% | noncoding (1/1195 nt) | IS5 | repeat region |
? | NC_000913 | 3393816 = | 298 (0.680) | coding (2443/3801 nt) | yhdP | DUF3971‑AsmA2 domains protein | |||||
* | ? | NC_000913 | = 2067353 | NA (NA) | 172 (0.390) | 48/300 | 20.7 | 38.3% | noncoding (1/1195 nt) | IS5 | repeat region |
? | NC_000913 | 2413638 = | 277 (0.630) | coding (169/1203 nt) | ackA | acetate kinase A and propionate kinase 2 |