| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,522,289 | C→G | 23.6% | G104G (GGG→GGC) | yicG ← | conserved inner membrane protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,522,289 | 0 | C | G | 23.6% | 116.0 / 33.5 | 72 | G104G (GGG→GGC) | yicG | conserved inner membrane protein |
| Reads supporting (aligned to +/- strand): ref base C (32/23); new base G (13/4); total (45/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.53e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01 | |||||||||||
ACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCAGGTAAGG > minE/2522233‑2522354 | 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:2888607/69‑1 (MQ=255) gcgACAACGGCAATAATTGGGCCGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGacc < 1:1973775/69‑1 (MQ=255) gcgACAACGGCAATAATTGGGCCGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGacc < 1:882605/69‑1 (MQ=255) gcgACAACGGCAATAATTGGGCCGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGacc < 1:1531079/69‑1 (MQ=255) ccGTGCCCCATATCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAg < 1:654687/68‑1 (MQ=255) ccGTGCCCCATATCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAg < 1:1224525/68‑1 (MQ=255) ccGTGCCCCATATCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAg < 1:448192/68‑1 (MQ=255) ccGTGCCCCATATCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc < 1:3359508/69‑1 (MQ=255) ccGTGCCCCATATCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc < 1:1733038/69‑1 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGaaaa > 1:2889432/1‑43 (MQ=38) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGCCCACCAGTCCGAGCGCATCGAGc > 1:2270878/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:1612999/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:1690697/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:736429/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:2446755/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:3098868/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:1198581/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:376349/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:2002678/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:209230/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:2896782/1‑68 (MQ=255) cGTGTCCCATATCCAGTGCAACCTGTGCGCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGc > 1:2270547/1‑68 (MQ=255) gTGCCCCATATCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCAc < 1:63577/69‑1 (MQ=255) tatCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAaca < 1:2875539/69‑1 (MQ=255) tatCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAaca < 1:1105191/69‑1 (MQ=255) tatCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAaca < 1:2390937/69‑1 (MQ=255) tatCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAaca < 1:1383332/69‑1 (MQ=255) tatCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAaca < 1:2282203/69‑1 (MQ=255) tatCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAaca < 1:926476/69‑1 (MQ=255) tatCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAaca < 1:3200086/69‑1 (MQ=255) tatCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAaca < 1:3279395/69‑1 (MQ=255) tatCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAaca < 1:3291794/69‑1 (MQ=255) tatCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAaca < 1:922542/69‑1 (MQ=255) tatCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAaca < 1:591093/69‑1 (MQ=255) atCCAGTGCAACCTGTGCCCCGATGATAGAAAAGac > 1:514465/1‑36 (MQ=255) atCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTc > 1:502154/1‑44 (MQ=255) atCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAacac > 1:681279/1‑69 (MQ=255) atCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAacac > 1:2409512/1‑69 (MQ=255) atCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAacac > 1:2405395/1‑69 (MQ=255) atCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAacac > 1:23946/1‑69 (MQ=255) atCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAacac > 1:235708/1‑69 (MQ=255) atCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAacac > 1:848975/1‑69 (MQ=255) atCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAacac > 1:464657/1‑69 (MQ=255) atCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAacac > 1:2000259/1‑69 (MQ=255) atCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAacac > 1:1915692/1‑69 (MQ=255) atCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAacac > 1:1879526/1‑69 (MQ=255) atCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAacac > 1:1861817/1‑69 (MQ=255) atCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAacac > 1:101450/1‑69 (MQ=255) atCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAacac > 1:1322709/1‑69 (MQ=255) tCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAacac > 1:2675763/1‑68 (MQ=255) tCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACt > 1:618973/1‑69 (MQ=255) tCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACt > 1:1728634/1‑69 (MQ=255) tCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACt > 1:2951207/1‑69 (MQ=255) cAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCa > 1:1432840/1‑69 (MQ=255) cAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCa > 1:1305677/1‑69 (MQ=255) cAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCa > 1:1735905/1‑69 (MQ=255) aCCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCAgg < 1:925402/69‑1 (MQ=255) aCCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCAgg < 1:1096313/69‑1 (MQ=255) aCCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCAgg < 1:2887300/69‑1 (MQ=255) aCCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCAgg < 1:2238478/69‑1 (MQ=255) aCCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCAgg < 1:1634081/69‑1 (MQ=255) gtgCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACCCTTTGCGCAGGTaa > 1:3191025/1‑68 (MQ=255) gtgCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCAGGTaa > 1:44750/1‑68 (MQ=255) gtgCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCAGGTaa > 1:3395034/1‑68 (MQ=255) gtgCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCAGGTaa > 1:2752143/1‑68 (MQ=255) gtgCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCAGGTaa > 1:2397359/1‑68 (MQ=255) gtgCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCAGGTaa > 1:1840747/1‑68 (MQ=255) gtgCCCCGATGATAGAAAAGACCA‑CAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCAGGTAAgg > 1:3149185/1‑69 (MQ=255) ccccGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCAGGTAAgg < 1:595705/67‑1 (MQ=255) ccccGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCAGGTAAgg < 1:2587660/67‑1 (MQ=255) ccccGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCAGGTAAgg < 1:1474227/67‑1 (MQ=255) ccccGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCAGGTAAgg < 1:1108963/67‑1 (MQ=255) ccccGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCAGGTAAgg < 1:94520/67‑1 (MQ=255) | ACCGCCGCGACAACGGCAATAATTGGGCCGTGCCCCATATCCAGTGCAACCTGTGCCCCGATGATAGAAAAGACCACCAGTCCGAGCGCATCGAGCACCAGAAACACTTTGCGCAGGTAAGG > minE/2522233‑2522354 |
| 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: ‑ |