Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,667,520 | 0 | A | C | 57.9% | 6.3 / 19.8 | 19 | G176G (GGT→GGG) | yphB | aldose 1‑epimerase |
Reads supporting (aligned to +/- strand): ref base A (6/2); new base C (0/11); total (6/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.03e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTGTTCACCCACTGGCGCGGCAACGGCGCGGGCTGGTTAAAATCCAGTTCCTGCGGTAGCTGCTCGCAAAACTCACCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAA > NZ_CP009273/2667445‑2667598 | aTTGTTCACCCACTGGCGCGGCAACGGCGCGGGCTGGTTAAAATCCAGTTCCTGCGGTAGCTGCTCGCAAAACTCACCCGCCAGCCACTg > 1:282447/1‑90 (MQ=255) tCACCCACTGGCGCGGCAACGGCGCGGGCTGGTTAAAATCCAGTTCCTGCGGTAGCTGCTCGCAAAACTCACCCGCCAGCCACTGCTccc > 1:470824/1‑90 (MQ=255) cattGGCGCGGCAGCGGCGCGGGAGGTTTAAAACCCCGTCCCGGCGGTATCTGCCGGCAAAACTCACCCGCCAGCCACTGCTCCCGCTcc < 1:134316/87‑1 (MQ=255) cACTGGCGCGGCAACGGCGCGGGCTGGTTAAAATCCAGTTCCTGCGGTAGCTGCTCGCAAAACTCCCCCGCCAGCCACTGCTCCCGCTcc < 2:1129218/90‑1 (MQ=255) cACGGGCGCGGCAACGGCGCGGGCTGGTTAAAACCCAGTTCCTGCGGTAGCTGCTCGCAAAACTCCCCCGCCAGCCACTGCTCCCGCTcc < 1:63143/90‑1 (MQ=255) gggCAGCGGCGGGGGCGGGTAAAAACCCAGTCCCTGCGGTCGCGGCCCGCAAACCCCCCCCGCCAGCCACTGCTCCCGCTCCAGCCAGTa < 1:928396/89‑1 (MQ=255) cgGCAACGGCGCGGGCTGGTTAAAATCCAGTTCCTGCGGTAGCTGCTCGCAAAACTCACCCGCCAGCCACTGCTCCCGCTCCAGCCAGTa > 2:943193/1‑90 (MQ=255) cGGCGCGGGCTGGTTAAAATCCAGTTCCTGCGGTAGCTGCTCGCAAAACTCCCCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCt < 2:454511/90‑1 (MQ=255) gggCTGGTAAAAACCCAGTCCCTGCGGTAGCTGCTCGAAAAACCCCCCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCtg < 1:436738/90‑1 (MQ=255) ggTTAAAATCCAGTTCCTGCGGTAGCTGCTCGCAAAACTCACCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCt > 2:600365/1‑90 (MQ=255) cAGTTCCTGCGGTAGCTGCTCGCAAAACTCACCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGt > 2:742752/1‑90 (MQ=255) cTGCGGTAGCTGCTCGCAAAACTCCCCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTgcgg < 2:218722/90‑1 (MQ=255) ggTAGCTGCTCGCAAAACTCCCCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGAc < 2:864775/90‑1 (MQ=255) ctgctCGCAAAACTCCCCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACgg < 1:1003995/90‑1 (MQ=255) ctgctCGCAAAACTCACCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACgg < 2:282447/90‑1 (MQ=255) ctgccCGAAAACCCCCCCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACgg < 1:228778/90‑1 (MQ=255) ctgccCGAAAACCCCCCCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACgg < 2:631607/90‑1 (MQ=255) gctCGCAAACCTCCCCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGaa < 2:636770/90‑1 (MQ=255) tCGCAAAACTCACCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGaaaa > 1:1056344/1‑90 (MQ=255) | ATTGTTCACCCACTGGCGCGGCAACGGCGCGGGCTGGTTAAAATCCAGTTCCTGCGGTAGCTGCTCGCAAAACTCACCCGCCAGCCACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAA > NZ_CP009273/2667445‑2667598 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |