| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 244,855 | C→A | 13.5% | G189G (GGG→GGT) | sbcD ← | exonuclease, dsDNA, ATP‑dependent |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 244,855 | 0 | C | A | 13.5% | 129.2 / 8.6 | 59 | G189G (GGG→GGT) | sbcD | exonuclease, dsDNA, ATP‑dependent |
| Reads supporting (aligned to +/- strand): ref base C (12/39); new base A (5/3); total (17/42) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.72e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGC > minE/244800‑244920 | gAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTa < 1:1134066/69‑1 (MQ=255) gtacAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTcc < 1:2621840/66‑1 (MQ=255) gTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTcc < 1:2508806/69‑1 (MQ=255) gTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCACCCACGGTCGTTAAATGTcc < 1:482626/69‑1 (MQ=255) gTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCACCCACGGTCGTTAAATGTcc < 1:2665473/69‑1 (MQ=255) gTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCACCCACGGTCGTTAAATGTcc < 1:352699/69‑1 (MQ=255) cGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGc < 1:2358031/69‑1 (MQ=255) cGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGc < 1:1221212/69‑1 (MQ=255) cGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGc < 1:2645406/69‑1 (MQ=255) cGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGc < 1:2253780/69‑1 (MQ=255) cGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGc < 1:2959293/69‑1 (MQ=255) cGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGc < 1:3389254/69‑1 (MQ=255) cGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGc < 1:1569623/69‑1 (MQ=255) gTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCg < 1:3120764/69‑1 (MQ=255) gCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCgat < 1:1979727/69‑1 (MQ=255) gCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCgat < 1:3037294/69‑1 (MQ=255) gCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCgat < 1:2162327/69‑1 (MQ=255) gCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCgat < 1:2201136/69‑1 (MQ=255) gCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCgat < 1:392761/69‑1 (MQ=255) gCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCgat < 1:2279141/69‑1 (MQ=255) gCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCgat < 1:1503754/69‑1 (MQ=255) gCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCgat < 1:1439670/69‑1 (MQ=255) ccAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTaa > 1:3168520/1‑53 (MQ=255) ccAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCgat < 1:312563/68‑1 (MQ=255) ataAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGt > 1:1772258/1‑45 (MQ=255) ataAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTaaa > 1:1306455/1‑49 (MQ=255) ataAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTa < 1:2534591/47‑1 (MQ=255) ataAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCgat < 1:242962/63‑1 (MQ=255) ataAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCgat < 1:870122/63‑1 (MQ=255) ataAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCgat < 1:670847/63‑1 (MQ=255) ataAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCgat < 1:637237/63‑1 (MQ=255) ataAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCgat < 1:1808267/63‑1 (MQ=255) ataAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCgat < 1:2104942/63‑1 (MQ=255) ataAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCgat < 1:292894/63‑1 (MQ=255) ataAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCgat < 1:246935/63‑1 (MQ=255) ataAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCgat < 1:2342865/63‑1 (MQ=255) taAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTg > 1:2108233/1‑57 (MQ=255) taAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTg > 1:2547552/1‑57 (MQ=255) gTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCagag > 1:747273/1‑68 (MQ=255) tCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTTAATGTCCCGTGGCGATGATGGGCAGAgg > 1:805202/1‑68 (MQ=255) tCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAgg > 1:2233780/1‑68 (MQ=255) aCTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGCTGATGGGCAGAgg > 1:2311794/1‑54 (MQ=255) aCTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAg > 1:1065596/1‑69 (MQ=255) aCTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATAGCCGCGCAg > 1:1651476/1‑69 (MQ=255) aCTTTTACTGGCCCCCACGGTCGTTAAATGTACCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAg > 1:1229124/1‑69 (MQ=255) cTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGc < 1:1867033/38‑1 (MQ=255) ttACTGGCACCCACGGTCGTTAAATGTCCAGTGGCgatga > 1:3181311/1‑40 (MQ=38) ttACTGGCACCCACGGTCGTTAAATGTCCAGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCa > 1:384444/1‑64 (MQ=255) ttACTGGCACCCACGGTCGTTAAATGTCCAGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCa > 1:387444/1‑64 (MQ=255) ttACTGGCACCCACGGTCGTTAAATGTCCAGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCa > 1:3104966/1‑64 (MQ=255) ttACTGGCACCCACGGTCGTTAAATGTCCAGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGtt > 1:1916487/1‑67 (MQ=255) ggCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGc < 1:2569534/69‑1 (MQ=255) ggCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGc < 1:2265744/69‑1 (MQ=255) ggCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGc < 1:1624038/69‑1 (MQ=255) ggCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGc < 1:3137999/69‑1 (MQ=255) ggCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGc < 1:3106905/69‑1 (MQ=255) ggCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGc < 1:843796/69‑1 (MQ=255) ggCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGc < 1:288550/69‑1 (MQ=255) ggCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGAGCAGAGGCTGATCGCCGCGCAGTTTGCAGGc < 1:2956011/69‑1 (MQ=255) | GAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGC > minE/244800‑244920 |
| 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: ‑ |