Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 760,187 | 0 | G | C | 55.6% | 8.1 / 26.6 | 18 | P184P (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.38e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGG > NZ_CP009273/760105‑760270 | 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) gATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGGGAGATcggcgg > 1:313221/1‑90 (MQ=255) | TGACCTATGAAGCGGTTAAACAGACCACGGATTACGGTTTCGGTCAGTCGACCTGTGTCGGTATCGGCGGTGACCCGATCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGG > NZ_CP009273/760105‑760270 |
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: ‑ |