| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 760,175 | 0 | T | G | 55.6% | 16.3 / 21.5 | 18 | G180G (GGT→GGG) | sucD | succinate‑‑CoA ligase subunit alpha |
| Reads supporting (aligned to +/- strand): ref base T (1/7); new base G (10/0); total (11/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.51e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTG > NZ_CP009273/760106‑760260 | gACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTaa < 1:142959/90‑1 (MQ=255)gACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGCCCCATCCCCGGcgctta > 1:341704/1‑88 (MQ=255)gACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCGTCCCCGGcgctta > 2:475865/1‑88 (MQ=255) gCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGCCCCATCCCCGGCTCTAAATTTTTTGAc > 1:311450/1‑90 (MQ=255) ggTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACAt < 2:351334/90‑1 (MQ=255) ggTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGCCCCCCCCCCGGCCCTAACTTTTTCTACAt > 2:277331/1‑90 (MQ=255) ggTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGAACCCACCCCCGGCTCTAACTTTATTGACAt > 2:354429/1‑90 (MQ=255) gACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGGCCCGTCCCCGGGGCTCACTTTTTTTACATTCTTGAAAt > 2:440624/1‑90 (MQ=255) gACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCAACCCCGGCTCTAACTTTTTTGACCATTTCGaaaa > 1:481179/1‑89 (MQ=255) tACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGGCCGCTCCCCGGCGCTTACTTTTTCAACATTCTCTTAAAGTTTGaaaaa > 2:205708/1‑90 (MQ=255) tACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCGACCCCGGCTCTTACTTTTTTGACATTTTCGAAATGTTTGaaaaa > 1:263682/1‑90 (MQ=255) aCGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAg < 1:266966/90‑1 (MQ=255) cGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCATCCCCGGCGCTTACTTTTTTTACATTTTCGAAAAGGTTGaaaaaaa > 1:285031/1‑88 (MQ=255) tCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGAc < 1:137430/90‑1 (MQ=255) aCCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCg < 1:277331/90‑1 (MQ=255) cTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGAt < 1:261514/90‑1 (MQ=255) gcggTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTg < 2:467449/90‑1 (MQ=255) gcggTGACCCCATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAAATCCGCAGACCGAAGCGATCGTGATGATCgggg > 1:88697/1‑88 (MQ=255) | GACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTG > NZ_CP009273/760106‑760260 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |