| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 349,785 | 0 | A | G | 25.4% | 83.9 / 29.6 | 55 | intergenic (+213/‑227) | prpB/prpC | 2‑methylisocitrate lyase/2‑methylcitrate synthase |
| Reads supporting (aligned to +/- strand): ref base A (19/22); new base G (5/9); total (24/31) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.47e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01 | |||||||||||
GATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATG > NC_000913/349695‑349884 | gttgtcTGATGCGTCGCTGACGTGACTCACCAGGCCTACGGGGTGCTCGGAACTGTAGGTCGCATAATTTGCATTGTTTCGCATTCGACAAGCGaaac < 1:860384/93‑4 (MQ=2) cgccgcTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCt < 2:564040/100‑1 (MQ=33) ggTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCAc < 2:1197976/100‑1 (MQ=32) gTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACt > 1:681410/1‑100 (MQ=33) tGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTg < 1:1652508/100‑1 (MQ=32) gCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCCCTGa > 2:817176/1‑100 (MQ=17) cACTGCCCTGTACGTCGGACAAGCCGCCTTTCGTCGCATCCGACAGCCGATGCCTGATGCGTCGCTGACGCTACTTATCAGGCCTACGGGGTGCACt < 1:1462172/97‑1 (MQ=11) cACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGaa > 2:1167592/1‑100 (MQ=32) aCTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAAc < 2:938773/100‑1 (MQ=33) aCTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAAc < 2:397711/100‑1 (MQ=32) aCTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATCCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAAc > 1:1034439/1‑100 (MQ=17) cTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACt < 1:629489/100‑1 (MQ=33) gAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGt > 2:484185/1‑100 (MQ=32) atcTGTAGGTCGGATAAGACGGATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGt < 2:1062363/97‑1 (MQ=18) aaCTGTAGGTCGGATAAGACGGATAGCGGCGCATCCGAGAACCGCTGCCTGATGTGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACt < 1:796465/97‑1 (MQ=11) aaCTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTa < 1:753586/100‑1 (MQ=32) aGGTCTGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCg < 1:257187/99‑1 (MQ=17) aGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCgg < 1:1439797/100‑1 (MQ=33) ggTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTAAACTGTAGGTCGGa < 2:979524/100‑1 (MQ=18) gTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATACGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGAt > 2:1047766/1‑100 (MQ=17) ggATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAg > 2:1264045/1‑100 (MQ=32) gATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa > 1:692249/1‑100 (MQ=32) gATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa < 1:1037747/100‑1 (MQ=33) gATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa > 1:391238/1‑100 (MQ=33) aTAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGAc < 2:136916/100‑1 (MQ=33) aGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCATGCCTACGGGGTGCACTCAACTGTAGGTCGAATAAGACGCa > 2:853038/1‑100 (MQ=255) aGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCa > 1:1383009/1‑100 (MQ=32) cgtatggcttCGCATCCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAg < 2:624313/91‑1 (MQ=1) cGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATa > 2:1697294/1‑99 (MQ=33) cATAGCGTCGCATCCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCg < 1:1793643/100‑1 (MQ=12) cATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCg > 2:1610603/1‑100 (MQ=33) cATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCg > 2:290866/1‑100 (MQ=32) cATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCg < 1:992293/100‑1 (MQ=33) atggCGTCGCATCCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGt < 2:1657052/97‑1 (MQ=255) aGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCg > 1:1672103/1‑100 (MQ=33) gCGTCGCATCCGATAGCCGATGCCGGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGc < 2:1285482/100‑1 (MQ=11) gCGTCGCATCCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGc > 1:237170/1‑100 (MQ=12) cGTCGCATTCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCg < 1:473477/98‑1 (MQ=12) gTCGCATCCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCAt < 1:1576093/100‑1 (MQ=12) gTCGCATCCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCAt < 2:420793/100‑1 (MQ=12) gTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCAt > 1:1205794/1‑100 (MQ=32) gTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCAt < 2:373196/100‑1 (MQ=32) tCGCATCCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATc < 2:1297948/100‑1 (MQ=255) cGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGGCGCATAGCGTCGCATcc > 1:415270/1‑100 (MQ=17) gCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCg > 2:292739/1‑100 (MQ=32) cATCCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCCCTGACCTGTAGGTCGGCTAAGACGCATAGCGTCGCATCCGa > 2:306588/1‑100 (MQ=11) cATCCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa > 2:713806/1‑100 (MQ=255) cATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa < 1:1734831/100‑1 (MQ=33) aTCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAc < 1:506479/100‑1 (MQ=33) cGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAc > 2:1037757/1‑100 (MQ=255) cGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCACATCCGACAAc > 2:531756/1‑100 (MQ=18) gACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGCTCGGATAAGACGCATAGCGTCGCATCCGACAAc > 1:958299/1‑99 (MQ=11) gACAAGCGATGTCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAcc < 2:188703/100‑1 (MQ=16) gACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAcc < 2:1260921/100‑1 (MQ=32) aCAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCg < 2:681414/100‑1 (MQ=32) aCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAc > 1:1621225/1‑95 (MQ=32) aCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAt > 2:1549197/1‑99 (MQ=33) aCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAGATg > 2:677727/1‑100 (MQ=18) | GATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATG > NC_000913/349695‑349884 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |