| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,111,517 | 0 | C | A | 36.9% | 95.7 / 57.8 | 65 | A54E (GCG→GAG) | agaI | galactosamine‑6‑phosphate isomerase |
| Reads supporting (aligned to +/- strand): ref base C (12/29); new base A (0/24); total (12/53) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.39e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAG > minE/2111456‑2111575 | cGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGc < 1:2773492/71‑1 (MQ=255)cGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGc < 1:580151/71‑1 (MQ=255)cGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGc < 1:4175852/71‑1 (MQ=255) gTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGc < 1:2023883/70‑1 (MQ=255) ccAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCAc < 1:4169638/70‑1 (MQ=255) aaTATTTATTGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCAtt > 1:1311414/1‑70 (MQ=255) tatttattGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACt > 1:3299494/1‑71 (MQ=255) tatttattGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACt > 1:588509/1‑71 (MQ=255) tatttattGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACt > 1:1931305/1‑71 (MQ=255) tatttattGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACt > 1:4212599/1‑71 (MQ=255) tattGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCa > 1:3150437/1‑62 (MQ=255) tattGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCa > 1:2994170/1‑62 (MQ=255) tattGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCa > 1:2299572/1‑62 (MQ=255) tattGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCa > 1:4675944/1‑62 (MQ=255) tattGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCa > 1:4849573/1‑62 (MQ=255) tattGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCa > 1:508985/1‑62 (MQ=255) tattGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCa > 1:1536746/1‑62 (MQ=255) ttGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGt < 1:2230393/70‑1 (MQ=255) ttGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGt < 1:4699148/70‑1 (MQ=255) ttGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGt < 1:2719336/70‑1 (MQ=255) ttGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGt < 1:3731486/70‑1 (MQ=255) ttGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACAGACGt < 1:845571/70‑1 (MQ=255) tGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGt < 1:933895/63‑1 (MQ=255) gATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGTCCGGAGCCTCGCCATTACTGACGTATCAtt < 1:3799941/68‑1 (MQ=255) gATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCAtt < 1:4256257/68‑1 (MQ=255) gATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCAtt < 1:4055568/68‑1 (MQ=255) gATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCAtt < 1:3938187/68‑1 (MQ=255) gATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCAtt < 1:3564088/68‑1 (MQ=255) gATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCAtt < 1:1256155/68‑1 (MQ=255) gATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCAtt < 1:3067324/68‑1 (MQ=255) gATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCAtt < 1:1034549/68‑1 (MQ=255) gATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCAtt < 1:295930/68‑1 (MQ=255) gATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCAtt < 1:2158347/68‑1 (MQ=255) gATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCAtt < 1:292732/68‑1 (MQ=255) gATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCAtt < 1:152638/68‑1 (MQ=255) gATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCAtt < 1:163906/68‑1 (MQ=255) aaCCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGaaaaa < 1:551734/71‑1 (MQ=255) aaCCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGaaaaa < 1:4859153/71‑1 (MQ=255) ccGAATGCCGTGATTTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:989243/54‑1 (MQ=255) ccGAATGCCGTGATTTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:585098/54‑1 (MQ=255) ccGAATGCCGTGATTTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:4925063/54‑1 (MQ=255) ccGAATGCCGTGATTTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:4898190/54‑1 (MQ=255) ccGAATGCCGTGATTTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:696279/54‑1 (MQ=255) ccGAATGCCGTGATTTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:1743194/54‑1 (MQ=255) ccGAATGCCGTGATTTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:2135318/54‑1 (MQ=255) ccGAATGCCGTGATTTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:4542465/54‑1 (MQ=255) ccGAATGCCGTGATTTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:1369889/54‑1 (MQ=255) ccGAATGCCGTGATTTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:4484204/54‑1 (MQ=255) ccGAATGCCGTGATTTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:4371621/54‑1 (MQ=255) ccGAATGCCGTGATTTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:3718574/54‑1 (MQ=255) ccGAATGCCGTGATTTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:235000/54‑1 (MQ=255) ccGAATGCCGTGATTTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:4099230/54‑1 (MQ=255) ccGAATGCCGTGATTTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:2354565/54‑1 (MQ=255) ccGAATGCCGTGATTTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:4030972/54‑1 (MQ=255) ccGAATGCCGTGATTTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:2538547/54‑1 (MQ=255) ccGAATGCCGTGATTTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:3744003/54‑1 (MQ=255) ccGAATGCCGTGATTTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:2553175/54‑1 (MQ=255) ccGAATGCCGTGATTTGCCTCGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:3674953/54‑1 (MQ=255) ccGAATGCCGTGATCTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:2206608/54‑1 (MQ=255) cGAATGCCGTGATTTGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:3852638/53‑1 (MQ=255) tGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCAtt < 1:35312/45‑1 (MQ=255) tGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACcagcag < 1:1605516/71‑1 (MQ=255) tGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACcagcag < 1:1854966/71‑1 (MQ=255) tttGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:3275197/41‑1 (MQ=255) tGCCTGGAGACCGGAGCCACGCCATTACTGACGTATCAt < 1:2325941/39‑1 (MQ=255) | CGTGCCAGCGAATATTTATTGGCCGTGATCCGTAGCAAACCGAATGCCGTGATTTGCCTGGCGACCGGAGCCACGCCATTACTGACGTATCATTATCTGGTAGAAAAAATCCACCAGCAG > minE/2111456‑2111575 |
| 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: ‑ |