breseq version 0.33.1 revision 8505477f25b3
mutation predictions | marginal predictions | summary statistics | genome diff | command line log |
Marginal read alignment evidence (highest frequency 20 of 209 shown, sorted by frequency from high to low) | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 417,533 | 0 | T | G | 54.0% | 8.0 / 30.7 | 38 | V131G (GTG→GGG) | phoB | response regulator in two‑component regulatory system with PhoR |
* | NC_000913 | 2,386,862 | 0 | G | A | 38.6% | 71.4 / 52.6 | 70 | intergenic (+31/‑72) | yfbL/yfbM | putative M28A family peptidase/DUF1877 family protein |
* | NC_000913 | 811,721 | 0 | T | G | 37.2% | 227.2 / 26.8 | 122 | V67G (GTG→GGG) | bioC | malonyl‑ACP O‑methyltransferase, SAM‑dependent |
* | NC_000913 | 3,862,664 | 0 | C | G | 37.0% | 23.0 / 13.2 | 50 | pseudogene (940/1617 nt) | glvC | pseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component |
* | NC_000913 | 2,684,779 | 0 | A | C | 34.7% | 65.4 / 11.1 | 50 | G243G (GGT→GGG) | glyA | serine hydroxymethyltransferase |
* | NC_000913 | 75,023 | 0 | A | C | 34.2% | 127.6 / 17.2 | 73 | V153G (GTG→GGG) | thiB | thiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein |
* | NC_000913 | 3,745,019 | 0 | C | G | 33.3% | 38.0 / 11.5 | 42 | G72G (GGC→GGG) | yiaN | 2,3‑diketo‑L‑gulonate TRAP transporter large permease protein |
* | NC_000913 | 540,558 | 0 | T | G | 33.0% | 192.0 / 37.7 | 100 | intergenic (+50/‑7) | allB/ybbY | allantoinase/putative uracil/xanthine transporter |
* | NC_000913 | 1,265,473 | 0 | T | G | 33.0% | 209.3 / 20.6 | 94 | G154G (GGT→GGG) | prfA | peptide chain release factor RF‑1 |
* | NC_000913 | 433,649 | 0 | T | G | 32.5% | 72.1 / 19.3 | 83 | G65G (GGT→GGG) | ribD | fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase |
* | NC_000913 | 354,288 | 0 | T | G | 31.6% | 132.3 / 16.6 | 57 | V528G (GTG→GGG) | prpE | propionate‑‑CoA ligase |
* | NC_000913 | 3,782,908 | 0 | T | G | 31.2% | 80.5 / 11.9 | 74 | T252P (ACC→CCC) | gpsA | glycerol‑3‑phosphate dehydrogenase (NAD+) |
* | NC_000913 | 2,509,880 | 0 | A | C | 31.2% | 302.5 / 29.0 | 128 | D84A (GAC→GCC) | yfeO | putative ion channel protein |
* | NC_000913 | 2,483,261 | 0 | C | A | 30.7% | 169.4 / 16.2 | 76 | A27S (GCG→TCG) | emrK | multidrug resistance efflux pump membrane fusion protein |
* | NC_000913 | 1,383,722 | 0 | G | A | 30.6% | 124.0 / 19.2 | 63 | F80F (TTC→TTT) | ycjW | LacI family putative transcriptional repressor |
* | NC_000913 | 753,125 | 0 | T | A | 30.6% | 213.2 / 46.2 | 98 | intergenic (‑330/+60) | ybgD/gltA | putative fimbrial‑like adhesin protein/citrate synthase |
* | NC_000913 | 3,760,384 | 0 | A | C | 30.3% | 120.3 / 19.8 | 76 | V289G (GTG→GGG) | selA | selenocysteine synthase |
* | NC_000913 | 2,878,947 | 0 | C | T | 30.2% | 212.3 / 25.6 | 106 | E276K (GAA→AAA) | ygbT | multifunctional endonuclease Cas1, CRISPR adaptation protein; DNA repair enzyme |
* | NC_000913 | 1,346,853 | 0 | G | T | 30.1% | 192.4 / 25.5 | 93 | V20L (GTA→TTA) | ymiB | uncharacterized protein |
* | NC_000913 | 29,496 | 0 | T | A | 30.0% | 107.3 / 20.4 | 50 | intergenic (+301/‑155) | dapB/carA | dihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase |
Marginal new junction evidence (lowest skew 10 of 30 shown) | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | NC_000913 | = 2713829 | 120 (0.890) | 10 (0.080) | 8/326 | 8.1 | 8.2% | coding (934/1335 nt) | srmB | ATP‑dependent RNA helicase |
? | NC_000913 | = 2713908 | 110 (0.870) | coding (1013/1335 nt) | srmB | ATP‑dependent RNA helicase | |||||
* | ? | NC_000913 | = 3698975 | 131 (0.970) | 8 (0.060) | 6/330 | 9.0 | 6.2% | coding (762/1680 nt) | bcsG | DUF3260 family cellulose production inner membrane protein |
? | NC_000913 | = 3699049 | 117 (0.910) | coding (836/1680 nt) | bcsG | DUF3260 family cellulose production inner membrane protein | |||||
* | ? | NC_000913 | 4542682 = | 72 (0.530) | 5 (0.040) | 5/330 | 9.6 | 4.9% | intergenic (+49/‑433) | fimE/fimA | tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin) |
? | NC_000913 | 4542996 = | 124 (0.960) | intergenic (+363/‑119) | fimE/fimA | tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin) | |||||
* | ? | NC_000913 | = 4498381 | 10 (0.070) | 16 (0.120) | 5/342 | 9.8 | 62.0% | noncoding (201/1331 nt) | IS2 | repeat region |
? | NC_000913 | = 4498773 | NA (NA) | noncoding (593/1331 nt) | IS2 | repeat region | |||||
* | ? | NC_000913 | = 4068629 | 183 (1.350) | 4 (0.030) | 4/324 | 10.0 | 2.6% | coding (648/2037 nt) | yihQ | putative sulpholipid alpha‑glucosidase; alpha‑glucosyl fluoride glucosidase |
? | NC_000913 | = 4068723 | 132 (1.050) | coding (554/2037 nt) | yihQ | putative sulpholipid alpha‑glucosidase; alpha‑glucosyl fluoride glucosidase | |||||
* | ? | NC_000913 | = 4095482 | 187 (1.380) | 4 (0.030) | 4/328 | 10.1 | 2.5% | coding (501/1260 nt) | rhaA | L‑rhamnose isomerase |
? | NC_000913 | = 4095710 | 132 (1.030) | coding (273/1260 nt) | rhaA | L‑rhamnose isomerase | |||||
* | ? | NC_000913 | = 3235030 | 95 (0.700) | 5 (0.040) | 4/338 | 10.2 | 4.6% | coding (846/1137 nt) | rlmG | 23S rRNA m(2)G1835 methyltransferase, SAM‑dependent |
? | NC_000913 | = 3235119 | 117 (0.890) | coding (757/1137 nt) | rlmG | 23S rRNA m(2)G1835 methyltransferase, SAM‑dependent | |||||
* | ? | NC_000913 | = 2283632 | 112 (0.830) | 4 (0.030) | 4/336 | 10.2 | 3.5% | coding (316/1008 nt) | yejK | spermidine nucleoid‑associated protein |
? | NC_000913 | = 2283824 | 110 (0.840) | coding (124/1008 nt) | yejK | spermidine nucleoid‑associated protein | |||||
* | ? | NC_000913 | 1060875 = | 130 (0.960) | 4 (0.030) | 4/336 | 10.2 | 2.8% | coding (1620/2547 nt) | torA | trimethylamine N‑oxide (TMAO) reductase I, catalytic subunit |
? | NC_000913 | 1061210 = | 149 (1.140) | coding (1955/2547 nt) | torA | trimethylamine N‑oxide (TMAO) reductase I, catalytic subunit | |||||
* | ? | NC_000913 | 992281 = | 129 (0.950) | 4 (0.030) | 4/336 | 10.2 | 3.3% | coding (1660/2613 nt) | pepN | aminopeptidase N |
? | NC_000913 | 992374 = | 109 (0.830) | coding (1753/2613 nt) | pepN | aminopeptidase N |