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.0% | 12.9 / 18.6 | 20 | Y293S (TAC→TCC) | ydgH | DUF1471 family protein YdgH |
Reads supporting (aligned to +/- strand): ref base A (6/3); new base C (0/11); total (6/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.17e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTAC > NZ_CP009273/1673479‑1673640 | aaaTCCTCGGGCACCGACGGCAACATGACCGAAGTCCCCTCCCAGGTTCCAAATCGTGCCGCGAAGGAGGGTGCTAAGTATTACcccacc < 1:607209/90‑3 (MQ=255) ctctGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCg > 2:1251300/1‑90 (MQ=255) ggCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAg > 2:1132613/1‑90 (MQ=255) gCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAAc > 2:332643/1‑90 (MQ=255) aaCATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAAAAAGGTGCTAAGTATTCCCACATCACCCGCCAGTGGCAGGAACGt < 2:1078138/90‑1 (MQ=255) tGACCGAAGTCCCCTATCAGGTTGCGAAACGTGCCGGGAAAAAAGGTGCAAAGTTTTACCCCAACCCCCGCCAGTGGCAGGAACGTGGta < 1:44600/90‑1 (MQ=255) tGACCGAAGCCCCCCATCCGTTTGCAAAAGGCGCCGAAAAGAAAGGTGCAAACTTTTCCAACATCCCCCGCCAGTGGCAGGAACGTGGta < 1:1251300/90‑1 (MQ=255) tctcCTATGAGGTTGAGAAGCGGGCAGAAAAAAAGGGTGCTAATTATTCCCCCACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGa < 1:1365695/90‑1 (MQ=255) tATCAGGTTGCGAAACGTGCCGAGAAAAAGGGTGATAAGTATTCCCCCACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc < 1:604361/90‑1 (MQ=255) aGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGGGGTAATAACCTGACCGTCAgcg > 1:861089/1‑90 (MQ=255) ggTTGCGAACCGTGCCGCGAAAAAAGGTGCTAAGTTTTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc < 1:396016/90‑1 (MQ=255) ggTTGCGAAACGTGCCGCGAAGAAAGGTGTTAAGTATTACCACACCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc < 1:114364/90‑1 (MQ=255) ggTTGAAAAAGGTCCCGCGAAAAAAGGTGATAATTTTTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc < 2:1033691/90‑1 (MQ=255) ggTGGCGAACCGTGCCGGGAAAAAAGGTGCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc < 2:258135/90‑1 (MQ=255) gTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGGACGGGGGAATAACCTGACCGTCAGCGCAGATCTGtata > 1:639819/1‑90 (MQ=255) gTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGtata > 2:1829/1‑90 (MQ=255) cgcgAAAAAAGGTGCTAATTATTCCCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataaata < 2:639819/90‑1 (MQ=255) agaaGGGTGCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTAt < 2:861089/90‑1 (MQ=255) tGCAAATTATTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTa < 2:43964/90‑1 (MQ=255) gCTAAGTTTTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTAc < 1:424705/90‑1 (MQ=255) | AAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTAC > NZ_CP009273/1673479‑1673640 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |