Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,017,067 | A→G | 71.6% | A267A (GCA→GCG) | manA → | mannose‑6‑phosphate isomerase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,017,067 | 0 | A | G | 71.6% | 57.5 / 21.1 | 56 | A267A (GCA→GCG) | manA | mannose‑6‑phosphate isomerase |
Reads supporting (aligned to +/- strand): ref base A (2/14); new base G (27/13); total (29/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.60e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01 |
TGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCAAACTCCGATAACGTGCTGCGTGCGGGTCTGACGCCTAAATACATTGA > minE/1017019‑1017114 | tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGtgctgc > 1:827285/1‑68 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGtgctgc > 1:2948451/1‑68 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg < 1:3220009/69‑1 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:2287138/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg < 1:2372813/69‑1 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg < 1:2387420/69‑1 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:2544982/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:254884/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg < 1:2692821/69‑1 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg < 1:2887910/69‑1 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:2907962/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:2992695/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:300387/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg < 1:3189809/69‑1 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:1073844/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg < 1:3223978/69‑1 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg < 1:335846/69‑1 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:3373618/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:3421099/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:365813/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:406059/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:447662/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:506082/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:563423/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg < 1:738647/69‑1 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:972618/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:220098/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:1867915/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:1152738/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg < 1:1805454/69‑1 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:1780788/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:1170815/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg < 1:145480/69‑1 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:149853/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg > 1:1688949/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGAGCtgcg > 1:2523315/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAAATCCGATAACGTGCtgcg > 1:982236/1‑69 (MQ=255) tGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGAGATGGCGAACTCCGATAACGTGCtgcg > 1:1919809/1‑69 (MQ=255) gCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCAAACTCCGATAACGTGCTGCGTGCGGGTCt < 1:2179594/69‑1 (MQ=255) gCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCAAACTCCGATAACGTGCTGCGTGCGGGTCt < 1:1345168/69‑1 (MQ=255) gCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCAAACTCCGATAACGTGCTGCGTGCGGGTCt < 1:1770361/69‑1 (MQ=255) cACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCGAACTCCGATAACGTGCtgcg < 1:101380/59‑1 (MQ=255) aCGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCAAACTCCGATAACGTGCTGCGTGCGGGTCTGa < 1:3171571/69‑1 (MQ=255) aCGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCAAACTCCGATAACGTGCTGCGTGCGGGTCTGa < 1:3047069/69‑1 (MQ=255) aCGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCAAACTCCGATAACGTGCTGCGTGCGGGTCTGa < 1:1346885/69‑1 (MQ=255) aCGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCAAACTCCGATAACGTGCTGCGTGCGGGTCTGa < 1:1259914/69‑1 (MQ=255) aCGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCAAACTCCGATAACGTGCTGCGTGCGGGTCTGa < 1:1216603/69‑1 (MQ=255) aCGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCAAACTCCGATAACGTGCTGCGTGCGGGTCTGa < 1:118075/69‑1 (MQ=255) aCGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCAAACTCCGATAACGTGCTGCGTGCGGGTCTGa < 1:911470/69‑1 (MQ=255) cGCTTACCTGCTAGGCGTGTCGCTGGAAGTGATGGCGAAGTGCGATAACGTGCtgcg < 1:1975152/57‑1 (MQ=255) cGTGGCGCTGGAAGTGATGGCAAACTCCGATAACGTGCTGCGTGCGGGTCTGACGCCTAAATACATTGa > 1:1685188/1‑69 (MQ=255) cGTGGCGCTGGAAGTGATGGCAAACTCCGATAACGTGCTGCGTGCGGGTCTGACGCCTAAATACATTGa > 1:1718250/1‑69 (MQ=255) ggCAAACTCCGATAACGTGCTGCGTGCGGGTCTGACGCCTAAATACATTGa < 1:1709103/51‑1 (MQ=255) ggCAAACTCCGATAACGTGCTGCGTGCGGGTCTGACGCCTAAATACATTGa < 1:1255064/51‑1 (MQ=255) ggCAAACTCCGATAACGTGCTGCGTGCGGGTCTGACGCCTAAATACATTGa < 1:1810462/51‑1 (MQ=255) gCAAACTCCGATAACGTGCTGCGTGCGGGTCTGACGCCTAAATACATTGa < 1:1313258/50‑1 (MQ=255) | TGAAACACCGCACGCTTACCTGCAAGGCGTGGCGCTGGAAGTGATGGCAAACTCCGATAACGTGCTGCGTGCGGGTCTGACGCCTAAATACATTGA > minE/1017019‑1017114 |
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: ‑ |