| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 2,939,281 | C→T | 19.4% | E207K (GAA→AAA) | ygdE ← | predicted methyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 2,939,281 | 0 | C | T | 19.4% | 77.2 / 13.0 | 36 | E207K (GAA→AAA) | ygdE | predicted methyltransferase |
| Reads supporting (aligned to +/- strand): ref base C (19/10); new base T (5/2); total (24/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01 | |||||||||||
AGGTCCAGCCGCCAGGGCAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAGCGTGGAACGACTCGGCG > W3110S.gb/2939221‑2939349 | aGGTCCAGCCGCCAGGGCAAGCGCCTAAATCTACCGCCCACATCCCGTTCGCCAGGCGTTCa > 1:3591508/1‑62 (MQ=255)aGGTCCAGCCGCCAGGGCAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCa > 1:2745768/1‑62 (MQ=255) ggTCCAGCCGCCAGGGCAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCAt > 1:49053/1‑71 (MQ=255) gTCCAGCCGCCAGGGCAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATc < 1:2750942/71‑1 (MQ=255) ccAGCCGCCAGGGCAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTTATCCCACTCATCCg > 1:1065707/1‑71 (MQ=255) ccgccAGGGCAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTTATCCCACTCATCCGCAgg < 1:1453129/71‑1 (MQ=255) gccAGGGCAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAg > 1:3763261/1‑68 (MQ=255) cAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAacac < 1:1826025/71‑1 (MQ=255) cAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAaca > 1:2285266/1‑70 (MQ=255) tAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTTATCCCACTCATCCGCAGGAATAAACACATGAAATg > 1:1298776/1‑71 (MQ=255) tAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCa > 1:746538/1‑37 (MQ=255) tAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCa > 1:2243517/1‑37 (MQ=255) tAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCATcc > 1:156823/1‑40 (MQ=255) aaaTCCACCGCCCACATCCCGTTAGCCAGGCGTTTATCCCACTCATCCGCAGGAATAAACACATGAAATGc > 1:1078584/1‑71 (MQ=255) aaTCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAAt > 1:3890413/1‑55 (MQ=255) aaTCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGc < 1:3956742/70‑1 (MQ=255) aTCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAgg < 1:637201/51‑1 (MQ=255) aCCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTc > 1:1962797/1‑71 (MQ=255) acctCCCGTTCGCCAGGCGTTTATGCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTg < 1:3716589/68‑1 (MQ=255) acaTCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGttt > 1:4086922/1‑70 (MQ=255) acaTCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTg < 1:3774830/71‑1 (MQ=255) caTCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGaaa > 1:3888177/1‑54 (MQ=255) tCCCGTTCGCCAGGCGTTTATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAGc > 1:1111793/1‑71 (MQ=255) cccGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGaaa > 1:1779247/1‑51 (MQ=255) cccGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAGCg > 1:2489209/1‑71 (MQ=255) ccGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAGCg < 1:4456292/70‑1 (MQ=255) ttCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAg > 1:12922/1‑65 (MQ=255) ttCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAGCGTGGa < 1:2725378/71‑1 (MQ=255) gCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTg < 1:3057537/55‑1 (MQ=255) gCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTg < 1:282538/55‑1 (MQ=255) gCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAGCGTGGAACGACTCg > 1:3956415/1‑71 (MQ=255) gTTTATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAGCGTgg > 1:4015154/1‑60 (MQ=255) gTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAGCGTGGAACGACTCgg > 1:4338399/1‑70 (MQ=255) gTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAGCGTGGAACGACTCGgc > 1:1597873/1‑71 (MQ=255) gTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAGCGTGGAACGACTCGgc > 1:2581023/1‑71 (MQ=255) ttCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAGCGTGGAACGACTCGgcg < 1:3757451/71‑1 (MQ=255) cATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAg > 1:315729/1‑46 (MQ=255) | AGGTCCAGCCGCCAGGGCAAGCGCCTAAATCCACCGCCCACATCCCGTTCGCCAGGCGTTCATCCCACTCATCCGCAGGAATAAACACATGAAATGCCTCTTCCAGTTTGAGCGTGGAACGACTCGGCG > W3110S.gb/2939221‑2939349 |
| 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: ‑ |