| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 835,956 | C→G | 24.2% | T225S (ACT→AGT) ‡ | 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,956 | 0 | C | G | 24.2% | 115.5 / 37.0 | 66 | T225S (ACT→AGT) ‡ | narH | nitrate reductase 1, beta (Fe‑S) subunit |
| Reads supporting (aligned to +/- strand): ref base C (19/31); new base G (7/9); total (26/40) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.72e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01 | |||||||||||
ACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTAT > minE/835887‑836026 | actaGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCAc < 1:1885904/67‑1 (MQ=255) cAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACt < 1:293862/70‑1 (MQ=255) gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGAt < 1:2449836/71‑1 (MQ=255) gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGAt < 1:2291124/71‑1 (MQ=255) gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGAt < 1:2115226/71‑1 (MQ=255) gCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGAt < 1:364089/71‑1 (MQ=255) gaagaaGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg < 1:1989575/64‑1 (MQ=255) gaagaaGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg < 1:2415818/64‑1 (MQ=255) gaagaaGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGccc > 1:1734274/1‑71 (MQ=255) gaagaaGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGccc > 1:3641783/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) 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) 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:100622/70‑1 (MQ=255) cTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACt < 1:219070/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:1843920/69‑1 (MQ=255) gcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaa < 1:3246154/69‑1 (MQ=255) gcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaa < 1:3241472/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) tgtgCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtc < 1:1561817/71‑1 (MQ=255) tgtgCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCAtc > 1:778851/1‑71 (MQ=255) gCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTgg > 1:934545/1‑41 (MQ=255) cATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGt > 1:207690/1‑53 (MQ=255) cATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtcttc < 1:1723363/70‑1 (MQ=255) tCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGc > 1:241643/1‑71 (MQ=255) aCTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTa < 1:501911/71‑1 (MQ=255) aCTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTa < 1:2942732/71‑1 (MQ=255) cTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTAt > 1:3086932/1‑71 (MQ=255) | ACAAGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTAT > minE/835887‑836026 |
| 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: ‑ |