breseq version 0.35.4 revision f352f80f4bc9
mutation predictions | marginal predictions | summary statistics | genome diff | command line log |
Marginal read alignment evidence (highest frequency 20 of 60 shown, sorted by frequency from high to low) | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,177,234 | 0 | C | A | 59.3% | ‑5.2 / 13.1 | 34 | intergenic (‑30/+278) | gatY/fbaB | tagatose‑1,6‑bisphosphate aldolase 2 subunit GatY/fructose‑bisphosphate aldolase class I |
* | NC_000913 | 729,372 | 0 | A | C | 57.7% | ‑5.5 / 11.4 | 26 | intergenic (+32/‑211) | ybfA/rhsC | DUF2517 domain‑containing protein YbfA/rhs element protein RhsC |
* | NC_000913 | 3,233,728 | 0 | T | G | 57.1% | 0.5 / 10.8 | 14 | *139Y (TAA→TAC) | higA | antitoxin/DNA‑binding transcriptional repressor HigA |
* | NC_000913 | 2,177,245 | 0 | G | T | 56.7% | ‑1.4 / 10.5 | 33 | intergenic (‑41/+267) | gatY/fbaB | tagatose‑1,6‑bisphosphate aldolase 2 subunit GatY/fructose‑bisphosphate aldolase class I |
* | NC_000913 | 2,470,427 | 0 | C | A | 55.0% | 9.2 / 21.1 | 23 | Q433K (CAA→AAA) | yfdI | CPS‑53 (KpLE1) prophage; serotype specific glucosyl transferase |
* | NC_000913 | 454,919 | 0 | G | T | 52.9% | 8.8 / 16.3 | 19 | intergenic (+130/‑214) | bolA/tig | DNA‑binding transcriptional dual regulator BolA/trigger factor |
* | NC_000913 | 708,004 | 0 | C | A | 52.9% | 8.0 / 15.5 | 18 | intergenic (+247/‑330) | glnS/chiP | glutamine‑‑tRNA ligase/chitobiose outer membrane channel |
* | NC_000913 | 454,873 | 0 | C | G | 52.0% | ‑4.8 / 11.7 | 27 | intergenic (+84/‑260) | bolA/tig | DNA‑binding transcriptional dual regulator BolA/trigger factor |
* | NC_000913 | 1,879,855 | 0 | C | A | 51.3% | 34.6 / 53.1 | 44 | D32Y (GAT→TAT) ‡ | yeaR | DUF1971 domain‑containing protein YeaR |
* | NC_000913 | 1,879,846 | 0 | T | G | 48.8% | 13.6 / 25.3 | 43 | T35P (ACG→CCG) | yeaR | DUF1971 domain‑containing protein YeaR |
* | NC_000913 | 1,879,848 | 0 | C | G | 48.8% | 31.0 / 19.5 | 44 | G34A (GGA→GCA) | yeaR | DUF1971 domain‑containing protein YeaR |
* | NC_000913 | 238,535 | 0 | G | T | 48.1% | 26.2 / 13.5 | 54 | H68N (CAT→AAT) | ykfM | uncharacterized protein YkfM |
* | NC_000913 | 2,177,255 | 0 | G | T | 47.4% | 9.7 / 22.7 | 38 | intergenic (‑51/+257) | gatY/fbaB | tagatose‑1,6‑bisphosphate aldolase 2 subunit GatY/fructose‑bisphosphate aldolase class I |
* | NC_000913 | 2,470,411 | 0 | A | T | 46.2% | 17.4 / 11.7 | 26 | L427F (TTA→TTT) | yfdI | CPS‑53 (KpLE1) prophage; serotype specific glucosyl transferase |
* | NC_000913 | 1,879,875 | 0 | C | A | 44.2% | 71.2 / 14.0 | 54 | G25V (GGA→GTA) | yeaR | DUF1971 domain‑containing protein YeaR |
* | NC_000913 | 1,879,853 | 0 | A | T | 38.6% | 42.7 / 15.1 | 44 | D32E (GAT→GAA) ‡ | yeaR | DUF1971 domain‑containing protein YeaR |
* | NC_000913 | 2,470,472 | 0 | G | T | 38.5% | 17.6 / 14.3 | 26 | intergenic (+10/+605) | yfdI/yfdK | CPS‑53 (KpLE1) prophage; serotype specific glucosyl transferase/CPS‑53 (KpLE1) prophage; putative tail fiber assembly protein YfdK |
* | NC_000913 | 2,731,445 | 0 | C | G | 38.5% | 36.6 / 10.8 | 27 | intergenic (‑288/+155) | rrsG/clpB | 16S ribosomal RNA/ClpB chaperone |
* | NC_000913 | 1,331,179 | 0 | C | A | 37.2% | 55.0 / 17.4 | 43 | A44A (GCC→GCA) | topA | DNA topoisomerase 1 |
* | NC_000913 | 1,636,452 | 0 | C | G | 37.1% | 43.4 / 10.2 | 40 | intergenic (‑85/‑304) | nohA/ynfO | Qin prophage; putative prophage DNA‑packaging protein NohA/Qin prophage; protein YnfO |
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 | 4498181 = | NA (NA) | 9 (0.150) | 5/256 | 5.2 | NA | noncoding (1/1331 nt) | IS2 | repeat region |
? | NC_000913 | = 4499512 | NA (NA) | intergenic (+12/+81) | insD‑6/yjgX | KpLE2 phage‑like element; IS2 insertion element protein InsB/KpLE2 phage‑like element; putative phosphoethanolamine transferase YjgX, N‑terminal fragment | |||||
* | ? | NC_000913 | 381258 = | NA (NA) | 8 (0.140) | 5/254 | 5.2 | NA | pseudogene (3/417 nt) | yaiX | putative acyltransferase, N‑terminal fragment |
? | NC_000913 | = 382589 | NA (NA) | noncoding (1330/1331 nt) | IS2 | repeat region | |||||
* | ? | NC_000913 | = 889313 | 53 (0.900) | 4 (0.070) | 4/234 | 5.3 | 7.9% | coding (507/1686 nt) | ybjL | putative transport protein YbjL |
? | NC_000913 | = 992846 | 45 (0.840) | coding (2225/2613 nt) | pepN | aminopeptidase N | |||||
* | ? | NC_000913 | 1467910 = | NA (NA) | 6 (0.100) +A |
4/254 | 5.7 | NA | noncoding (1331/1331 nt) | IS2 | repeat region |
? | NC_000913 | = 1469240 | NA (NA) | noncoding (1/1331 nt) | IS2 | repeat region | |||||
* | ? | NC_000913 | = 1179385 | 53 (0.900) | 4 (0.070) | 4/252 | 5.7 | 6.8% | coding (793/912 nt) | nagK | N‑acetyl‑D‑glucosamine kinase |
? | NC_000913 | 1248808 = | 58 (1.010) | coding (307/1419 nt) | dhaM | dihydroxyacetone kinase subunit M | |||||
* | ? | NC_000913 | 453046 = | 57 (0.970) | 4 (0.070) | 4/252 | 5.7 | 6.6% | coding (500/1476 nt) | ampG | muropeptide:H(+) symporter |
? | NC_000913 | 3672936 = | 58 (1.010) | coding (523/900 nt) | yhjC | putative DNA‑binding transcriptional regulator YhjC | |||||
* | ? | NC_000913 | 381260 = | NA (NA) | 4 (0.070) | 4/254 | 5.7 | NA | noncoding (1/1331 nt) | IS2 | repeat region |
? | NC_000913 | = 382591 | NA (NA) | pseudogene (282/282 nt) | yaiX | putative acyltransferase, N‑terminal fragment |