| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,255,725 | 0 | A | C | 53.3% | 16.1 / 19.3 | 15 | G40G (GGT→GGG) | fruK | 1‑phosphofructokinase |
| Reads supporting (aligned to +/- strand): ref base A (6/1); new base C (0/8); total (6/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.40e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTGTCTTTACCCAGGAAGCCGCCAACGGTGACATCAATTCCCAGGTCTTTTAATACTTTGGCCACGTTGATGCCTTTACCCGCCGCATGCAGACCGGTGGTTTTCACCAGGTTCACTTCGCCGCGTTCAATTTCCGGGCAGAAACCAACAAGGTCATAAGCC > NZ_CP009273/2255648‑2255809 | ttGTCTTTCCCCAGGAAGCCCCCAACGGTGACATCAATCCCCAGGTTTTTTAATATTTTGGCCACGTTGATGCCTTTCCCCGCCGCATGc < 1:121411/90‑1 (MQ=255) aGGAAGCCGCCAACGGTGACATCAATTCCCAGGTCTTTTAATACTTTGGCCACGTTGATGCCTTTACCCGCCGCATGCAGACCGGTGGtt > 2:370048/1‑90 (MQ=255) aGGAAGCCGCCAACGGTGACATCAATTCCCAGGTCTTTTAATACTTTGGCCACGTTGATGCCTTTACCCGCCGCATGCAGACCGGTGGtt > 2:85087/1‑90 (MQ=255) cgccAACGGTGACATCAATTCCCAGGTTTTTTAATATTTTGGCCACGTTGATGCCTTTCCCCGCCGCATGCAGACCGGTGGTTTTCACCa < 2:361366/90‑1 (MQ=255) cTTTTAATCTTTTGCCCCGGTTGATGCCTTTCCCCGCCGCATGCAGACCGGTGGTTTTCACCAGGTTCACTTCGCCGCGTTCAATTTCCg < 1:50165/90‑1 (MQ=255) cTTTTAATACTTTGGCCACGTTGATGCCTTTACCCGCCGCATGCAGACCGGTGGTTTTCACCAGGTTCACTTCGCCGCGTTCAATTTCCg > 2:65591/1‑90 (MQ=255) ttttttGGCCATGTTGATCCTTTTCCCCGCCGCATGCAGACCGGTGGTTTTCACCAGGTTCACTTCGCCGCGTTCAATTTCCGGGCAGaa < 1:224248/87‑1 (MQ=255) ttGGCCACGTTGATGCCTTTACCCGCCGCATGCAGACCGGTGGTTTTCACCAGGTTCACTTCGCCGCGTTCAATTTCCGGGCAGAAACca > 1:440627/1‑90 (MQ=255) cACGTTGATGCTTTTCCCCGCCGCATGCAGACCGGTGGTTTTCACCAGGTTCACTTCGCCGCGTTCAATTTCCGGGCAGAAACCAACAAg < 2:47633/90‑1 (MQ=255) cttttCCCCGCCGCATGCAGACCGGTGGTTTTCACCAGGTTCACTTCGCCGCGTTCAATTTCCGGGCAGAAACCAACAAGGTCATAAGcc < 1:108600/88‑1 (MQ=255) cttttCCCCGCCGCATGCAGACCGGTGGTTTTCACCAGGTTCACTTCGCCGCGTTCAATTTCCGGGCAGAAACCAACAAGGTCATAAGcc < 1:204041/88‑1 (MQ=255) ccTTTCCCCGCCGCATGCAGACCGGTGGTTTTCACCAGGTTCACTTCGCCGCGTTCAATTTCCGGGCAGAAACCAACAAGGTCATAAGcc < 2:520582/90‑1 (MQ=255) tttACCCGCCGCATGCAGACCGGTGGTTTTCACCAGGTTCACTTCGCCGCGTTCAATTTCCGGGCAGAAACCAACAAGGTCATAAGcc > 1:19217/1‑88 (MQ=255) tttACCCGCCGCATGCAGACCGGTGGTTTTCACCAGGTTCACTTCGCCGCGTTCAATTTCCGGGCAGAAACCAACAAGGTCATAAGcc < 2:19217/88‑1 (MQ=255) tACCCGCCGCATGCAGACCGGTGGTTTTCACCAGGTTCACTTCGCCGCGTTCAATTTCCGGGCAGAAACCAACAAGGTCATAAGcc > 2:85438/1‑86 (MQ=255) | TTGTCTTTACCCAGGAAGCCGCCAACGGTGACATCAATTCCCAGGTCTTTTAATACTTTGGCCACGTTGATGCCTTTACCCGCCGCATGCAGACCGGTGGTTTTCACCAGGTTCACTTCGCCGCGTTCAATTTCCGGGCAGAAACCAACAAGGTCATAAGCC > NZ_CP009273/2255648‑2255809 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |