| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 1,931,125 | A→G | E82G (GAG→GGG) | purT → | phosphoribosylglycinamide formyltransferase 2 |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,931,125 | 0 | A | G | 94.3% | 105.5 / 0.1 | 36 | coding (245/1179 nt) | purT | phosphoribosylglycinamide formyltransferase 2 |
| Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (18/15); minor base . (1/1); total (20/16) | |||||||||||
| 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.44e-01 | |||||||||||
TGCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGAGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGAGCTGCAG > NC_000913/1930997‑1931258 | tGCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAGGCt < 1:524446/139‑1 (MQ=255) gTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCCAAGCTATTg < 1:516751/139‑1 (MQ=255) tCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGc < 1:540619/139‑1 (MQ=255) cGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCg < 2:750839/139‑1 (MQ=255) cGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCa > 1:482729/1‑125 (MQ=255) cGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCa < 2:482729/125‑1 (MQ=255) cgcgCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGc > 2:97152/1‑105 (MQ=255) cgcgCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGc < 1:97152/105‑1 (MQ=255) gcgcATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGAc > 2:568819/1‑139 (MQ=255) cATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGa < 1:665265/139‑1 (MQ=255) tGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTcccc > 2:399647/1‑139 (MQ=255) ttAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTgcgc > 2:671320/1‑139 (MQ=255) ttAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTgcgc > 2:178364/1‑139 (MQ=255) tatGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTcgc > 2:659546/1‑139 (MQ=255) gCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGc < 2:734494/78‑1 (MQ=255) gCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGc > 1:734494/1‑78 (MQ=255) gCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGc < 2:128931/78‑1 (MQ=255) gCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGc > 1:128931/1‑78 (MQ=255) tGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCGCCTGCGCTCGCGCAACGAAATTa > 2:317477/1‑139 (MQ=255) aTGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGAt > 2:262217/1‑82 (MQ=255) aTGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGAt < 1:262217/82‑1 (MQ=255) aTGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATcc < 2:286173/84‑1 (MQ=255) aTGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATcc > 1:286173/1‑84 (MQ=255) aCGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCC‑‑‑GATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAAtgtt > 1:734105/1‑100 (MQ=255) aCGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCC‑‑‑GATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAAtgtt < 2:734105/100‑1 (MQ=255) ccGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAAt < 2:310492/125‑1 (MQ=255) ccGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAAt > 1:310492/1‑125 (MQ=255) cGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATcgcg < 2:404097/139‑1 (MQ=255) ggAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGc < 1:317477/139‑1 (MQ=255) ggAAAAACCACACTATATCGTACCGGCGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGc > 2:369934/1‑139 (MQ=255) gAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGCAGGAATGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGcc > 2:660744/1‑139 (MQ=255) gAAAAACCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTCAAGAGGAAGGACTGAATGTTGTCCCCTGCGCGCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGcc > 2:501365/1‑139 (MQ=255) aaaaaCCACATTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCt < 1:24719/139‑1 (MQ=255) aTTATATCGTGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCaga > 2:346192/1‑139 (MQ=255) tGCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGagctgca > 1:251188/1‑139 (MQ=255) gCCGGGGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGagctgcag > 1:435647/1‑139 (MQ=255) | TGCCGTCGATCGCTATGCCGACGCACCAGCCATGCATGTCGCGCATCGCTCCCATGTCATTAATATGCTTGATGGTGATGCATTACGCCGTGTGGTTGAACTGGAAAAACCACATTATATCGTGCCGGAGATCGAAGCTATTGCCACCGATATGCTGATCCAACTTGAAGAGGAAGGACTGAATGTTGTCCCCTGCGCTCGCGCAACGAAATTAACGATGAATCGCGAGGGTATCCGTCGCCTGGCGGCAGAAGAGCTGCAG > NC_000913/1930997‑1931258 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |