| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 835,798 | C→T | 13.9% | F172F (TTC→TTT) | 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,798 | 0 | C | T | 13.9% | 222.1 / 22.7 | 86 | F172F (TTC→TTT) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
| Reads supporting (aligned to +/- strand): ref base C (35/39); new base T (8/4); total (43/43) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.51e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01 | |||||||||||
GGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCT > minE/835730‑835868 | ggTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGa < 1:3582048/71‑1 (MQ=255) gTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa < 1:2308116/71‑1 (MQ=255) gTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa < 1:456674/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) cAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCa < 1:1209796/71‑1 (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: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) gAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTg < 1:3467917/64‑1 (MQ=255) caaCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGcc > 1:2392205/2‑71 (MQ=255) caaCTTCGAGGCCATGCAAAAGGAGATGTACGCGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGcc > 1:2376532/2‑71 (MQ=255) aaCTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCt < 1:2457881/71‑1 (MQ=255) aaCTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCt < 1:1312012/71‑1 (MQ=255) aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGACAACACCTTCATGATGTACCTGCCGCGCCtg < 1:3120492/71‑1 (MQ=255) aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCtg > 1:1546133/1‑71 (MQ=255) aCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCtg < 1:2015105/71‑1 (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:1714013/1‑71 (MQ=255) tCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGc > 1:2269706/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:945051/71‑1 (MQ=255) aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg < 1:940361/71‑1 (MQ=255) aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg < 1:489827/71‑1 (MQ=255) aGAAGGCGAGGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt > 1:3215519/1‑70 (MQ=255) gAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGa < 1:374529/71‑1 (MQ=255) aGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAc < 1:135266/71‑1 (MQ=255) aGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAc > 1:1848602/1‑71 (MQ=255) gCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAccc < 1:3876236/71‑1 (MQ=255) cGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCg > 1:3273268/1‑71 (MQ=255) gATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAaca > 1:3508944/1‑57 (MQ=255) tGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGc > 1:937686/1‑59 (MQ=255) tGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGc > 1:2309667/1‑59 (MQ=255) tGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCa > 1:2081803/1‑71 (MQ=255) tGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCa > 1:2271441/1‑71 (MQ=255) tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAac > 1:2661842/1‑52 (MQ=255) tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtg > 1:2041892/1‑71 (MQ=255) tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAt < 1:1593906/70‑1 (MQ=255) taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg > 1:3336679/1‑71 (MQ=255) taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg < 1:3181309/71‑1 (MQ=255) taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg < 1:1472040/71‑1 (MQ=255) taGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgt < 1:3621570/70‑1 (MQ=255) aGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtg < 1:3689200/70‑1 (MQ=255) gCCAGTTCGAAAACACCTTCATGATGTATTTGCcg > 1:1766361/1‑35 (MQ=255) gCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtgt < 1:3018493/70‑1 (MQ=255) gCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCAtgtgt < 1:3867308/70‑1 (MQ=255) ccAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc < 1:3935445/37‑1 (MQ=255) ccAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc < 1:2860980/37‑1 (MQ=255) ccAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTgg < 1:3039361/71‑1 (MQ=255) cAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGc < 1:3913462/71‑1 (MQ=255) aGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAAtccca > 1:1090451/1‑58 (MQ=255) aGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTTCGAACACTGCCTGAACCCGGCATGTGTGGc > 1:349795/1‑70 (MQ=255) aGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAcc > 1:3218393/1‑57 (MQ=255) aGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGc > 1:2971718/1‑70 (MQ=255) aGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGc > 1:289658/1‑70 (MQ=255) gTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACt < 1:1907587/47‑1 (MQ=38) gTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACt < 1:2882346/47‑1 (MQ=38) ttCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGAc < 1:3779713/71‑1 (MQ=255) ttCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGAc < 1:1829041/71‑1 (MQ=255) tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCg < 1:3163158/39‑1 (MQ=255) tCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg > 1:2069830/1‑57 (MQ=255) cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGAcc < 1:3935453/70‑1 (MQ=255) cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGAcc < 1:2028352/70‑1 (MQ=255) cGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCt > 1:2973485/1‑71 (MQ=255) | GGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGGCATGTGTGGCGACCT > minE/835730‑835868 |
| 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: ‑ |