| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,834,094 | A→G | 32.3% | D65G (GAT→GGT) | ygbM → | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,834,094 | 0 | A | G | 32.3% | 79.1 / 58.0 | 65 | D65G (GAT→GGT) | ygbM | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base A (20/24); new base G (14/7); total (34/31) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.22e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01 | |||||||||||
ACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCAGACATCGA > minE/1834024‑1834163 | aCCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGgaga < 1:3315725/71‑1 (MQ=255) cTGCAAATCCAAAAGCAACTGGAGCAAACCCATCTGACACTGGCGCTATTCAATACCGCTCCCGGAGatat > 1:2331865/1‑71 (MQ=255) cAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTAtt < 1:3611625/69‑1 (MQ=255) cAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTAtt < 1:914718/69‑1 (MQ=255) cAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATAtt < 1:2046742/69‑1 (MQ=255) cAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATAtt < 1:2792581/69‑1 (MQ=255) cAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATAtt < 1:1543667/69‑1 (MQ=255) cAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATAtt < 1:1895877/69‑1 (MQ=255) tCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCg < 1:2055442/71‑1 (MQ=255) tCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCg < 1:2535731/71‑1 (MQ=255) ccAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCgg > 1:2255739/1‑71 (MQ=255) aaaGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCggggg < 1:3521309/71‑1 (MQ=255) aaaGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCggggg < 1:2273538/71‑1 (MQ=255) cAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCGGGGGAATg < 1:1836460/71‑1 (MQ=255) cAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAAt < 1:3823521/70‑1 (MQ=255) tGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGa > 1:3292893/1‑71 (MQ=255) tGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGa > 1:3034428/1‑71 (MQ=255) tGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGa > 1:2038446/1‑71 (MQ=255) tGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGa > 1:959227/1‑71 (MQ=255) tGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGa > 1:1243599/1‑71 (MQ=255) ggAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGAt > 1:397009/1‑71 (MQ=255) gAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGAtt < 1:3777754/71‑1 (MQ=255) gAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGAtt < 1:1716496/71‑1 (MQ=255) gAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGAtt < 1:2517228/71‑1 (MQ=255) aGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTa < 1:2558303/71‑1 (MQ=255) aGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTa < 1:1425215/71‑1 (MQ=255) cAAAACCATCTGACACTGGCGCTGTTCAATACGGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATc > 1:1696644/1‑71 (MQ=255) aaaaCCATCTGACACTGGCGCAGTTCAATACCGCTCCCCGAGGTATAAATCCGGGGGAATGGGGATTATc > 1:1724596/1‑70 (MQ=255) aaCCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGAGGGGGAATGGGGATTATCCGc < 1:2416559/71‑1 (MQ=255) tCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCGGGGGAATg < 1:952563/54‑1 (MQ=255) tCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAAt < 1:692068/53‑1 (MQ=255) tCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAAt < 1:2728382/53‑1 (MQ=255) gACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTaa > 1:2141198/1‑39 (MQ=255) acTGGCGCTGTTCAATACCGCTCCCGGAGATATTaa > 1:1999277/1‑36 (MQ=255) acTGGCGCTGTTCAATACCGCTCCCGGAGATATTaa > 1:1757218/1‑36 (MQ=255) acTGGCGCTGTTCAATACCGCTCCCGGAGATATTaa > 1:1356681/1‑36 (MQ=255) acTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGa > 1:3159756/1‑52 (MQ=255) acTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGa > 1:197672/1‑70 (MQ=255) ggCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGt > 1:1322505/1‑70 (MQ=255) ttCAATACCGCTCCCGGAGATATTAATGCGGGGGa < 1:2140339/35‑1 (MQ=255) ttCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGccc < 1:561831/52‑1 (MQ=255) ttCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGccc < 1:1538587/52‑1 (MQ=255) ttCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCAc < 1:182876/69‑1 (MQ=255) ttCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCAc < 1:1528281/69‑1 (MQ=255) ttCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCAc < 1:2897877/69‑1 (MQ=255) ttCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGa < 1:947178/71‑1 (MQ=255) ttCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGa < 1:2494029/71‑1 (MQ=255) cAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCTGACGTGAGCACGAAg > 1:2415753/1‑71 (MQ=255) cAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACg > 1:3206158/1‑60 (MQ=255) cAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAg > 1:3251547/1‑71 (MQ=255) cAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAg > 1:3477268/1‑71 (MQ=255) cAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGATTATCCGCCCTTCACGGACGTGAGCACGAAg > 1:3147602/1‑71 (MQ=255) cAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAg > 1:1726625/1‑71 (MQ=255) tACCGCTCCCGGAGATTTTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCAc < 1:2447540/64‑1 (MQ=255) ccGCTCCCGGAGGTATTAATGCGGGGGAATGGGGATTATc > 1:2117047/1‑40 (MQ=255) ccGCTCCCGGAGGTATTAATGCGGGGGAATGGGGATTATc > 1:718418/1‑40 (MQ=255) ccGCTCCCGGAGGTATTAATGCGGGGGAATGGGGATTATc > 1:1424494/1‑40 (MQ=255) ccGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGcac > 1:388448/1‑69 (MQ=255) ccGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGcac > 1:2137130/1‑69 (MQ=255) ccGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCacac > 1:876619/1‑71 (MQ=255) ccGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCacac > 1:2753752/1‑71 (MQ=255) cccGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAg > 1:1840929/1‑54 (MQ=255) cccGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACg > 1:2954791/1‑58 (MQ=255) cccGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGa > 1:1890503/1‑59 (MQ=255) cccGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGa > 1:2419710/1‑59 (MQ=255) cccGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGa > 1:2834609/1‑59 (MQ=255) cGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCAc < 1:617069/55‑1 (MQ=255) gaTATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCAGACATCGa < 1:203511/71‑1 (MQ=255) | ACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCAGACATCGA > minE/1834024‑1834163 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |