| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 835,993 | G→A | 10.0% | K237K (AAG→AAA) | 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,993 | 0 | G | A | 10.0% | 208.8 / 11.0 | 70 | K237K (AAG→AAA) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
| Reads supporting (aligned to +/- strand): ref base G (31/32); new base A (3/4); total (34/36) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01 | |||||||||||
GATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCT > minE/835925‑836060 | gATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAagag > 1:3268113/1‑71 (MQ=255) tAAATGCCGTGGCTGGCGTTTGTGCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCg > 1:3606226/1‑71 (MQ=255) tAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTAAAAAAAAATCTACTTCAACTGGAAGAGc < 1:3777339/70‑1 (MQ=255) aaTGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCgg > 1:3173533/1‑70 (MQ=255) tggctggcGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTg < 1:2854133/71‑1 (MQ=255) gctggcGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTgag < 1:1387852/71‑1 (MQ=255) tATGTGCATCACTGTATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCa < 1:1657722/71‑1 (MQ=255) tATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCa < 1:1182694/71‑1 (MQ=255) aTGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCa > 1:1844103/1‑70 (MQ=255) tgtgCATAAGCGGTTGTCCGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCAtc > 1:1746518/1‑71 (MQ=255) tgCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtctt < 1:1243368/71‑1 (MQ=255) cATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGt < 1:2695539/53‑1 (MQ=255) cATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtcttct < 1:1351853/71‑1 (MQ=255) tCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAg < 1:250896/50‑1 (MQ=255) tCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGc > 1:1081727/1‑71 (MQ=255) cACTGGATGCCCGAACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt > 1:3821169/1‑63 (MQ=255) aCTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTc > 1:854190/1‑50 (MQ=255) aCTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCt < 1:760497/70‑1 (MQ=255) aCTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTa > 1:1611484/1‑71 (MQ=255) aCTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTa > 1:1308476/1‑71 (MQ=255) gATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATcc < 1:3891101/70‑1 (MQ=255) gATGCCAGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa > 1:540667/1‑52 (MQ=255) tGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCcgcg < 1:2440273/71‑1 (MQ=255) cGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTgaga < 1:1620997/45‑1 (MQ=255) cGTACAAAAAAATCTACTTCAACTGGAAAAGCGGCAAGt < 1:207541/39‑1 (MQ=38) gTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTgag < 1:1735364/43‑1 (MQ=255) tACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTg < 1:3202829/40‑1 (MQ=255) tACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGa < 1:975383/71‑1 (MQ=255) aCAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTg < 1:2922962/39‑1 (MQ=255) aCAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa < 1:3135431/43‑1 (MQ=255) aCAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa < 1:1913945/43‑1 (MQ=255) cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGt > 1:893888/1‑35 (MQ=255) cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTgaga > 1:992809/1‑41 (MQ=255) cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTga > 1:1926320/1‑39 (MQ=255) cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTga > 1:3792626/1‑39 (MQ=255) cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTga < 1:3459567/39‑1 (MQ=255) cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa > 1:913501/1‑42 (MQ=255) cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa > 1:4024130/1‑42 (MQ=255) cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa > 1:2902138/1‑42 (MQ=255) cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa < 1:2943451/42‑1 (MQ=255) cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg > 1:206487/1‑71 (MQ=255) cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg > 1:3513274/1‑71 (MQ=255) aaaaaaTCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGagtc < 1:3543412/70‑4 (MQ=255) aaaaaaTCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGagtc < 1:1046301/70‑4 (MQ=255) aaaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtct < 1:1759857/47‑1 (MQ=255) aaaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt > 1:1905121/1‑45 (MQ=255) aaaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt > 1:591480/1‑45 (MQ=255) aaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCggg < 1:2924591/71‑1 (MQ=255) aaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa > 1:3983243/1‑36 (MQ=255) aTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa > 1:2057001/1‑35 (MQ=255) aTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGaa < 1:1196987/63‑1 (MQ=255) aTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCa < 1:1158317/71‑1 (MQ=255) aTCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGa > 1:3249947/1‑62 (MQ=255) aTCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGAGTCCGGTCa < 1:3802249/71‑1 (MQ=255) tACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTg > 1:3051319/1‑58 (MQ=255) cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAg < 1:3433341/67‑1 (MQ=255) cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAg < 1:1827491/67‑1 (MQ=255) cTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGa > 1:3220856/1‑71 (MQ=255) cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtctt > 1:3751946/1‑36 (MQ=255) cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCg > 1:2901892/1‑58 (MQ=255) cAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg > 1:1894805/1‑71 (MQ=255) aaCTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg > 1:922271/1‑55 (MQ=255) aaCTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGAcc < 1:3493792/69‑1 (MQ=255) aaCTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg < 1:1790614/70‑1 (MQ=255) aCTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgtg < 1:3979604/71‑1 (MQ=255) cTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTa < 1:2609985/39‑1 (MQ=255) cTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgt > 1:1707050/1‑69 (MQ=255) cTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgt > 1:3961544/1‑69 (MQ=255) ggAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCcg > 1:1260889/1‑41 (MQ=255) gAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGc < 1:2191634/70‑1 (MQ=255) aaGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCt > 1:3657457/1‑70 (MQ=255) agagCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCg < 1:1445157/63‑1 (MQ=255) gagCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCt > 1:754678/1‑68 (MQ=255) | GATAAATGCCGTGGCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCT > minE/835925‑836060 |
| 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: ‑ |