| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 835,774 | C→G | 15.1% | I164M (ATC→ATG) | 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,774 | 0 | C | G | 15.1% | 128.2 / 13.5 | 53 | I164M (ATC→ATG) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
| Reads supporting (aligned to +/- strand): ref base C (21/24); new base G (6/2); total (27/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.50e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01 | |||||||||||
AGATGACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCC > minE/835717‑835844 | aGATGACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATc > 1:3223780/1‑58 (MQ=255) gACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAg > 1:3178455/1‑60 (MQ=255) gACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAg > 1:2417978/1‑71 (MQ=255) gggCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCg > 1:2033963/1‑58 (MQ=255) gggCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCg > 1:159614/1‑58 (MQ=255) ggCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGc < 1:3930149/56‑1 (MQ=255) ggTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGa < 1:2630266/49‑1 (MQ=255) ggTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGa < 1:3582048/71‑1 (MQ=255) gTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATcc < 1:1136780/45‑1 (MQ=255) gTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa < 1:456674/71‑1 (MQ=255) gTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa < 1:2308116/71‑1 (MQ=255) aGTTTGACAAACTGGCGAAAGACAAGAGCTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa > 1:2126997/1‑70 (MQ=255) aGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa > 1:2025126/1‑70 (MQ=255) aGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa > 1:1919247/1‑68 (MQ=255) tttGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAacac > 1:2167966/1‑71 (MQ=255) tGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGt < 1:375901/59‑1 (MQ=255) cAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGta < 1:3410205/48‑1 (MQ=255) cAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCa < 1:1209796/71‑1 (MQ=255) aaaCTGGCGAAAGACAAGAACTTCGCCAACATCCAGAAGGCGATGtat > 1:3596208/1‑48 (MQ=255) ggCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCatgatg > 1:27783/1‑70 (MQ=255) gCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaa < 1:981039/55‑1 (MQ=255) gCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaa < 1:2872116/55‑1 (MQ=255) gCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTa < 1:1794791/71‑1 (MQ=255) cGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAt < 1:2822654/71‑1 (MQ=255) aaaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCa < 1:2760474/61‑1 (MQ=255) aaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa < 1:724718/50‑1 (MQ=255) aaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTg > 1:99782/1‑71 (MQ=255) aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGc > 1:1473338/1‑71 (MQ=255) cAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCcgc < 1:3433131/71‑1 (MQ=255) aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGcc > 1:1396351/1‑67 (MQ=255) aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGcc > 1:2299302/1‑67 (MQ=255) aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc > 1:3482164/1‑71 (MQ=255) aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc > 1:1582637/1‑71 (MQ=255) gAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGt < 1:1399271/38‑1 (MQ=255) gAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTg < 1:3467917/64‑1 (MQ=255) caaCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGcc > 1:2392205/2‑71 (MQ=255) caaCTTCGAGGCCATGCAAAAGGAGATGTACGCGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGcc > 1:2376532/2‑71 (MQ=255) aaCTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCt < 1:1312012/71‑1 (MQ=255) aaCTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCt < 1:2457881/71‑1 (MQ=255) aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGACAACACCTTCATGATGTACCTGCCGCGCCtg < 1:3120492/71‑1 (MQ=255) aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCtg < 1:2015105/71‑1 (MQ=255) aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCtg > 1:1546133/1‑71 (MQ=255) aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCtg > 1:614351/1‑71 (MQ=255) ttCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCcgc < 1:1348976/64‑1 (MQ=255) cGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGc > 1:960591/1‑70 (MQ=255) gACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCg < 1:3804234/70‑1 (MQ=255) acaacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGc < 1:3075465/68‑1 (MQ=255) ggccaTGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGaa > 1:403939/4‑70 (MQ=255) caacaTCCAGAAGGCGATGTATAGCCAGTTCGaaa > 1:1676077/1‑35 (MQ=255) caacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAac < 1:3166509/71‑1 (MQ=255) caTGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCatga > 1:1306037/1‑44 (MQ=12) tCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc < 1:1159264/56‑1 (MQ=255) tCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGc > 1:2269706/1‑71 (MQ=255) tCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGc > 1:3526955/1‑71 (MQ=255) tCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGc > 1:1714013/1‑71 (MQ=255) ccAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGcc > 1:903867/1‑71 (MQ=255) ccAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGcc > 1:304097/1‑71 (MQ=255) | AGATGACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCC > minE/835717‑835844 |
| 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: ‑ |