Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,443,798 | 0 | G | C | 35.9% | 79.2 / 80.6 | 78 | L446L (CTC→CTG) | ilvD | dihydroxyacid dehydratase |
Reads supporting (aligned to +/- strand): ref base G (29/21); new base C (0/28); total (29/49) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.06e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TACCGCCGAGAATCGCTTCTACCGCATCGTCCTGGCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTA > minE/2443731‑2443860 | tACCGCCGAGAATCGCTTCTACCGCATCGTCCTGGCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGa < 1:4022479/69‑1 (MQ=255) tACCGCCGAGAATCGCTTCTACCGCATCGTCCTGGCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGa < 1:2322991/69‑1 (MQ=255) tACCGCCGAGAATCGCTTCTACCGCATCGTCCTGGCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAg < 1:1578005/70‑1 (MQ=255) tACCGCCGAGAATCGCTTCTACCGCATCGTCCTGGCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAg < 1:2177242/70‑1 (MQ=255) gagaATCGCTTCTACCGCATCGTCCTGGCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGt > 1:4383172/1‑71 (MQ=255) gaATCGCTTCTACCGCATCGTCCTGGCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCa > 1:521045/1‑71 (MQ=255) gaATCGCTTCTACCGCATCGTCCTGGCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCa > 1:192341/1‑71 (MQ=255) cGCTTCTACCGCATCGTCCTGGCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGa > 1:80546/1‑71 (MQ=255) cGCTTCTACCGCATCGTCCTGGCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGa > 1:4371271/1‑71 (MQ=255) cTTCTACCGCATCGTCCTGGCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAg < 1:3556088/55‑1 (MQ=255) ggCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATg > 1:995224/1‑71 (MQ=255) ggCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATg > 1:4726078/1‑71 (MQ=255) ggCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATg > 1:468400/1‑71 (MQ=255) ggCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATg > 1:3686239/1‑71 (MQ=255) ggCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATg > 1:3595796/1‑71 (MQ=255) ggCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATg > 1:2231432/1‑71 (MQ=255) tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAg < 1:77511/69‑1 (MQ=255) tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAg < 1:1697054/69‑1 (MQ=255) tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAg < 1:58512/69‑1 (MQ=255) tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAg < 1:3789686/69‑1 (MQ=255) tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAg < 1:2452915/69‑1 (MQ=255) tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAg < 1:3438269/69‑1 (MQ=255) tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGcc > 1:1041811/1‑71 (MQ=255) tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGcc > 1:4661111/1‑71 (MQ=255) tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGcc > 1:3748545/1‑71 (MQ=255) tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGcc > 1:3066121/1‑71 (MQ=255) tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGcc > 1:2963267/1‑71 (MQ=255) tCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGcc > 1:2940503/1‑71 (MQ=255) tgccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTa < 1:2620903/70‑1 (MQ=255) ggccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTa < 1:4454948/71‑1 (MQ=255) ggccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTa < 1:4505832/71‑1 (MQ=255) ggccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTa < 1:2638649/71‑1 (MQ=255) ggccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTa < 1:4649532/71‑1 (MQ=255) gccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCg > 1:210122/1‑38 (MQ=255) gccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAg > 1:4498068/1‑67 (MQ=255) gccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAg > 1:4290689/1‑67 (MQ=255) gccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAg > 1:588065/1‑67 (MQ=255) gccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAg > 1:3572857/1‑67 (MQ=255) gccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGtag > 1:4685860/1‑68 (MQ=255) gccggAGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCa > 1:1696116/1‑51 (MQ=255) ccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTAcc < 1:2960280/71‑1 (MQ=255) ccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTAcc < 1:2706261/71‑1 (MQ=255) ccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTAcc < 1:4613248/71‑1 (MQ=255) ccggTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTAcc < 1:3093379/71‑1 (MQ=255) ccggTGAATTTGAGGAAGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTAcc < 1:1682618/71‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:1157673/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:525205/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:82374/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:2313096/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:68449/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:4625977/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:826293/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:1097462/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:2199953/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:78479/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:883823/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:1277433/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:4323995/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:2871520/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:2652028/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:3015535/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:3178736/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:3282539/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:2481977/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:3570949/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:3599102/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:1587670/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:2329256/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:3789365/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:3834846/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:1545348/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:4087833/69‑1 (MQ=255) gagAACTTCAGGATGCTTTCGTCCACGCCGGCGGTCTTGACGATGCAGCCCTTCTCGGCGATGTTGCCGTa < 1:1257426/69‑1 (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:798064/1‑69 (MQ=255) gAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTa > 1:1169459/1‑69 (MQ=255) | TACCGCCGAGAATCGCTTCTACCGCATCGTCCTGGCTTTCGTACACTTTCGCCGGGCCGGTGAATTTGAGGATGCTGTCATCGACGCCTGCCGTTTTCACGATGCAGCCGTTTTCCGCAAAGTTACCGTA > minE/2443731‑2443860 |
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: ‑ |