| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 4,607,212 | 0 | A | C | 53.8% | 13.0 / 16.7 | 13 | T25P (ACC→CCC) | deoC | deoxyribose‑phosphate aldolase |
| Reads supporting (aligned to +/- strand): ref base A (4/2); new base C (0/7); total (4/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.10e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATC > NZ_CP009273/4607135‑4607298 | atgaaatgaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAg > 1:49999/1‑90 (MQ=255)atgaaatgaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAg > 2:105965/1‑90 (MQ=255)atgaaatgaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAg > 2:229967/1‑90 (MQ=255) aaatgaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGa < 2:305516/90‑1 (MQ=255) cTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGccc < 2:100900/90‑1 (MQ=255) aaaGTGATGGACCTGACCACCCTGAATGCCGACGACCCCGCCGAGAAAGTGATCGCCCTGTGTCACCAGGCCAAACCCCCGGTCGGCAAt < 1:11236/90‑1 (MQ=255) ggcccTGCCCCCCCTGAATGACGACGACCCCGACGAAAAAGTGATCGCCCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACcgccgc < 1:29093/87‑1 (MQ=255) ggcccTGACCCCCCTGAATGACGCCGACCCCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACcgccgc < 2:245061/87‑1 (MQ=255) cTGCCCCCCCTGAATGACGCCGACCCCGACGAGAAAGTGATCCCCTTGTGTCATCAGGCAAAACCCCCGGTCGGCAATACCGCCGCTATc < 1:353517/90‑1 (MQ=255) cTGACCCCCCTGAATGACGACGACCCCGACGAAAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATc < 2:171557/90‑1 (MQ=255) cccTAAATGACGACGCCCCCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAACCCCCGGCCGGCAATCCCGCCGCTATCTGtatct < 1:236434/90‑1 (MQ=255) aaTGCCGACGCCCCCGACGAGAAAGTGATCCCCCTGTGTCATCAGGCCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCt < 1:29313/90‑1 (MQ=255) gacACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATc > 1:191399/1‑90 (MQ=255) | ATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATC > NZ_CP009273/4607135‑4607298 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |