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% | 8.6 / 21.7 | 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.11e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTG > NZ_CP009273/760106‑760260 | gACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTaa < 1:142959/90‑1 (MQ=255) gACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGCCCCATCCCCGGcgctta > 1:341704/1‑88 (MQ=255) gACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCGTCCCCGGcgctta > 2:475865/1‑88 (MQ=255) gCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGCCCCATCCCCGGCTCTAAATTTTTTGAc > 1:311450/1‑90 (MQ=255) ggTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACAt < 2:351334/90‑1 (MQ=255) ggTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGCCCCCCCCCCGGCCCTAACTTTTTCTACAt > 2:277331/1‑90 (MQ=255) ggTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGAACCCACCCCCGGCTCTAACTTTATTGACAt > 2:354429/1‑90 (MQ=255) gACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGGCCCGTCCCCGGGGCTCACTTTTTTTACATTCTTGAAAt > 2:440624/1‑90 (MQ=255) gACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCAACCCCGGCTCTAACTTTTTTGACCATTTCGaaaa > 1:481179/1‑89 (MQ=255) tACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGGGCCGCTCCCCGGCGCTTACTTTTTCAACATTCTCTTAAAGTTTGaaaaa > 2:205708/1‑90 (MQ=255) tACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCGACCCCGGCTCTTACTTTTTTGACATTTTCGAAATGTTTGaaaaa > 1:263682/1‑90 (MQ=255) aCGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAg < 1:266966/90‑1 (MQ=255) cGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGGGACCCCATCCCCGGCGCTTACTTTTTTTACATTTTCGAAAAGGTTGaaaaaaa > 1:285031/1‑88 (MQ=255) tCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGAc < 1:137430/90‑1 (MQ=255) aCCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCg < 1:277331/90‑1 (MQ=255) cTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGAt < 1:261514/90‑1 (MQ=255) gcggTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTg < 2:467449/90‑1 (MQ=255) gcggTGACCCCATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAAATCCGCAGACCGAAGCGATCGTGATGATCgggg > 1:88697/1‑88 (MQ=255) | GACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTG > NZ_CP009273/760106‑760260 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |