| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,238,977 | C→T | 37.2% | A418A (GCG→GCA) | zwf ← | glucose‑6‑phosphate dehydrogenase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,238,977 | 0 | C | T | 37.2% | 33.6 / 34.1 | 59 | A418A (GCG→GCA) | zwf | glucose‑6‑phosphate dehydrogenase |
| Reads supporting (aligned to +/- strand): ref base C (21/16); new base T (8/14); total (29/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.80e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01 | |||||||||||
TACAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCAGATGCGTCTGATTAAAGGTTTCTGAATAGCTCAGATCCAGCTTGGTGATTTGCAGGTTATGTTT > minE/1238920‑1239045 | tACAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCa > 1:1740511/1‑62 (MQ=255)tACAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCa > 1:2078740/1‑62 (MQ=255)tACAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCa > 1:2128033/1‑62 (MQ=255)tACAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCa > 1:863655/1‑62 (MQ=255)tACAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCa > 1:2711822/1‑62 (MQ=255)tACAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCa > 1:2816142/1‑62 (MQ=255)tACAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCa > 1:1603918/1‑62 (MQ=255)tACAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCa > 1:2956402/1‑62 (MQ=255)tACAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCa > 1:314780/1‑62 (MQ=255)tACAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCa > 1:1184933/1‑62 (MQ=255) cAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCAGATGCGTCt > 1:2824773/1‑69 (MQ=255) cAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCAGATGCGTCt > 1:2485864/1‑69 (MQ=255) cAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCAGATGCGTCt < 1:3041113/69‑1 (MQ=255) cAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCAGATGCGTCt < 1:2668771/69‑1 (MQ=255) cAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCAGATGCGTCt < 1:2548941/69‑1 (MQ=255) cAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCAGATGCGTCt > 1:2483829/1‑69 (MQ=255) cAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCAGATGCGTCt < 1:561870/69‑1 (MQ=255) cAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCAGATGCGTCt < 1:895784/69‑1 (MQ=255) cAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCAGATGCGTCt > 1:1837413/1‑69 (MQ=255) gAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCAGATGCGTCt < 1:934418/57‑1 (MQ=255) aaaCGTTCATAGGCATCCGCCAGATGCGTCTGATTAAAGGTTTCTGAATAGCTCAGATCCAGCTTGGTg < 1:2789143/69‑1 (MQ=255) aaaCGTTCATAGGCATCCGCCAGATGCGTCTGATTAAAGGTTTCTGAATAGCTCAGATCCAGCTTGGTg < 1:371882/69‑1 (MQ=255) aaaCGTTCATAGGCATCCGCCAGATGCGTCTGATTAAAGGTTTCTGAATAGCTCAGATCCAGCTTGGTg < 1:3296989/69‑1 (MQ=255) aaaCGTTCATAGGCATCCGCCAGATGCGTCTGATTAAAGGTTTCTGAATAGCTCAGATCCAGCTTGGTg < 1:836692/69‑1 (MQ=255) aaaCGTTCATAGGCATCCGCCAGATGCGTCTGATTAAAGGTTTCTGAATAGCTCAGATCCAGCTTGGTg < 1:671885/69‑1 (MQ=255) gTTCATAGGCATGTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGAt < 1:812224/54‑1 (MQ=255) gTTCATAGGCATCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGAt < 1:2949980/54‑1 (MQ=38) gTTCATAGGCATCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGAt < 1:2194741/54‑1 (MQ=38) gTTCATAGGCATCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGAt < 1:950885/54‑1 (MQ=38) gTTCATAGGCATCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGAt < 1:824563/54‑1 (MQ=38) gTTCATAGGCATCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGAt < 1:59431/54‑1 (MQ=38) gTTCATAGGCATCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGAt < 1:3106221/54‑1 (MQ=38) gTTCATAGGCATCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGAt < 1:3099524/54‑1 (MQ=38) gTTCATAGGCATCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGAt < 1:2409029/54‑1 (MQ=38) gTTCATAGGCATCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGAt < 1:2942287/54‑1 (MQ=38) gTTCATAGGCATCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGAt < 1:1735233/54‑1 (MQ=38) gTTCATAGGCATCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGAt < 1:1755411/54‑1 (MQ=38) gTTCATAGGCATCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGAt < 1:1957294/54‑1 (MQ=38) gTTCATAGGCATCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGAt < 1:1889050/54‑1 (MQ=38) aGGCATCCGCCAGATGCGTCTGATTAAAGGTTTCTGAATAGCTCAGATCCAGCTTGGTGATTTGCAGGt < 1:204245/69‑1 (MQ=255) aGGCATCCGCCAGATGCGTCTGATTAAAGGTTTCTGAATAGCTCAGATCCAGCTTGGTGATTTGCAGGt < 1:27087/69‑1 (MQ=255) aGGCATCCGCCAGATGCGTCTGATTAAAGGTTTCTGAATAGCTCAGATCCAGCTTGGTGATTTGCAGGt < 1:355574/69‑1 (MQ=255) aGGCATCCGCCAGATGCGTCTGATTAAAGGTTTCTGAATAGCTCAGATCCAGCTTGGTGATTTGCAGGt < 1:2538543/69‑1 (MQ=255) ggCATCCGCCAGATGCGTCTGATTAAAGGTTTCTGAATAGCTCAGATCCAGCTTGGTGATTTGCAGGtt > 1:868844/1‑69 (MQ=255) ggCATCCGCCAGATGCGTCTGATTAAAGGTTTCTGAATAGCTCAGATCCAGCTTGGTGATTTGCAGGtt > 1:989889/1‑69 (MQ=255) ggCATCCGCCAGATGCGTCTGATTAAAGGTTTCTGAATAGCTCAGATCCAGCTTGGTGATTTGCAGGtt > 1:649619/1‑69 (MQ=255) ggCATCCGCCAGATGCGTCTGATTAAAGGTTTCTGAATAGCTCAGATCCAGCTTGGTGATTTGCAGGtt > 1:3312071/1‑69 (MQ=255) ggCATCCGCCAGATGCGTCTGATTAAAGGTTTCTGAATAGCTCAGATCCAGCTTGGTGATTTGCAGGtt > 1:1540686/1‑69 (MQ=255) ggCATCCGCCAGATGCGTCTGATTAAAGGTTTCTGAATAGCTCAGATCCAGCTTGGTGATTTGCAGGt < 1:2542920/68‑1 (MQ=255) ggCATCCGCCAGATGCGTCTGATTAAAGGTTTCTGAATAGCTCAGATCCAGCTTGGTGATTTGCAGGt > 1:1348102/1‑68 (MQ=255) ggCATCCGCCAGATGCGTCTGATTAAAGGTTTCTGAATAGCTCAGATCCAGCTTGGTGATTTGCAGGt > 1:1325732/1‑68 (MQ=255) aTCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGATCCAGCTTGGTGAt > 1:3111647/1‑57 (MQ=255) aTCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGATCCAGCTTGGTGAt > 1:2520459/1‑57 (MQ=255) aTCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGATCCAGCTTGGTGAt > 1:2377604/1‑57 (MQ=255) aTCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGATCCAGCTTGGTGAt > 1:2207402/1‑57 (MQ=255) aTCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGATCCAGCTTGGTGAt > 1:1893416/1‑57 (MQ=255) aTCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGATCCAGCTTGGTGAt > 1:843067/1‑57 (MQ=255) aTCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGATCCAGCTTGGTGAt > 1:1865542/1‑57 (MQ=255) aTCTGCCAGATGCGTTTGATTAAAAGTTTCTGAATAGCTCAGATCCAGCTTGGTGAt > 1:1285861/1‑57 (MQ=255) cGCCAGATGCGTCTGATTAAAGGTTTCTGAATGGCTCAGATCCAGCTTGGTGATTTGCAGGTTATGttt > 1:94693/1‑69 (MQ=255) cGCCAGATGCGTCTGATTAAAGGTTTCTGAATAGCTCAGATCCAGCTTGGTGATTTGCAGGTTATGttt > 1:1591124/1‑69 (MQ=255) | TACAAACAGTGCCTGAATACCACGCATGGTTTCCAGCAGCAAACGTTCATAGGCATCCGCCAGATGCGTCTGATTAAAGGTTTCTGAATAGCTCAGATCCAGCTTGGTGATTTGCAGGTTATGTTT > minE/1238920‑1239045 |
| 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: ‑ |