breseq version 0.33.1 revision 8505477f25b3
mutation predictions | marginal predictions | summary statistics | genome diff | command line log |
Marginal read alignment evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,226,550 | 0 | C | A | 65.0% | 8.8 / 20.3 | 20 | A908E (GCG→GAG) | metH | homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent |
* | NC_000913 | 4,296,060 | 0 | C | T | 30.9% | 61.2 / 31.2 | 42 | intergenic (+266/+376) | gltP/yjcO | glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein |
* | NC_000913 | 2,801,056 | 0 | A | T | 24.0% | 43.5 / 10.9 | 25 | N31I (AAT→ATT) | nrdI | NrdEF cluster assembly flavodoxin |
* | NC_000913 | 4,512,350 | 1 | . | G | 21.7% | 58.0 / 19.1 | 23 | coding (65/957 nt) | fecD | ferric citrate ABC transporter permease |
* | NC_000913 | 4,512,344 | 0 | C | . | 20.8% | 63.3 / 12.7 | 24 | coding (71/957 nt) | fecD | ferric citrate ABC transporter permease |
* | NC_000913 | 2,587,131 | 0 | G | C | 20.7% | 62.3 / 13.8 | 29 | Q467H (CAG→CAC) | narQ | sensory histidine kinase in two‑component regulatory system with NarP |
* | NC_000913 | 2,587,132 | 0 | G | C | 20.7% | 65.0 / 13.6 | 29 | V468L (GTG→CTG) | narQ | sensory histidine kinase in two‑component regulatory system with NarP |
* | NC_000913 | 4,498,223 | 0 | A | T | 20.7% | 59.9 / 11.5 | 29 | intergenic (+283/‑49) | intB/insC1 | pseudogene, integrase homology;IS, phage, Tn; Phage‑related functions and prophages; KpLE2 phage‑like element; P4‑like integrase/IS2 repressor TnpA |
* | NC_000913 | 1,932,095 | 0 | T | A | 20.0% | 75.0 / 14.3 | 36 | intergenic (+36/+20) | purT/eda | phosphoribosylglycinamide formyltransferase 2/KHG/KDPG aldolase; 2‑dehydro‑3‑deoxy‑phosphogluconate/4‑hydroxy‑2‑ oxoglutarate aldolase |
* | NC_000913 | 2,561,205 | 0 | A | T | 20.0% | 86.9 / 14.6 | 40 | intergenic (+307/‑163) | yffL/yffM | CPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein |
* | NC_000913 | 2,561,206 | 0 | A | T | 20.0% | 85.1 / 14.8 | 40 | intergenic (+308/‑162) | yffL/yffM | CPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein |
* | NC_000913 | 2,906,608 | 0 | G | C | 20.0% | 72.4 / 12.3 | 30 | intergenic (+25/+35) | ygcG/eno | TPM domain protein, putative phosphatase/enolase |
* | NC_000913 | 2,906,609 | 0 | G | C | 20.0% | 69.9 / 14.5 | 30 | intergenic (+26/+34) | ygcG/eno | TPM domain protein, putative phosphatase/enolase |
* | NC_000913 | 2,906,610 | 0 | A | G | 20.0% | 60.3 / 10.5 | 30 | intergenic (+27/+33) | ygcG/eno | TPM domain protein, putative phosphatase/enolase |
* | NC_000913 | 3,101,541 | 0 | G | C | 20.0% | 73.8 / 10.6 | 30 | P28R (CCG→CGG) | yggN | DUF2884 family putative periplasmic protein |
Marginal new junction evidence (lowest skew 10 of 70 shown) | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | NC_000913 | = 3798990 | 179 (4.740) | 7 (0.200) | 7/238 | 2.5 | 3.8% | coding (323/1074 nt) | waaU | lipopolysaccharide core biosynthesis |
? | NC_000913 | = 3798995 | 190 (5.360) | coding (318/1074 nt) | waaU | lipopolysaccharide core biosynthesis | |||||
* | ? | NC_000913 | = 3722005 | 127 (3.360) | 7 (0.200) | 7/238 | 2.5 | 5.4% | noncoding (1373/1443 nt) | IS150 | repeat region |
? | NC_000913 | = 3722024 | 128 (3.610) | noncoding (1392/1443 nt) | IS150 | repeat region | |||||
* | ? | NC_000913 | 4335640 = | 48 (1.270) | 6 (0.180) | 6/226 | 2.6 | 13.1% | intergenic (‑50/+54) | eptA/adiC | lipid A phosphoethanolamine transferase/arginine:agmatine antiporter |
? | NC_000913 | 4335670 = | 37 (1.100) | intergenic (‑80/+24) | eptA/adiC | lipid A phosphoethanolamine transferase/arginine:agmatine antiporter | |||||
* | ? | NC_000913 | 1411899 = | 61 (1.610) | 5 (0.170) | 5/202 | 2.6 | 10.1% | coding (50/936 nt) | ttcA | tRNA s(2)C32 thioltransferase, iron sulfur cluster protein |
? | NC_000913 | = 1434958 | 41 (1.360) | pseudogene (51/51 nt) | ttcC | pseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related | |||||
* | ? | NC_000913 | 546854 = | 39 (1.030) | 6 (0.170) | 6/240 | 2.7 | 13.8% | coding (175/1668 nt) | fdrA | putative NAD(P)‑binding acyl‑CoA synthetase |
? | NC_000913 | 546880 = | 38 (1.060) | coding (201/1668 nt) | fdrA | putative NAD(P)‑binding acyl‑CoA synthetase | |||||
* | ? | NC_000913 | = 4623082 | 33 (0.870) | 6 (0.160) | 6/246 | 2.8 | 15.6% | coding (1314/1332 nt) | yjjJ | putative protein kinase |
? | NC_000913 | = 4623098 | 33 (0.900) | coding (1330/1332 nt) | yjjJ | putative protein kinase | |||||
* | ? | NC_000913 | 3844903 = | 109 (2.880) | 6 (0.170) | 6/244 | 2.8 | 5.5% | coding (940/1767 nt) | adeD | cryptic adenine deaminase |
? | NC_000913 | 3844942 = | 102 (2.810) | coding (979/1767 nt) | adeD | cryptic adenine deaminase | |||||
* | ? | NC_000913 | 3805341 = | 154 (4.070) | 6 (0.160) | 6/250 | 2.8 | 3.9% | coding (610/798 nt) | waaP | kinase that phosphorylates core heptose of lipopolysaccharide |
? | NC_000913 | 3805362 = | 144 (3.870) | coding (589/798 nt) | waaP | kinase that phosphorylates core heptose of lipopolysaccharide | |||||
* | ? | NC_000913 | 382237 = | NA (NA) | 5 (0.140) | 5/242 | 3.0 | NA | noncoding (978/1331 nt) | IS2 | repeat region |
? | NC_000913 | 382267 = | NA (NA) | noncoding (1008/1331 nt) | IS2 | repeat region | |||||
* | ? | NC_000913 | 4285004 = | 31 (0.820) | 5 (0.140) | 5/244 | 3.1 | 13.7% | coding (210/315 nt) | yjcH | DUF485 family inner membrane protein |
? | NC_000913 | 4285020 = | 33 (0.910) | coding (194/315 nt) | yjcH | DUF485 family inner membrane protein |