| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,606,282 | T→G | 57.6% | intergenic (+3/‑1) | amiA → / → hemF | N‑acetylmuramoyl‑l‑alanine amidase I/coproporphyrinogen III oxidase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,606,282 | 0 | T | G | 57.6% | 35.7 / 57.8 | 66 | intergenic (+3/‑1) | amiA/hemF | N‑acetylmuramoyl‑l‑alanine amidase I/coproporphyrinogen III oxidase |
| Reads supporting (aligned to +/- strand): ref base T (2/26); new base G (22/16); total (24/42) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.01e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01 | |||||||||||
GGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCT > minE/1606220‑1606350 | ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa < 1:1881221/69‑1 (MQ=255)ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa < 1:2051070/69‑1 (MQ=255)ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa < 1:2134287/69‑1 (MQ=255)ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa < 1:1748171/69‑1 (MQ=255)ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa < 1:2135826/69‑1 (MQ=255)ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa < 1:1533701/69‑1 (MQ=255)ggCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa < 1:2358139/69‑1 (MQ=255) cGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGaaa < 1:1301254/67‑1 (MQ=255) gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc < 1:1104405/69‑1 (MQ=255) gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc < 1:2373057/69‑1 (MQ=255) gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc < 1:2807089/69‑1 (MQ=255) gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc < 1:1072763/69‑1 (MQ=255) tGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc < 1:2874272/68‑1 (MQ=255) aGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGca < 1:1524868/69‑1 (MQ=255) aGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGca < 1:2061367/69‑1 (MQ=255) ttATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGca < 1:494087/67‑1 (MQ=255) tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACg < 1:673727/64‑1 (MQ=255) tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACg < 1:2049228/64‑1 (MQ=255) tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACg < 1:1976084/64‑1 (MQ=255) tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCacac < 1:2186453/69‑1 (MQ=255) tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCacac < 1:7288/69‑1 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGtgatga > 1:965430/1‑46 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAGACCCGACGCACACCAGGTTaa > 1:1631841/1‑68 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGcaca > 1:1276619/1‑59 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa > 1:249986/1‑68 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa > 1:3328248/1‑68 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa > 1:3350073/1‑68 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa > 1:3365461/1‑68 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa > 1:3428131/1‑68 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa > 1:1989972/1‑68 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa > 1:283026/1‑68 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa > 1:2815214/1‑68 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa > 1:2734313/1‑68 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa > 1:2703380/1‑68 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa > 1:2401693/1‑68 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa > 1:666148/1‑68 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa > 1:1529633/1‑68 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa > 1:1175700/1‑68 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa > 1:1765700/1‑68 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa > 1:2272388/1‑68 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa > 1:1611271/1‑68 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTaa > 1:2203143/1‑68 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTa > 1:2378217/1‑67 (MQ=255) tGGTTTGACAACCAGAAAGCACACTCGAGAAAGCGCTAAGTGATGAAACCCGACGCACACCAGGTTaa > 1:878926/1‑68 (MQ=255) ggTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAAc < 1:2904586/69‑1 (MQ=255) cGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCg > 1:2788395/1‑47 (MQ=255) cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc < 1:364641/68‑1 (MQ=255) cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc < 1:124013/68‑1 (MQ=255) cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc < 1:1178808/68‑1 (MQ=255) cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc < 1:3433238/68‑1 (MQ=255) cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc < 1:784282/68‑1 (MQ=255) cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc < 1:612957/68‑1 (MQ=255) cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc < 1:932091/68‑1 (MQ=255) cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc < 1:1315863/68‑1 (MQ=255) cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc < 1:1951301/68‑1 (MQ=255) cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc < 1:3000247/68‑1 (MQ=255) cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc < 1:135215/68‑1 (MQ=255) cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc < 1:1844939/68‑1 (MQ=255) cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc < 1:2521342/68‑1 (MQ=255) cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc < 1:1764234/68‑1 (MQ=255) cAACCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc < 1:1413607/68‑1 (MQ=255) aaCCAGAAAGCACACTCGAGAAAGCGGTAAGTGATGAAACCCGACGCACACCAGGTTAAACAGTTTc < 1:2040196/67‑1 (MQ=255) aaaGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAt > 1:2735172/1‑69 (MQ=255) aGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagca < 1:633066/69‑1 (MQ=255) aGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagca < 1:901962/69‑1 (MQ=255) aGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagca < 1:885422/69‑1 (MQ=255) aGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagca < 1:3422448/69‑1 (MQ=255) tATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt < 1:2264820/69‑1 (MQ=255) tATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt < 1:2167927/69‑1 (MQ=255) tATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt < 1:3012266/69‑1 (MQ=255) tATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt < 1:535753/69‑1 (MQ=255) tATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt < 1:420518/69‑1 (MQ=255) tATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt < 1:2831933/69‑1 (MQ=255) tATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt < 1:3001603/69‑1 (MQ=255) | GGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCT > minE/1606220‑1606350 |
| 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: ‑ |