| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 835,957 | T→C | 47.0% | T225T (ACT→ACC) | 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,957 | 0 | T | C | 47.0% | 4.2 / 72.0 | 83 | T225T (ACT→ACC) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
| Reads supporting (aligned to +/- strand): ref base T (28/16); new base C (20/19); total (48/35) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.75e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01 | |||||||||||
AGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGC > minE/835890‑836023 | aGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg < 1:2387346/69‑1 (MQ=255)aGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg < 1:632221/69‑1 (MQ=255)aGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg < 1:614342/69‑1 (MQ=255)aGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg < 1:2303718/69‑1 (MQ=255)aGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg < 1:1006738/69‑1 (MQ=255)aGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg < 1:2230782/69‑1 (MQ=255) aagaagATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTg < 1:3034073/63‑1 (MQ=255) gaTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTa > 1:1260928/1‑68 (MQ=255) gaTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTa > 1:1268773/1‑68 (MQ=255) gaTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTa > 1:591325/1‑68 (MQ=255) gaTGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTa > 1:1859093/1‑68 (MQ=255) gTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAAt > 1:3072649/1‑68 (MQ=255) gTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATc > 1:41865/1‑69 (MQ=255) gTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATc > 1:1678988/1‑69 (MQ=255) cTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTgg > 1:539026/1‑47 (MQ=255) cTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAAt > 1:3239597/1‑65 (MQ=255) cTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATc > 1:3210532/1‑66 (MQ=255) cTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATc > 1:2205119/1‑66 (MQ=255) cTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATc > 1:2399534/1‑66 (MQ=255) cTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATc > 1:2635207/1‑66 (MQ=255) cTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATc > 1:2474735/1‑66 (MQ=255) cTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATc > 1:445598/1‑66 (MQ=255) cTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATc > 1:996096/1‑66 (MQ=255) cTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATc > 1:1745873/1‑66 (MQ=255) cTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATc > 1:1661322/1‑66 (MQ=255) cTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTAc > 1:2621855/1‑69 (MQ=255) gATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACt > 1:2118148/1‑68 (MQ=255) gATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACt > 1:233324/1‑68 (MQ=255) gATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACt > 1:2415787/1‑68 (MQ=255) cGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATTAGCTGTTGTCCGTACAAAAAAATCTACTTCaa > 1:2503906/1‑69 (MQ=255) cGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATTAGCGGTTGTCCGTACAAAAAAATCTACTTCaa > 1:3377566/1‑69 (MQ=255) cGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATTAGCGGTTGTCCGTACAAAAAAATCTACTTCaa > 1:2863268/1‑69 (MQ=255) cGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATTAGCGGTTGTCCGTACAAAAAAATCTACTTCaa > 1:268636/1‑69 (MQ=255) cGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATTAGCGGTTGTCCGTACAAAAAAATCTACTTCaa > 1:2914853/1‑69 (MQ=255) cGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATTAGCGGTTGTCCGTACAAAAAAATCTACTTCaa > 1:296109/1‑69 (MQ=255) cGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATTAGCGGTTGTCCGTACAAAAAAATCTACTTCaa > 1:3110189/1‑69 (MQ=255) cGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATTAGCGGTTGTCCGTACAAAAAAATCTACTTCaa > 1:3195430/1‑69 (MQ=255) cGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATTAGCGGTTGTCCGTACAAAAAAATCTACTTCaa > 1:3217106/1‑69 (MQ=255) cGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATTAGCGGTTGTCCGTACAAAAAAATCTACTTCaa > 1:1615436/1‑69 (MQ=255) cGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATTAGCGGTTGTCCGTACAAAAAAATCTACTTCaa > 1:1655519/1‑69 (MQ=255) cGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATTAGCGGTTGTCCGTACAAAAAAATCTACTTCaa > 1:2347448/1‑69 (MQ=255) cGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATTAGCGGTTGTCCGTACAAAAAAATCTACTTCaa > 1:3422146/1‑69 (MQ=255) cGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATTAGCGGTTGTCCGTACAAAAAAATCTACTTCaa > 1:2257433/1‑69 (MQ=255) cGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATTAGCGGTTGTCCGTACAAAAAAATCTACTTCaa > 1:485639/1‑69 (MQ=255) cGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATTAGCGGTTGTCCGTACAAAAAAATCTACTTCaa > 1:2203306/1‑69 (MQ=255) cGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATTAGCGGTTGTCCGTACAAAAAAATCTACTTCaa > 1:1346399/1‑69 (MQ=255) cGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATTAGCGGTTGTCCGTACAAAAAAATCTACTTCaa > 1:623630/1‑69 (MQ=255) cGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATTAGCGGTTGTCCGTACAAAAAAATCTACTTCaa > 1:1916758/1‑69 (MQ=255) cGATCAGGATAAATGCCGTGGCTGGCGTTTGTGCATTAGCGGTTGTCCGTACAAAAAAATCTACTTCa > 1:2042869/1‑68 (MQ=255) aGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGa < 1:1110335/69‑1 (MQ=255) aGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGa < 1:2637181/69‑1 (MQ=255) tggcGTATGTGCATCACCGGATGCCCGTACAAAAAAATCTACTTCAACTgg < 1:1399420/51‑1 (MQ=255) tggcGTATGTGCATCACCGGATGCCCGTACAAAAAAATCTACTTCAACTgg < 1:1268685/51‑1 (MQ=255) tggcGTATGTGCATCACCGGATGCCCGTACAAAAAAATCTACTTCAACTgg < 1:206763/51‑1 (MQ=255) tggcGTATGTGCATCACCGGATGCCCGTACAAAAAAATCTACTTCAACTgg < 1:534812/51‑1 (MQ=255) tggcGTATGTGCATCACCGGATGCCCGTACAAAAAAATCTACTTCAACTgg < 1:509049/51‑1 (MQ=255) tggcGTATGTGCATCACCGGATGCCCGTACAAAAAAATCTACTTCAACTgg < 1:1902704/51‑1 (MQ=255) tggcGTATGTGCATCACCGGATGCCCGTACAAAAAAATCTACTTCAACTgg < 1:1352304/51‑1 (MQ=255) tggcGTATGTGCATCACCGGATGCCCGTACAAAAAAATCTACTTCAACTgg < 1:433195/51‑1 (MQ=255) tggcGTATGTGCATCACCGGATGCCCGTACAAAAAAATCTACTTCAACTgg < 1:2891493/51‑1 (MQ=255) tggcGTATGTGCATCACCGGATGCCCGTACAAAAAAATCTACTTCAACTgg < 1:2343997/51‑1 (MQ=255) tggcGTATGTGCATCACCGGATGCCCGTACAAAAAAATCTACTTCAACTgg < 1:3262679/51‑1 (MQ=255) tggcGTATGTGCATCACCGGATGCCCGTACAAAAAAATCTACTTCAACTgg < 1:1577879/51‑1 (MQ=255) tggcGTATGTGCATCACCGGATGCCCGTACAAAAAAATCTACTTCAACTgg < 1:1413475/51‑1 (MQ=255) tggcGTATGTGCATCACCGGATGCCCGTACAAAAAAATCTACTTCAACTgg < 1:1224244/51‑1 (MQ=255) tggcGTATGTGCATCACCGGATGCCCGTACAAAAAAATCTACTTCAACTgg < 1:731635/51‑1 (MQ=255) tggcGTATGTGCATCACCGGATGCCCGTACAAAAAAATCTACTTCAACTgg < 1:1453819/51‑1 (MQ=255) tggcGTATGTGCATCACCGGATGCCCGTACAAAAAAATCTACTTCAACTgg < 1:1489207/51‑1 (MQ=255) tggcGTATGTGCATCACCGGATGCCCGTACAAAAAAATCTACTTCAACTgg < 1:2914771/51‑1 (MQ=255) tggcGTATGTGCATCAACGGATGCCCGTACAAAAAAATCTACTTCAACTgg < 1:1378344/51‑1 (MQ=255) tgtgCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTgg > 1:2591720/1‑44 (MQ=255) tgtgCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTgg > 1:1462942/1‑44 (MQ=255) tgtgCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTgg < 1:2968074/44‑1 (MQ=255) tgtgCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTgg > 1:217486/1‑44 (MQ=255) tgCATCACTGGATGCCCGTACAAAAAAATCTACTTc < 1:2296006/36‑1 (MQ=255) tgCATCACTGGATGCCCGTACAAAAAAATCTACTTc < 1:848541/36‑1 (MQ=255) tgCATCACTGGATGCCCGTACAAAAAAATCTACTTc < 1:621931/36‑1 (MQ=255) tgCATCACTGGATGCCCGTACAAAAAAATCTACTTc < 1:1332468/36‑1 (MQ=255) tgCATCACTGGATGCCCGTACAAAAAAATCTACTTc < 1:3208100/36‑1 (MQ=255) tgCATCACTGGATGCCCGTACAAAAAAATCTACTTc < 1:2504409/36‑1 (MQ=255) aCTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGc > 1:569719/1‑69 (MQ=255) aCTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGc > 1:3209292/1‑69 (MQ=255) aCTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGc > 1:2017355/1‑69 (MQ=255) | AGCGTGAAGAAGATGGCATCGTCCTGATCGACCAGGATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGC > minE/835890‑836023 |
| 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: ‑ |