| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 835,351 | A→G | 14.3% | S23S (TCA→TCG) | 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,351 | 0 | A | G | 14.3% | 196.9 / 23.6 | 77 | S23S (TCA→TCG) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
| Reads supporting (aligned to +/- strand): ref base A (36/30); new base G (3/8); total (39/38) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.14e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01 | |||||||||||
AAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAAC > minE/835281‑835419 | aaaTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCa < 1:3013980/71‑1 (MQ=255) aaTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCTGCTGCCACACCTGTTCa < 1:1130553/70‑1 (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) gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa < 1:1209399/71‑1 (MQ=255) gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa < 1:584646/71‑1 (MQ=255) gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa < 1:3949171/71‑1 (MQ=255) gTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa < 1:2609513/71‑1 (MQ=255) ttCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaa < 1:2759346/70‑1 (MQ=255) tCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGt > 1:959615/1‑59 (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:3414327/71‑1 (MQ=255) gTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTg < 1:142995/71‑1 (MQ=255) tCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTgg < 1:1611548/71‑1 (MQ=255) cGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTgg > 1:435241/1‑70 (MQ=255) ggCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGAc < 1:3251052/71‑1 (MQ=255) gTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGcc < 1:377194/70‑1 (MQ=255) gTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCg > 1:829858/1‑71 (MQ=255) gCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTaaa < 1:2487693/52‑1 (MQ=255) cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg > 1:696405/1‑70 (MQ=255) cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg > 1:701179/1‑70 (MQ=255) cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg < 1:2688763/70‑1 (MQ=255) tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCa < 1:2111351/38‑1 (MQ=255) tGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCa < 1:3554844/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:776056/69‑1 (MQ=255) ctcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGc < 1:1376149/69‑1 (MQ=255) tcGATAAGTGCATCGGCTGCCACACCTGTTCAGtt < 1:2048618/35‑1 (MQ=255) tcGATAAGTGCATCGGCTGCCACACCTGTTCAGtt > 1:739812/1‑35 (MQ=255) tcGATAAGTGCATCGGCTGCCACACCTGTTCAGtt > 1:3667292/1‑35 (MQ=255) tcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTa > 1:535111/1‑36 (MQ=255) tcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTaaaaa > 1:1745641/1‑46 (MQ=255) gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg < 1:687478/69‑1 (MQ=255) gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg < 1:1161257/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) aTCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGAcc > 1:654250/1‑45 (MQ=255) tCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata > 1:2318982/1‑64 (MQ=255) ggCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGtt > 1:3936105/1‑71 (MQ=255) ggCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGtt > 1:3887534/1‑71 (MQ=255) cTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTca < 1:796317/71‑1 (MQ=255) gCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg < 1:771015/45‑1 (MQ=255) ccATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac < 1:3895725/70‑1 (MQ=255) ccATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac < 1:815090/70‑1 (MQ=255) ccACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGaa < 1:1636282/55‑1 (MQ=255) acacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTgg < 1:3749559/51‑1 (MQ=255) acacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata > 1:247609/1‑55 (MQ=255) acacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata > 1:3281131/1‑55 (MQ=255) acacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCaacaac > 1:2683398/1‑71 (MQ=255) tacCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACg > 1:244932/2‑71 (MQ=255) cacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACg > 1:2469016/1‑71 (MQ=255) ccTGTTCGGTGACCTGTAAAAACGTCTGTACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGtc > 1:2059072/1‑70 (MQ=255) ccTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata > 1:1167491/1‑52 (MQ=255) ccTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg < 1:3461557/71‑1 (MQ=255) cTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGaa < 1:3328891/49‑1 (MQ=255) tGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGt < 1:2760596/58‑1 (MQ=255) tGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaacaa > 1:2640198/1‑65 (MQ=255) gTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaacaa > 1:3949634/1‑64 (MQ=255) ttCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATg < 1:940882/43‑1 (MQ=16) ttCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGtt < 1:3319294/57‑1 (MQ=255) ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg < 1:512008/35‑1 (MQ=255) ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata > 1:2868468/1‑48 (MQ=255) ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTc > 1:2777811/1‑58 (MQ=255) ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTc > 1:197828/1‑58 (MQ=255) ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGaa < 1:1224112/70‑1 (MQ=255) tCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCaacaac > 1:2321343/1‑63 (MQ=255) cAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg > 1:2559532/1‑65 (MQ=255) cAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAAc > 1:1769383/1‑70 (MQ=255) | AAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAAC > minE/835281‑835419 |
| 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: ‑ |