| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 728,532 | 0 | A | . | 66.0% | 54.2 / 79.7 | 50 | coding (519/954 nt) | fabH | 3‑oxoacyl‑[acyl‑carrier‑protein] synthase III |
| Reads supporting (aligned to +/- strand): ref base A (0/17); new base . (33/0); total (33/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.02e-13 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGGGCGCTGCGGTGCTGGCTGCCTCTGAAGAGCCGGGAATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGC > minE/728495‑728594 | cgGGCGCTGCGGTGCTGGCTGCCTCTGAAGAGCCGGGAATCATTTCCACCCATCTGCATGCCGACGGTAGt < 1:1090543/71‑1 (MQ=255)cgGGCGCTGCGGTGCTGGCTGCCTCTGAAGAGCCGGGAATCATTTCCACCCATCTGCATGCCGACGGTAGt < 1:280649/71‑1 (MQ=255)cgGGCGCTGCGGTGCTGGCTGCCTCTGAAGAGCCGGGAATCATTTCCACCCATCTGCATGCCGACGGTAGt < 1:288943/71‑1 (MQ=255)cgGGCGCTGCGGTGCTGGCTGCCTCTGAAGAGCCGGGAATCATTTCCACCCATCTGCATGCCGACGGTAGt < 1:3015983/71‑1 (MQ=255)cgGGCGCTGCGGTGCTGGCTGCCTCTGAAGAGCCGGGAATCATTTCCACCCATCTGCATGCCGACGGTAGt < 1:2033656/71‑1 (MQ=255)cgGGCGCTGCGGTGCTGGCTGCCTCTGAAGAGCCGGGAATCATTTCCACCCATCTGCATGCCGACGGTAGt < 1:317542/71‑1 (MQ=255)cgGGCGCTGCGGTGCTGGCTGCCTCTGAAGAGCCGGGAATCATTTCCACCCATCTGCATGCCGACGGTAGt < 1:393992/71‑1 (MQ=255)cgGGCGCTGCGGTGCTGGCTGCCTCTGAAGAGCCGGGAATCATTTCCACCCATCTGCATGCCGACGGTAGt < 1:734051/71‑1 (MQ=255)cgGGCGCTGCGGTGCTGGCTGCCTCTGAAGAGCCGGGAATCATTTCCACCCATCTGCATGCCGACGGTAGt < 1:1507574/71‑1 (MQ=255)cgGGCGCTGCGGTGCTGGCTGCCTCTGAAGAGCCGGGAATCATTTCCACCCATCTGCATGCCGACGGTAGt < 1:877098/71‑1 (MQ=255) gggCGCTGCGGTGCTGGCTGCCTCTGAAGAGCCGGGAATCATTTCCACCCATCTGCATGCCGACGGTAGt < 1:2466502/70‑1 (MQ=255) gggCGCTGCGGTGCTGGCTGCCTCTGAAGAGCCGGGAATCATTTCCACCCATCTGCATGCCGACGGTAGt < 1:957018/70‑1 (MQ=255) gggCGCTGCGGTGCTGGCTGCCTCTGAAGAGCCGGGAATCATTTCCACCCATCTGCATGCCGACGGTAGt < 1:2038195/70‑1 (MQ=255) gggCGCTGCGGTGCTGGCTGCCTCTGAAGAGCCGGGAATCATTTCCACCCATCTGCATGCCGACGGTAGt < 1:1983650/70‑1 (MQ=255) ggCGCTGCGGTGCTGGCTGCCTCTGAAGAGCCGGGAATCATTTCCACCCATCTGCATGCCGACGGTAGt < 1:528842/69‑1 (MQ=255) ggTGCTGGCTGCCTCTGAAGAGCCGGGAATCATTTCCACCCATCTGCATGCCGACGGTAGt < 1:2637384/61‑1 (MQ=255) ggTGCTGGCTGCCTCTGAAGAGCCGGGAATCATTTCCACCCATCTGCATGCCGACGGTAGt < 1:1418168/61‑1 (MQ=255) aaGAGCCGG‑‑ATCATTTCCGCCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:1901492/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGa > 1:2136510/1‑35 (MQ=25) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:894747/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:3059961/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:923554/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:2990145/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:2964322/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:2934657/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:3232854/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:383241/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:991279/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:2693798/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:2566979/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:1115911/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:1135063/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:1265061/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:1406492/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:1673667/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:177790/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:1868852/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:187718/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:1936727/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:1973394/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:2140216/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:2249074/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:2292757/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:2314211/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:2334540/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:2359298/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:2479767/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:2531877/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGc > 1:2627000/1‑71 (MQ=255) aaGAGCCGG‑‑ATCATTTCCAACCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCaaa > 1:1977902/1‑68 (MQ=255) | CGGGCGCTGCGGTGCTGGCTGCCTCTGAAGAGCCGGGAATCATTTCCACCCATCTGCATGCCGACGGTAGTTATGGTGAATTGCTGACGCTGCCAAACGC > minE/728495‑728594 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |