| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 835,330 | C→T | 10.0% | C16C (TGC→TGT) | 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,330 | 0 | C | T | 10.0% | 240.1 / 13.4 | 80 | C16C (TGC→TGT) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
| Reads supporting (aligned to +/- strand): ref base C (31/41); new base T (2/6); total (33/47) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.59e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01 | |||||||||||
ACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGT > minE/835261‑835400 | aCCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa < 1:271334/71‑1 (MQ=255)aCCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa > 1:3291216/1‑71 (MQ=255)aCCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa > 1:3325566/1‑71 (MQ=255)aCCAGGTAAAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa > 1:3377333/1‑71 (MQ=255) cAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCCCAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATc < 1:412281/71‑1 (MQ=255) cAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATc < 1:3986898/71‑1 (MQ=255) aGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCg < 1:3659142/71‑1 (MQ=255) aGGTACAGGAG‑GCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATc < 1:3340659/70‑1 (MQ=255) ggTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCg > 1:3315964/1‑70 (MQ=255) aCAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCAt > 1:3697308/1‑65 (MQ=255) aCAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTg > 1:650800/1‑71 (MQ=255) aCAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTg > 1:1720348/1‑71 (MQ=255) cAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGc < 1:1878940/71‑1 (MQ=255) gagaGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa > 1:1536728/1‑60 (MQ=255) gagaGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa > 1:2403855/1‑60 (MQ=255) agagCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCcaca < 1:2470985/71‑1 (MQ=255) cG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTg < 1:456784/71‑1 (MQ=255) tAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGtt > 1:1315777/1‑71 (MQ=255) aaaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGtt > 1:1565500/2‑71 (MQ=255) aaaTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCa < 1:3013980/71‑1 (MQ=255) aaaTGAAAATTCGTTCACAAATCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc < 1:2036912/70‑1 (MQ=255) aaTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCTGCTGCCACACCTGTTCa < 1:1130553/70‑1 (MQ=255) aTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCAt > 1:1942934/1‑50 (MQ=255) tGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGtt < 1:2563524/71‑1 (MQ=255) tGAAAATCCGTACACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGt > 1:760565/1‑70 (MQ=255) aaaaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGTCTGCCACACCTGTTCAGTTAc > 1:1714717/1‑71 (MQ=255) aaaaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTAc > 1:525186/1‑71 (MQ=255) atccGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTc < 1:1220197/59‑1 (MQ=255) gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa < 1:1209399/71‑1 (MQ=255) gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa < 1:584646/71‑1 (MQ=255) gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa < 1:2609513/71‑1 (MQ=255) gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa < 1:3949171/71‑1 (MQ=255) ttCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa < 1:2759346/70‑1 (MQ=255) tCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCAt > 1:4119775/1‑38 (MQ=255) tCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGt > 1:959615/1‑59 (MQ=255) cacaAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc < 1:2091827/55‑1 (MQ=255) cacaAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc < 1:1258608/55‑1 (MQ=255) cacaAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAAc > 1:218643/1‑71 (MQ=255) acaAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAAc < 1:1864425/70‑1 (MQ=255) acaAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAAc < 1:1398341/70‑1 (MQ=255) caAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGt < 1:2737055/71‑1 (MQ=255) aaGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCt > 1:3289879/1‑60 (MQ=255) gTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg > 1:2017075/1‑71 (MQ=255) gTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg < 1:142995/71‑1 (MQ=255) gTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg < 1:3414327/71‑1 (MQ=255) tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc < 1:2673847/49‑1 (MQ=255) tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc < 1:555095/49‑1 (MQ=255) tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc < 1:3542954/49‑1 (MQ=255) tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTgg < 1:1611548/71‑1 (MQ=255) cGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTgg > 1:435241/1‑70 (MQ=255) ggCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGAc < 1:3251052/71‑1 (MQ=255) tGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc < 1:1698773/43‑1 (MQ=255) gTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGcc < 1:377194/70‑1 (MQ=255) gTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCg > 1:829858/1‑71 (MQ=255) tGCTGAATCTCGATAAGTGCATCGGCTGCCAAACCTg > 1:620591/1‑37 (MQ=255) gCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTaaa < 1:2487693/52‑1 (MQ=255) cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg > 1:701179/1‑70 (MQ=255) cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg > 1:696405/1‑70 (MQ=255) cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg < 1:2688763/70‑1 (MQ=255) tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc < 1:883458/37‑1 (MQ=255) tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc < 1:2638249/37‑1 (MQ=255) tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCa < 1:3554844/38‑1 (MQ=255) tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCa < 1:2111351/38‑1 (MQ=255) gAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAg > 1:2275248/1‑71 (MQ=255) acctcGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGc < 1:3752733/69‑1 (MQ=255) tctcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTAc < 1:1352567/39‑1 (MQ=255) ctcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGc < 1:1376149/69‑1 (MQ=255) ctcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGc < 1:776056/69‑1 (MQ=255) tcGATAAGTGCATCGGCTGCCACACCTGTTCAGtt < 1:2048618/35‑1 (MQ=255) tcGATAAGTGCATCGGCTGCCACACCTGTTCAGtt > 1:3667292/1‑35 (MQ=255) tcGATAAGTGCATCGGCTGCCACACCTGTTCAGtt > 1:739812/1‑35 (MQ=255) tcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTa > 1:535111/1‑36 (MQ=255) tcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaaaa > 1:1745641/1‑46 (MQ=255) gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg < 1:1161257/69‑1 (MQ=255) gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg < 1:687478/69‑1 (MQ=255) gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg < 1:612248/69‑1 (MQ=255) gTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATg < 1:372211/71‑1 (MQ=255) tGCATCGGCTGCCACACCTGTTCAGTTACCTGTATAAACGTCTGGACCAGcc > 1:2302690/1‑52 (MQ=255) tGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg > 1:3829178/1‑44 (MQ=255) tGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg > 1:908627/1‑63 (MQ=255) cATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGt > 1:947244/1‑71 (MQ=255) | ACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGT > minE/835261‑835400 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |