| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,673,569 | 0 | A | C | 57.9% | 16.2 / 23.5 | 19 | T296P (ACC→CCC) | ydgH | DUF1471 family protein YdgH |
| Reads supporting (aligned to +/- strand): ref base A (7/1); new base C (0/11); total (7/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.59e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-05 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCT > NZ_CP009273/1673486‑1673655 | ctGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGcc > 2:571862/1‑90 (MQ=255) tGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCa < 1:571862/90‑1 (MQ=255) aaCTACGGCAACATGACCGAAGTCTCCTGTCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGTTAAGTATTACCACATCACCCGCCAGTgg > 2:508256/1‑90 (MQ=255) cTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCCCCCGCCAGTGGCa < 1:114604/90‑1 (MQ=255) acaaTGACCGAATCCCCCTATCAGGTGCCGAAACGGCCGGAGAAGAAGGGTGGTAATTTTTACCACATCCCCCGCCAGTGGCAGGAACGt < 2:687810/87‑1 (MQ=255) aaCATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCCCCCGCCAGTGGCAGGAACGt < 1:234416/90‑1 (MQ=255) tGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCCCCCGCCAGTGGCAGGAACGTGGta < 1:271016/90‑1 (MQ=255) tATCAGGTTGCAAACCGGGCCGCAAAAAAAGGTGCAAAGTATTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc < 2:180909/90‑1 (MQ=255) tCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAg > 1:167419/1‑90 (MQ=255) ggTTGCAAACGGCGCCGCGAAGAAGGGTGCTAATTTTTCCCCCATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc < 2:726531/90‑1 (MQ=255) ggTTGCAAACCGTGCCGAAAAAAAGGGTGCTAAGTATTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc < 2:167419/90‑1 (MQ=255) cgcaAAAAAGGGTGTAAAGTTTTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataaata < 2:211928/90‑1 (MQ=255) aaaaGGTGCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTAtt < 1:385953/89‑1 (MQ=255) gCTAAGTATTCCCCCACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTAc < 2:191482/90‑1 (MQ=255) ttACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATaa > 1:860100/1‑90 (MQ=255) ttACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATaa > 1:353630/1‑90 (MQ=255) ccACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAggc < 1:651949/90‑1 (MQ=255) acaTCACCCGCCAGTGGCAGGAACGTGGGAATAACCTGACCGTCAGCGCAGATCTGTATAAAAAACTATTCAGTTACTCCTATAAggcgg > 2:584268/1‑90 (MQ=255) aTCACCCGCCAGTGGCAGGAACGGGGTAATAACCTGACCGTCAGCCCAGAACTGTATAAATAACTATTCAGTTTCTCCTAAAAggcgggt > 1:480256/1‑88 (MQ=255) | CTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCT > NZ_CP009273/1673486‑1673655 |
| 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: ‑ |