| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 836,008 | G→A | 8.6% | E242E (GAG→GAA) | 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 | 836,008 | 0 | G | A | 8.6% | 179.1 / 5.9 | 58 | E242E (GAG→GAA) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
| Reads supporting (aligned to +/- strand): ref base G (28/25); new base A (2/3); total (30/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.65e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01 | |||||||||||
GCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGT > minE/835938‑836077 | 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) cATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtcttct < 1:1351853/71‑1 (MQ=255) tCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGc > 1:1081727/1‑71 (MQ=255) cACTGGATGCCCGAACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt > 1:3821169/1‑63 (MQ=255) aCTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCt < 1:760497/70‑1 (MQ=255) aCTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTa > 1:1308476/1‑71 (MQ=255) aCTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTa > 1:1611484/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) gTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTgag < 1:1735364/43‑1 (MQ=255) tACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGa < 1:975383/71‑1 (MQ=255) aCAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa < 1:3135431/43‑1 (MQ=255) aCAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa < 1:1913945/43‑1 (MQ=255) cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTgaga > 1:992809/1‑41 (MQ=255) cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa > 1:913501/1‑42 (MQ=255) cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa > 1:4024130/1‑42 (MQ=255) cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa < 1:2943451/42‑1 (MQ=255) cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGaa > 1:2902138/1‑42 (MQ=255) cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg > 1:3513274/1‑71 (MQ=255) cAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAg > 1:206487/1‑71 (MQ=255) aaaaaaTCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGagtc < 1:1046301/70‑4 (MQ=255) aaaaaaTCTACTTCAACTGGAAAAGCGGCAAGTCAGAAAAATGCATCTTCTGTTACCCACGAATTGagtc < 1:3543412/70‑4 (MQ=255) aaaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAtct < 1:1759857/47‑1 (MQ=255) aaaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt > 1:591480/1‑45 (MQ=255) aaaaaTCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCAt > 1:1905121/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:3961544/1‑69 (MQ=255) cTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCgt > 1:1707050/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) ggTAAGGCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAAc < 1:3216790/71‑1 (MQ=255) aaGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCg > 1:3771073/1‑42 (MQ=255) aaGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCtg > 1:1992359/1‑71 (MQ=255) tCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTgtgt < 1:4030856/71‑1 (MQ=255) agaAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtcggt < 1:2091282/71‑1 (MQ=255) agaAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTgtcgg < 1:588547/70‑1 (MQ=255) | GCTGGCGTATGTGCATCACTGGATGCCCGTACAAAAAAATCTACTTCAACTGGAAGAGCGGTAAGTCTGAGAAGTGCATCTTCTGCTATCCGCGTATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGT > minE/835938‑836077 |
| 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: ‑ |