| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,773,670 | Δ1 bp | 14.0% | intergenic (+59/‑124) | aceA → / → aceK | isocitrate lyase/isocitrate dehydrogenase kinase/phosphatase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,773,670 | 0 | T | . | 14.0% | 127.7 / 20.1 | 50 | intergenic (+59/‑124) | aceA/aceK | isocitrate lyase/isocitrate dehydrogenase kinase/phosphatase |
| Reads supporting (aligned to +/- strand): ref base T (26/17); new base . (5/2); total (31/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.95e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01 | |||||||||||
TCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGG > minE/2773600‑2773736 | tCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAt > 1:2486932/1‑71 (MQ=255) gCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATc < 1:442572/70‑1 (MQ=255) gCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCg > 1:1850798/1‑71 (MQ=255) gTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTg < 1:957818/71‑1 (MQ=255) tCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCa > 1:3790738/1‑71 (MQ=255) aGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCAc < 1:226811/68‑1 (MQ=255) aGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCAc < 1:1614326/68‑1 (MQ=255) aGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑t > 1:1865811/1‑71 (MQ=255) gCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCa > 1:1667529/1‑66 (MQ=255) caacaacaaCCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTTCACGA‑Tg < 1:252253/70‑1 (MQ=255) acaacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCt > 1:1136191/1‑71 (MQ=255) acaacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCt < 1:2673142/71‑1 (MQ=255) aacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGa > 1:1055052/1‑71 (MQ=255) aacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGa > 1:3113033/1‑71 (MQ=255) acaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGAt < 1:1211463/71‑1 (MQ=255) aacCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATGc > 1:670578/1‑71 (MQ=255) acCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTg > 1:2195773/1‑66 (MQ=255) ccGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGa > 1:2646798/1‑71 (MQ=255) cGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGAc > 1:413405/1‑71 (MQ=255) cTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGcgcg > 1:3192539/1‑71 (MQ=255) tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGa < 1:2402179/55‑1 (MQ=255) cgTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAA‑CG‑GTGC‑CGACTGCCTGATGCGACGC‑TTGCGCGTc > 1:1854368/2‑71 (MQ=11) tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATg > 1:1408607/1‑55 (MQ=255) tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATg > 1:2317835/1‑55 (MQ=255) tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCt > 1:675486/1‑71 (MQ=255) aGGCCGGATAAGGCGTTTACGCCGCATCCGGCAA‑CG‑GTGC‑CGACTGCCTGATGCGACGC‑TTGCGCGTCTTa > 1:1395403/1‑71 (MQ=11) ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCtt < 1:2546433/70‑1 (MQ=255) ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGACGA‑TTGCGCGTCTTa < 1:512364/71‑1 (MQ=255) gCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTAt < 1:732176/71‑1 (MQ=255) accGGATAAGGCGTTCACGCCGCATCCGGCAGTC‑‑GTGCGCCA‑TGCCTGATGCGACGCTTTGCGCGTCTTAt < 1:954391/70‑1 (MQ=11) cGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCa < 1:1140473/71‑1 (MQ=255) cGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCa > 1:3720013/1‑71 (MQ=255) ggATAAGGCGTTTACGCCGCATCCGGCAA‑CG‑GTGC‑CGACTGCCTGATGCGACGC‑TTGCGCGTCTTATCagg < 1:306651/71‑3 (MQ=11) aTAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGc > 1:735708/1‑71 (MQ=255) tAAGGCGTTTACGCCGCATCCGGCAA‑CG‑GTGC‑CGACTGCCTGATGCGACGC‑TTGCGCGTCTTATCAGGcc > 1:1940358/1‑70 (MQ=11) tAAGGCGCTCGCGCCGCATCCGGCAA‑CG‑GTGCCAGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCAGGcc > 1:158527/1‑70 (MQ=14) aaGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGcc > 1:2629731/1‑70 (MQ=255) aGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTa > 1:1611301/1‑71 (MQ=255) ggCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTa > 1:951994/1‑70 (MQ=255) ggCGCTCGCGCCGCATCCGGCAA‑CG‑GTGCCAGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCAGGCCTACa < 1:3578553/71‑1 (MQ=11) gCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTa > 1:3285718/1‑69 (MQ=255) cACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTc > 1:1701766/1‑52 (MQ=255) cACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCt > 1:1239757/1‑53 (MQ=255) cACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGt > 1:3047258/1‑71 (MQ=255) gccgcATCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGc < 1:2061280/71‑1 (MQ=255) aTCCGGCAGTCGTGTGCA‑GA‑TGCCTGATGCGACGC‑t < 1:3825615/36‑1 (MQ=0) tCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACg < 1:3985695/71‑1 (MQ=255) tCCGGCAA‑CG‑GTGC‑CGACTGCCTGATGCGACGC‑TTGCGCGTCTTATCa > 1:1573553/1‑48 (MQ=11) cGGCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTa < 1:2009784/71‑1 (MQ=255) gCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAgg > 1:3179210/1‑71 (MQ=255) gCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAgg > 1:2921912/1‑71 (MQ=255) tCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCtt < 1:1073959/37‑1 (MQ=255) | TCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCG‑GTGCACGA‑TGCCTGATGCGACGC‑TTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGG > minE/2773600‑2773736 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |