| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 4,219,725 | A→G | N377S (AAC→AGC) | aceK → | isocitrate dehydrogenase kinase/phosphatase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,219,725 | 0 | A | G | 80.7% | 62.6 / 8.6 | 31 | N377S (AAC→AGC) | aceK | isocitrate dehydrogenase kinase/phosphatase |
| Reads supporting (aligned to +/- strand): ref base A (5/1); new base G (10/15); total (15/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.29e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACC > NC_000913/4219647‑4219801 | cTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAGCt < 1:793749/81‑1 (MQ=255) gTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAGCTTTGTGCTgg < 2:793749/81‑1 (MQ=255) gTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTg > 2:793790/1‑80 (MQ=255) ccTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAGCTTTGTGCTGgag < 1:793751/80‑1 (MQ=255) tGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAGCTTTGTGCTGGAGAAg > 1:793782/1‑81 (MQ=255) cTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAGCTTTGTGCTGGAGAAGCg < 2:793752/81‑1 (MQ=255) cTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAGCTTTGTGCTGGAGAAGCg < 2:793751/81‑1 (MQ=255) aaCTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAGCTTTGTGCTGGAGAAGCGGCat < 2:793753/80‑1 (MQ=255) gTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAGCTTTGTGCTGGAGAAGCGGCATATTTc < 2:793754/80‑1 (MQ=255) aaaGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAGCTTTGTGCTGGAGAAGCGGCATATTTcccc < 1:793756/80‑1 (MQ=255) aaaGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAGCTTTGTGCTGGAGAAGCGGCATATTTCCCCg > 1:793788/1‑81 (MQ=255) aaaGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCg > 1:793795/1‑81 (MQ=255) agCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCAt > 2:793792/1‑81 (MQ=255) gATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAGCTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATg > 1:793798/1‑81 (MQ=255) gATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATg < 2:793756/81‑1 (MQ=255) aTCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAGCTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATgg < 2:793757/81‑1 (MQ=255) aTCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAGCTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATg > 2:793798/1‑80 (MQ=255) gcctGGGCCGAATGGCGGACACCCAGGAGTTTGAAAGCTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAAt < 2:793758/78‑1 (MQ=255) tGGGCCGAATGGCGGACACCCAGGAGTTTGAAAGCTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTAc > 2:793796/1‑81 (MQ=255) gggCCGAATGGCGGACACCCAGGAGTTTGAAAGCTTTGTGCTGGAGAAGTGGCATATTTCCCCGGCATTAATGGAATTACt > 2:793799/1‑81 (MQ=255) aTGGCGGACACCCAGGAGTTTGAAAGCTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAg > 2:793804/1‑81 (MQ=255) gCGGACACCCAGGAGTTTGAAAGCTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGaa < 1:793760/81‑1 (MQ=255) cGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAag > 1:793787/1‑81 (MQ=255) ggACACCCAGGAGTTTGAAAGCTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAagc < 2:793761/81‑1 (MQ=255) ggACACCCAGGAGTTTGAAAGCTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAagc < 1:793761/81‑1 (MQ=255) ggACACCCAGGAGTTTGAAAGCTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAag > 2:793808/1‑80 (MQ=255) aGGAGTTTGAAAGCTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGaaa < 1:793763/81‑1 (MQ=255) ggAGTTTTAAAACTTTGTTCTCGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGtaa > 2:793806/1‑80 (MQ=255) gTTTGAAAGCTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAAt > 2:793793/1‑81 (MQ=255) tttGAAAGCTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATc < 2:793763/81‑1 (MQ=255) tGAAAGCTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCAc > 2:793811/1‑81 (MQ=255) aaagCTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCAcc < 1:793765/80‑1 (MQ=255) | CTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACC > NC_000913/4219647‑4219801 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |