Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,727,747 | A→T | 46.3% | V161V (GTT→GTA) | purL ← | phosphoribosylformyl‑glycineamide synthetase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,727,747 | 0 | A | T | 46.3% | 3.6 / 43.4 | 41 | V161V (GTT→GTA) | purL | phosphoribosylformyl‑glycineamide synthetase |
Reads supporting (aligned to +/- strand): ref base A (13/9); new base T (3/16); total (16/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.25e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-01 |
AGCCAAGACGCAGGTTAGCGTCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCAAAAAA > minE/1727692‑1727812 | ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc < 1:1638554/66‑1 (MQ=255) ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc < 1:2687414/66‑1 (MQ=255) ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc < 1:2881304/66‑1 (MQ=255) ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc < 1:3073610/66‑1 (MQ=255) ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc < 1:3368557/66‑1 (MQ=255) ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc < 1:200277/66‑1 (MQ=255) ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc < 1:395604/66‑1 (MQ=255) ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc < 1:14080/66‑1 (MQ=255) ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc < 1:1398471/66‑1 (MQ=255) ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc < 1:788299/66‑1 (MQ=255) ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc < 1:1340355/66‑1 (MQ=255) ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc < 1:1312948/66‑1 (MQ=255) ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc < 1:101536/66‑1 (MQ=255) ggccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCCGCCCCAGCAAATCTACGCTGGTGAc < 1:3005027/66‑1 (MQ=255) gccaagACGCAGGTTAGCGTCAATTAGCGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc < 1:3208896/66‑1 (MQ=255) gCAGGTTAGCGTCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGtt < 1:3388933/69‑1 (MQ=255) cAGGTTAGCGTCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGtt < 1:2194700/68‑1 (MQ=255) aGGTTAGCGTCAATTAGAGCCTGACGGCCCTGCCCCAGCAAATCTACGCTGGTGAc < 1:1227277/56‑1 (MQ=255) tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa > 1:772448/1‑69 (MQ=255) tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa > 1:83341/1‑69 (MQ=255) tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa > 1:3355250/1‑69 (MQ=255) tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa > 1:2104785/1‑69 (MQ=255) tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa > 1:1233591/1‑69 (MQ=255) tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa > 1:1580735/1‑69 (MQ=255) tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa > 1:2621650/1‑69 (MQ=255) tCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCa > 1:2236198/1‑69 (MQ=255) gCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTg < 1:231309/55‑1 (MQ=255) gCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTg < 1:2036553/55‑1 (MQ=255) gCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTg < 1:1719437/55‑1 (MQ=255) gCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTg < 1:694959/55‑1 (MQ=255) gCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTg < 1:2778854/55‑1 (MQ=255) gCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTg < 1:1356538/55‑1 (MQ=255) cAAATCTACGCTGGTGACCGGAGTCGGTTGATGGTGTGCAAACAACTGCTCTGCATCATCTAAAGCaa > 1:3026406/1‑68 (MQ=255) cAAATCTACGCTGGTGACCGGAGTCGGTTGATGGTGTGCAAACAACTGCTCTGCATCATCTAAAGCaa > 1:3030001/1‑68 (MQ=255) cAAATCTACGCTGGTGACCGGAGTCGGTTGATGGTGTGCAAACAACTGCTCTGCATCATCTAAAGCaa > 1:1594409/1‑68 (MQ=255) aaaTCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTg > 1:2578648/1‑52 (MQ=255) aaaTCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCaa > 1:2254551/1‑67 (MQ=255) aaaTCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCaa > 1:2845279/1‑67 (MQ=255) aaaTCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCaa > 1:465013/1‑67 (MQ=255) aaaTCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCaa > 1:2740670/1‑67 (MQ=255) aTCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCaaaaaa < 1:558475/69‑1 (MQ=255) | AGCCAAGACGCAGGTTAGCGTCGATCAGCGCCTGACGGCCCTGCCCCAGCAAATCAACGCTGGTAACCGGAGTCGGTTGATGGTGGGCAAACAACTGCTCTGCATCATCTAAAGCAAAAAA > minE/1727692‑1727812 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |