| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,023,923 | A→G | 45.0% | intergenic (‑17/+94) | ygiQ ← / ← sufI | conserved hypothetical protein/repressor protein for FtsI |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,023,923 | 0 | A | G | 45.0% | 6.9 / 48.4 | 60 | intergenic (‑17/+94) | ygiQ/sufI | conserved hypothetical protein/repressor protein for FtsI |
| Reads supporting (aligned to +/- strand): ref base A (10/23); new base G (17/10); total (27/33) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.85e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCA > minE/2023871‑2023987 | gTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTAAAAAAATAATCaa > 1:288187/1‑68 (MQ=255)gTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTAAAAAAATAATCaa > 1:2160371/1‑68 (MQ=255)gTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTAAAAAAATAATCaa > 1:496829/1‑68 (MQ=255)gTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTAAAAAAATAATCaa > 1:2237908/1‑68 (MQ=255)gTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTAAAAAAATAATCaa > 1:3375875/1‑68 (MQ=255)gTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTAAAAAAATAATCaa > 1:3371811/1‑68 (MQ=255)gTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTAAAAAAATAATCaa > 1:3147489/1‑68 (MQ=255)gTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTAAAAAAATAATCaa > 1:1601196/1‑68 (MQ=255)gTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTAAAAAAATAATCaa > 1:3009029/1‑68 (MQ=255)gTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTAAAAAAATAATCaa > 1:2668501/1‑68 (MQ=255)gTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTAAAAAAATAATCaa > 1:834945/1‑68 (MQ=255)gTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTAAAAAAATAATCaa > 1:954215/1‑68 (MQ=255)gTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTAAAAAAATAATCaa > 1:956943/1‑68 (MQ=255)gTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTAAAAAAATAATCaa > 1:1151678/1‑68 (MQ=255)gTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTAAAAAAATAATCaa > 1:2783094/1‑68 (MQ=255)gTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTAAAAAAATAATCaa > 1:1033280/1‑68 (MQ=255)gTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTAAAAAAATAATCa > 1:1845018/1‑67 (MQ=255) cgcgATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTaaaaaaa < 1:2828435/59‑1 (MQ=255) cgcgATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTaaaaaaa < 1:2748606/59‑1 (MQ=255) cgcgATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTaaaaaaa < 1:3111723/59‑1 (MQ=255) cgcgATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTaaaaaaa < 1:275891/59‑1 (MQ=255) cgcgATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTaaaaaaa < 1:2021501/59‑1 (MQ=255) cgcgATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTaaaaaaa < 1:199814/59‑1 (MQ=255) cgcgATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTaaaaaaa < 1:727696/59‑1 (MQ=255) cgcgATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTaaaaaaa < 1:1290061/59‑1 (MQ=255) gcgATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTaaaaaaa < 1:536519/58‑1 (MQ=255) tCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCggg < 1:1607128/69‑1 (MQ=255) tCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCggg < 1:1540958/69‑1 (MQ=255) tCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCggg < 1:1446738/69‑1 (MQ=255) ccGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCgggg > 1:431229/1‑69 (MQ=255) ccGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCgggg > 1:1146386/1‑69 (MQ=255) ggATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGGTaaaaaaa < 1:1429854/47‑1 (MQ=255) agagCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGAtt < 1:1874539/53‑1 (MQ=255) agagCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGAtt < 1:1393028/53‑1 (MQ=255) agagCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGAtt < 1:2753646/53‑1 (MQ=255) gagCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGa > 1:1129796/1‑66 (MQ=255) ttGCTTCCAGAAATGATAAAAAAATAATCAAAGTCCGGGGATTATAAGACGGAaagaaaga > 1:3249767/1‑61 (MQ=255) ttGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATc > 1:961977/1‑67 (MQ=255) ttGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATc > 1:3350127/1‑67 (MQ=255) tCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAgg < 1:2582445/68‑1 (MQ=255) tCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAgg < 1:942849/68‑1 (MQ=255) tCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAgg < 1:3208635/68‑1 (MQ=255) tCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAgg < 1:3216362/68‑1 (MQ=255) tCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAgg < 1:3299475/68‑1 (MQ=255) tCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAgg < 1:249543/68‑1 (MQ=255) tCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAgg < 1:2378292/68‑1 (MQ=255) aGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTa < 1:2596941/67‑1 (MQ=255) aGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTa < 1:1349137/67‑1 (MQ=255) aGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTa < 1:1898312/67‑1 (MQ=255) aGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTa < 1:2498383/67‑1 (MQ=255) aGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTa < 1:2634744/67‑1 (MQ=255) gAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTc < 1:2217546/69‑1 (MQ=255) gAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTc < 1:858465/69‑1 (MQ=255) gAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTc < 1:3165353/69‑1 (MQ=255) gAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATAGAAAGGTATTc < 1:859882/69‑1 (MQ=255) aaTGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTcc > 1:1184209/1‑68 (MQ=255) aaTGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCa > 1:2162872/1‑69 (MQ=255) aaTGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCa < 1:1599639/69‑1 (MQ=255) aTGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTa > 1:1380908/1‑63 (MQ=255) aTGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTca > 1:701791/1‑66 (MQ=255) | GTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCA > minE/2023871‑2023987 |
| 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: ‑ |