| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 217,643 | T→A | 100% | Q89L (CAG→CTG) | yafU ← | predicted inner membrane protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 217,643 | 0 | T | A | 88.9% | 127.1 / 8.3 | 45 | Q89L (CAG→CTG) | yafU | predicted inner membrane protein |
| Reads supporting (aligned to +/- strand): ref base T (4/1); new base A (40/0); total (44/1) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.11e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
GTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCC > minE/217581‑217659 | ggggAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGtt < 1:646129/64‑1 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGtt > 1:1544035/1‑54 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGtt > 1:2139809/1‑54 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGtttt > 1:792790/1‑62 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGtttt > 1:2100579/1‑62 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGttt > 1:1265519/1‑61 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGttt > 1:2040137/1‑61 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTcc > 1:1937900/1‑67 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTcc > 1:1795929/1‑67 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTcc > 1:1942200/1‑67 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTcc > 1:2055607/1‑67 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTcc > 1:2118137/1‑67 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTcc > 1:509073/1‑67 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTcc > 1:467756/1‑67 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTcc > 1:412208/1‑67 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTcc > 1:409393/1‑67 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTcc > 1:410997/1‑67 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTcc > 1:1853675/1‑67 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTcc > 1:1774118/1‑67 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTcc > 1:655566/1‑67 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTcc > 1:755779/1‑67 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTcc > 1:1682520/1‑67 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTcc > 1:15951/1‑67 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTcc > 1:1496180/1‑67 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTcc > 1:756853/1‑67 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTca > 1:646759/1‑66 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTc > 1:2096703/1‑66 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTc > 1:561614/1‑66 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTc > 1:941219/1‑66 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTc > 1:961240/1‑66 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTc > 1:991320/1‑66 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTc > 1:349389/1‑66 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTc > 1:33576/1‑66 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTc > 1:1079306/1‑66 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTc > 1:1987918/1‑66 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTc > 1:173488/1‑66 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTc > 1:1411584/1‑66 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTc > 1:1176170/1‑66 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTc > 1:1111302/1‑66 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTc > 1:1108833/1‑66 (MQ=255) gAACATACATAAATGGAGTCATGTTTTCCCTTCTCCATTTATCAAGTTCCAGTTGCCGTTTTAGTc > 1:286152/1‑66 (MQ=255) taAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTcc > 1:1754975/1‑58 (MQ=255) taAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTcc > 1:1720005/1‑58 (MQ=255) taAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTcc > 1:1102111/1‑58 (MQ=255) taAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTc > 1:751402/1‑57 (MQ=255) | GTGGAGTGTTAAGAACATACATAAATGGAGTCATGTTTTCCCTTTTCCATTTATCAAGTTCCTGTTGCCGTTTTAGTCC > minE/217581‑217659 |
| 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: ‑ |