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 228 shown, sorted by frequency from high to low) | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,515,909 | 0 | A | C | 54.9% | 59.5 / 106.1 | 71 | V260G (GTG→GGG) | ypfI | predicted hydrolase |
* | REL606 | 720,464 | 0 | C | G | 50.0% | 3.0 / 23.9 | 20 | A22G (GCC→GGC) | phrB | deoxyribodipyrimidine photolyase, FAD‑binding |
* | REL606 | 1,854,401 | 0 | A | C | 45.9% | 74.7 / 53.3 | 85 | K88N (AAA→AAC) | yeaO | hypothetical protein |
* | REL606 | 433 | 0 | T | G | 44.3% | 71.5 / 50.5 | 70 | V33G (GTG→GGG) | thrA | bifunctional aspartokinase I/homeserine dehydrogenase I |
* | REL606 | 2,437,406 | 0 | G | T | 40.7% | 35.0 / 17.7 | 27 | D435E (GAC→GAA) | ypdD | fused predicted PTS enzymes: Hpr component/enzyme I component/enzyme IIA component |
* | REL606 | 2,376,348 | 0 | T | G | 39.5% | 211.6 / 22.8 | 131 | T63P (ACC→CCC) | folC | bifunctional folylpolyglutamate synthase/ dihydrofolate synthase |
* | REL606 | 895,856 | 0 | G | C | 39.5% | 51.9 / 31.1 | 38 | A26P (GCC→CCC) | ECB_00836 | putative Regulatory protein |
* | REL606 | 1,906,430 | 0 | A | C | 39.1% | 251.7 / 52.7 | 133 | V329G (GTG→GGG) | ptrB | protease II |
* | REL606 | 4,111,429 | 0 | T | G | 38.8% | 105.1 / 18.9 | 67 | G496G (GGT→GGG) | metL | bifunctional aspartate kinase II/homoserine dehydrogenase II |
* | REL606 | 4,180,249 | 0 | T | G | 38.2% | 198.0 / 24.3 | 111 | G24G (GGT→GGG) | yjaH | hypothetical protein |
* | REL606 | 3,378,885 | 0 | T | G | 36.6% | 463.2 / 66.0 | 224 | T129P (ACC→CCC) | rplB | 50S ribosomal protein L2 |
* | REL606 | 157,994 | 0 | G | T | 36.4% | 254.0 / 163.9 | 140 | Q93K (CAA→AAA) | htrE | predicted outer membrane usher protein |
* | REL606 | 1,485,450 | 0 | A | G | 36.1% | 348.0 / 120.0 | 194 | G193G (GGA→GGG) | ydcV | predicted spermidine/putrescine transporter subunit |
* | REL606 | 2,523,349 | 0 | A | C | 35.9% | 226.6 / 44.1 | 117 | H235P (CAC→CCC) | hyfB | NADH dehydrogenase subunit N |
* | REL606 | 1,118,209 | 0 | G | T | 35.3% | 38.5 / 25.9 | 34 | intergenic (‑444/‑311) | csgD/csgB | DNA‑binding transcriptional activator in two‑component regulatory system/curlin nucleator protein, minor subunit in curli complex |
* | REL606 | 3,783,565 | 0 | A | C | 34.6% | 178.1 / 45.0 | 101 | T744P (ACA→CCA) | ECB_03521 | conserved hypothetical protein |
* | REL606 | 3,217,434 | 0 | G | C | 34.1% | 41.6 / 23.6 | 44 | A76P (GCC→CCC) | agaS | tagatose‑6‑phosphate ketose/aldose isomerase |
* | REL606 | 2,685,638 | 0 | T | G | 33.6% | 235.4 / 38.1 | 110 | V326G (GTG→GGG) | ygaF | predicted enzyme |
* | REL606 | 851,606 | 0 | G | A | 33.6% | 154.2 / 80.1 | 113 | S3F (TCT→TTT) | ybiU | hypothetical protein |
* | REL606 | 4,595,070 | 0 | A | G | 33.6% | 278.0 / 80.6 | 149 | G11G (GGA→GGG) | yjjZ | hypothetical protein |
Marginal new junction evidence (sorted from low to high skew) | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | REL606 | 1532002 = | 227 (1.050) | 5 (0.030) | 3/130 | 19.4 | 2.3% | coding (954/987 nt) | yddP | D‑ala‑D‑ala transporter subunit |
? | REL606 | = 3014081 | 228 (1.220) | coding (662/1137 nt) | ECB_02815 | Capsule polysaccharide export inner‑membrane protein kpsE | |||||
* | ? | REL606 | 984122 = | 261 (1.210) | 3 (0.020) | 3/130 | 19.4 | 1.3% | coding (410/1749 nt) | msbA | fused lipid transporter subunits of ABC superfamily: membrane component/ATP‑binding component |
? | REL606 | = 984118 | 219 (1.180) | coding (406/1749 nt) | msbA | fused lipid transporter subunits of ABC superfamily: membrane component/ATP‑binding component | |||||
* | ? | REL606 | = 2525354 | 173 (0.800) | 3 (0.020) | 3/132 | 19.5 | 1.6% | pseudogene (453/720 nt) | hyfC | hydrogenase 4 membrane subunit; b2483_1 |
? | REL606 | = 2718592 | 210 (1.110) | coding (37/1062 nt) | recA | recombinase A |