| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 835,376 | A→G | 11.4% | S32G (AGC→GGC) ‡ | 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,376 | 0 | A | G | 11.4% | 185.0 / 15.1 | 70 | S32G (AGC→GGC) ‡ | narH | nitrate reductase 1, beta (Fe‑S) subunit |
| Reads supporting (aligned to +/- strand): ref base A (37/25); new base G (3/5); total (40/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.75e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01 | |||||||||||
GTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGA > minE/835310‑835446 | gTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGcc < 1:377194/70‑1 (MQ=255)gTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCg > 1:829858/1‑71 (MQ=255) cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg < 1:2688763/70‑1 (MQ=255) cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg > 1:696405/1‑70 (MQ=255) cTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTg > 1:701179/1‑70 (MQ=255) gAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAg > 1:2275248/1‑71 (MQ=255) acctcGACAAATGTATCGGCTGCCATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGc < 1:3752733/69‑1 (MQ=255) ctcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGc < 1:1376149/69‑1 (MQ=255) ctcGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGc < 1:776056/69‑1 (MQ=255) gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg < 1:1161257/69‑1 (MQ=255) gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg < 1:612248/69‑1 (MQ=255) gATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg < 1:687478/69‑1 (MQ=255) gTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATg < 1:372211/71‑1 (MQ=255) tGCATCGGCTGCCACACCTGTTCAGTTACCTGTATAAACGTCTGGACCAGcc > 1:2302690/1‑52 (MQ=255) tGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTg > 1:908627/1‑63 (MQ=255) cATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGt > 1:947244/1‑71 (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:815090/70‑1 (MQ=255) ccATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac < 1:3895725/70‑1 (MQ=255) ccACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGaa < 1:1636282/55‑1 (MQ=255) acacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTgg < 1:3749559/51‑1 (MQ=255) acacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata > 1:3281131/1‑55 (MQ=255) acacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAata > 1:247609/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:197828/1‑58 (MQ=255) ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTc > 1:2777811/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) gTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTgg < 1:1442120/40‑1 (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:3292367/1‑65 (MQ=255) aaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCg > 1:2804611/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:818788/44‑1 (MQ=255) ggACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg < 1:1358102/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) | GTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGA > minE/835310‑835446 |
| 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: ‑ |