| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 835,946 | A→T | 22.4% | M222L (ATG→TTG) | narH → | nitrate reductase 1, beta (Fe‑S) subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 835,946 | 0 | A | T | 22.4% | 135.0 / 37.5 | 67 | M222L (ATG→TTG) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
| Reads supporting (aligned to +/- strand): ref base A (23/29); new base T (6/9); total (29/38) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01 | |||||||||||
TGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAT > minE/835879‑836016 | tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTa < 1:2532702/68‑1 (MQ=255)tGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTa < 1:1792878/68‑1 (MQ=255) gCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgtg > 1:3399541/1‑71 (MQ=255) cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgtg < 1:2331074/70‑1 (MQ=255) cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtgtg > 1:3524678/1‑70 (MQ=255) cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtg > 1:2288057/1‑68 (MQ=255) cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAtg > 1:2190091/1‑68 (MQ=255) cGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAt > 1:943696/1‑67 (MQ=255) aTTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCAt > 1:163782/1‑71 (MQ=255) ttACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCa > 1:2083247/1‑71 (MQ=255) actaGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCAc < 1:1885904/67‑1 (MQ=255) cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACt < 1:293862/70‑1 (MQ=255) gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGAt < 1:364089/71‑1 (MQ=255) gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGAt < 1:2291124/71‑1 (MQ=255) gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGAt < 1:2449836/71‑1 (MQ=255) gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGAt < 1:2115226/71‑1 (MQ=255) gaagaaGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg < 1:2415818/64‑1 (MQ=255) gaagaaGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg < 1:1989575/64‑1 (MQ=255) gaagaaGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGccc > 1:3641783/1‑71 (MQ=255) gaagaaGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGccc > 1:1734274/1‑71 (MQ=255) aagaagATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCg < 1:986271/71‑1 (MQ=255) aagaagATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCg < 1:3582159/71‑1 (MQ=255) agaagaTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCg > 1:3266908/1‑70 (MQ=255) aagaTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTgg < 1:1780459/61‑1 (MQ=255) agaTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGa > 1:933712/1‑61 (MQ=255) gacggCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACa > 1:2196522/4‑70 (MQ=255) gaTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACaa < 1:781931/71‑1 (MQ=255) tGGCATCGTCCTGAACGACCAGGATAAATGCCGTGGCTGGCGTAtgt > 1:2608690/1‑47 (MQ=255) cggCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACaaa > 1:155852/2‑70 (MQ=255) ggCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACaaaa > 1:1867789/1‑70 (MQ=255) gCATTGTGCTGATTGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACaaaaaa < 1:71648/71‑1 (MQ=255) gCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACaaaaaa < 1:635731/71‑1 (MQ=255) gCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACaaaaaa < 1:3279261/71‑1 (MQ=255) cATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTAt > 1:3091753/1‑42 (MQ=255) cATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACaaaaaaa < 1:177164/71‑1 (MQ=255) tCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATc > 1:3691146/1‑71 (MQ=255) gTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCt < 1:2782806/70‑1 (MQ=255) tCCTGATCGACCAGGATAAATGCCGTGGCTGGCTTATGTGCATCACTGGATGCCCGTACAAAAAAATCTAc < 1:925842/71‑1 (MQ=255) cTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACt < 1:219070/70‑1 (MQ=255) cTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACt < 1:100622/70‑1 (MQ=255) tGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACtt > 1:2159777/1‑70 (MQ=255) gATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGccc > 1:256887/1‑51 (MQ=255) gACCAGGATAAATGCCGTGGGTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTg < 1:2925096/71‑1 (MQ=255) tcAGGATAAATGCCGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGa < 1:3108239/70‑1 (MQ=255) gATAAATGCCGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAaaag > 1:3892619/1‑71 (MQ=255) gATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAagag > 1:3515872/1‑71 (MQ=255) tAAATGCCGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCg < 1:271886/71‑1 (MQ=255) tAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCg < 1:3585669/71‑1 (MQ=255) aTGCCGTGGCTGGCGTATGTGCATCACTGGATGccc > 1:2835075/1‑36 (MQ=255) aTGCCGTGGCTGGCGTATGTGCATCACTGGATGccc > 1:1655903/1‑36 (MQ=255) tGCCGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCaa < 1:838168/71‑1 (MQ=255) ccGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGt < 1:2965818/71‑1 (MQ=255) ccGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGt < 1:604585/71‑1 (MQ=255) gTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCt < 1:239832/71‑1 (MQ=255) tggctggcGTATGTGCATCTCTGGATGCCCGTACAAAAAAATCTACTTCAACTGGaa < 1:3299187/57‑1 (MQ=255) tggctggcGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTg < 1:2128510/71‑1 (MQ=255) ggctggcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTcag > 1:1195984/1‑68 (MQ=255) ggctggcGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCt > 1:2222021/1‑69 (MQ=255) gctggcGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTgag < 1:203394/71‑1 (MQ=255) tggcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAgaaa > 1:1152407/1‑68 (MQ=255) gcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaaat > 1:1719179/1‑69 (MQ=255) gcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaa < 1:3241472/69‑1 (MQ=255) gcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaa < 1:1843920/69‑1 (MQ=255) gcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaa < 1:3246154/69‑1 (MQ=255) gcGTATGTTCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAg > 1:2915481/1‑70 (MQ=255) gcGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGt > 1:3747510/1‑71 (MQ=255) cGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTgag < 1:934189/66‑1 (MQ=255) tATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCa < 1:1710470/71‑1 (MQ=255) tATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCa < 1:1503231/71‑1 (MQ=255) aTGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt > 1:1890133/1‑71 (MQ=255) aTGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt > 1:309516/1‑71 (MQ=255) | TGCGATTTACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAT > minE/835879‑836016 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |