| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,286,937 | C→G | 66.1% | A150A (GCG→GCC) | argC ← | N‑acetyl‑gamma‑glutamylphosphate reductase, NAD(P)‑binding |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,286,937 | 0 | C | G | 66.1% | 54.3 / 39.4 | 56 | A150A (GCG→GCC) | argC | N‑acetyl‑gamma‑glutamylphosphate reductase, NAD(P)‑binding |
| Reads supporting (aligned to +/- strand): ref base C (14/5); new base G (22/15); total (36/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.82e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAATTTATTACCGCACCACTCCGCCAGACCGTAG > minE/2286878‑2286991 | tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATaaag > 1:1424902/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag > 1:1713827/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag > 1:872402/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag > 1:718995/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag > 1:2460240/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag > 1:2497330/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag > 1:2571436/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag > 1:1872846/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag > 1:268451/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag > 1:1805164/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag > 1:641308/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag > 1:1647529/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag > 1:1572715/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag > 1:1551126/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag > 1:281663/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag > 1:2984344/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaag > 1:3428406/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaa > 1:482846/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaa > 1:258202/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaa > 1:3358785/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaa > 1:1112898/1‑68 (MQ=255)tCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCaa > 1:1377326/1‑68 (MQ=255) atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc < 1:245206/69‑1 (MQ=255) atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc < 1:532208/69‑1 (MQ=255) atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc < 1:544889/69‑1 (MQ=255) atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc < 1:3234805/69‑1 (MQ=255) atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc < 1:3117192/69‑1 (MQ=255) atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc < 1:3095463/69‑1 (MQ=255) atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc < 1:3078080/69‑1 (MQ=255) atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc < 1:1176309/69‑1 (MQ=255) atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc < 1:2929654/69‑1 (MQ=255) atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc < 1:1068179/69‑1 (MQ=255) atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc < 1:688170/69‑1 (MQ=255) atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc < 1:271149/69‑1 (MQ=255) atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc < 1:950053/69‑1 (MQ=255) atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc < 1:1886862/69‑1 (MQ=255) atcaatcaACGGTTTCAGCGCCAGCTGTGCTGCCGTCGGATAACAGCCTGGCACGGCAATCAAGTTCGc < 1:2115201/69‑1 (MQ=255) aatcaaCGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGcttc > 1:1972490/1‑69 (MQ=255) aatcaaCGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGcttc > 1:3244726/1‑69 (MQ=255) aatcaaCGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGcttc > 1:1071325/1‑69 (MQ=255) caaCGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCa > 1:918793/1‑56 (MQ=255) caaCGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCa > 1:2194924/1‑56 (MQ=255) ttCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAAttt > 1:621768/1‑68 (MQ=255) ttCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAAttt > 1:995849/1‑68 (MQ=255) ttCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAAttt > 1:396870/1‑68 (MQ=255) atCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAAttt > 1:2722005/2‑68 (MQ=255) tCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAATttat < 1:2814048/69‑1 (MQ=255) tCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAATttat < 1:2787722/69‑1 (MQ=255) cgcCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCtttt < 1:534484/58‑1 (MQ=255) cgcCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCtttt < 1:2001921/58‑1 (MQ=255) cGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCtttt < 1:1813406/40‑1 (MQ=255) tAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAATTTATTACCGCACCACTCCGCCAGAcc > 1:38202/1‑65 (MQ=255) tAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAATTTATTACCGCACCACTCCGCCAGACCGTAg > 1:223440/1‑69 (MQ=255) tAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAATTTATTACCGCACCACTCCGCCAGACCGTAg > 1:940097/1‑69 (MQ=255) tAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAATTTATTACCGCACCACTCCGCCAGACCGTAg > 1:1699649/1‑69 (MQ=255) tAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAATTTATTACCGCACCACTCAGCCAGACCGTAg > 1:2363374/1‑69 (MQ=255) | TCGGCATCAATCAACGGTTTCAGCGCCAGCTGTGCCGCCGTCGGATAACAGCCCGGCACCGCAATCAAATTCGCTTCTTTTAATTTATTACCGCACCACTCCGCCAGACCGTAG > minE/2286878‑2286991 |
| 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: ‑ |