| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 555,000 | G→A | 10.7% | D148N (GAC→AAC) | fsaA → | fructose‑6‑phosphate aldolase 1 |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 555,000 | 0 | G | A | 10.7% | 163.1 / 9.4 | 56 | D148N (GAC→AAC) | fsaA | fructose‑6‑phosphate aldolase 1 |
| Reads supporting (aligned to +/- strand): ref base G (26/24); new base A (4/2); total (30/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.75e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01 | |||||||||||
CGGAATATGTTGCGCCTTACGTTAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCAAAAC > minE/554932‑555064 | cGGAATATGTTGCGCCTTACGTTAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGAc > 1:1180218/1‑71 (MQ=255)cGGAATATGTTGCGCCTTACGTTAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGAc > 1:3921690/1‑71 (MQ=255)cGGAATATGTTGCGCCTTACGTTAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGAc > 1:1787975/1‑71 (MQ=255) ggAATATGTTGCGCCTTACGTTAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACt < 1:1158869/71‑1 (MQ=255) ggAATATGTTGCGCCTTACGTTAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACt < 1:2693464/71‑1 (MQ=255) ggAATATGTTGCGCCTTACGTTAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACt < 1:1746343/71‑1 (MQ=255) tGCGCCTTACGTTAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGt > 1:1879095/1‑71 (MQ=255) tACGTTAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGa > 1:3909980/1‑70 (MQ=255) aCGTTAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTa > 1:523348/1‑56 (MQ=255) aCGTTAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGaa > 1:2278373/1‑70 (MQ=255) aCGTTAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGaa > 1:2417041/1‑70 (MQ=255) aCGTTAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGaa > 1:2782146/1‑70 (MQ=255) gTTAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGaa < 1:654137/68‑1 (MQ=255) ttAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAatg < 1:69650/71‑1 (MQ=255) tAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGa > 1:1521196/1‑48 (MQ=255) tAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGa > 1:3132115/1‑48 (MQ=255) tAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGa > 1:2781110/1‑48 (MQ=255) aaTCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAatgc < 1:3469600/70‑1 (MQ=255) aaTCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAatgc < 1:521784/70‑1 (MQ=255) gTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAatgcatgc < 1:2718729/70‑1 (MQ=255) ttGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGcgcc < 1:1968081/70‑1 (MQ=255) tGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGcgccg > 1:3941721/1‑70 (MQ=255) gATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGcgccg > 1:2725363/1‑69 (MQ=255) aTGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAg < 1:633444/71‑1 (MQ=255) tGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAgg > 1:1059207/1‑71 (MQ=255) tGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAgg > 1:3201857/1‑71 (MQ=255) aGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAatg < 1:1412425/53‑1 (MQ=255) aGGGCGGTACCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGa > 1:3664508/1‑69 (MQ=255) gCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTg < 1:1478486/71‑1 (MQ=255) tAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGa < 1:1478291/62‑1 (MQ=255) tAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGa < 1:3979952/62‑1 (MQ=255) tAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGa < 1:2013446/62‑1 (MQ=255) aGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAg > 1:3611598/1‑64 (MQ=255) aTTCAGACTGTGACCAACTTACACCAGTTATTGaa < 1:2672476/35‑1 (MQ=255) aTTCAGACTGTGACCAACTTACACCAGTTATTGaa < 1:2350769/35‑1 (MQ=255) ttCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGa < 1:857114/54‑1 (MQ=255) ttCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGc < 1:3516034/60‑1 (MQ=255) tCAGACTGTGACCAACTTACACCAGTTATTGaaaa > 1:1528927/1‑35 (MQ=255) tCAGACTGTGACCAACTTACACCAGTTATTGaaaa > 1:852695/1‑35 (MQ=255) tCAGACTGTGACCAACTTACACCAGTTATTGaaaa > 1:2160254/1‑35 (MQ=255) tCAGACTGTGACCAACTTACACCAGTTATTGaaaa > 1:322815/1‑35 (MQ=255) cAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCa > 1:2805975/1‑63 (MQ=255) cAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGttt > 1:3870545/1‑71 (MQ=255) cAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGttt > 1:3589906/1‑71 (MQ=255) cAGACTGTGACCGACTTACACCAGTCATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCa > 1:2711096/1‑63 (MQ=255) aCTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCaa < 1:3905651/71‑1 (MQ=255) aCTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCaa < 1:2584650/71‑1 (MQ=255) cTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCaaa < 1:464268/71‑1 (MQ=255) cTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCaa < 1:3016148/70‑1 (MQ=255) tgtgACCGACTTACACCAGTTATTGAAAATGCATGcgcc < 1:2693611/39‑1 (MQ=255) tgtgACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCaa < 1:3261725/69‑1 (MQ=255) tgtgACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCaa < 1:575848/69‑1 (MQ=255) tgACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGc > 1:3446832/1‑51 (MQ=255) tgACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCAAAAc > 1:3040523/1‑70 (MQ=255) cGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCaaaa > 1:400755/1‑65 (MQ=255) cGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCaaaa > 1:4157797/1‑65 (MQ=255) | CGGAATATGTTGCGCCTTACGTTAATCGTATTGATGCTCAGGGCGGTAGCGGCATTCAGACTGTGACCGACTTACACCAGTTATTGAAAATGCATGCGCCGCAGGCGAAAGTGCTGGCAGCGAGTTTCAAAAC > minE/554932‑555064 |
| 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: ‑ |