| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 835,789 | T→C | 12.5% | Y169Y (TAT→TAC) | 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,789 | 0 | T | C | 12.5% | 157.3 / 12.7 | 64 | Y169Y (TAT→TAC) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
| Reads supporting (aligned to +/- strand): ref base T (25/31); new base C (6/2); total (31/33) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.42e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01 | |||||||||||
GACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGG > minE/835721‑835853 | gACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAg > 1:4795866/1‑71 (MQ=255)gACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAg > 1:795753/1‑71 (MQ=255) aCCTGGGCGGTGTGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGc < 1:1143967/71‑1 (MQ=255) gTGAGTTTGACAAACTGGCGAAATACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa < 1:2148505/71‑1 (MQ=255) tGAGTTTGACAAACTGGCGAAAGACAATAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGa < 1:5368627/69‑1 (MQ=255) gAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGc < 1:2461712/60‑1 (MQ=255) gAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa > 1:1527763/1‑71 (MQ=255) gAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa > 1:3329932/1‑71 (MQ=255) tGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCt > 1:832129/1‑71 (MQ=255) tGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCt > 1:2327926/1‑71 (MQ=255) aCAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTc < 1:1505982/71‑1 (MQ=255) aaaCTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCat < 1:4135534/71‑1 (MQ=255) aaaCTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCat < 1:4764202/71‑1 (MQ=255) ggCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCatgatg > 1:4345582/1‑70 (MQ=255) gCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCatgatg < 1:194374/69‑1 (MQ=255) gCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTa < 1:3107848/71‑1 (MQ=255) aaaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa < 1:3098997/53‑1 (MQ=255) aaaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaaaa < 1:693211/53‑1 (MQ=255) aaaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAttt < 1:3531289/71‑1 (MQ=255) aaaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAttt < 1:945759/71‑1 (MQ=255) aaaGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAttt > 1:44376/1‑71 (MQ=255) aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGaa < 1:91621/49‑1 (MQ=255) aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGc > 1:4756403/1‑71 (MQ=255) aGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGc > 1:2123703/1‑71 (MQ=255) gAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTa > 1:2749590/1‑60 (MQ=255) gAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTa > 1:2371877/1‑60 (MQ=255) caaCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCgcgc > 1:594367/2‑70 (MQ=255) caaCTTCGAGGCCATGCAAAAGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGcc > 1:3923764/2‑71 (MQ=255) aCTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAt > 1:125370/1‑59 (MQ=255) aCTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAt > 1:1884637/1‑59 (MQ=255) cTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCcgc < 1:4191217/65‑1 (MQ=255) tCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTc > 1:1805003/1‑47 (MQ=255) tCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTAt > 1:2531262/1‑56 (MQ=255) tCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCtgtg > 1:773307/1‑70 (MQ=255) tCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGc > 1:2789212/1‑71 (MQ=255) cGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCgcgc < 1:3201256/64‑1 (MQ=255) cGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGc > 1:631081/1‑70 (MQ=255) acaacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGag < 1:686382/71‑2 (MQ=255) acaacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGaa < 1:1707990/71‑1 (MQ=255) acaacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGaa < 1:5528516/71‑1 (MQ=255) caacaTGCAGAAGGCGATGTATAGCCAGTTCGAAAacac > 1:2496/1‑39 (MQ=255) caacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAacac > 1:3657043/1‑39 (MQ=255) aacaTCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAaca < 1:4207407/71‑1 (MQ=255) caaaaGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCt > 1:2853531/4‑71 (MQ=255) caaaaGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCt > 1:432562/4‑71 (MQ=255) cAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGcc < 1:1291312/70‑1 (MQ=255) cAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt < 1:5149149/71‑1 (MQ=255) cAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt > 1:5462064/1‑71 (MQ=255) aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt > 1:4348567/1‑70 (MQ=255) aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg < 1:1663632/71‑1 (MQ=255) aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg < 1:1781324/71‑1 (MQ=255) aGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTg < 1:3186099/71‑1 (MQ=255) aaaGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCtca < 1:3666103/70‑3 (MQ=255) aGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCaat > 1:3935238/1‑70 (MQ=255) aGGAGATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCaat > 1:3079316/1‑70 (MQ=255) gagATGTACGGGCAGTTTGAAAACACCTTCATGATGTACCTGCCGCGCCTGTGCGAACACTGCCTCAAtcc < 1:5372620/69‑1 (MQ=255) gCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAccc < 1:560082/71‑1 (MQ=255) gCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAAcc < 1:691822/70‑1 (MQ=255) cGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCt > 1:3961053/1‑64 (MQ=255) gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCtg > 1:2120516/1‑49 (MQ=255) gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg < 1:5643594/71‑1 (MQ=255) gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg < 1:5176742/71‑1 (MQ=255) gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg < 1:2563743/71‑1 (MQ=255) gATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCgg < 1:1687524/71‑1 (MQ=255) tataGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCtg < 1:698975/45‑1 (MQ=255) | GACCTGGGCGGTGAGTTTGACAAACTGGCGAAAGACAAGAACTTCGACAACATCCAGAAGGCGATGTATAGCCAGTTCGAAAACACCTTCATGATGTATTTGCCGCGCCTGTGCGAACACTGCCTGAACCCGG > minE/835721‑835853 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |