| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,833,787 | A→G | 64.5% | E367E (GAA→GAG) | groL → | Cpn60 chaperonin GroEL, large subunit of GroESL |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,833,787 | 0 | A | G | 64.5% | 39.5 / 31.8 | 59 | E367E (GAA→GAG) | groL | Cpn60 chaperonin GroEL, large subunit of GroESL |
| Reads supporting (aligned to +/- strand): ref base A (15/6); new base G (28/10); total (43/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01 | |||||||||||
TTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAGGCGGCGTTGCAGTTATCAAAGTGGGTGCTGCTACCGAAGTTGAAATG > minE/2833723‑2833853 | ttGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGcg > 1:1165711/1‑67 (MQ=255)ttGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGcg > 1:254793/1‑67 (MQ=255)ttGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGcg > 1:2214113/1‑67 (MQ=255)ttGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGcg > 1:2003713/1‑67 (MQ=255)ttGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGcg > 1:1976394/1‑67 (MQ=255)ttGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGcg > 1:2810542/1‑67 (MQ=255)ttGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGcg > 1:339685/1‑67 (MQ=255)ttGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGcg > 1:474411/1‑67 (MQ=255)ttGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGcg > 1:565235/1‑67 (MQ=255) gCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTa > 1:3234524/1‑69 (MQ=255) aTCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGa < 1:1935393/67‑1 (MQ=255) cggTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACt > 1:3385208/3‑69 (MQ=255) ccGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACt > 1:2451392/1‑69 (MQ=255) ccGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACt > 1:248588/1‑69 (MQ=255) ccGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACt > 1:2513735/1‑69 (MQ=255) ccGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACt > 1:2522357/1‑69 (MQ=255) ccGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACt > 1:130523/1‑69 (MQ=255) ccGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACt > 1:60474/1‑69 (MQ=255) ccGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACt > 1:2569030/1‑69 (MQ=255) ccGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACt > 1:2573031/1‑69 (MQ=255) ccGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACt > 1:2600784/1‑69 (MQ=255) ccGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACt > 1:2795902/1‑69 (MQ=255) ccGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACt > 1:2803589/1‑69 (MQ=255) ccGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACt > 1:3400626/1‑69 (MQ=255) ccGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACt > 1:2892465/1‑69 (MQ=255) ccGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACt > 1:1982832/1‑69 (MQ=255) ccGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACt > 1:2017978/1‑69 (MQ=255) ccGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACt > 1:1620796/1‑69 (MQ=255) ccGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACt > 1:1618208/1‑69 (MQ=255) ccGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACt > 1:1439983/1‑69 (MQ=255) tGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGa > 1:2471018/1‑53 (MQ=255) tGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGa > 1:3246170/1‑53 (MQ=255) tGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGa > 1:1908013/1‑53 (MQ=255) tGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGa > 1:1721041/1‑53 (MQ=255) tGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGa > 1:2150638/1‑53 (MQ=255) tGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGa > 1:2205341/1‑53 (MQ=255) tGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGa > 1:2221633/1‑53 (MQ=255) tGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGa > 1:612064/1‑53 (MQ=255) tGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGa > 1:2248292/1‑53 (MQ=255) tGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTg < 1:2334492/58‑1 (MQ=255) tGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTg < 1:3269688/58‑1 (MQ=255) tGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAGGCGGCGt > 1:787337/1‑69 (MQ=255) tGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAGGCGGCGt > 1:2992397/1‑69 (MQ=255) tGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAGGCGGCGt > 1:2972672/1‑69 (MQ=255) gaagaaGCAACTTCTGACTACGACAGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAGGCGGCGt > 1:261931/1‑68 (MQ=255) agCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAGGCGGCGTTGCAGt > 1:3405179/1‑69 (MQ=255) cGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACTGGCAGGCGGCGTTGCAGtt < 1:2946368/55‑1 (MQ=255) cGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACTGGCAGGCGGCGTTGCAGtt < 1:2294366/55‑1 (MQ=255) cGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACTGGCAGGCGGCGTTGCAGtt < 1:622859/55‑1 (MQ=255) cGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACTGGCAGGCGGCGTTGCAGtt < 1:1343606/55‑1 (MQ=255) cGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACTGGCAGGCGGCGTTGCAGtt < 1:1613554/55‑1 (MQ=255) cGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACTGGCAGGCGGCGTTGCAGtt < 1:2704980/55‑1 (MQ=255) cGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACTGGCAGGCGGCGTTGCAGtt < 1:3284938/55‑1 (MQ=255) cGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACTGGCAGGCGGCGTTGCAGtt < 1:2899654/55‑1 (MQ=255) cGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACTGGCAGGCGGCGTTGCAGtt < 1:1781684/55‑1 (MQ=255) cGACCGTGAAAAACTGCAGGAGCGCGTAGCGAAACTGGCAGGCGGCGTTGCAGtt < 1:3117967/55‑1 (MQ=255) gCAGGAACGCGTAGCGAAACTGGCAGGCGGCGTTGCAGTTATCAAAGTGGGTGCTGCTACCGAAGTtga < 1:2107622/69‑1 (MQ=255) gCAGGAACGCGTAGCGAAACTGGCAGGCGGCGTTGCAGTTATCAAAGTGGGTGCTGCTACCGAAGTtga < 1:223499/69‑1 (MQ=255) gAACGCGTAGCGAAACTGGCAGGCGGCGTTGCAGTTATCAAAGTGGGTGCTGCTACCGAAGTtgaaatg < 1:2688925/69‑1 (MQ=255) aCGCGTAGCGAAACTGGCAGGCGGCGTTGCAGTTATCAAAGTGGGTGCTGCTACCGAAGTtgaaat > 1:3354759/1‑66 (MQ=255) aCGCGTAGCGAAACTGGCAGGCGGCGTTGCAGTTATCAAAGTGGGTGCTGCTACCGAAGTtgaaa > 1:1897090/1‑65 (MQ=255) | TTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAACGCGTAGCGAAACTGGCAGGCGGCGTTGCAGTTATCAAAGTGGGTGCTGCTACCGAAGTTGAAATG > minE/2833723‑2833853 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |