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 1418 shown) | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,386,862 | 0 | G | A | 66.1% | 92.2 / 112.4 | 121 | intergenic (+31/‑72) | yfbL/yfbM | putative M28A family peptidase/DUF1877 family protein |
* | NC_000913 | 357,890 | 0 | C | A | 61.2% | 4.0 / 48.1 | 50 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
* | NC_000913 | 2,332,823 | 0 | A | G | 56.7% | 54.3 / 107.6 | 93 | pseudogene (1580/4605 nt) | yfaS | pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein |
* | NC_000913 | 354,288 | 0 | T | G | 56.2% | 69.4 / 133.8 | 91 | V528G (GTG→GGG) | prpE | propionate‑‑CoA ligase |
* | NC_000913 | 3,133,302 | 0 | T | G | 54.4% | 129.8 / 191.6 | 125 | H219P (CAC→CCC) | yghS | putative ATP‑binding protein |
* | NC_000913 | 4,604,344 | 0 | C | T | 54.4% | 23.5 / 147.6 | 125 | S62F (TCT→TTT) | bglJ | bgl operon transcriptional activator |
* | NC_000913 | 3,665,702 | 0 | C | A | 54.4% | 67.8 / 201.2 | 149 | A37S (GCG→TCG) | gadX | acid resistance regulon transcriptional activator; autoactivator |
* | NC_000913 | 3,665,698 | 0 | T | C | 52.4% | 72.6 / 199.7 | 145 | D38G (GAT→GGT) | gadX | acid resistance regulon transcriptional activator; autoactivator |
* | NC_000913 | 993,917 | 0 | A | C | 51.9% | 54.3 / 81.1 | 52 | V43G (GTG→GGG) | ssuB | aliphatic sulfonate ABC transporter ATPase |
* | NC_000913 | 1,222,233 | 0 | G | A | 51.5% | 38.5 / 91.0 | 68 | pseudogene (60/75 nt) | ycgH | pseudogene; putative ATP‑binding component of a transport system |
* | NC_000913 | 1,554,000 | 0 | T | G | 51.2% | 81.9 / 105.5 | 87 | Y557S (TAC→TCC) | maeA | malate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme |
* | NC_000913 | 1,356,201 | 0 | T | G | 50.7% | 99.2 / 123.4 | 70 | T304P (ACC→CCC) | sapA | antimicrobial peptide transport ABC transporter periplasmic binding protein |
* | NC_000913 | 29,496 | 0 | T | A | 50.0% | 74.2 / 72.2 | 86 | intergenic (+301/‑155) | dapB/carA | dihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase |
* | NC_000913 | 2,386,868 | 0 | T | A | 49.6% | 78.3 / 86.9 | 120 | intergenic (+37/‑66) | yfbL/yfbM | putative M28A family peptidase/DUF1877 family protein |
* | NC_000913 | 2,571,036 | 0 | T | G | 49.2% | 75.0 / 22.1 | 67 | N239T (AAC→ACC) | eutE | aldehyde oxidoreductase, ethanolamine utilization protein |
* | NC_000913 | 2,386,863 | 0 | C | A | 49.1% | 105.2 / 65.9 | 121 | intergenic (+32/‑71) | yfbL/yfbM | putative M28A family peptidase/DUF1877 family protein |
* | NC_000913 | 2,386,883 | 0 | C | T | 48.9% | 74.5 / 120.0 | 131 | intergenic (+52/‑51) | yfbL/yfbM | putative M28A family peptidase/DUF1877 family protein |
* | NC_000913 | 4,622,250 | 0 | A | G | 48.0% | 105.6 / 57.1 | 103 | E161G (GAG→GGG) | yjjJ | putative protein kinase |
* | NC_000913 | 650,458 | 0 | A | G | 47.7% | 114.1 / 24.0 | 90 | T11T (ACT→ACC) | citE | citrate lyase, citryl‑ACP lyase (beta) subunit |
* | NC_000913 | 2,903,715 | 0 | G | A | 47.4% | 36.9 / 84.8 | 78 | intergenic (+341/+1032) | ygcE/queE | putative kinase/7‑carboxy‑7‑deazaguanine synthase; queosine biosynthesis |
Marginal new junction evidence (lowest skew 10 of 217 shown) | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | NC_000913 | = 1207805 | 79 (0.700) | 17 (0.170) | 14/266 | 5.0 | 19.8% | coding (305/630 nt) | stfP | e14 prophage; uncharacterized protein |
? | NC_000913 | = 1209602 | 67 (0.670) | pseudogene (18/501 nt) | stfE | pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related | |||||
* | ? | NC_000913 | = 381689 | NA (NA) | 18 (0.170) | 15/280 | 5.0 | NA | noncoding (430/1331 nt) | IS2 | repeat region |
? | NC_000913 | = 381743 | NA (NA) | noncoding (484/1331 nt) | IS2 | repeat region | |||||
* | ? | NC_000913 | 4409174 = | NA (NA) | 10 (0.100) | 10/268 | 6.4 | 6.9% | noncoding (22/98 nt) | RIP331 (repetitive extragenic palindromic) element; contains 2 REP sequences and 1 IHF site | RIP331 (repetitive extragenic palindromic) element; contains 2 REP sequences and 1 IHF site |
? | NC_000913 | 4409235 = | 135 (1.340) | noncoding (83/98 nt) | RIP331 (repetitive extragenic palindromic) element; contains 2 REP sequences and 1 IHF site | RIP331 (repetitive extragenic palindromic) element; contains 2 REP sequences and 1 IHF site | |||||
* | ? | NC_000913 | 1207790 = | 73 (0.650) | 11 (0.110) | 10/266 | 6.4 | 14.1% | coding (290/630 nt) | stfP | e14 prophage; uncharacterized protein |
? | NC_000913 | 1209619 = | 69 (0.690) | pseudogene (1/501 nt) | stfE | pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related | |||||
* | ? | NC_000913 | = 3880655 | 124 (1.110) | 13 (0.120) | 10/278 | 6.6 | 10.0% | coding (567/1074 nt) | recF | gap repair protein |
? | NC_000913 | = 3880717 | 118 (1.130) | coding (505/1074 nt) | recF | gap repair protein | |||||
* | ? | NC_000913 | 4005815 = | 122 (1.090) | 9 (0.080) | 9/282 | 7.1 | 7.3% | intergenic (+84/‑49) | pldA/recQ | outer membrane phospholipase A/ATP‑dependent DNA helicase |
? | NC_000913 | 4005837 = | 114 (1.070) | intergenic (+106/‑27) | pldA/recQ | outer membrane phospholipase A/ATP‑dependent DNA helicase | |||||
* | ? | NC_000913 | = 3529405 | 133 (1.190) | 11 (0.110) | 9/278 | 7.1 | 7.6% | coding (58/402 nt) | hslR | ribosome‑associated heat shock protein Hsp15 |
? | NC_000913 | = 3529461 | 142 (1.360) | coding (114/402 nt) | hslR | ribosome‑associated heat shock protein Hsp15 | |||||
* | ? | NC_000913 | = 373065 | 120 (1.070) | 10 (0.100) | 8/262 | 7.1 | 8.6% | coding (145/951 nt) | mhpF | acetaldehyde‑CoA dehydrogenase II, NAD‑binding |
? | NC_000913 | = 373047 | 107 (1.090) | coding (127/951 nt) | mhpF | acetaldehyde‑CoA dehydrogenase II, NAD‑binding | |||||
* | ? | NC_000913 | 326549 = | 115 (1.030) | 10 (0.100) | 9/278 | 7.1 | 8.3% | coding (699/1671 nt) | betA | choline dehydrogenase, a flavoprotein |
? | NC_000913 | 326652 = | 115 (1.100) | coding (596/1671 nt) | betA | choline dehydrogenase, a flavoprotein | |||||
* | ? | NC_000913 | = 4062810 | 119 (1.060) | 9 (0.090) | 8/280 | 7.5 | 7.4% | coding (564/1266 nt) | yihN | MFS transporter family protein |
? | NC_000913 | = 4062816 | 114 (1.080) | coding (570/1266 nt) | yihN | MFS transporter family protein |