| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 4,297,122 | 0 | A | C | 52.2% | 26.3 / 31.3 | 23 | G165G (GGT→GGG) | alsK | allose kinase |
| Reads supporting (aligned to +/- strand): ref base A (9/2); new base C (0/12); total (9/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.73e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.04e-05 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGA > NZ_CP009273/4297042‑4297201 | gcCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCg < 1:45070/90‑1 (MQ=255) agCAATTGGTTTCCAGGCACCCAGGATTCCCACACGCGCAGTGTTGGGTCATATCCCCCGGGGGGATATGACCCAGTTCGCCTGCCacac < 1:246915/90‑1 (MQ=39) ggTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGCCCCAGTTCGCCTGCCACACCGTGTGc < 1:73499/90‑1 (MQ=255) aGGCCCCCAGGTTCCCCCCCCGCGCGGTGTGGGGTCATACCCCCGGGGGGGATATGCCCCAGTTCCCCTGCCACACCGTGTGCACCCGTc < 2:138771/90‑1 (MQ=255) aTTGCCACACGCGCAGTGTTGGGTCATATCCCCCAGGGGGATATGCCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCAcc < 1:124273/90‑1 (MQ=255) tcccACCCGCGCGGTGTTGGTTCATACCCCCCGGGGGGATATGCCCCGGTCCCCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGt < 2:29742/88‑1 (MQ=255) tcccACACGCGCAGTGTGGGGTCATACCCCCCGGGGGGATATGCCCCAGTCCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGt < 1:63378/88‑1 (MQ=255) ccgcacGCGGGGTTTGGGGTCATCCCCCCCAGGGGCATAGGCCCAAGTTCCCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTc < 1:241937/87‑1 (MQ=255) gcgcAGTGTTGGGTAATATCTCCCGGGGGGATATGCCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCcac < 1:146982/90‑1 (MQ=255) cgcgGTGTGGGGTCATATCCCCCGGGGGGATATGCCCCAGTCCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCcaca < 1:166926/90‑1 (MQ=255) acgggtgtTGGTTCATACCCCCCGGGGGAATATCCCCCATTCCCCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCcaca < 2:118205/86‑1 (MQ=255) atatCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAcccctgt > 1:118205/1‑87 (MQ=255) atatCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATc > 2:85684/1‑90 (MQ=255) atatCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATc > 2:184467/1‑90 (MQ=255) atatCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATc > 1:17742/1‑90 (MQ=255) atatCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATc > 1:135025/1‑90 (MQ=255) tcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGt > 1:191835/1‑90 (MQ=255) cccAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTa > 1:310568/1‑90 (MQ=255) cgggggggATATGCCCCTGTCCCCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTAc < 2:135025/87‑1 (MQ=255) ccggggggATATACCCCTGTTCCCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTAc < 1:224787/87‑1 (MQ=255) gggggggATATGCCCCAGTTCCCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTAcc < 1:280025/88‑1 (MQ=255) gggggATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCga > 2:131349/1‑90 (MQ=255) gggggATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCga > 1:220392/1‑90 (MQ=255) | GCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGA > NZ_CP009273/4297042‑4297201 |
| 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: ‑ |