| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,673,561 | 0 | A | C | 57.1% | 11.6 / 14.7 | 14 | Y293S (TAC→TCC) | ydgH | DUF1471 family protein YdgH |
| Reads supporting (aligned to +/- strand): ref base A (5/1); new base C (0/8); total (5/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.00e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAA > NZ_CP009273/1673472‑1673648 | cAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTa < 2:211238/90‑1 (MQ=255) ctctGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCg > 1:447005/1‑90 (MQ=255) ctctGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCg > 2:159823/1‑90 (MQ=255) tctGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGc < 2:447005/90‑1 (MQ=255) cGGCAACATGACCGAAGTCTCCTATCAGGTTGCAAAACGTGCCGCGAAAAAAGGTGCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGa < 1:136182/90‑1 (MQ=255) ggCAACATGACCGAAGCCCCCTATCAGGTTGCCAAACGGCCCGCGAAGAAGGGTGCTAATTATTCCCACACCCCCCGCCAGTGGCAGGaa < 1:217007/90‑1 (MQ=255) aaaaTGCCCAAAGTCTCCTATAAGGTTGCGAAACGTGCCGCGAAAAAGGGTGTAAAGTTTTCCCACACCCCCCGCCAGTGGCAGGAACGt < 2:579343/87‑1 (MQ=255) aaCATGACCGAAGTCCCCTTTCGGGTTGCGAAAGGTGCCGAGAAGAAAGGTGTTAAGTATTCCCCCACCCCCCGCCAGTGGCAGGAACGt < 2:354427/90‑1 (MQ=255) ggTTGCGAACCGTGCCGCGAAAAAGGGTGCTAAGTTTTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc < 1:159823/90‑1 (MQ=255) tGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGGGGTAATAACCTGACCGTCAGCGCAGa > 1:129562/1‑90 (MQ=255) tGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGGGGGAATAACCTGACCGTCAGCGCAGa > 1:422428/1‑90 (MQ=255) cGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATc < 2:239879/90‑1 (MQ=255) ccGCAAAAAAAGGTGCTAAGTTTTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataaat < 2:184277/90‑1 (MQ=255) atAAGTTTTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACt < 1:356714/89‑1 (MQ=255) ttACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATaa > 1:259524/1‑90 (MQ=255) | CAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAA > NZ_CP009273/1673472‑1673648 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |