| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,517,987 | 0 | T | A | 43.8% | 23.8 / 60.8 | 48 | I70F (ATT→TTT) | folC | bifunctional folylpolyglutamate synthase and dihydrofolate synthase |
| Reads supporting (aligned to +/- strand): ref base T (20/7); new base A (21/0); total (41/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.36e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGCGCTCGGTATAACGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAATCGACTCCAGCGTACGGCAGGTGGTGCCTTTGCCATTCGTAC > minE/1517918‑1518028 | cgcgCTCGGTATAACGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAATc < 1:578453/71‑1 (MQ=255) aCGCACGAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:560838/1‑59 (MQ=255) aCGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:2466710/1‑59 (MQ=255) aCGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:3280380/1‑59 (MQ=255) aCGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:4175816/1‑59 (MQ=255) aCGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:2782976/1‑59 (MQ=255) aCGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:2699893/1‑59 (MQ=255) aCGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:265393/1‑59 (MQ=255) aCGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:446048/1‑59 (MQ=255) aCGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:47759/1‑59 (MQ=255) aCGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:2595616/1‑59 (MQ=255) aCGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:2509875/1‑59 (MQ=255) aCGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:2425213/1‑59 (MQ=255) aCGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:2325855/1‑59 (MQ=255) aCGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:2098189/1‑59 (MQ=255) aCGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:1753175/1‑59 (MQ=255) aCGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:858337/1‑59 (MQ=255) aCGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:884750/1‑59 (MQ=255) aCGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:1534029/1‑59 (MQ=255) aCGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:1533033/1‑59 (MQ=255) aCGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:889677/1‑59 (MQ=255) aCGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAAACgc > 1:926744/1‑59 (MQ=255) gCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAATCGACTCCAGCGTACGGCAGGTGGTGcc > 1:813993/1‑70 (MQ=255) gCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAATCGACTCCAGCGTACGGCAGGTGGTGCCt > 1:3175769/1‑71 (MQ=255) gCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAATCGACTCCAGCGTACGGCAGGTGGTGCCt > 1:4020169/1‑71 (MQ=255) gCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAATCGACTCCAGCGTACGGCAGGTGGTGCCt > 1:3799101/1‑71 (MQ=255) gCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAATCGACTCCAGCGTACGGCAGGTGGTGCCt > 1:1268058/1‑71 (MQ=255) gCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAATCGACTCCAGCGTACGGCAGGTGGTGCCt > 1:1706367/1‑71 (MQ=255) gCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAATCGACTCCAGCGTACGGCAGGTGGTGCCt > 1:2574400/1‑71 (MQ=255) gCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAATCGACTCCAGCGTACGGCAGGTGGTGCCt > 1:2599927/1‑71 (MQ=255) gTAGACGCCCACTTTGTACCCTGCCGCCATCAGAATCGACTCCAGCGTACGGCAggtggt < 1:3758046/60‑1 (MQ=255) tCAGAATCGACTCCAGCGTACGGCAGGTGGTGCCTTTGCCATTCGt < 1:3721213/46‑1 (MQ=255) tCAGAATCGACTCCAGCGTACGGCAGGTGGTGCCTTTGCCATTCGt < 1:2220135/46‑1 (MQ=255) tCAGAATCGACTCCAGCGTACGGCAGGTGGTGCCTTTGCCATTCGt < 1:1057930/46‑1 (MQ=255) cAGAATCGACTCCAGCGTACGGCAGGTGGTGCCTTTGCCATTCGt < 1:3526357/45‑1 (MQ=255) aGAATCGACTCCAGCGTACGGCAGGTGGTGCCTTTGCCATTCGt < 1:3531649/44‑1 (MQ=255) aGAATCGACTCCAGCGTACGGCAGGTGGTGCCTTTGCCATTCGTa > 1:205582/1‑45 (MQ=255) aGAATCGACTCCAGCGTACGGCAGGTGGTGCCTTTGCCATTCGTa > 1:3267300/1‑45 (MQ=255) aGAATCGACTCCAGCGTACGGCAGGTGGTGCCTTTGCCATTCGTAc > 1:2631235/1‑46 (MQ=255) aGAATCGACTCCAGCGTACGGCAGGTGGTGCCTTTGCCATTCGTAc > 1:2818769/1‑46 (MQ=255) aGAATCGACTCCAGCGTACGGCAGGTGGTGCCTTTGCCATTCGTAc > 1:550013/1‑46 (MQ=255) aGAATCGACTCCAGCGTACGGCAGGTGGTGCCTTTGCCATTCGTAc > 1:2836961/1‑46 (MQ=255) aGAATCGACTCCAGCGTACGGCAGGTGGTGCCTTTGCCATTCGTAc > 1:372165/1‑46 (MQ=255) aGAATCGACTCCAGCGTACGGCAGGTGGTGCCTTTGCCATTCGTAc > 1:646145/1‑46 (MQ=255) aGAATCGACTCCAGCGTACGGCAGGTGGTGCCTTTGCCATTCGTAc > 1:65837/1‑46 (MQ=255) aGAATCGACTCCAGCGTACGGCAGGTGGTGCCTTTGCCATTCGTAc > 1:861893/1‑46 (MQ=255) aGAATCGACTCCAGCGTACGGCAGGTGGTGCCTTTGCCATTCGTAc > 1:1606248/1‑46 (MQ=255) aGAATCGACTCCAGCGTACGGCAGGTGGTGCCTTTGCCATTCGTAc > 1:1332964/1‑46 (MQ=255) | CGCGCTCGGTATAACGCACCAGATGAGGCGAACTGTAGACGCCCACTTTGTACCCTGCCGCCATCAGAATCGACTCCAGCGTACGGCAGGTGGTGCCTTTGCCATTCGTAC > minE/1517918‑1518028 |
| 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: ‑ |