| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,474,863 | C→G | 13.4% | intergenic (+161/‑73) | nrdA → / → nrdB | ribonucleoside diphosphate reductase 1, alpha subunit/ribonucleoside diphosphate reductase 1, beta subunit, ferritin‑like |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,474,863 | 0 | C | G | 13.4% | 118.4 / 6.8 | 43 | intergenic (+161/‑73) | nrdA/nrdB | ribonucleoside diphosphate reductase 1, alpha subunit/ribonucleoside diphosphate reductase 1, beta subunit, ferritin‑like |
| Reads supporting (aligned to +/- strand): ref base C (15/21); new base G (2/3); total (18/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01 | |||||||||||
GTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAG > minE/1474797‑1474926 | gTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGc > 1:139728/1‑71 (MQ=255)gTGCCGGATGCAGCGTGAACGCATTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGc > 1:1755111/1‑71 (MQ=255) ggATGCAGCGTGAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCg > 1:2861943/1‑70 (MQ=255) gATGCAGCGTGAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa > 1:3642708/1‑71 (MQ=255) gATGCAGCGTGAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCa < 1:2643893/71‑1 (MQ=255) aTGCAGCGTGAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCAt < 1:1398706/71‑1 (MQ=255) aGCGTAAACGCCTGATCCGTCCAAC‑‑‑‑GGTTCG‑G‑‑ATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAg > 1:3457726/1‑71 (MQ=2) gCGTGAACGCCTTATCCGGCCTACGGATGGCGCG‑GGAA‑TTTGTAGGCCTGATAAGACGCGTTAGCGTCGCa < 1:2546502/71‑1 (MQ=0) cGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGA‑A‑TTTGTAGGCCTGATAAGACGCGTTAGCGTCGCAt > 1:3206446/1‑71 (MQ=0) taaaCGCCTGATCCGTCCAAC‑‑‑‑GGTTCG‑G‑‑ATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGcat < 1:3519389/69‑2 (MQ=2) aaCGCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAgg < 1:1702689/71‑1 (MQ=255) aaCGCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAg < 1:2833867/70‑1 (MQ=255) aaCGCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAg < 1:3225702/70‑1 (MQ=255) aaCGCCTGATCCGTCCAAC‑‑‑‑GGTTCG‑G‑‑ATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGcatct < 1:1790819/71‑4 (MQ=2) aaCGCCTGATCCGTCCAAC‑‑‑‑GGTTCG‑G‑‑ATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCATCAGGcatc < 1:3415668/70‑3 (MQ=2) aCGCCTGATCCGTCCAAC‑‑‑‑GGTTCG‑G‑‑ATTTTGTAGGCCTGATAAGACGCGGTAGCGTCGCAt > 1:2284845/1‑61 (MQ=2) gCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATg < 1:3217889/71‑1 (MQ=255) gCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATg < 1:2401244/71‑1 (MQ=255) tATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATgcggc > 1:3370307/1‑71 (MQ=255) tATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATgcggc > 1:1032332/1‑71 (MQ=255) gCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATg < 1:2240725/61‑1 (MQ=255) ccTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAAt > 1:237384/1‑71 (MQ=255) cTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATg < 1:2909558/71‑1 (MQ=255) tAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTaa < 1:1571474/66‑1 (MQ=255) aC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAg < 1:3692041/53‑1 (MQ=255) cTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTAt < 1:2901458/71‑1 (MQ=255) tCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATc > 1:788138/1‑71 (MQ=255) g‑g‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATcc < 1:669431/70‑1 (MQ=255) g‑g‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCg < 1:3873256/71‑1 (MQ=255) g‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCgg < 1:294594/71‑1 (MQ=255) g‑‑A‑TTTGTAGGCCTGATAAGACGAGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCgg < 1:604547/71‑1 (MQ=255) ttGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCAt > 1:475180/1‑71 (MQ=255) tAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTa < 1:901173/70‑1 (MQ=255) ggCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAAc > 1:816354/1‑71 (MQ=255) gCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACt < 1:910614/71‑1 (MQ=255) ccTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTc > 1:1699658/1‑71 (MQ=255) tAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCaa > 1:1346502/1‑70 (MQ=255) aGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAg < 1:1082570/71‑1 (MQ=255) aCGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCa < 1:545070/54‑1 (MQ=255) cgcgCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGc > 1:2492379/1‑42 (MQ=255) cgcgCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTa > 1:409811/1‑46 (MQ=255) cgcgCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTa > 1:2663659/1‑46 (MQ=255) cgcgCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACt < 1:2138498/60‑1 (MQ=255) cgcCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCaa < 1:3782427/63‑1 (MQ=255) gcCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAAt > 1:1750236/1‑37 (MQ=255) | GTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCG‑G‑‑A‑TTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCACAGGATGCGGCGTAAAATGCCTTATCCGGCATTAAACTCCCAACAG > minE/1474797‑1474926 |
| 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: ‑ |