| 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% | 21.2 / 24.4 | 18 | G180G (GGT→GGG) | sucD | succinate‑‑CoA ligase subunit alpha |
| Reads supporting (aligned to +/- strand): ref base T (0/8); new base G (10/0); total (10/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.29e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCGTTCCGGTACACTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGA > NZ_CP009273/760091‑760263 | ccGTTCCGGTACACTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGAccc < 1:498165/90‑1 (MQ=255) tGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTa < 1:502988/90‑1 (MQ=255) gACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGCCCGATCCCCGGGTCTaa > 1:68845/1‑90 (MQ=255) gACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGCCCCATCCCCGGCTCTaa > 1:47083/1‑90 (MQ=255) aTGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGCCCCCTCCCCGGGCCTTACtttt > 1:478190/1‑89 (MQ=255) aaGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCCCCCCCGGCGCTAACTTTTTCg > 2:248192/1‑90 (MQ=255) ggTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCGTCCCCGGCGCTAACTTTTTCTACAt > 2:582655/1‑90 (MQ=255) aGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGaaa < 2:516714/90‑1 (MQ=255) gACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCGACCCCCGGCCCTAACTTTATCGACATTCTCGaaaa > 2:470596/1‑89 (MQ=255) gACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCATCCCCGGCTCTAACTTTTTCTACATTCTCCAAAt > 1:525852/1‑90 (MQ=255) tACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGGCCCCTCCCCGGGCCTAAATTTTTTGTCATTTTTGTAATGTTCGaaaaa > 2:564718/1‑90 (MQ=255) tACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGCCCCATCCCCGGCCCTAACTTTTTTGACATTCTCGAAAAGTTTGaaaaa > 2:200408/1‑90 (MQ=255) tCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGc < 1:458602/90‑1 (MQ=255) tCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGAc < 1:250940/90‑1 (MQ=255) cAGTCGACCTGTGTCGGTATCGGCGGGGGCCCCATCCCCGGCTCTAACTTTTTTGACATTTTCGAAAAGTTTGAAAAAAAACCCCCGAcc > 2:112686/1‑90 (MQ=255) gACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGc < 2:47083/90‑1 (MQ=255) ccTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGa < 2:525852/90‑1 (MQ=255) gTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTgaga < 1:582655/90‑1 (MQ=255) | CCGTTCCGGTACACTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGA > NZ_CP009273/760091‑760263 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |