| 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.2% | 19.1 / 23.4 | 23 | V785G (GTT→GGT) | mutS | DNA mismatch repair protein MutS |
| Reads supporting (aligned to +/- strand): ref base T (5/6); new base G (12/0); total (17/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.58e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAAT > NZ_CP009273/2852723‑2852889 | tttATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAg > 2:104303/1‑90 (MQ=255) tATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGGTATTAAacg > 2:208953/1‑90 (MQ=255) cacaGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCa < 1:127387/90‑1 (MQ=255) acaGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGGTATTAAGCGCGCAc > 2:93398/1‑90 (MQ=255) aGCGTGCAGGATGGCGCGTCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAAAGGTTATTAAACGCGCACgg > 2:405918/1‑90 (MQ=255) aGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGGGCCAAAAGAGGGTATTAAACGCGCACgg > 2:471583/1‑90 (MQ=255) aGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCCCTGGCAGGGGGGCCAAAAGAGGGTATTAAACGCGCACgg > 1:449699/1‑90 (MQ=255) gCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGc > 2:439115/1‑90 (MQ=255) gCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGGGCCAAAAGAGGGTATTAAGCGCGCACGGc > 1:116227/1‑90 (MQ=255) gCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGGGCCAAAAGAGGGTATTAAACGCGCACGGc > 2:97724/1‑90 (MQ=255) cAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAg < 1:182764/90‑1 (MQ=255) tGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCg > 2:42204/1‑90 (MQ=255) tGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGGTATTAAGCGCGCACGGCAAAAGCTGCg > 2:221237/1‑90 (MQ=255) tGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGGTATTAAACGCGCACGGCAAAAACTGCg > 2:405927/1‑90 (MQ=255) gcgcGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAgcgc > 2:262216/1‑71 (MQ=255) gcgcGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAgcgc < 1:262216/71‑1 (MQ=255) gcgcGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGGTATTAAGCGCGCACGGCAAAAGCTGCggg > 2:374371/1‑88 (MQ=255) gcgcGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGGGCCAAAAGAGGGTATTAAGCGCGCACGGCAAAAGCTGCggg > 1:17782/1‑88 (MQ=255) gcggcgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGGCGGGGGGCCAAAAAAGGGTATTAAACGCGCACGGCAAAAGCTGCggggg > 1:314282/1‑86 (MQ=255) gcggcgAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGGGCCAAAAAAGGGTATTAAGCGCGCACGGCAAAAGCTGCGGGAg > 1:444148/1‑90 (MQ=255) cGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATgg < 1:471583/90‑1 (MQ=255) agagGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAAt < 2:444148/90‑1 (MQ=255) agagGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAAt < 1:42204/90‑1 (MQ=255) | TTTATGCACAGCGTGCAGGATGGCGCGGCGAGCAAAAGCTACGGCCTGGCGGTTGCAGCTCTGGCAGGCGTGCCAAAAGAGGTTATTAAGCGCGCACGGCAAAAGCTGCGTGAGCTGGAAAGCATTTCGCCGAACGCCGCCGCTACGCAAGTGGATGGTACGCAAAT > NZ_CP009273/2852723‑2852889 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |