Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,673,561 | 0 | A | C | 56.2% | 14.0 / 17.5 | 16 | Y293S (TAC→TCC) | ydgH | DUF1471 family protein YdgH |
Reads supporting (aligned to +/- strand): ref base A (4/3); new base C (0/9); total (4/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.92e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.98e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACT > NZ_CP009273/1673476‑1673641 | aTCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACcac < 2:666200/90‑1 (MQ=255) aaCATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGt < 2:272714/90‑1 (MQ=255) aCCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGGACGTGGtaat > 2:667503/1‑90 (MQ=255) ggTTGCGAAACGTGCCGGGGAGAAAGGTGCTAAGTTTTCCCCCATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc < 2:187749/90‑1 (MQ=255) ggTTGCGAAACGTGCCGCAAAGAAAGGTGAAAATTTTTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc < 2:184456/90‑1 (MQ=255) gCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGGAATAACCTGACCGTCAGCGCAGAt > 2:68058/1‑90 (MQ=255) cGAAACGTGCCGCAAAGAAAGGTGCTAATTATTCCCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATc < 1:295948/90‑1 (MQ=255) gAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCt > 1:791677/1‑90 (MQ=255) ggCCGAGAAGAAAGGTGAAAAGTTTTCCCCCACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataa < 2:844886/89‑1 (MQ=255) gcgAAGAAAGGTGCTAAGTATTCCCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataaataa < 2:333806/90‑1 (MQ=255) gTGCTAAGTATTCCCCCACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGtt < 2:366456/90‑1 (MQ=255) gTGCTAAGTATTCCCACCTCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGtt < 1:667503/90‑1 (MQ=255) gTGCTAAGTATTCCCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGtt < 1:304294/90‑1 (MQ=255) cTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACt < 1:32157/90‑1 (MQ=255) ttACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc < 1:187136/49‑1 (MQ=255) ttACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc > 2:187136/1‑49 (MQ=255) | ATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACT > NZ_CP009273/1673476‑1673641 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |