Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,237,256 | C→T | 48.1% | T267T (ACC→ACT) | purD → | phosphoribosylglycinamide synthetase phosphoribosylamine‑glycine ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,237,256 | 0 | C | T | 48.1% | 4.7 / 68.1 | 79 | T267T (ACC→ACT) | purD | phosphoribosylglycinamide synthetase phosphoribosylamine‑glycine ligase |
Reads supporting (aligned to +/- strand): ref base C (26/15); new base T (9/29); total (35/44) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.97e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01 |
CCATGGAACGTATTATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTATTGAATTTA > minE/2237195‑2237317 | cgaTGGAACTTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc < 1:3413730/67‑1 (MQ=255) cgaTGGAACGTATCATCTTGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc < 1:1767025/67‑1 (MQ=255) cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc < 1:316764/67‑1 (MQ=255) cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc < 1:2296628/67‑1 (MQ=255) cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc < 1:625602/67‑1 (MQ=255) cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc < 1:3009562/67‑1 (MQ=255) cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc < 1:1543870/67‑1 (MQ=255) cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc < 1:1539981/67‑1 (MQ=255) cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc < 1:1466786/67‑1 (MQ=255) cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc < 1:2381697/67‑1 (MQ=255) cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc < 1:2747145/67‑1 (MQ=255) cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc < 1:1392844/67‑1 (MQ=255) cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc < 1:3169456/67‑1 (MQ=255) cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc < 1:1300613/67‑1 (MQ=255) cgaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc < 1:2724498/67‑1 (MQ=255) ccATGGAACGTATTATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTtc < 1:2328584/69‑1 (MQ=255) ccATGGAACGTATTATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTtc < 1:1742077/69‑1 (MQ=255) gaTGGAACGTATCATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc < 1:2374016/67‑1 (MQ=255) gaTGGAACGTATCATCTGGCCAACCGGGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTtc < 1:2691266/67‑1 (MQ=255) aTGGAACGTATTATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTtctc < 1:2224677/69‑1 (MQ=255) aTGGAACGTATTATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTtctc < 1:2516787/69‑1 (MQ=255) tGGAACGTATTATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTtctc < 1:803330/68‑1 (MQ=255) ggAACGTATTATCTTGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTa < 1:1573690/69‑1 (MQ=255) ggAACGTATTATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTa < 1:3174764/69‑1 (MQ=255) atCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCtg > 1:177560/1‑69 (MQ=255) atCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCtg > 1:3138806/1‑69 (MQ=255) atCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCtg > 1:1821432/1‑69 (MQ=255) atCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCtg > 1:402404/1‑69 (MQ=255) atCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCtg > 1:697330/1‑69 (MQ=255) atCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCtg > 1:1365191/1‑69 (MQ=255) atCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCtg > 1:922085/1‑69 (MQ=255) atCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCt < 1:1195699/68‑1 (MQ=255) atCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACAACGGTTTTCTCTACGCGGGCCtg > 1:612252/1‑69 (MQ=255) cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCTGGCCTGATGATc < 1:727317/68‑1 (MQ=255) cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc < 1:3161769/68‑1 (MQ=255) cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc < 1:3116437/68‑1 (MQ=255) cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc < 1:97953/68‑1 (MQ=255) cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc < 1:3767/68‑1 (MQ=255) cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc < 1:2989194/68‑1 (MQ=255) cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc < 1:72294/68‑1 (MQ=255) cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc < 1:1737031/68‑1 (MQ=255) cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc < 1:232747/68‑1 (MQ=255) cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc < 1:2087831/68‑1 (MQ=255) cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc < 1:2029630/68‑1 (MQ=255) cAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATc < 1:2012588/68‑1 (MQ=255) aaCCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCg > 1:1453894/1‑68 (MQ=255) aaCCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCg > 1:3331799/1‑68 (MQ=255) gTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACaaa < 1:920067/69‑1 (MQ=255) gTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACaaa < 1:2686332/69‑1 (MQ=255) gTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACaaa < 1:3278146/69‑1 (MQ=255) gTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACaaa < 1:585971/69‑1 (MQ=255) ggCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCtgatg > 1:945704/1‑54 (MQ=255) ggCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCtgatg > 1:1938698/1‑54 (MQ=255) ggcggcGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAAt > 1:3433055/1‑67 (MQ=255) ggcggcGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATc > 1:385004/1‑68 (MQ=255) ggcggcGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATc > 1:1609250/1‑68 (MQ=255) ggcggcGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATc > 1:1607183/1‑68 (MQ=255) ggcggcGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATc > 1:1602593/1‑68 (MQ=255) ggcggcGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATc > 1:2948436/1‑68 (MQ=255) ggcggcGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATc > 1:2845646/1‑68 (MQ=255) ggcggcGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATc > 1:1238327/1‑68 (MQ=255) ggcggcGGAAGGCAACACCTACACTGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATc > 1:410654/1‑68 (MQ=255) ggAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGt < 1:2083197/69‑1 (MQ=255) aGGCAACACCTACACCGTTTTTCTCTACGCGGGCCTGATTATCGACaaa < 1:109093/49‑1 (MQ=255) ggCAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTAt > 1:457574/1‑68 (MQ=255) cAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACa > 1:629919/1‑44 (MQ=255) cAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACa > 1:1780586/1‑44 (MQ=255) cAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACa > 1:1621190/1‑44 (MQ=255) cAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACa > 1:362252/1‑44 (MQ=255) cAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACa > 1:2528163/1‑44 (MQ=255) cAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACa > 1:2704454/1‑44 (MQ=255) acCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTATTGAAttt > 1:2074734/1‑69 (MQ=255) acCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTATTGAAttt > 1:3076071/1‑69 (MQ=255) acCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTATTGAAttt > 1:1305898/1‑69 (MQ=255) acCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTATTGAAttt < 1:2766999/69‑1 (MQ=255) acCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTATTGAAttt > 1:2458426/1‑69 (MQ=255) acCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTATTGAAtt > 1:203328/1‑68 (MQ=255) ccTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTATTGAATTTa > 1:3011119/1‑69 (MQ=255) ccTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTATTGAATTTa > 1:1407553/1‑69 (MQ=255) | CCATGGAACGTATTATCTGGCCAACCGTGAAAGGCATGGCGGCGGAAGGCAACACCTACACCGGTTTTCTCTACGCGGGCCTGATGATCGACAAACAGGGCAATCCGAAGGTTATTGAATTTA > minE/2237195‑2237317 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |