| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,522,265 | C→T | 36.9% | G112G (GGG→GGA) | yicG ← | conserved inner membrane protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,522,265 | 0 | C | T | 36.9% | 25.6 / 25.3 | 46 | G112G (GGG→GGA) | yicG | conserved inner membrane protein |
| Reads supporting (aligned to +/- strand): ref base C (17/12); new base T (13/4); total (30/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.38e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01 | |||||||||||
GCGCAAAACGCCGCCAAACACGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCAC > minE/2522202‑2522330 | gcgcAAAACGCCGCCAAACCCGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCat > 1:3427853/1‑69 (MQ=255)gcgcAAAACGCCGCCAAACACGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCat > 1:2443052/1‑69 (MQ=255)gcgcAAAACGCCGCCAAACACGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCat > 1:3360506/1‑69 (MQ=255)gcgcAAAACGCCGCCAAACACGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCat > 1:694488/1‑69 (MQ=255)gcgcAAAACGCCGCCAAACACGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCAGTGCCCCat > 1:506334/1‑69 (MQ=255) ccgccAAACACGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCa < 1:985455/69‑1 (MQ=255) ccgccAAACACGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCa < 1:2641891/69‑1 (MQ=255) cgccAACCACGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCa < 1:475373/68‑1 (MQ=255) cgccAAACACGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCaa < 1:25748/69‑1 (MQ=255) cgccAAACACGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCa < 1:3340595/68‑1 (MQ=255) cacGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCAAc > 1:2199289/1‑63 (MQ=255) cacGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCAACCTGTGc > 1:180955/1‑69 (MQ=255) cacGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCAACCTGTGc > 1:2404236/1‑69 (MQ=255) cacGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCAACCTGTGc > 1:3331406/1‑69 (MQ=255) cacGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCAACCTGTGc > 1:3071270/1‑69 (MQ=255) cacGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCAACCTGTGc > 1:2602994/1‑69 (MQ=255) cacGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCAACCTGTGc > 1:2660885/1‑69 (MQ=255) cacGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCAACCTGTGc > 1:2996379/1‑69 (MQ=255) cacGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCAACCTGTGc > 1:393131/1‑69 (MQ=255) aCCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCAACCTGTGCCCCGATGATa > 1:636723/1‑66 (MQ=255) aCCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCAACCTGTGCCCCGATGATAg > 1:3026965/1‑67 (MQ=255) aCCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCAACCTGTGCCCCGATGATAg > 1:1175152/1‑67 (MQ=255) gccgcGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCAACCTGTGCCCCGATGATAGaaaa < 1:1343284/68‑1 (MQ=255) gcgACAACGGCAATAATTGGGCCGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGacc < 1:882605/69‑1 (MQ=255) gcgACAACGGCAATAATTGGGCCGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGacc < 1:1973775/69‑1 (MQ=255) gcgACAACGGCAATAATTGGGCCGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGacc < 1:2888607/69‑1 (MQ=255) gcgACAACGGCAATAATTGGGCCGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGacc < 1:1531079/69‑1 (MQ=255) ccGTGCCCCATATCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAg < 1:448192/68‑1 (MQ=255) ccGTGCCCCATATCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAg < 1:1224525/68‑1 (MQ=255) ccGTGCCCCATATCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAg < 1:654687/68‑1 (MQ=255) ccGTGCCCCATATCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc < 1:1733038/69‑1 (MQ=255) ccGTGCCCCATATCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc < 1:3359508/69‑1 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGaaaa > 1:2889432/1‑43 (MQ=38) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGCCCACCAGTCCGAGCGCATCGAGc > 1:2270878/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:3098868/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:376349/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:2896782/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:2446755/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:2270547/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:209230/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:2002678/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:1690697/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:736429/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:1612999/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:1198581/1‑68 (MQ=255) gTGCCCCATATCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCAc < 1:63577/69‑1 (MQ=255) | GCGCAAAACGCCGCCAAACACGCCGGTGGTTACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCAC > minE/2522202‑2522330 |
| 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: ‑ |