| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,893,859 | G→T | 31.4% | T119T (ACC→ACA) | 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,859 | 0 | G | T | 31.4% | 38.4 / 18.2 | 35 | T119T (ACC→ACA) | mltA | membrane‑bound lytic murein transglycosylase A |
| Reads supporting (aligned to +/- strand): ref base G (20/4); new base T (3/8); total (23/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.24e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-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: ‑ |