| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,039,561 | C→T | 41.2% | G116G (GGC→GGT) | nth → | DNA glycosylase and apyrimidinic (AP) lyase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,039,561 | 0 | C | T | 41.2% | 24.1 / 42.9 | 63 | G116G (GGC→GGT) | nth | DNA glycosylase and apyrimidinic (AP) lyase |
| Reads supporting (aligned to +/- strand): ref base C (13/24); new base T (14/12); total (27/36) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.97e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01 | |||||||||||
TCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCT > minE/1039497‑1039621 | tCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAg > 1:2877045/1‑69 (MQ=255)tCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAg > 1:3272381/1‑69 (MQ=255) gCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTc > 1:2731509/1‑68 (MQ=255) gCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTc > 1:1113102/1‑68 (MQ=255) gCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTc > 1:3243870/1‑68 (MQ=255) gCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTc > 1:1004949/1‑68 (MQ=255) gCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTc > 1:95574/1‑68 (MQ=255) ggAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTaa < 1:992508/69‑1 (MQ=255) ggAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTaa < 1:2365172/69‑1 (MQ=255) ggAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTaa < 1:2361242/69‑1 (MQ=255) ggAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTaa < 1:22935/69‑1 (MQ=255) ggAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTaa < 1:2244559/69‑1 (MQ=255) tAATGGCGAGGTTCCTGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCa > 1:791961/1‑68 (MQ=255) tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACa > 1:1672572/1‑64 (MQ=255) tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa > 1:1537942/1‑69 (MQ=255) tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa > 1:3093530/1‑69 (MQ=255) tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa > 1:2955234/1‑69 (MQ=255) tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa > 1:2653480/1‑69 (MQ=255) tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa > 1:2598168/1‑69 (MQ=255) tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa > 1:2498331/1‑69 (MQ=255) tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa > 1:2469014/1‑69 (MQ=255) tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa > 1:913401/1‑69 (MQ=255) tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa > 1:2132180/1‑69 (MQ=255) tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa > 1:1608366/1‑69 (MQ=255) tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa > 1:2013940/1‑69 (MQ=255) tAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCaa > 1:1804106/1‑69 (MQ=255) aTGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAcg > 1:1443645/1‑69 (MQ=255) aTGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAcg > 1:2537885/1‑69 (MQ=255) aTGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAcg > 1:417172/1‑69 (MQ=255) ggCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAAc < 1:312947/59‑1 (MQ=255) ggCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAAc < 1:3426769/59‑1 (MQ=255) ggCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAAc < 1:1811971/59‑1 (MQ=255) gAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa < 1:1687090/69‑1 (MQ=255) gAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa < 1:513989/69‑1 (MQ=255) gAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa < 1:1131766/69‑1 (MQ=255) gAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa < 1:2525071/69‑1 (MQ=255) gAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa < 1:2457228/69‑1 (MQ=255) gAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa < 1:1264407/69‑1 (MQ=255) gAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa < 1:237546/69‑1 (MQ=255) gAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa < 1:2319572/69‑1 (MQ=255) gAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa < 1:22482/69‑1 (MQ=255) gAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa < 1:86150/69‑1 (MQ=255) gAGGTTCCGGAAGATCGTGCTGCCCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa < 1:981471/69‑1 (MQ=255) aGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCTGGTGTAGGTCGTAAAACAGCCAACGTCGTa < 1:1917544/68‑1 (MQ=255) gTTCCGGGAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa < 1:2453568/69‑1 (MQ=255) gTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa < 1:3014795/69‑1 (MQ=255) gTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa < 1:1095858/69‑1 (MQ=255) gTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa < 1:1518378/69‑1 (MQ=255) gTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa < 1:18220/69‑1 (MQ=255) gTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa < 1:806535/69‑1 (MQ=255) ttCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa > 1:316670/1‑68 (MQ=255) ttCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa < 1:2025755/68‑1 (MQ=255) ttCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa < 1:871516/68‑1 (MQ=255) ccGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAac > 1:3419881/1‑69 (MQ=255) ctgcGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCt > 1:3201873/1‑69 (MQ=255) gcgcTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTg < 1:1090486/68‑1 (MQ=255) tGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGAc < 1:159355/69‑1 (MQ=255) tGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGAc < 1:176119/69‑1 (MQ=255) tGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGAc < 1:1544071/69‑1 (MQ=255) tGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGAc < 1:353748/69‑1 (MQ=255) tGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGAc < 1:2652711/69‑1 (MQ=255) cTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCt < 1:1757263/69‑1 (MQ=255) tGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCt < 1:2030818/68‑1 (MQ=255) | TCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCT > minE/1039497‑1039621 |
| 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: ‑ |