| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,607,584 | G→C | 58.3% | A219A (GCC→GCG) | yfeG ← | predicted DNA‑binding transcriptional regulator |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,607,584 | 0 | G | C | 58.3% | ‑0.5 / 7.6 | 12 | A219A (GCC→GCG) | yfeG | predicted DNA‑binding transcriptional regulator |
| Reads supporting (aligned to +/- strand): ref base G (3/2); new base C (5/2); total (8/4) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01 | |||||||||||
AGCAATCGACGGTAACTCTGATGACTGATGCTTTCCGCCGTCACCATTGGTTGCGCTTCTTCCAGCATGGCCCCCATCGCCATTAGCAAATTATCCCCCAGCACTTTTCGCACTGCAGGCTGATGGAGATTT > minE/1607516‑1607647 | aGCAATCGACGGTAACTCTGATGACTGATGCTTTCCGCCGTCACCATTGGTTGCGCTTCTTCCAGCATg < 1:886970/69‑1 (MQ=255) gACGGTAACTCTGATGACTGATGCTTTCCGCCGTCACCATTGGTTGCGCTTCTTCCAGCATGGCCCCCa > 1:2361997/1‑69 (MQ=255) aaCTCTGATGACTGATGCTTTCCGCCGTCACCATTGGTTGCGCTTCTTCCAGCATg > 1:2227706/1‑56 (MQ=255) aaCTCTGATGACTGATGCTTTCCGCCGTCACCATTGGTTGCGCTTCTTCCAGCATg > 1:3166549/1‑56 (MQ=255) aaCTCTGATGACTGATGCTTTCCGCCGTCACCATTGGTTGCGCTTCTTCCAGCATg > 1:2221374/1‑56 (MQ=255) aaCTCTGATGACTGATGCTTTCCGCCGTCACCATTGGTTGCGCTTCTTCCAGCATg > 1:2163251/1‑56 (MQ=255) aaCTCTGATGACTGATGCTTTCCGCCGTCACCATTGGTTGCGCTTCTTCCAGCATg > 1:1758366/1‑56 (MQ=255) aaCTCTGATGACTGATGCTTTCCGCCGTCACCATTGGATGCGCTTCTTCCAGCATg > 1:336987/1‑56 (MQ=255) gaCTGATGCTTTCCGCCGTCACCATTGGTTGCGCTTCTTCCTGCATGGCCCCCATCGCCATTAGCAAAt < 1:2399661/69‑1 (MQ=255) aCTGATGCTTTCCGCCGTCACCATTGGTTGCGCTTCTTCCTGCATGGCCCCCATCGCCATTAGCAAAt < 1:2633095/68‑1 (MQ=255) gTCACCATTGGTTGCGCTTCTTCCAGCATGGCCCCCATCGCCATTAGCAAATTATCCCCCAGCACttt > 1:2472383/1‑68 (MQ=255) cATTGGTTGCGCTTCTTCCAGCATGGCCCCCATCGCCATTAGCAAATTATCCCCCAGCACTTTTCGCa > 1:1449858/1‑68 (MQ=255) cttcCAGCATCGCCCCCATCGCCATTAGCAAATTATCCCCCAGCACTTTTCGCACTGCAGGCTGa > 1:1842108/1‑65 (MQ=255) cttcCAGCATCGCCCCCATCGCCATTAGCAAATTATCCCCCAGCACTTTTCGCACTGCAGGCTGa > 1:1023545/1‑65 (MQ=255) cttcCAGCATCGCCCCCATCGCCATTAGCAAATTATCCCCCAGCACTTTTCGCACTGCAGGCTGa > 1:1479166/1‑65 (MQ=255) cttcCAGCATCGCCCCCATCGCCATTAGCAAATTATCCCCCAGCACTTTTCGCACTGCAGGCTGa > 1:43869/1‑65 (MQ=255) cttcCAGCATCGCCCCCATCGCCATTAGCAAATTATCCACCAGCACtt > 1:1425215/1‑48 (MQ=38) gCATCGCCCCCATCGCCATTAGCAAATTATCCCCCAGCACTTTTCGCACTGCAGGCTGATGGAGAttt < 1:2034960/68‑1 (MQ=255) gCATCGCCCCCATCGCCATTAGCAAATTATCCCCCAGCACTTTTCGCACTGCAGGCTGATGGAGAttt < 1:351729/68‑1 (MQ=255) | AGCAATCGACGGTAACTCTGATGACTGATGCTTTCCGCCGTCACCATTGGTTGCGCTTCTTCCAGCATGGCCCCCATCGCCATTAGCAAATTATCCCCCAGCACTTTTCGCACTGCAGGCTGATGGAGATTT > minE/1607516‑1607647 |
| 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: ‑ |