| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 760,187 | 0 | G | C | 55.2% | 20.9 / 43.3 | 29 | P184P (CCG→CCC) | sucD | succinate‑‑CoA ligase subunit alpha |
| Reads supporting (aligned to +/- strand): ref base G (4/9); new base C (16/0); total (20/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.14e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-07 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGTACACTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCG‑ATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTA > NZ_CP009273/760098‑760272 | ggTACACTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCG‑ATCCCg < 1:1308555/90‑1 (MQ=255) gACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGAGCCG‑ATCCCCGGCTCTaa > 1:289328/1‑90 (MQ=255) gACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCC‑ATCCCCGGccctta > 1:108339/1‑88 (MQ=255) gACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCC‑ATCCCCGGccctta > 1:254983/1‑88 (MQ=255) gACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGGCCC‑GTCCCCGGCGCTaa > 2:991796/1‑90 (MQ=255) gACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGCCCC‑ATCCCCGGCGCTaa > 2:503510/1‑90 (MQ=255) aTGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCC‑GTCCCCGGCTCTAACcttt > 1:756076/1‑89 (MQ=255) gCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCC‑ATCCCCGGGTCTTACTTTTTTGAc > 1:340314/1‑90 (MQ=255) gCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGAGCCC‑GTCCCCGGCGCTAACTTTTTTGAc > 2:675600/1‑90 (MQ=255) ggTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCG‑ATCCCGGGCTCTAACTTTATCGACAt < 2:820631/90‑1 (MQ=255) ggTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCC‑CTCCCCGGCTCTTAATTTTTCGACAt > 1:1032359/1‑90 (MQ=255) ggTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGAACCC‑CCCCCCGGGGCTAACCTTTTCGAcca > 1:264172/1‑88 (MQ=255) aGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCG‑ATCCCGGGCTCTAACTTTATCGACATTCTCGaaa < 1:976630/90‑1 (MQ=255) gACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCC‑CTCCCCGGCGCTAACTTTTTTGACATTTTCCaaaa > 1:1229987/1‑89 (MQ=255) ccACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCC‑ATCCCCGGCTCTAACTTTATCGACATTTTCGAAAAGt > 1:14553/1‑90 (MQ=255) ccACGGATTACAGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCC‑ATCCCCGGCTCTAACTTTTTTCACATTTTCGAAAAGt > 1:473333/1‑90 (MQ=255) tACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCC‑ATCCCCGGCTCTAACTTTTTTGACATTTTCCAAAAGTTTGaaaaa > 1:1277573/1‑90 (MQ=255) cGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCG‑ATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGa < 2:1212171/90‑1 (MQ=255) tCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCGCACCCCCGGC‑CTGACTTTTTCGACATTATCGCAAAGTTTCAAAAAAATcccc > 1:1212171/1‑88 (MQ=255) gACCTGTGTCGGTATCGGCGGTGACCCG‑ATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGc < 2:1228185/90‑1 (MQ=255) gACCTGTGTCGGTATCGGCGGTGACCCG‑ATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGc < 2:476207/90‑1 (MQ=255) cTGTGTCGGTATCGGCGGTGACCCG‑ATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGAt < 1:554632/90‑1 (MQ=255) gtgtCGGTATCGGCGGTGACCCG‑ATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCg < 2:1202377/90‑1 (MQ=255) gtCGGTATCGGCGGTGACCCG‑ATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGtg < 1:181648/90‑1 (MQ=255) ggTATCGGCGGTGACCCC‑ATCCCCGGCTCTTACTTTTTTGACATTTTCTAAAAGTTTGAAAAAAAACCCCCGACCGAAACGGACgtggtg > 2:159189/1‑90 (MQ=255) g‑ATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATcggcgg > 1:872449/1‑90 (MQ=255) g‑ATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATcggcgg > 2:138225/1‑90 (MQ=255) g‑ATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATcggcgg > 2:253856/1‑90 (MQ=255) g‑ATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATcggcgg > 2:467935/1‑90 (MQ=255) tCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTa < 2:795927/90‑1 (MQ=255) | GGTACACTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCG‑ATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTA > NZ_CP009273/760098‑760272 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |