| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 835,393 | A→G | 12.5% | E37E (GAA→GAG) | 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,393 | 0 | A | G | 12.5% | 123.5 / 10.1 | 48 | E37E (GAA→GAG) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
| Reads supporting (aligned to +/- strand): ref base A (28/14); new base G (3/3); total (31/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.51e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01 | |||||||||||
GTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAA > minE/835327‑835455 | gTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATg < 1:372211/71‑1 (MQ=255) cATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGt > 1:947244/1‑71 (MQ=255) tCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata > 1:2318982/1‑64 (MQ=255) ggCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGtt > 1:3887534/1‑71 (MQ=255) ggCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGtt > 1:3936105/1‑71 (MQ=255) cTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTca < 1:796317/71‑1 (MQ=255) ccATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac < 1:3895725/70‑1 (MQ=255) ccATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac < 1:815090/70‑1 (MQ=255) ccACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGaa < 1:1636282/55‑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) ttCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGtt < 1:3319294/57‑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) gtgaCCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGAAACCa > 1:3987689/4‑70 (MQ=255) gTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata > 1:2612771/1‑44 (MQ=255) tACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAat > 1:3478885/1‑41 (MQ=255) tACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTgg > 1:236573/1‑49 (MQ=255) aCCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata > 1:53289/1‑41 (MQ=255) tGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAatat < 1:195268/39‑1 (MQ=255) tGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGc < 1:3215261/67‑1 (MQ=255) gTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaa > 1:2514953/1‑49 (MQ=255) tAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGc < 1:3336516/70‑1 (MQ=255) aaaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGc > 1:3875650/1‑69 (MQ=255) aaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCg > 1:2804611/1‑65 (MQ=255) aaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCg > 1:3292367/1‑65 (MQ=255) aaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAg < 1:2807869/71‑1 (MQ=255) tGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGaaa > 1:3995356/1‑49 (MQ=255) tGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGaa > 1:3169442/1‑48 (MQ=255) ggACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg < 1:1358102/44‑1 (MQ=255) ggACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg < 1:818788/44‑1 (MQ=255) gACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGcc > 1:2079717/1‑54 (MQ=255) gACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGc > 1:2084506/1‑53 (MQ=255) aCCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGa > 1:233150/1‑44 (MQ=255) aCCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCga > 1:3414097/1‑70 (MQ=255) cAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAggg < 1:3489602/60‑1 (MQ=255) aGCCGTAAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactga < 1:2591958/71‑1 (MQ=255) cGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTc > 1:1178231/1‑60 (MQ=255) gTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGaaa > 1:3321462/1‑39 (MQ=255) tGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGc < 1:1629368/49‑1 (MQ=255) aaGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCaa > 1:1061547/1‑40 (MQ=255) aaGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGaa < 1:4137788/71‑1 (MQ=255) ggCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGaaaa > 1:1071799/1‑71 (MQ=255) gAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTccc > 1:2147852/1‑50 (MQ=255) | GTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAA > minE/835327‑835455 |
| 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: ‑ |