| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 423,211 | A→G | 100% | R139R (CGT→CGC) | nagD ← | UMP phosphatase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 423,211 | 0 | A | G | 78.0% | 46.7 / 6.6 | 36 | R139R (CGT→CGC) | nagD | UMP phosphatase |
| Reads supporting (aligned to +/- strand): ref base A (7/1); new base G (9/19); total (16/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.21e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
CAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTC > minE/423151‑423276 | caACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCTGATTGCTGGCGATAAAACGTg < 1:901136/65‑1 (MQ=255) aGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGcgacga > 1:1362888/1‑68 (MQ=255) aGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGcgacga > 1:2372855/1‑68 (MQ=255) aGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGcgacga > 1:2151048/1‑68 (MQ=255) aGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGagacga > 1:1392934/1‑68 (MQ=255) ggTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGt > 1:1405713/1‑68 (MQ=255) ggTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGt > 1:815268/1‑68 (MQ=255) ggTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGt > 1:2179084/1‑68 (MQ=255) gtatCCGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat > 1:3115491/4‑62 (MQ=255) gtatCCGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat > 1:3129417/4‑62 (MQ=255) gtatCCGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat > 1:2644389/4‑62 (MQ=255) gtatCCGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat > 1:2437252/4‑62 (MQ=255) gtatCCGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat > 1:864941/4‑62 (MQ=255) gtatCCGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat > 1:2050943/4‑62 (MQ=255) gtatCCGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat > 1:1615966/4‑62 (MQ=255) gtatCCGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat > 1:916785/4‑62 (MQ=255) gtatCCGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcat > 1:1313863/4‑62 (MQ=255) cGGGTTGGTTGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat < 1:3157537/60‑1 (MQ=255) cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat < 1:926658/60‑1 (MQ=255) cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat < 1:699879/60‑1 (MQ=255) cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat < 1:568845/60‑1 (MQ=255) cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat < 1:3210528/60‑1 (MQ=255) cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat < 1:998982/60‑1 (MQ=255) cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat < 1:2651558/60‑1 (MQ=255) cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat < 1:28295/60‑1 (MQ=255) cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat < 1:2791785/60‑1 (MQ=255) cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat < 1:1084625/60‑1 (MQ=255) cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat < 1:2287086/60‑1 (MQ=255) cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat < 1:2069954/60‑1 (MQ=255) cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat < 1:1689854/60‑1 (MQ=255) cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat < 1:1592616/60‑1 (MQ=255) cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat < 1:1544053/60‑1 (MQ=255) cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat < 1:109947/60‑1 (MQ=255) cGGGTTGGTGGCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCCTTATGcatcat < 1:3096913/60‑1 (MQ=255) cGGGTTGGT‑GCGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat < 1:3297173/59‑1 (MQ=38) cGATAAAGCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat < 1:2539170/49‑1 (MQ=255) aaaaCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTc > 1:3158055/1‑69 (MQ=255) aaaaCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTc > 1:929371/1‑69 (MQ=255) | CAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTC > minE/423151‑423276 |
| 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: ‑ |