| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,738,361 | 0 | T | G | 52.5% | 26.6 / 40.5 | 40 | G66G (GGT→GGG) | yiaN | 2,3‑diketo‑L‑gulonate transporter large permease YiaN |
| Reads supporting (aligned to +/- strand): ref base T (4/15); new base G (21/0); total (25/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.64e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.22e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGG > NZ_CP009273/3738274‑3738439 | gTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGTg > 2:615772/1‑89 (MQ=255)gTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGTg < 1:615772/89‑1 (MQ=255)gTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGGGa > 2:191025/1‑90 (MQ=255) gATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGGGAAATc > 2:684407/1‑90 (MQ=255) aTGTCCAGATCATGGCACAAACGCTGGTGAACGGGGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGGGGAATCa > 1:385827/1‑90 (MQ=255) ccAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGTGaaa < 1:1438453/83‑1 (MQ=255) ccAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGTGaaa > 2:1438453/1‑83 (MQ=255) ccAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGGGGAATCAAGaa > 1:225901/1‑90 (MQ=255) tGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGGGGAATCAAAAATgcgggcg > 1:1003657/1‑90 (MQ=255) tGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGGGAAAACATGAATgcggggg > 1:1435681/1‑88 (MQ=255) acaAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGTGAAATCATGAATGCGGGCGGCtt < 2:291617/90‑1 (MQ=255) caAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGGGGAAAAATGAAAgcgggggggctt > 1:1525597/1‑86 (MQ=255) caAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGGGAAATCATGAATGCGGGGGGCTTg > 2:630996/1‑90 (MQ=255) caAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGGGAAAACATGAATGCGGGGGGCCTg > 2:466015/1‑90 (MQ=255) aaCGCTGGTGAACGGATCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGGGGAATTATTAAAGCGGGGGGCCTGTc > 1:1308250/1‑90 (MQ=255) aaCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGGGGAAACATGAAAGCGGGGGGCCTggc > 1:80154/1‑88 (MQ=255) cGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGTGAAATCATGAATGCGGGCGGCTTGTCaa < 1:732838/90‑1 (MQ=255) cGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGTGAAATCATGAATGCGGGCGGCTTGTCaa < 1:758797/90‑1 (MQ=255) gCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGGGAAATCATGAATGCGGGCGGCTTGTCaaa > 2:1288609/1‑90 (MQ=255) cTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGTGAAATCATGAATGCGGGCGGCTTGTCaaaa < 2:1212617/90‑1 (MQ=255) tGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGTGAAATCATGAATgcgggc < 2:1343450/77‑1 (MQ=255) tGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGGGAAATCATGAATgcgggc > 1:1343450/1‑77 (MQ=255) gAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGGGGAAACATGAATGCGGGGGGCTTGTCAAAACGCCttgttg > 1:911121/1‑90 (MQ=255) gAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGGGAAATCATGAATGCGGGGGGCTTGTCAAAACGCAttgttt > 2:581385/1‑89 (MQ=255) gCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAATCTTGATTG‑GGGGGGGTTGTCAAAAAGCAttttttgt > 1:982186/1‑87 (MQ=255) ccGATAGCTTCTCCCTGCGGGCGATTCCTTTTTTTGTGCTGGC‑GGGTAAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGccc < 1:921907/90‑1 (MQ=255) cTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGGGAAATCATGAATGCGGGGGGCTTGTCAAAAAGCATTGTTGACCTGCCggt > 2:1029755/1‑88 (MQ=255) tctcCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGGGAAATCATGAATGCGGGGGGCTTGTCAAAACGCATTGTTGACCTGCCGATGa > 2:1300219/1‑90 (MQ=255) tctcCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGGGAAATCATGAATGCGGGGGGCTTGTCAAAAAGCATTTTTGACCCGCCGGTGa > 2:179833/1‑90 (MQ=255) ctcCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGaa < 2:175855/90‑1 (MQ=255) tcCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGGGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTTTTTACCTGCCGATGaaa > 2:170851/1‑90 (MQ=255) cccTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGGGGAAACATTAATGCGGGCGGCCTGTCAAAAAGCATTGTTGACCTGCCGATGgaaa > 1:175855/1‑89 (MQ=255) tgGCGTTTCCGTTTTTTGTGCGGGG‑GGGTAAAATCATGAAGGGGGGGGGCTTTAAAAAACGCATTGTTGACCTGCCGATGAAACtggtgg < 1:1020179/90‑1 (MQ=255) gATTCCTTTCTTTGTGCTGGG‑GGGTGAAATCATGAAGGCGGGCGGTTTGTAAAAACGCATTTTTGACCTGCCGATGAAACTGGTGGGGCa < 2:526599/90‑1 (MQ=255) gATTCCGTTCTTTGTGCTGGC‑GGGTGAAATCATGAAAGCGGGCGGCTTGTCAAAACGCATTGTTTACCTGCCGATGAAACTGGTGGGGCa > 2:797900/1‑90 (MQ=255) tttCCGTTTTTTGTGGGGGG‑GGGTAAAATCATGAATGGGGGGGGTTTGTCAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCAt < 1:466015/89‑1 (MQ=255) aTTCCGTTCTTTGGGCTGGG‑GGGTAAAATCATGAAGGGGGGGGGCTTGAAAAAACGCTTTGTTGCCCTGCCGATGAAACTGGTGGGGCAt < 1:1604566/90‑1 (MQ=255) tttGTGCGGGT‑GGGTGAAATCATGAATGGGGGCGGCTTGTCAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGc < 2:694286/90‑1 (MQ=255) ttGTGCTGGC‑GGGTGGAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGGCGATGAAAATGGTGGGGGATAAACCAGGCg > 2:968896/1‑90 (MQ=255) gtgCGGGC‑GGGTAAAATCATGAATGGGGGGGGTTTGTCAAAACGCTTTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCggg < 1:191025/90‑1 (MQ=255) | GTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGC‑GGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGG > NZ_CP009273/3738274‑3738439 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |