| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,852,805 | 0 | T | G | 52.0% | 19.0 / 28.9 | 25 | V785G (GTT→GGT) | mutS | DNA mismatch repair protein MutS |
| Reads supporting (aligned to +/- strand): ref base T (4/8); new base G (13/0); total (17/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.58e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGG > NZ_CP009273/2852719‑2852880 | tGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGttat > 2:549413/1‑90 (MQ=255) tttATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGGCGGCGGGCCAAAAAAGGGTATTaaa > 1:674299/1‑89 (MQ=255) cacaGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCa < 1:513386/90‑1 (MQ=255) cacaGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCa < 1:549413/90‑1 (MQ=255) cacaGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCa < 1:572973/90‑1 (MQ=255) cacaGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCa < 1:180264/90‑1 (MQ=255) acaGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGGGGGCCAAAAGAGGGTATTAAACGCGCAc > 2:572973/1‑90 (MQ=255) acaGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGGTGCAGCTCTGGCAGGCGGGGCCAAAGAAGGTATTAAacgcgccc > 1:539284/1‑88 (MQ=255) aGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGGGGCAAAAAAGGGTATTTAGCGCGCCCgg > 1:75533/1‑90 (MQ=255) gCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGc > 1:34613/1‑90 (MQ=255) gCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGGTATTAAGCGCGCACGGc > 2:150142/1‑90 (MQ=255) gCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGGTATTAAACGCGCACGGc > 2:757387/1‑90 (MQ=255) tGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCaaa > 2:419205/1‑90 (MQ=255) tGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAAAGGGTATTAAGCGCGCACGGCaaa > 1:945778/1‑90 (MQ=255) cAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGGGCCAAAAGAGGGTATTAAGCGCGCACGGCAAAAg > 2:858292/1‑90 (MQ=255) aGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGc < 2:584301/90‑1 (MQ=255) tGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTTGCCGGCGTGCCAAAAGAGGTTATTAAGCGCCCACGGCCAAAGCTGCg > 2:111238/1‑90 (MQ=255) tGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGGGGGCCAAAAAAGGGTATTAAACGCGCCCGGCAAAAGCTGCg > 1:916873/1‑90 (MQ=255) tGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGGTATTAAGCGCGCACGGCAAAAGCTGCg > 2:733595/1‑90 (MQ=255) gcgcGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGGAGGCGGGCCAAAAGAGGGTATTAAACGCGCCCCGCAAAAGCTGCggg > 2:636971/1‑88 (MQ=255) gcgcGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGGGGGCCCAAAAAGGGTATTAAACGCGCCCGGGAAAAGCTGCGTg > 1:448569/1‑90 (MQ=255) gcgcGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGGGCCAAAAGAGGGTATTAAGCGCGCACGGCAAAAACTGCGTg > 1:515962/1‑90 (MQ=255) cggcgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGc < 2:202808/90‑1 (MQ=255) tGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAg > 1:195857/1‑90 (MQ=255) ggCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGAt < 2:34613/90‑1 (MQ=255) cGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATgg < 2:715011/90‑1 (MQ=255) | TGCCTTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGG > NZ_CP009273/2852719‑2852880 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |