| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,378,968 | 0 | C | T | 55.9% | 2.8 / 30.0 | 34 | E98K (GAA→AAA) | menD | 2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylate synthase; SEPHCHC synthase |
| Reads supporting (aligned to +/- strand): ref base C (4/11); new base T (19/0); total (23/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.77e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.40e-06 | |||||||||||
GGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCA‑GTTTTTCTCC‑GGTTAA‑CCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAAT > NC_000913/2378828‑2379105 | ggggATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGGTTAGAGGAACA‑GTGTTTTTCC‑GGTTaa > 1:271310/1‑150 (MQ=255) gggATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTTAGGGAATCA‑GTTTTTTTCC‑GGGTAA‑c > 2:642238/1‑150 (MQ=255) gggATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCC‑GTTTTTCTCC‑GGTTAA‑c > 1:537476/1‑150 (MQ=255) gggATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGCTTAGAGAATCAGGTTTTTTCCC‑GGTTac > 1:615631/1‑149 (MQ=255) ggATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCTCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCA‑GTTTTTCTCC‑GGTTAT‑cc < 2:490667/150‑1 (MQ=255) ccTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGGGAATCA‑GTTTTTTTCC‑CGGTAA‑CCCCGcct > 2:510918/1‑148 (MQ=255) gcgcGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCA‑GTTTTTCTCC‑GGTTAA‑CCCGGCTTCAATCAGTGc < 2:565728/150‑1 (MQ=255) gCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCA‑GTTTTTCTCC‑GGTTAA‑CCCGGCTTCAATCAGTGCCGGAt < 1:530601/150‑1 (MQ=255) aTGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTTAGAGAATCA‑AGTTTTTTCC‑GGTTAA‑CCCCGGTTCAATCAGTGCCCGGTagag > 1:528949/1‑151 (MQ=255) aTGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTTAGAGAAACC‑GTTTTTTTCC‑CGGTAA‑ACCCGCCTCAATCAGTGCCGGGTagag > 1:346816/1‑151 (MQ=255) aTGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGGGAATCA‑GGTTTTTTCC‑GGGTAA‑CCCCGCTTCAATCAGTGCCgggaaaag > 2:526223/1‑145 (MQ=255) aTGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCA‑GTTTTTTTCC‑GGTTAA‑CCCGGCTTCAATCAATGGCGGGAagag > 2:557448/1‑151 (MQ=255) aTGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAAACAGGTTTTTTTCC‑GGTTAA‑CCGGGCTT‑AAACAATGCCGGGTagag > 2:610932/1‑151 (MQ=255) aTGAAATACTGTGCGTGGGATGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTTAGAGAATCA‑GGTTTTCTCC‑CGTTAA‑CCCCGCTTTAATCAAAGCCGGGTaaag > 2:68095/1‑151 (MQ=255) tGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCA‑GTTTTTCTCC‑GGTTAA‑CCCGGCTTCAATCAGTGCCGGATAgaga < 1:220949/151‑1 (MQ=255) gAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTTAGAGAATCA‑GGTTTTTTCC‑CGTTAA‑CCCGGCCTCACTCAGTGTCCGGTAgagag > 1:608300/1‑150 (MQ=255) gAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCA‑GTTTTTTTCC‑CGGTAA‑CCCGGCTTTCATCACGGCCGGGTAGAGAt > 2:122507/1‑151 (MQ=255) gAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGGTTAGAGAATCC‑GGTTTTTTCC‑GGTTAA‑CCCCGGTTCCATCAGCGGCGGATAAAGAt > 2:173622/1‑151 (MQ=255) gAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGGTAAAAGAATCA‑GTTTTTTTTC‑GGGTAA‑ACCCGCTTTAATCAGTGCCGGGGAGAGAt > 1:644280/1‑151 (MQ=255) gtgCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAAGGAATTA‑GGTTTTTTCCTGGTTAA‑CCCCGCTTC‑ATCAGGTCCGGGTAGAGAttttgcac > 2:345840/1‑145 (MQ=255) cGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCA‑GTTTTTCTCC‑GGTTAA‑CCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCAccgc < 2:34529/150‑1 (MQ=255) aCATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCA‑GTTTTTCTCC‑GGTTAA‑CCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAAt < 1:208376/151‑1 (MQ=255) cATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTTAGAGAATCA‑GTTTTTTTCC‑GGCTAC‑CCCCGCTTCCCTCAGTGCCCGGTAAAGATTTTCCCCCCCCCGGCCCGCGGTcccaca > 2:272863/1‑148 (MQ=255) ttCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCA‑GTTTTTCTCC‑GGTTAA‑CCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCAc < 2:41978/151‑1 (MQ=255) cGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCA‑GTTTTTCTCC‑GGTTAC‑CCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGgc < 2:350268/150‑1 (MQ=255) gATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCC‑CTTTTTTTCC‑GGTTAACCCCGGCTT‑AATCCAGGCCGGGTAGAGATTTCCCCCCCCCGCGCCGGCGGCCACAATCACCGCCCCCCGCTGCTTTCTgtt > 2:472136/1‑149 (MQ=255) ccgcAGTCAATTAGCTCCGGCGGGCGATCGGCGGGTAANANAANCA‑GGTTTTTTCC‑GGGTAA‑CCCCGGTTCAATCAGTGCCCGGTAGAGATTTTCCACCGCCGTGCCGGGGGGCCAAATCACCGCCCCCCGCTGCTTTCTGTCTTTCTCCa > 1:627400/1‑151 (MQ=255) gTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGGATCA‑GTTTTTTTCC‑CGTTAA‑CCCGGCTTCAATCAGTGCCGGGTAGAGATTTTCCCCCGCCGTGCCGGAGGGCACAATCACCCCCCCCCGCTGCCTGCTTCCTTTCCCCCGccc > 2:95532/1‑150 (MQ=255) gCTCCGGCGGGCGATCGGCGGGTAAGAGGATTA‑GTTTTTTTCC‑GGTTAA‑CCCCGCTTCAATCAGTGCCGGGTAGAGATTTGCCACCCCCCTGCCGGAGGGCACAATCACCGCCCCCCGCCGCTTTCTGTCCTTTGCCAGCCCCCGCGccca > 1:626171/1‑149 (MQ=255) cccGGCGGGCGATCGGCGGTTAAGAGAATCA‑GTTTTTCTCC‑GGTTAA‑CCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGAt < 1:122482/150‑1 (MQ=255) cgATCGGCGGTTAAGAGAATCA‑GTTTTTCTCC‑GGTTAA‑CCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAAcc < 1:305831/151‑1 (MQ=255) gcggTTAAGAGAATCA‑GGTTTTCTCC‑GGTTAA‑CCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCCCCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGc > 2:285669/1‑150 (MQ=255) gagaATCA‑GTTTTTCTCC‑GGTTAA‑CCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGaaa > 2:458533/1‑151 (MQ=255) aaTCA‑TTTTTTCCCC‑GGTTAA‑CCCGGCTTCAATCAGGGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAAt < 2:515231/149‑1 (MQ=255) | GGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGTTAAGAGAATCA‑GTTTTTCTCC‑GGTTAA‑CCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGCTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAAT > NC_000913/2378828‑2379105 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |