| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | Exported | 833,993 | C→G | 43.6% | A817A (GCC→GCG) | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | Exported | 833,993 | 0 | C | G | 43.6% | 18.0 / 64.6 | 55 | A817A (GCC→GCG) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base C (13/18); new base G (8/16); total (21/34) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.84e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01 | |||||||||||
AAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGT > Exported/833923‑834058 | aaaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGcc > 1:159557/1‑71 (MQ=255) aaaaGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGcc < 1:2971644/70‑1 (MQ=255) aaaGACGATATGAACACCTCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTg < 1:527492/71‑1 (MQ=255) gacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg < 1:2757281/70‑1 (MQ=255) gacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg < 1:673802/70‑1 (MQ=255) gacATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg < 1:5605469/70‑1 (MQ=255) tGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGa < 1:2070046/69‑1 (MQ=255) gAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGaa > 1:4454508/1‑69 (MQ=255) aaTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAg > 1:1726496/1‑64 (MQ=255) aTACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc < 1:3035158/71‑1 (MQ=255) ctCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGaa < 1:2503919/57‑1 (MQ=255) ctCGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCAGCGTGGGaa < 1:2231905/57‑1 (MQ=255) cGGATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg > 1:5187583/1‑71 (MQ=255) ggATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTggg < 1:1846389/71‑1 (MQ=255) ggATATGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg > 1:3689816/1‑70 (MQ=255) tatGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGTGaaa < 1:1481207/71‑1 (MQ=255) tatGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGaaa < 1:2915590/71‑1 (MQ=255) tatGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGa < 1:3322576/69‑1 (MQ=255) tGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTggg < 1:4323523/48‑1 (MQ=37) tGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTggg < 1:2517544/48‑1 (MQ=37) tGCATCCGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTggg < 1:1551834/48‑1 (MQ=37) tGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGc < 1:2124322/52‑1 (MQ=255) gCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGa > 1:2439283/1‑60 (MQ=255) cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCt < 1:3401193/70‑1 (MQ=255) cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTa < 1:4377260/71‑1 (MQ=255) cATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTa < 1:771430/71‑1 (MQ=255) tCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc < 1:78578/55‑1 (MQ=255) cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACaaa < 1:1469837/71‑1 (MQ=255) cGTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACa > 1:4287315/1‑69 (MQ=255) cGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCt > 1:2225820/1‑39 (MQ=255) cGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACa > 1:3522696/1‑69 (MQ=255) gTTTATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAAt < 1:4112380/63‑1 (MQ=255) tttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGa < 1:2210751/46‑1 (MQ=255) tttATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGc < 1:5515961/51‑1 (MQ=255) ttATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAgg > 1:4407565/1‑69 (MQ=255) tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg > 1:5369519/1‑55 (MQ=255) tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTgg > 1:5665279/1‑55 (MQ=255) tATTCATCCGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAggta > 1:4662936/1‑68 (MQ=255) tATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCa > 1:5064259/1‑71 (MQ=255) tATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCa > 1:3288003/1‑71 (MQ=255) tccGCTTTCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATcgcca < 1:5645485/70‑3 (MQ=255) gCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCgaag > 1:4232274/1‑70 (MQ=255) ctttCTGCGGCGGTCGATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATCGCcaa > 1:2376203/4‑68 (MQ=255) gTCGATCCGGCGGGGGAATCACGCAGCGACTGGGAAATCTACaa < 1:2095785/44‑1 (MQ=255) gTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAgt < 1:593408/71‑1 (MQ=255) cGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGaaa > 1:3149122/1‑58 (MQ=255) gATCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATCGCCAAAGCATTTTCGCAAGtgtg < 1:4110679/71‑1 (MQ=255) aTCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACAAAGGTATCGCCAAAGCATTTTCGCAAGTGTGc < 1:1316784/71‑1 (MQ=255) aTCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGc < 1:3609045/43‑1 (MQ=255) aTCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAgtgt < 1:2696619/69‑1 (MQ=255) aTCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGc > 1:217339/1‑71 (MQ=255) tCCGGCGTGGGAATCACGCAGCGACTGGGAAATCTACaa < 1:2132342/39‑1 (MQ=16) tCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCg > 1:1893072/1‑71 (MQ=255) tCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCg > 1:1222129/1‑71 (MQ=255) ccGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGt > 1:3690035/1‑71 (MQ=255) cGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGt < 1:5064572/70‑1 (MQ=255) gCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGcc < 1:349508/39‑1 (MQ=255) | AAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCACCCGCTGTCTGCGGCGGTCGATCCGGCCTGGGAAGCGAAAAGCGACTGGGAAATCTACAAAGCCATCGCGAAGAAATTCTCCGAAGTGTGCGT > Exported/833923‑834058 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |