| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 837,722 | T→G | 12.8% | G65G (GGT→GGG) | narI → | nitrate reductase 1, gamma (cytochrome b(NR)) subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 837,722 | 0 | T | G | 12.8% | 78.6 / 9.3 | 39 | G65G (GGT→GGG) | narI | nitrate reductase 1, gamma (cytochrome b(NR)) subunit |
| Reads supporting (aligned to +/- strand): ref base T (22/12); new base G (2/3); total (24/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.54e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01 | |||||||||||
TCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAA > minE/837659‑837785 | tCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCActtct < 1:3077969/71‑1 (MQ=255)tCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCActtct > 1:2943459/1‑71 (MQ=255) aaaGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCgg > 1:2214863/1‑71 (MQ=255) aaGGGATGAACCTGGCGTCGAACCTGTTCCATATAGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCgg > 1:3945069/1‑70 (MQ=255) gggATGAACCTGGCGTCGCACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCgg > 1:338379/1‑68 (MQ=255) gATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGc < 1:2548172/71‑1 (MQ=255) aTGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt > 1:498578/1‑71 (MQ=255) aTGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCt > 1:174865/1‑71 (MQ=255) aCCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCa > 1:3327032/1‑52 (MQ=255) ccTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgc > 1:2300660/1‑71 (MQ=255) cTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgc < 1:3764346/70‑1 (MQ=255) ggCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgccg > 1:1505350/1‑70 (MQ=255) gCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGAcgccg < 1:2446114/69‑1 (MQ=255) tCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTg < 1:2020350/71‑1 (MQ=255) gAATCTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTgg < 1:999754/70‑1 (MQ=255) gAATCTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTGCGGCATGTTAACGCCGCACTg < 1:1437546/69‑1 (MQ=255) gAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGa > 1:1273085/1‑71 (MQ=255) aCCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATg > 1:1463209/1‑71 (MQ=255) tcTGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGATg > 1:2560582/2‑70 (MQ=255) tGTTCCATATCGGCATTTTGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTgg < 1:1092000/65‑1 (MQ=255) ttCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGa < 1:1064298/71‑1 (MQ=255) tCCATTTCGGGATTCTGGGGATTTTTGTCGGTCAGTTCTTCGGTATGCTGACGCCGCACTGGATGTATGa < 1:2305298/70‑1 (MQ=255) tCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCg > 1:3939225/1‑42 (MQ=255) tCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATg > 1:3013147/1‑69 (MQ=255) atCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGa < 1:2855573/46‑1 (MQ=255) tGGGGATTTTCTTCGGGCACTTGTTCGGCATGTTAACGCCGCACTGGATGTa > 1:723555/1‑52 (MQ=255) gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCt > 1:1730880/1‑62 (MQ=255) gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCt > 1:210283/1‑62 (MQ=255) gggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGa < 1:766058/71‑1 (MQ=255) ggATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGc > 1:941868/1‑55 (MQ=255) aTTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGt < 1:1054019/71‑1 (MQ=255) aTTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGt < 1:994830/71‑1 (MQ=255) tttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGt > 1:3873386/1‑70 (MQ=255) tttttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGt > 1:3052113/1‑70 (MQ=255) ttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGaaa > 1:2970484/1‑71 (MQ=255) ttGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGaa < 1:2830759/70‑1 (MQ=255) tGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCt > 1:3248328/1‑50 (MQ=255) tGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCt > 1:3079998/1‑50 (MQ=255) tCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGaaa > 1:2326049/1‑68 (MQ=255) tCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGaa < 1:403390/40‑1 (MQ=255) tCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGaa < 1:304793/40‑1 (MQ=255) | TCGCAAAGGGATGAACCTGGCGTCGAACCTGTTCCATATCGGGATTCTGGGGATTTTTGTCGGTCACTTCTTCGGTATGCTGACGCCGCACTGGATGTATGAAGCCTGGCTGCCGATTGAAGTGAAA > minE/837659‑837785 |
| 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: ‑ |