| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 760,181 | 0 | G | C | 55.6% | 13.0 / 24.5 | 18 | P182P (CCG→CCC) | sucD | succinate‑‑CoA ligase subunit alpha |
| Reads supporting (aligned to +/- strand): ref base G (1/7); new base C (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. | |||||||||||
GGTACACTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGG > NZ_CP009273/760098‑760270 | ggTACACTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCg < 2:1019309/90‑1 (MQ=255) cTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCGTCCCCGGcccc > 2:989966/1‑87 (MQ=255) cTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCATCCCCGGcgct > 2:813238/1‑90 (MQ=255) gACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCCATCCCCGGCTCTaa > 2:171385/1‑90 (MQ=255) gACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGCCCGGTCCCGGGGGCTaa > 1:1143609/1‑90 (MQ=255) ggTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGGCCCCATCCCGGGGTCTTAATTTTTTGACAt > 2:237138/1‑90 (MQ=255) ggTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTACCGACAt < 2:940183/90‑1 (MQ=255) tAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCCATCCCCGGGCCCAACTTTTTTTACAtttt > 2:1228385/1‑88 (MQ=255) tAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCAACCCCGGCTCTAACTTTATCGACATtct > 2:22369/1‑90 (MQ=255) ccACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGt < 1:813238/90‑1 (MQ=255) ccACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCCCTCCCCGGCTCTAACTTTTTTGACATTTTCGAAATGt > 2:362658/1‑90 (MQ=255) gATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGAGCCCCGCCCCGGCTCCGATTTTTTTGACATTCTCGGAATGttttaa > 1:997402/1‑86 (MQ=255) aCGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAg < 1:1129251/90‑1 (MQ=255) ttCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGGCCCCTCCCCGGGCCTTACTTTTTTTACATTTTCGAAAAGTTTGAAAAAAAtacc > 1:852263/1‑87 (MQ=255) tCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGc < 1:1077022/90‑1 (MQ=255) tCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGCCCCCTCCCCGGCCCTAACTTTTTTTACATTCTCGAAAAGTTCGAAAAAAAACCGc > 1:196020/1‑90 (MQ=255) tATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTgatgat < 2:686751/90‑1 (MQ=255) gTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTgaga < 2:885195/90‑1 (MQ=255) gATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATcggcgg > 1:675148/1‑90 (MQ=255) | GGTACACTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGG > NZ_CP009273/760098‑760270 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |