Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 760,187 | 0 | G | C | 59.3% | 18.4 / 38.8 | 27 | P184P (CCG→CCC) | sucD | succinate‑‑CoA ligase subunit alpha |
Reads supporting (aligned to +/- strand): ref base G (2/9); new base C (16/0); total (18/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.17e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGTACACTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGC > NZ_CP009273/760098‑760276 | 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:1379787/90‑1 (MQ=255) cTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGAt < 2:192930/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) cGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAgcg > 2:119342/1‑90 (MQ=255) cGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAgcg > 1:505319/1‑90 (MQ=255) gggCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTgaga > 2:286249/1‑77 (MQ=255) gggCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTgaga < 1:286249/77‑1 (MQ=255) gggCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAgcgc < 2:100214/90‑1 (MQ=255) | GGTACACTGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGC > NZ_CP009273/760098‑760276 |
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: ‑ |