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% | 22.0 / 39.2 | 27 | P182P (CCG→CCC) | sucD | succinate‑‑CoA ligase subunit alpha |
Reads supporting (aligned to +/- strand): ref base G (1/11); new base C (15/0); total (16/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.20e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.36e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGGTACACTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAG > NZ_CP009273/760097‑760262 | cGGTACACTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATccc < 1:1440684/90‑1 (MQ=255) ggTACACTGACCTAAGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCg < 2:276529/90‑1 (MQ=255) cTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGGCCCCCTCCCCGGccct > 2:711969/1‑90 (MQ=255) cTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCCCTCCCCGGctct > 1:1379787/1‑90 (MQ=255) cTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGCCCGCTCCCCGGgtct > 1:192930/1‑90 (MQ=255) gACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCATCCCCGGctctta > 1:248123/1‑88 (MQ=255) gACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCAGCCCCGGCTCTaa > 2:87042/1‑90 (MQ=255) cTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGCCCCAGCCCCGGCCCTTACtt > 2:12116/1‑90 (MQ=255) aTGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCCCTCCCCGGCTCTTACtttt > 1:1383387/1‑89 (MQ=255) gCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGAc < 1:850421/90‑1 (MQ=255) gCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCCACCCCCGGGCCTTACTTTTTTGAc > 1:440706/1‑90 (MQ=255) gCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGCCCCCTCCCCGGGCCTCACTTTATctcc > 2:725591/1‑87 (MQ=255) ggTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCCATCCCCGGGTCTAACTTTTTTGACAt > 1:520151/1‑90 (MQ=255) ggACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGaaa < 1:562940/89‑1 (MQ=255) tACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCCATCCCCGGCTCTCACTTTTTTTTCATTTTCGGAATGTTTGgaaaa > 2:562940/1‑90 (MQ=255) tACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGCCCCCTCCCCGGGTCTTACTTTTTTTACATTTTCGAAAAGTTGGgaaaa > 2:1017908/1‑90 (MQ=255) aCGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAg < 1:591407/90‑1 (MQ=255) cGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGa < 1:87042/90‑1 (MQ=255) gACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGc < 2:889662/90‑1 (MQ=255) cTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGAt < 2:192930/90‑1 (MQ=255) cTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGAt < 2:1379787/90‑1 (MQ=255) gtgtCGGTATCGGCGGGGGACCCCTCCCCCGCCCTAACTTTTTTGACATTCTCCCAATGTTTGGAAAAAAAACGCAGACCGAAGGGATCg > 2:860351/1‑90 (MQ=255) gtCGGTATCGGCGGGGACCCCCTCCCCGGGTCTTACTTTTTTGACATTTTCGAAAAGTTTTGAAAAAATCCCCGCACCCAAGCGACCggg > 1:813964/1‑88 (MQ=255) gtCGGTATCGGCGGGGACCCCATCCCCGGCGCTAAATTTTTCTACATTTTCGAAATGTTTGAAAAAAAACCGCAGAACGAaacgaacgcg > 1:927389/1‑88 (MQ=255) ggTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGtgatg < 1:989825/90‑1 (MQ=255) gcggTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTg < 1:515803/90‑1 (MQ=255) ggTGGCCCCATCCCCGGCCCTAAATTTTTTCACCTTTTCCAAATGTTTGAAAAAAAACCCCCGACCGAAGCGATCGTGGTGGTCGGGgag > 2:1280420/1‑90 (MQ=255) | CGGTACACTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAG > NZ_CP009273/760097‑760262 |
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: ‑ |