Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 835,315 | G→T | 10.9% | L11L (CTG→CTT) | 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,315 | 0 | G | T | 10.9% | 295.9 / 12.2 | 91 | L11L (CTG→CTT) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
Reads supporting (aligned to +/- strand): ref base G (35/46); new base T (5/5); total (40/51) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.44e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01 |
GATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAG > minE/835245‑835385 | gATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTg < 1:1145921/71‑1 (MQ=255) ggCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt > 1:1561820/1‑71 (MQ=255) ggCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt > 1:1741163/1‑71 (MQ=255) gCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtc < 1:2442591/71‑1 (MQ=255) cGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtct < 1:1533646/71‑1 (MQ=255) aaGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtc < 1:2869373/68‑1 (MQ=255) ggCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGAt < 1:3187346/71‑1 (MQ=255) ggCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGAt < 1:538969/71‑1 (MQ=255) aaTGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt < 1:3007371/62‑1 (MQ=255) aTGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGt < 1:3008300/71‑1 (MQ=255) gACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTg > 1:2091815/1‑70 (MQ=255) aCCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa > 1:3325566/1‑71 (MQ=255) aCCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa > 1:3291216/1‑71 (MQ=255) aCCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa < 1:271334/71‑1 (MQ=255) aCCAGGTAAAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa > 1:3377333/1‑71 (MQ=255) tcAGGTACAGGAG‑GCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTaa > 1:1715360/2‑56 (MQ=255) tcAGGTACAGGAG‑GCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTaa > 1:1778164/2‑56 (MQ=255) tcAGGTACAGGAG‑GCGAAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTaa > 1:2259667/2‑56 (MQ=255) cAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCCCAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATc < 1:412281/71‑1 (MQ=255) cAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCg < 1:1992091/60‑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) gTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTg < 1:1424238/50‑1 (MQ=255) gTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTg > 1:1295622/1‑64 (MQ=255) tACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtct > 1:2213868/1‑54 (MQ=255) aCAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCAt > 1:3697308/1‑65 (MQ=255) aCAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTg > 1:1720348/1‑71 (MQ=255) aCAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTg > 1:650800/1‑71 (MQ=255) cAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGaa > 1:3334287/1‑49 (MQ=255) cAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt > 1:3197154/1‑50 (MQ=255) cAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGc < 1:1878940/71‑1 (MQ=255) cAGGAGAGCG‑TAAAATGAAAATTCGATCACACGTCGTCATGGTGCTGAAtc > 1:2101013/1‑51 (MQ=255) aGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtctc < 1:2356822/52‑1 (MQ=255) aGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGAt < 1:2451042/55‑1 (MQ=255) gagaGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATaa > 1:56236/1‑55 (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) gagCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGa > 1:742718/1‑50 (MQ=255) gagCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGa > 1:53865/1‑50 (MQ=255) cG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTg < 1:456784/71‑1 (MQ=255) tAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGaa > 1:1742712/1‑39 (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) aTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATa < 1:2149539/37‑1 (MQ=255) ttCGTTCACAAGTCGGCATGGTGCTGAATCTCGAt < 1:4516/35‑1 (MQ=255) gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa < 1:3949171/71‑1 (MQ=255) gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa < 1:2609513/71‑1 (MQ=255) gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa < 1:1209399/71‑1 (MQ=255) gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa < 1:584646/71‑1 (MQ=255) ttCACAAGTCGGCATGGTGCTGAATCTCGATAAGt < 1:3827467/35‑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:1258608/55‑1 (MQ=255) cacaAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc < 1:2091827/55‑1 (MQ=255) cacaAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAAc > 1:218643/1‑71 (MQ=255) acaAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAAc < 1:1398341/70‑1 (MQ=255) acaAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAAc < 1:1864425/70‑1 (MQ=255) caAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGt < 1:2737055/71‑1 (MQ=255) aaGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCt > 1:3289879/1‑60 (MQ=255) gTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg < 1:3414327/71‑1 (MQ=255) gTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg < 1:142995/71‑1 (MQ=255) gTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg > 1:2017075/1‑71 (MQ=255) tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc < 1:3542954/49‑1 (MQ=255) tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc < 1:555095/49‑1 (MQ=255) tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc < 1:2673847/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:2688763/70‑1 (MQ=255) cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg > 1:696405/1‑70 (MQ=255) cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg > 1:701179/1‑70 (MQ=255) tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc < 1:2638249/37‑1 (MQ=255) tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTc < 1:883458/37‑1 (MQ=255) tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCa < 1:2111351/38‑1 (MQ=255) tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCa < 1:3554844/38‑1 (MQ=255) gAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAg > 1:2275248/1‑71 (MQ=255) | GATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑TAAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAG > minE/835245‑835385 |
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: ‑ |