| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 835,963 | C→T | 27.6% | C227C (TGC→TGT) | 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,963 | 0 | C | T | 27.6% | 90.0 / 37.4 | 58 | C227C (TGC→TGT) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
| Reads supporting (aligned to +/- strand): ref base C (15/27); new base T (7/9); total (22/36) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.63e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCC > minE/835895‑836028 | gaagaaGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGccc > 1:1734274/1‑71 (MQ=255)gaagaaGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGccc > 1:3641783/1‑71 (MQ=255) aagaagATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCg < 1:3582159/71‑1 (MQ=255) aagaagATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCg < 1:986271/71‑1 (MQ=255) agaagaTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCg > 1:3266908/1‑70 (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:3279261/71‑1 (MQ=255) gCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACaaaaaa < 1:635731/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: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:604585/71‑1 (MQ=255) ccGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGt < 1:2965818/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:3246154/69‑1 (MQ=255) gcGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGaaaa < 1:1843920/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:1503231/71‑1 (MQ=255) tATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCa < 1:1710470/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) ggATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATcc > 1:2096601/1‑71 (MQ=255) gATGCCCGTACAAAAAAATCTACTTCAATTGGAAGAGCGGTAAGTc < 1:2849429/46‑1 (MQ=255) gATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTc < 1:583075/46‑1 (MQ=255) aTGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa < 1:2546401/51‑1 (MQ=255) tGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTga < 1:136854/47‑1 (MQ=255) | GAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCC > minE/835895‑836028 |
| 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: ‑ |