Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,872,275 | T→C | 77.9% | C176R (TGC→CGC) | ulaF → | L‑ribulose 5‑phosphate 4‑epimerase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,872,275 | 0 | T | C | 77.9% | 110.5 / 22.0 | 81 | C176R (TGC→CGC) | ulaF | L‑ribulose 5‑phosphate 4‑epimerase |
Reads supporting (aligned to +/- strand): ref base T (16/2); new base C (30/33); total (46/35) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.38e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01 |
TTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATT > minE/2872226‑2872338 | ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGATAGATGCCCACGACGCGGTGCATAACGCggtg > 1:74650/1‑68 (MQ=255) ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg > 1:2533167/1‑68 (MQ=255) ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg > 1:984315/1‑68 (MQ=255) ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg > 1:681360/1‑68 (MQ=255) ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg > 1:520097/1‑68 (MQ=255) ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg > 1:506326/1‑68 (MQ=255) ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg > 1:340682/1‑68 (MQ=255) ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg > 1:3059951/1‑68 (MQ=255) ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg > 1:3034215/1‑68 (MQ=255) ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg > 1:2885186/1‑68 (MQ=255) ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg > 1:2370091/1‑68 (MQ=255) ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg > 1:2156777/1‑68 (MQ=255) ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg > 1:1881764/1‑68 (MQ=255) ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg > 1:1708479/1‑68 (MQ=255) ttGTGGTGTATCAGCACGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCggtg > 1:1614997/1‑68 (MQ=255) tgtgGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCggtggt < 1:3303790/69‑1 (MQ=255) gtggtgTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCggtgg > 1:912962/1‑67 (MQ=255) gtggtgTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGgtggtg > 1:2115442/1‑69 (MQ=255) tCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaaga > 1:1289280/1‑69 (MQ=255) tCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaaga > 1:3089497/1‑69 (MQ=255) tCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaaga > 1:2184702/1‑69 (MQ=255) tCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaaga > 1:2710180/1‑69 (MQ=255) tCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaaga > 1:2077280/1‑69 (MQ=255) aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGga > 1:1180736/1‑64 (MQ=255) aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGaagaag > 1:540956/1‑69 (MQ=255) aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGaagaag > 1:2089168/1‑69 (MQ=255) aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGaagaag > 1:2721466/1‑69 (MQ=255) aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc > 1:1647915/1‑37 (MQ=25) aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc > 1:3197161/1‑37 (MQ=25) aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc > 1:3186756/1‑37 (MQ=25) aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc > 1:357796/1‑37 (MQ=25) aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc > 1:2956298/1‑37 (MQ=25) aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc > 1:3020342/1‑37 (MQ=25) aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc > 1:2468780/1‑37 (MQ=25) aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc > 1:1326801/1‑37 (MQ=25) aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc > 1:1401308/1‑37 (MQ=25) aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc > 1:554411/1‑37 (MQ=25) aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc > 1:1790361/1‑37 (MQ=25) aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc > 1:745873/1‑37 (MQ=25) aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc > 1:1838338/1‑37 (MQ=25) aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc > 1:2177427/1‑37 (MQ=25) aCGGACCGTTCGCCTGGGGGAAAGATGCCCACGACGc > 1:2061916/1‑37 (MQ=25) gCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaa < 1:902955/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGGCGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:3265477/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:2735877/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:830459/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:77709/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:747051/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:328907/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:1258623/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:3397796/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:1379368/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:627055/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:383387/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:1327710/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:2976774/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:172296/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:1758086/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:1782504/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:1633490/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:2049078/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:2217703/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:2252857/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:1476252/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:2375466/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:244404/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:2601116/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:26156/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:279014/69‑1 (MQ=255) cGTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:1404523/69‑1 (MQ=255) gTTCGCCTGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:284224/68‑1 (MQ=255) tGGGGGAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:114505/61‑1 (MQ=255) gggtgAAAGATGCCCACGACGCGGTGCATAACGCGGTGGTGATGTAAGAAGTGGCGaaaa < 1:2198662/60‑1 (MQ=255) cccACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:265950/46‑1 (MQ=38) cccACGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:2160215/46‑1 (MQ=38) cACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCAtt > 1:2332336/1‑69 (MQ=255) cACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCAtt > 1:2180039/1‑69 (MQ=255) cACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCAtt > 1:2730504/1‑69 (MQ=255) cACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCAtt > 1:75472/1‑69 (MQ=255) cACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCAtt > 1:1406114/1‑69 (MQ=255) cGACGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa < 1:2024953/44‑1 (MQ=255) | TTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATT > minE/2872226‑2872338 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |