| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 835,782 | C→A | 13.0% | A167E (GCG→GAG) | 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,782 | 0 | C | A | 13.0% | 185.5 / 9.2 | 61 | A167E (GCG→GAG) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
| Reads supporting (aligned to +/- strand): ref base C (25/28); new base A (6/2); total (31/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.55e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01 | |||||||||||
GACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCG > minE/835721‑835852 | 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) ggTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGa < 1:3582048/71‑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:2872116/55‑1 (MQ=255) gCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaa < 1:981039/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:2299302/1‑67 (MQ=255) aGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGcc > 1:1396351/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:1546133/1‑71 (MQ=255) aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCtg > 1:614351/1‑71 (MQ=255) aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCtg < 1:2015105/71‑1 (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:1714013/1‑71 (MQ=255) tCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGc > 1:3526955/1‑71 (MQ=255) ccAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGcc > 1:304097/1‑71 (MQ=255) ccAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGcc > 1:903867/1‑71 (MQ=255) cAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt > 1:608740/1‑71 (MQ=255) aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg < 1:489827/71‑1 (MQ=255) aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg < 1:940361/71‑1 (MQ=255) aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg < 1:945051/71‑1 (MQ=255) aGAAGGCGAGGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt > 1:3215519/1‑70 (MQ=255) gAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGa < 1:374529/71‑1 (MQ=255) aGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAc > 1:1848602/1‑71 (MQ=255) aGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAc < 1:135266/71‑1 (MQ=255) gCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAccc < 1:3876236/71‑1 (MQ=255) cGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCg > 1:3273268/1‑71 (MQ=255) | GACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCG > minE/835721‑835852 |
| 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: ‑ |