Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,733,864 | G→A | 45.9% | S184S (AGC→AGT) | pdxJ ← | pyridoxine 5'‑phosphate synthase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,733,864 | 0 | G | A | 45.9% | 3.7 / 26.4 | 37 | S184S (AGC→AGT) | pdxJ | pyridoxine 5'‑phosphate synthase |
Reads supporting (aligned to +/- strand): ref base G (6/14); new base A (14/3); total (20/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.57e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GGCGGCAATGGCTTTCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGGCTTGCGGCAAAGGTCGCGGCTTTGGCGATACGCGCCAGCTCTTGCGCCTGTTCGGCGTCAGTTTTGGC > minE/1733801‑1733932 | ggcggcAATGGCTTTCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGACtag > 1:1065702/1‑66 (MQ=255) ggcggcAATGGCTTTCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGACtag > 1:786683/1‑66 (MQ=255) ggcggcAATGGCTTTCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGACtag > 1:597175/1‑66 (MQ=255) ggcggcAATGGCTTTCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGACtag > 1:447709/1‑66 (MQ=255) ggcggcAATGGCTTTCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGACtag > 1:1285570/1‑66 (MQ=255) ggcggcAATGGCTTTCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGACtag > 1:3159655/1‑66 (MQ=255) ggcggcAATGGCTTTCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGACtag > 1:2558802/1‑66 (MQ=255) ggcggcAATGGCTTTCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGACtag > 1:2494032/1‑66 (MQ=255) ggcggcAATGGCTTTCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGACtag > 1:170665/1‑66 (MQ=255) ggcggcAATGGCTTTCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGACtag > 1:2491486/1‑66 (MQ=255) ggcggcAATGGCTTTCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGACtag > 1:1873513/1‑66 (MQ=255) ggcggcAATGGCTTTCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGACtag > 1:1000414/1‑66 (MQ=255) ggcggcAATGGCTTTCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGACtag > 1:2235628/1‑66 (MQ=255) ggcggcAATGGCTTTCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGAACGAGACtag > 1:3365982/1‑66 (MQ=255) gcggcAATGGCTTTCACGTTGTGATATGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGGCTTGCg < 1:2490936/69‑1 (MQ=255) gcggcAATGGCTTTCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGGCTTGCg < 1:2044174/69‑1 (MQ=255) gcggcAATGGCTTTCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGGCTTGCg < 1:1280498/69‑1 (MQ=255) gcggcAATGGCTTTCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGGCTTGCg < 1:1283429/69‑1 (MQ=255) gcggcAATGGCTTTCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGGCTTGCg < 1:2824772/69‑1 (MQ=255) ttCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGGCTTGCGGCAAAGGTcgcg < 1:955171/69‑1 (MQ=255) ttCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGGCTTGCGGCAAAGGTcgcg < 1:1205117/69‑1 (MQ=255) ttCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGGCTTGCGGCAAAGGTcgcg < 1:129055/69‑1 (MQ=255) ttCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGGCTTGCGGCAAAGGTcgcg < 1:2960095/69‑1 (MQ=255) cGTGTCCGGCGTTAACTTTCAGACCGAGGCTTGCGGCAAAGGTCGCGGCTTTGGCGATACGCGCCAGct < 1:2617863/69‑1 (MQ=255) cGTGTCCGGCGTTAACTTTCAGACCGAGGCTTGCGGCAAAGGTCGCGGCTTTGGCGATACGCGCCAGct < 1:2652494/69‑1 (MQ=255) cGTGTCCGGCGTTAACTTTCAGACCGAGGCTTGCGGCAAAGGTCGCGGCTTTGGCGATACGCGCCAGct < 1:2477044/69‑1 (MQ=255) gtgtCCGGCGTTAACTTTCAGACCGAGGCTTGCGGCAAAGGTCGCGGCTTTGGCGATACGCGCCAGct < 1:2355328/68‑1 (MQ=255) gtgtCCGGCGTTAACTTTCAGACCGAGGCTTGCGGCAAAGGTCGCGGCTTTGGCGATACGCGCCAGct < 1:347949/68‑1 (MQ=255) cGGCGTTAACTTTCAGACCGAGGCTTGCGGCAAAGGTCGCGGCTTTGGCGATACGCGCCAGCTCTTgcg > 1:242214/1‑69 (MQ=255) tCAGACCGAGGCTTGCGGCAAAGGTCGCGGCTTTGGCGATACGCGCCAGCTCTTGCGCCTGTTCGGc > 1:3294934/1‑67 (MQ=255) tCAGACCGAGGCTTGCGGCAAAGGTCGCGGCTTTGGCGATACGCGCCAGCTCTTGCGCCTGTTCGGc > 1:2325202/1‑67 (MQ=255) tCAGACCGAGGCTTGCGGCAAAGGTCGCGGCTTTGGCGATACGCGCCAGCTCTTGCGCCTGTTCGGc > 1:1908548/1‑67 (MQ=255) tCAGACCGAGGCTTGCGGCAAAGGTCGCGGCTTTGGCGATACGCGCCAGCTCTTGCGCCTGTTCGGc > 1:1582007/1‑67 (MQ=255) tCAGACCGAGGCTTGCGGCAAAGGTCGCGGCTTTGGCGATACGCGCCAGCTCTTGCGCCTGTTCGGc > 1:1717008/1‑67 (MQ=255) tCAGACCGAGACTAGTGGCAAAGGTCGAGGCTTTGGCGATa < 1:3290805/41‑1 (MQ=25) tCAGACCGAGACTAGTGGCAAAGGTCGAGGCTTTGGCGATa < 1:61041/41‑1 (MQ=25) tCAGACCGAGACTAGTGGCAAAGGTCGAGGCTTTGGCGATa < 1:888857/41‑1 (MQ=25) ggCTTGCGGCAAAGGTCGCGGCTTTGGCGATACGCGCCAGCTCTTGCGCCTGTTCGGCGTCAGTTTTgg < 1:1685848/69‑1 (MQ=255) gCTTGCGGCAAAGGTCGCGGCTTTGGCGATACGCGCCAGCTCTTGCGCCTGTTCGGCGTCAGTTTTGGc < 1:1382542/69‑1 (MQ=255) gCTTGCGGCAAAGGTCGCGGCTTTGGCGATACGCGCCAGCTCTTGCGCCTGTTCGGCGTCAGTTTTGGc < 1:2570186/69‑1 (MQ=255) | GGCGGCAATGGCTTTCACGTTGTGATAGGTCAGACCGTGTCCGGCGTTAACTTTCAGACCGAGGCTTGCGGCAAAGGTCGCGGCTTTGGCGATACGCGCCAGCTCTTGCGCCTGTTCGGCGTCAGTTTTGGC > minE/1733801‑1733932 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |