Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 349,847 | T→C | intergenic (+275/‑165) | prpB → / → prpC | 2‑methylisocitrate lyase/2‑methylcitrate synthase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 349,847 | 0 | T | C | 93.9% | 112.8 / ‑0.4 | 50 | intergenic (+275/‑165) | prpB/prpC | 2‑methylisocitrate lyase/2‑methylcitrate synthase |
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (24/22); total (26/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01 |
ACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAC > NC_000913/349719‑349969 | aCTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGAc > 1:1660374/1‑139 (MQ=17) tATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCa < 2:1675730/139‑1 (MQ=17) cTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCg < 2:1004737/139‑1 (MQ=17) cACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACaa > 2:71996/1‑139 (MQ=17) cGGATAAGACGTATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgcg > 1:224165/1‑139 (MQ=25) cGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgcg > 1:981400/1‑139 (MQ=255) cGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgcg > 2:293401/1‑139 (MQ=255) aaGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCt < 2:1235317/139‑1 (MQ=37) gACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa < 2:1277048/139‑1 (MQ=37) gACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa < 2:171203/139‑1 (MQ=37) cGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTa > 1:1022685/1‑139 (MQ=18) gTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTAt < 1:293401/139‑1 (MQ=17) gTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTAt < 2:732474/139‑1 (MQ=17) gTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTAt < 2:1660374/139‑1 (MQ=18) tCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGc < 1:23614/92‑1 (MQ=18) tCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGc > 2:23614/1‑92 (MQ=18) tCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATc > 2:847390/1‑139 (MQ=17) cATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAgg > 1:1496851/1‑139 (MQ=18) cAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAg > 1:890851/1‑73 (MQ=255) cAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAg < 2:890851/73‑1 (MQ=19) aaCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGATCTGTAGGACGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACggg < 1:847390/139‑1 (MQ=17) aaCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTAcagg < 2:1496851/139‑4 (MQ=16) ccGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGt < 2:72705/139‑1 (MQ=17) cGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtg < 2:224165/139‑1 (MQ=18) aTGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccg < 2:643651/108‑1 (MQ=33) aTGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccg > 1:643651/1‑108 (MQ=33) aTGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAg > 2:930584/1‑125 (MQ=18) aTGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAg < 1:930584/125‑1 (MQ=17) gCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtgaact < 2:1309566/139‑1 (MQ=37) gCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGCGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATACGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtgaact > 1:1083996/1‑139 (MQ=255) ccTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtgaactg > 2:1151988/1‑139 (MQ=25) gccgcTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAg > 2:1495054/1‑139 (MQ=21) ccgcTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAgg > 2:1793851/1‑139 (MQ=37) ccgcTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAgg > 1:1650106/1‑139 (MQ=25) gcTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTc < 2:1022685/139‑1 (MQ=21) gcTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTc > 2:332911/1‑139 (MQ=37) tGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCgg < 2:214406/139‑1 (MQ=37) cgACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAg < 1:332911/139‑1 (MQ=18) cGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtgaact < 2:1652171/104‑1 (MQ=18) cGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtgaact > 1:1652171/1‑104 (MQ=18) cGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCAt > 2:349835/1‑139 (MQ=18) ggggTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATc > 2:533639/1‑139 (MQ=17) tGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGac > 1:1716625/1‑139 (MQ=18) tGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGac > 1:549415/1‑139 (MQ=17) gCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCACATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtgaactgaa > 1:1437009/1‑101 (MQ=255) gCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCACATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtgaactgaa < 2:1437009/101‑1 (MQ=12) tGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGt < 2:660019/96‑1 (MQ=21) tGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGt > 1:660019/1‑96 (MQ=18) tAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAgg > 2:108401/1‑97 (MQ=18) tAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAgg < 1:108401/97‑1 (MQ=21) | ACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAC > NC_000913/349719‑349969 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |