| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,222,231 | G→C | 14.0% | G269A (GGT→GCT) | yebT → | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,222,231 | 0 | G | C | 14.0% | 140.9 / 14.1 | 57 | G269A (GGT→GCT) | yebT | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base G (29/20); new base C (5/3); total (34/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01 | |||||||||||
GGTTAACGGTGCGATTGCCTTCGATTCACCAGAAGAGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGTAATAATAAAACTGGAACTGCCGAGTGGGGCCGGA > minE/1222163‑1222301 | ggTTAACGGTGCGATTGCCTTCGATTCACCAGAAGAGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTc > 1:89793/1‑71 (MQ=255) gTTAACGGTGCGATTGCCTTCGATTCACCAGAAGAGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGt < 1:3012401/69‑1 (MQ=255) gTTAACGGTGCGATTGCCTTCGATTCACCAGAAGAGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCt < 1:601374/71‑1 (MQ=255) aaCGGTGCGATTGCCTTCGATTCACCAGAAGAGTCGAAACCTGCCGAGGCGGAAGATACCTTTGCTCTg < 1:3701210/69‑1 (MQ=255) aaCGGTGCGATTGCCTTCCATTCACCAGAAGAGTCGAAACCTGCCGAGGCGGAAGATACCTTTGCTCTg < 1:3174101/69‑1 (MQ=255) cGGTGCGATTGCCTTCGATTCACCAGAAGAGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGt < 1:2609971/64‑1 (MQ=255) cGGTGCGATTGCCTTCGATTCACCAGAAGAGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCTGTATg < 1:2582780/71‑1 (MQ=255) gATTGCCTTCGATTCACCAGAAGAGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGt > 1:3531365/1‑58 (MQ=255) ttGCCTTCGATTCACCAGAAGAGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTg > 1:3343197/1‑71 (MQ=255) ttGCCTTCGATTCACCAGAAGAGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTg > 1:2421395/1‑71 (MQ=255) gCCTTCGATTCACCAGAAGAGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCt > 1:204769/1‑56 (MQ=255) tCGATTCACCAGAAGAGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCca > 1:1615051/1‑70 (MQ=255) cGATTCACCAGAAGAGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCca > 1:1784352/1‑69 (MQ=255) cACCAGAAGAGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCt > 1:3149798/1‑58 (MQ=255) cACCAGAAGAGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCAcagccag > 1:2012804/1‑71 (MQ=255) ccAGAAGAGTCGAAACCTGCCGAGGCGGAAGATACCTTTGCTCTGTATGAAGATCTGGCCCACAGCCAGCg < 1:3645716/71‑1 (MQ=255) cAGAAGAGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGt > 1:3923610/1‑71 (MQ=255) cAGAAGAGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGt > 1:2609296/1‑71 (MQ=255) agaagaGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCcaca > 1:3363986/1‑62 (MQ=255) agaagaGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCcaca > 1:3095276/1‑62 (MQ=255) agaagaGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCcaca > 1:2365330/1‑62 (MQ=255) agaagaGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTg < 1:3627540/71‑1 (MQ=255) agaagaGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTg > 1:1759276/1‑71 (MQ=255) agaagaGTCGAAACCTGCCGAGGCGGAAGATACCTTTGCTCTGTATGAAGAt > 1:1460574/1‑52 (MQ=255) gaagaGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTgg < 1:627105/71‑1 (MQ=255) gaagaGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTgg < 1:535561/71‑1 (MQ=255) agagTCGAAACCTGCCGAGGCGGAAGATACCTTTgg < 1:2405093/36‑1 (MQ=255) gTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGtaa > 1:355122/1‑71 (MQ=255) gTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGtaa < 1:1734466/71‑1 (MQ=255) gAAACCTGCCGAGGCGGAAGATACCTTTGCTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGtaataa > 1:2645244/1‑71 (MQ=255) aaaCCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGaa > 1:2365724/1‑39 (MQ=255) aCCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGtaataata > 1:2712740/1‑70 (MQ=255) ccTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGtaataataa > 1:2255141/1‑70 (MQ=255) ccTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGTAATAATaaa < 1:2935928/71‑1 (MQ=255) gCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGcc < 1:1837208/42‑1 (MQ=255) ccGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGTAATAATAAAAc < 1:1800679/69‑1 (MQ=255) ccGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGTAATAATAAAAc < 1:2850219/69‑1 (MQ=255) ccGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGTAATAATAAAACt < 1:1841716/70‑1 (MQ=255) cGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCca > 1:3874167/1‑42 (MQ=255) cGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCca > 1:1863155/1‑42 (MQ=255) cGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGTAATAATAAAACTgg < 1:2958736/71‑1 (MQ=255) cGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGTAATAATAAAACTgg < 1:1529651/71‑1 (MQ=255) aGGCGGAAGATACCTTTGGTCTTTATGAAGATCTGGCCCACAGCCAGCGTGGCGTAATAATAAAACTGGa > 1:2531300/1‑70 (MQ=255) aGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGTAATAATAAAACTGGa > 1:3919568/1‑70 (MQ=255) cGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGtaataataa < 1:18716/60‑1 (MQ=255) gAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGTAATAATAAAACTGGAACTGcc > 1:346391/1‑71 (MQ=255) aTACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGtaa < 1:894799/48‑1 (MQ=255) aTACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGtaa < 1:3729817/48‑1 (MQ=255) tACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGtaa < 1:749300/47‑1 (MQ=255) tACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGTAATAATAAAACTGGAACTGCCGAGTg > 1:858644/1‑71 (MQ=255) tACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGTAATAATAAAACTGGAACTGCCGAGTg > 1:1841631/1‑71 (MQ=255) ccTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGTAATAATAAAACTGGAACTGCCGAGTgg < 1:1016456/70‑1 (MQ=255) cTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTg > 1:3199924/1‑38 (MQ=255) ttGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCg > 1:1870787/1‑39 (MQ=255) ttGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGTAATAATAAAACTGGAACt > 1:1110109/1‑58 (MQ=255) ttGCTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGt > 1:1966714/1‑40 (MQ=255) ttGCTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGt > 1:199555/1‑40 (MQ=255) ttGCTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGt > 1:3228839/1‑40 (MQ=255) gTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGTAATAATAAAACTGGAACTGCCGAGTGGGGCCGGa < 1:3135126/71‑1 (MQ=255) | GGTTAACGGTGCGATTGCCTTCGATTCACCAGAAGAGTCGAAACCTGCCGAGGCGGAAGATACCTTTGGTCTGTATGAAGATCTGGCCCACAGCCAGCGTGGCGTAATAATAAAACTGGAACTGCCGAGTGGGGCCGGA > minE/1222163‑1222301 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |