| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,016,905 | G→A | 32.3% | P213P (CCG→CCA) | manA → | mannose‑6‑phosphate isomerase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,016,905 | 0 | G | A | 32.3% | 27.8 / 16.0 | 37 | P213P (CCG→CCA) | manA | mannose‑6‑phosphate isomerase |
| Reads supporting (aligned to +/- strand): ref base G (6/19); new base A (6/6); total (12/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.46e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01 | |||||||||||
TGAAGAAAAATCCCGCGCGCTGGCGATTTTAAAATCGGCCCTCGATAGCCAGCAGGGTGAACCGTGGC‑AAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCGGTCTGTTCTCCCCGCT > minE/1016842‑1016967 | tGAAGAAAAATCCCGCGCGCTGGCGATTTTAAAATCGGCCCTCGATAGCCAGCAGGGTGAACCGTGGC‑a > 1:263585/1‑69 (MQ=255) aaaaaTCCCGCGCGCTGGCGATTTTAAAATCGGCCCTCGATAGCCAGCAGGGTGAACCGTGGC‑AAACGa < 1:3108541/69‑1 (MQ=255) aaaaaTCCCGCGCGCTGGCGATTTTAAAATCGGCCCTCGATAGCCAGCAGGGTGAACCGTGGC‑AAACGa < 1:7274/69‑1 (MQ=255) cTGGCGATTTTAAAATCGGCCCTCGATAGCCAGCAGGGTGAACCGTGGC‑AAACgt < 1:833981/55‑2 (MQ=255) gCGATTTTAAAATCGGCCCTCGACAGCCAGCAGGGTGAACCATGGC‑AAACGATTCGTTTAATTTCTg < 1:1177144/67‑1 (MQ=255) gCGATTTTAAAATCGGCCCTCGACAGCCAGCAGGGTGAACCATGGC‑AAACGATTCGTTTAATTTCTg < 1:2755133/67‑1 (MQ=255) gCGATTTTAAAATCGGCCCTCGACAGCCAGCAGGGTGAACCATGGC‑AAACGATTCGTTTAATTTCTg < 1:1447953/67‑1 (MQ=255) gCGATTTTAAAATCGGCCCTCGACAGCCAGCAGGGTGAACCATGGC‑AAACGATTCGTTTAATTTCTg < 1:1091704/67‑1 (MQ=255) gCGATTTTAAAATCGGCCCTCGACAGCCAGCAGGGTGAACCATGGC‑AAACGATTCGTTTAATTTCTg < 1:1377963/67‑1 (MQ=255) tttAAAATCGGCCCTCGATAGCCAGCAGGGTGAACCGTGGC‑AAACGATTCGTTTAATTTCTGAATTTTa < 1:1800306/69‑1 (MQ=255) tttAAAATCGGCCCTCGATAGCCAGCAGGGTGAACCGTGGC‑AAACGATTCGTTTAATTTCTGAATTTTa < 1:1930464/69‑1 (MQ=255) tttAAAATCGGCCCTCGATAGCCAGCAGGGTGAACCGTGGC‑AAACGATTCGTTTAATTTCTGAATTTTa < 1:3371217/69‑1 (MQ=255) aaaTCGGCCCTCGATAGCCAGCAGGGTGAACCGTGGC‑AAACGATTCGTTTAATTTCTGAATTTTAccc < 1:2238006/68‑1 (MQ=255) ccTCGACAGCCAGCAGGGTGAAGCATGGC‑AAACGATTCGTTTAATTtgtg < 1:3217766/50‑1 (MQ=255) gCCAGCAGGGTGAACCGTGGCAAAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCGGTCTg < 1:1536560/69‑1 (MQ=255) gCCAGCAGGGTGAACCGTGGCAAAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCGGTCTg < 1:2376680/69‑1 (MQ=255) gCCAGCAGGGTGAACCGTGGCAAAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCGGTCTg < 1:2080217/69‑1 (MQ=255) agcagGGTGAACCGTGGC‑AAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCGGTCTGTtct > 1:1409169/1‑69 (MQ=255) agcagGGTGAACCGTGGC‑AAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCGGTCTGTtct > 1:811867/1‑69 (MQ=255) agcagGGTGAACCGTGGC‑AAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCGGTCTGTtct > 1:1866334/1‑69 (MQ=255) gTGAACCGTGGC‑AAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCGGTCTGTTCTCCCCGc < 1:433602/69‑1 (MQ=255) gTGAACCGTGGC‑AAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCGGTCTGTTCTCCCCGc < 1:3221359/69‑1 (MQ=255) gTGAACCGTGGC‑AAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCGGTCTGTTCTCCCCGc < 1:816868/69‑1 (MQ=255) gTGAACCGTGGC‑AAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCGGTCTGTTCTCCCCGc < 1:893984/69‑1 (MQ=255) gTGAACCGTGGC‑AAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCGGTCTGTTCTCCCCGc < 1:2916831/69‑1 (MQ=255) gTGAACCGTGGC‑AAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCGGTCTGTTCTCCCCGc < 1:2041204/69‑1 (MQ=255) gTGAACCGTGGC‑AAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCGGTCTGTTCTCCCCGc < 1:184268/69‑1 (MQ=255) gTGAACCGTGGC‑AAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCGGTCTGTTCTCCCCGc < 1:1258339/69‑1 (MQ=255) tGAACCGTGGC‑AAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCGGTCTGTTCTCCCCGCt < 1:2381947/69‑1 (MQ=255) tGAACCGTGGC‑AAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCGGTCTGTTCTCCCCGCt > 1:225369/1‑69 (MQ=255) tGAACCGTGGC‑AAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCGGTCTGTTCTCCCCGCt > 1:1937706/1‑69 (MQ=255) gAACCATGGC‑AAACGATTCGTTTAATTTCTGAATTTTAccc > 1:1494272/1‑41 (MQ=255) gAACCATGGC‑AAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCggcc > 1:2882948/1‑53 (MQ=255) gAACCATGGC‑AAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCggcc > 1:312035/1‑53 (MQ=255) gAACCATGGC‑AAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCggcc > 1:2606266/1‑53 (MQ=255) gAACCATGGC‑AAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCggcc > 1:576608/1‑53 (MQ=255) gAACCATGGC‑AAACGATTCGTTTAATTTCAGAATTTTACCCGGAAGACAGCggcc > 1:284010/1‑53 (MQ=255) | TGAAGAAAAATCCCGCGCGCTGGCGATTTTAAAATCGGCCCTCGATAGCCAGCAGGGTGAACCGTGGC‑AAACGATTCGTTTAATTTCTGAATTTTACCCGGAAGACAGCGGTCTGTTCTCCCCGCT > minE/1016842‑1016967 |
| 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: ‑ |