| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,673,561 | 0 | A | C | 55.6% | 13.9 / 17.3 | 18 | Y293S (TAC→TCC) | ydgH | DUF1471 family protein YdgH |
| Reads supporting (aligned to +/- strand): ref base A (6/2); new base C (0/10); total (6/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.51e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCC > NZ_CP009273/1673486‑1673643 | ctGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGcc > 2:459513/1‑90 (MQ=255) gCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGt > 2:626665/1‑90 (MQ=255) gCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAAc > 1:629950/1‑90 (MQ=255) aTGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGt > 1:745073/1‑90 (MQ=255) taccAGTTTGCAAAACGTGCCGCGAAAAAAGGTGCTAAGTTTTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc < 2:199197/87‑1 (MQ=255) tATCGGGTTGCGAACCGTGCCGCGAAGAAGGGTGCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc < 2:345557/90‑1 (MQ=255) tATCAGGTTGCGAAACGTGCCGAGAAAAAAGGTGCTAAGTATCCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc < 1:626665/90‑1 (MQ=255) tATCAGGTTGCAAAACGTGCCGCAAAAAAAGGTGCAAAGTATTACCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc < 2:268150/90‑1 (MQ=255) tATCAGGTGGCGAAACGGCCCGGGAAAAAAGGTGATAATTATTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc < 1:656618/90‑1 (MQ=255) gtttGCAAAACGTGCCGCGAAAAAAGGTCCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc < 2:745073/88‑1 (MQ=255) ggTTGCGAAACGTGCCGCGAAAAAAGGTGCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc < 2:629950/90‑1 (MQ=255) ggTTGCGAAACGTGCCGCGAAAAAAGGTGATAATTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc < 1:116350/90‑1 (MQ=255) ggTTGCGAAACGGGCCGGAAAAAAAGGTGCTAAGTTTTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc < 1:130341/90‑1 (MQ=255) ggTTGCGAAACGCGCCGCGAAAAAAGGTCATAAGTATTCCCCCCCCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc < 1:47879/90‑1 (MQ=255) cGAACCGTGCCGGAAAAAAAGCTGCTAAGTATTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATc < 2:316673/90‑1 (MQ=255) gCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGGACGAGGGAAAAAACTGACCGGCAGCGCAGATCTGTataaa > 1:652964/1‑90 (MQ=255) caaaGTATAACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACt < 2:652964/88‑1 (MQ=255) aaGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTcc > 2:306754/1‑90 (MQ=255) | CTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCC > NZ_CP009273/1673486‑1673643 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |