| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,443,843 | 0 | T | C | 40.6% | 42.5 / 82.4 | 69 | E431E (GAA→GAG) | ilvD | dihydroxyacid dehydratase |
| Reads supporting (aligned to +/- strand): ref base T (20/21); new base C (0/28); total (20/49) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.40e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCTGTAGGCGTGTTCCAGCGAGCGGA > minE/2443785‑2443910 | tgccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTa < 1:2620903/70‑1 (MQ=255)ggccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTa < 1:4649532/71‑1 (MQ=255)ggccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTa < 1:4505832/71‑1 (MQ=255)ggccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTa < 1:4454948/71‑1 (MQ=255)ggccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTa < 1:2638649/71‑1 (MQ=255) gccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAg > 1:3572857/1‑67 (MQ=255) gccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAg > 1:4498068/1‑67 (MQ=255) gccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAg > 1:588065/1‑67 (MQ=255) gccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAg > 1:4290689/1‑67 (MQ=255) gccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGtag > 1:4685860/1‑68 (MQ=255) ccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTAcc < 1:2706261/71‑1 (MQ=255) ccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTAcc < 1:2960280/71‑1 (MQ=255) ccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTAcc < 1:3093379/71‑1 (MQ=255) ccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTAcc < 1:4613248/71‑1 (MQ=255) ccggTGAATTTGAGGAAGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTAcc < 1:1682618/71‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:2481977/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:2871520/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:3015535/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:2652028/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:3178736/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:3282539/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:3570949/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:3599102/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:4625977/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:3789365/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:3834846/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:4087833/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:4323995/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:525205/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:68449/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:1257426/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:1277433/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:826293/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:1157673/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:1545348/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:1587670/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:883823/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:82374/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:2199953/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:1097462/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:2313096/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:2329256/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:78479/69‑1 (MQ=255) gAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTa > 1:798064/1‑69 (MQ=255) gAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTa > 1:229448/1‑69 (MQ=255) gAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTa > 1:113447/1‑69 (MQ=255) gAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTa > 1:1626976/1‑69 (MQ=255) gAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTa > 1:1169459/1‑69 (MQ=255) ggATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTcc > 1:3387701/1‑47 (MQ=255) ggATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGc > 1:2227309/1‑70 (MQ=255) ggATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGc > 1:4266210/1‑70 (MQ=255) ggATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGc > 1:4187558/1‑70 (MQ=255) ggATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGc > 1:4081607/1‑70 (MQ=255) ggATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGc > 1:3587132/1‑70 (MQ=255) ggATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGc > 1:3388268/1‑70 (MQ=255) ggATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACC‑TAGAGCACCGc > 1:788576/1‑69 (MQ=255) ggATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAACTTACCg > 1:4585782/1‑59 (MQ=255) ccGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCTGTAg > 1:1516587/1‑71 (MQ=255) cAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCTGTAGGCGTGTTCCagcg < 1:4906801/70‑1 (MQ=255) cAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCTGTAGGCGTGTTCCagcg < 1:1266543/70‑1 (MQ=255) cAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCTGTAGGCGTGTTCCagcg < 1:1489368/70‑1 (MQ=255) cAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCTGTAGGCGTGTTCCagcg < 1:600829/70‑1 (MQ=255) cAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCTGTAGGCGTGTTCCagcg < 1:3488544/70‑1 (MQ=255) cAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCTGTAGGCGTGTTCCagcg < 1:344489/70‑1 (MQ=255) cAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCTGTAGGCGTGTTCCagcg < 1:2946359/70‑1 (MQ=255) cAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCTGTAGGCGTGTTCCagcg < 1:254618/70‑1 (MQ=255) gTTTTCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCTGTAGGCGTGTTCCAGCGAGCGGa < 1:4816233/71‑1 (MQ=255) gTTTTCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCTGTAGGCGTGTTCCAGCGAGCGGa < 1:3758132/71‑1 (MQ=255) gTTTTCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCTGTAGGCGTGTTCCAGCGAGCGGa < 1:1397113/71‑1 (MQ=255) ttttCCGCAAAGTTACCGTAGAGCACCGCCAGGCCTCCGTCTTTGCt < 1:2100939/47‑1 (MQ=255) ttttCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCt < 1:1038444/47‑1 (MQ=255) ttttCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCt < 1:951516/47‑1 (MQ=255) ttttCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCt < 1:4698225/47‑1 (MQ=255) ttttCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCt < 1:4552703/47‑1 (MQ=255) ttttCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCt < 1:4436592/47‑1 (MQ=255) ttttCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCt < 1:3730824/47‑1 (MQ=255) ttttCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCt < 1:2467746/47‑1 (MQ=255) ttttCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCt < 1:2443447/47‑1 (MQ=255) ttttCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCt < 1:2362251/47‑1 (MQ=255) ttttCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCt < 1:2212931/47‑1 (MQ=255) ttttCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCt < 1:1985375/47‑1 (MQ=255) ttttCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCt < 1:1450750/47‑1 (MQ=255) ttttCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCt < 1:1141328/47‑1 (MQ=255) ttttCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCt < 1:1052108/47‑1 (MQ=255) | GGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTAGAGCACCGCCAGGCCGCCGTCTTTGCTGTAGGCGTGTTCCAGCGAGCGGA > minE/2443785‑2443910 |
| 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: ‑ |