| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,893,856 | A→G | 32.3% | G120G (GGT→GGC) | mltA ← | membrane‑bound lytic murein transglycosylase A |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,893,856 | 0 | A | G | 32.3% | 22.7 / 14.2 | 34 | G120G (GGT→GGC) | mltA | membrane‑bound lytic murein transglycosylase A |
| Reads supporting (aligned to +/- strand): ref base A (19/4); new base G (3/8); total (22/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.31e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01 | |||||||||||
AAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCGGCACCTTCCATCTGCCAGGCATCAATGCCGAACTGGCG > minE/1893795‑1893924 | aaaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTa > 1:2764166/1‑68 (MQ=255)aaaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTa > 1:2534312/1‑68 (MQ=255)aaaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTa > 1:2043328/1‑68 (MQ=255)aaaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTa > 1:2025315/1‑68 (MQ=255)aaaTAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTa > 1:1852046/1‑68 (MQ=255) ggATACTGGAACTCGCCCTTACGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAGCCTGTAAACTGc > 1:2734274/1‑69 (MQ=255) ggATACTGGAACTCGCCCTGACGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAGCCTGTAAACTGc > 1:2575965/1‑69 (MQ=255) ggATACTGGAACTCGCCCTGACGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAGCCTGTAAACTGc > 1:1785175/1‑69 (MQ=255) aCTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGt > 1:2245052/1‑69 (MQ=255) cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAAc > 1:121105/1‑61 (MQ=255) cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt > 1:1094074/1‑69 (MQ=255) cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt > 1:891671/1‑69 (MQ=255) cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt > 1:819391/1‑69 (MQ=255) cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt > 1:2457449/1‑69 (MQ=255) cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt > 1:1006853/1‑69 (MQ=255) cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt > 1:1806833/1‑69 (MQ=255) cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt > 1:180677/1‑69 (MQ=255) cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt > 1:1121640/1‑69 (MQ=255) cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt > 1:1776822/1‑69 (MQ=255) cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt > 1:1644289/1‑69 (MQ=255) cTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt > 1:1494682/1‑69 (MQ=255) cTGGAACTAGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt > 1:3346090/1‑69 (MQ=255) cTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATa < 1:3060947/69‑1 (MQ=255) gcgcgcCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCGGCACCTTCCAt < 1:1499678/69‑1 (MQ=255) ccTGAATTACCGGCGTGTAATAGCCTGTAAACTGCACGTTACCATAGTTGTCGACACCTTCCATCTGcc < 1:1829411/69‑1 (MQ=255) ccTGAATTACCGGCGTGTAATAGCCTGTAAACTGCACGTTACCATAGTTGTCGACACCTTCCATCTGcc < 1:714763/69‑1 (MQ=255) ccTGAATTACCGGCGTGTAATAGCCTGTAAACTGCACGTTACCATAGTTGTCGACACCTTCCATCTGcc < 1:623765/69‑1 (MQ=255) ccTGAATTACCGGCGTGTAATAGCCTGTAAACTGCACGTTACCATAGTTGTCGACACCTTCCATCTGcc < 1:610011/69‑1 (MQ=255) ccTGAATTACCGGCGTGTAATAGCCTGTAAACTGCACGTTACCATAGTTGTCGACACCTTCCATCTGcc < 1:2982644/69‑1 (MQ=255) ccTGAATTACCGGCGTGTAATAGCCTGTAAACTGCACGTTACCATAGTTGTCGACACCTTCCATCTGcc < 1:1749320/69‑1 (MQ=255) ccTGAATTACCGGCGTGTAATAGCCTGTAAACTGCACGTTACCATAGTTGTCGACACCTTCCATCTGcc < 1:2057087/69‑1 (MQ=255) ccTGAATTACCGGCGTGTAATAGCCTGTAAACTGCACGTTACCATAGTTGTCGACACCTTCCATCTGcc < 1:1830551/69‑1 (MQ=255) ttACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCGGCACCTTCCATCTGCCAGGCAt < 1:1495160/69‑1 (MQ=255) ttACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCGGCACCTTCCATCTGCCAGGCAt < 1:90547/69‑1 (MQ=255) aCCGGTAAACTGCACGTTACCATAGTTGTCGGCACCTTCCATCTGCCAGGCATCAATGCCGAACTGgcg > 1:2771094/1‑69 (MQ=255) | AAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCGGCACCTTCCATCTGCCAGGCATCAATGCCGAACTGGCG > minE/1893795‑1893924 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |