| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 611,667 | C→T | 53.1% | P555P (CCG→CCA) | ycaO ← | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 611,667 | 0 | C | T | 53.1% | ‑2.4 / 24.6 | 32 | P555P (CCG→CCA) | ycaO | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base C (9/6); new base T (8/9); total (17/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.02e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CGCTGCAGCTTTTCGTAGGCCTTCAACAACGACTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCTCATTGCCGCACTGGCGGCTTCTACGG > minE/611599‑611725 | cgcTGCAGCTTTTCGTAGGCCTTCAACAACGACTGATGTGCAGCAAACGCGTGCAGATCGCTATCTAcc < 1:67275/69‑1 (MQ=255) acGACTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTcgc > 1:1365801/1‑69 (MQ=255) acGACTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTcgc > 1:2979467/1‑69 (MQ=255) acGACTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTcgc > 1:287714/1‑69 (MQ=255) acGACTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTcgc > 1:2544832/1‑69 (MQ=255) cTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCt < 1:1033378/69‑1 (MQ=255) cTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCt < 1:790856/69‑1 (MQ=255) cTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCt < 1:507658/69‑1 (MQ=255) cTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCt < 1:426740/69‑1 (MQ=255) cTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCt < 1:2340131/69‑1 (MQ=255) cagcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCg < 1:1092155/69‑1 (MQ=255) cagcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCg < 1:985842/69‑1 (MQ=255) cagcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCg < 1:314063/69‑1 (MQ=255) cagcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCg < 1:3128333/69‑1 (MQ=255) cagcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCg < 1:2656571/69‑1 (MQ=255) cagcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCg < 1:2256487/69‑1 (MQ=255) cagcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCg < 1:1860650/69‑1 (MQ=255) cagcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCg < 1:1635585/69‑1 (MQ=255) cagcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCg < 1:1539667/69‑1 (MQ=255) gcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCGCa > 1:2516240/1‑69 (MQ=255) gcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCGCa > 1:2680260/1‑69 (MQ=255) gcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCGCa > 1:2399793/1‑69 (MQ=255) gcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCGCa > 1:2524026/1‑69 (MQ=255) gcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCGCa > 1:386582/1‑69 (MQ=255) gcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCGCa > 1:645264/1‑69 (MQ=255) gcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCGCa > 1:1490521/1‑69 (MQ=255) gcaAACGCGTGCAGATCGCTATCTACTGGTTGCAAGCCATAAAACGCCGCTTCGCCGCTCATTGCCGCa > 1:1385783/1‑69 (MQ=255) cGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCTCATTGCCGCACTGGCGGCTTCt > 1:886684/1‑65 (MQ=255) cGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCTCATTGCCGCACTGGCGGCTTCTa > 1:110650/1‑66 (MQ=255) cGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCTCATTGCCGCACTGGCGGCTTCTa > 1:3070453/1‑66 (MQ=255) cGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCTCATTGCCGCACTGGCGGCTTCTACgg > 1:3350109/1‑69 (MQ=255) cGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCTCATTGCCGCACTGGCGGCTTCTACgg > 1:2630909/1‑69 (MQ=255) cGCTAACTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCTCATTGCCGCACTGGCGGCTTCTa < 1:2249221/66‑1 (MQ=255) | CGCTGCAGCTTTTCGTAGGCCTTCAACAACGACTGATGTGCAGCAAACGCGTGCAGATCGCTATCTACCGGTTGCAGGCCGTAAAACGCCGCTTCGCCGCTCATTGCCGCACTGGCGGCTTCTACGG > minE/611599‑611725 |
| 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: ‑ |