| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,515,074 | C→T | 66.7% | E139E (GAG→GAA) | purF ← | amidophosphoribosyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,515,074 | 0 | C | T | 66.7% | 63.8 / 48.7 | 90 | E139E (GAG→GAA) | purF | amidophosphoribosyltransferase |
| Reads supporting (aligned to +/- strand): ref base C (25/5); new base T (6/54); total (31/59) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.59e-12 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01 | |||||||||||
TGTGGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTT > minE/1515011‑1515135 | tgtgGCAGCAATGGCAGCGAAACTATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCt > 1:506017/1‑69 (MQ=255)tgtgGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGc > 1:1692178/1‑68 (MQ=255)tgtgGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGc > 1:400334/1‑68 (MQ=255)tgtgGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGc > 1:2568088/1‑68 (MQ=255)tgtgGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCt > 1:2002063/1‑69 (MQ=255)tgtgGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCt > 1:2717117/1‑69 (MQ=255)tgtgGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCt > 1:3323577/1‑69 (MQ=255)tgtgGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCt > 1:2119379/1‑69 (MQ=255)tgtgGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCt > 1:1070412/1‑69 (MQ=255) aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGGTTGCGAAAata < 1:2650322/69‑1 (MQ=255) aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata < 1:1893354/69‑1 (MQ=255) aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata < 1:1894572/69‑1 (MQ=255) aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata < 1:3344128/69‑1 (MQ=255) aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata < 1:1831299/69‑1 (MQ=255) aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata < 1:3386327/69‑1 (MQ=255) aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata < 1:194461/69‑1 (MQ=255) aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata < 1:341039/69‑1 (MQ=255) aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata < 1:1027692/69‑1 (MQ=255) aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata < 1:2918623/69‑1 (MQ=255) aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata < 1:360076/69‑1 (MQ=255) aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata < 1:1406928/69‑1 (MQ=255) aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata < 1:1375410/69‑1 (MQ=255) aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata < 1:1373223/69‑1 (MQ=255) aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata < 1:532387/69‑1 (MQ=255) aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata < 1:818529/69‑1 (MQ=255) aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata < 1:1260100/69‑1 (MQ=255) aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata < 1:2537015/69‑1 (MQ=255) aTGGCGGCGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata < 1:1234445/69‑1 (MQ=255) ggcggcGAAAATGTTGTCGGCTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata < 1:3091806/67‑1 (MQ=255) atTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGa > 1:2771308/1‑53 (MQ=255) atTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGa > 1:2332282/1‑53 (MQ=255) atTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGa > 1:3145195/1‑53 (MQ=255) atTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGa > 1:1853129/1‑53 (MQ=255) atTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGa > 1:2091478/1‑53 (MQ=255) atTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGa > 1:918423/1‑53 (MQ=255) atTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGa > 1:1656683/1‑53 (MQ=255) cTTCCAGCGGGTAATGGCGGAAGTTGTCCAGTTCGCTTGCGAAAata < 1:2092724/47‑1 (MQ=25) ggTACTGGCGGAAGTTGTCCAGTTCTCTTGCGAAAata < 1:84760/38‑1 (MQ=16) aGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGa < 1:176064/69‑1 (MQ=255) aGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGa < 1:365576/69‑1 (MQ=255) aGTGGCGGAAGGTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGa < 1:495147/69‑1 (MQ=255) gTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGa < 1:2142539/68‑1 (MQ=255) gTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGa < 1:2167487/68‑1 (MQ=255) ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:2561580/68‑1 (MQ=255) ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:24598/68‑1 (MQ=255) ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:126527/68‑1 (MQ=255) ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:2344193/68‑1 (MQ=255) ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:1626131/68‑1 (MQ=255) ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:2282202/68‑1 (MQ=255) ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:3066869/68‑1 (MQ=255) ggCGGAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:1885209/68‑1 (MQ=255) gCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAgt > 1:10641/1‑57 (MQ=255) gCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGa > 1:1239033/1‑65 (MQ=255) gCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGa > 1:1271811/1‑65 (MQ=255) gCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGa > 1:1566859/1‑65 (MQ=255) gCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGa > 1:2983261/1‑65 (MQ=255) gCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtgtg > 1:2544573/1‑69 (MQ=255) gCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtgtg > 1:236374/1‑69 (MQ=255) gCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtgtg > 1:1428367/1‑69 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:3173778/65‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:3229136/65‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:3299601/65‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:1882832/65‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:1846310/65‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:1804484/65‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:1504802/65‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:1264711/65‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:1171863/65‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:964651/65‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:2172000/65‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:2191977/65‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:2203646/65‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:2567625/65‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:2448392/65‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGAtg < 1:2888177/65‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTggc > 1:25598/1‑69 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTggc > 1:164612/1‑69 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTggc > 1:597161/1‑69 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg > 1:2565512/1‑68 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg < 1:2567770/68‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg < 1:783440/68‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg > 1:2717476/1‑68 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg < 1:1371423/68‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg < 1:2319720/68‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg > 1:2813559/1‑68 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg < 1:398712/68‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg < 1:1803357/68‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg < 1:2220680/68‑1 (MQ=255) ggAAGTTGTCCAGTTCGCTTGCGAAAATATTGAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTgg < 1:3040862/68‑1 (MQ=255) tGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGtt > 1:3154759/1‑69 (MQ=255) | TGTGGCAGCAATGGCAGCGAAAATATTGTCGGCTTCCAGCGGGTAGTGGCGGAAGTTGTCCAGCTCGCTGGCGAAGATATTAAGCAGAATTTCCGAGTCGGAAGTGGTGTTGATGTGGCGGCGTT > minE/1515011‑1515135 |
| 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: ‑ |