| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,653,313 | A→C | 100% | G82G (GGT→GGG) | upp ← | uracil phosphoribosyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,653,313 | 0 | A | C | 88.0% | 70.6 / 1.5 | 33 | G82G (GGT→GGG) | upp | uracil phosphoribosyltransferase |
| Reads supporting (aligned to +/- strand): ref base A (2/2); new base C (29/0); total (31/2) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.14e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
CCGACAACGCTGATGCGCGCGCTCGGAACGTTTTCCAGCACACCGTCCATCATACCAAGACCCGCACGCAGAATTGGCACAACGGTAATTTTCT > minE/1653260‑1653353 | ccGACAACGCTGATGCGCGCGCTCGGAACGTTTTCCAGCACACCGTCCATCATACCAAGACCcgca > 1:2684427/1‑66 (MQ=255)ccGACAACGCTGATGCGCGCGCTCGGAACGTTTTCCAGCACACCGTCCATCATACCAAGACCcgca > 1:1608107/1‑66 (MQ=255) cGGAATGTTTTCCAGCACACCGTCCATCATCCCaa > 1:3185897/1‑35 (MQ=25) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:2560815/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:2504186/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:2981544/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:2985572/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:3059056/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:3089254/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:3143434/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:3164667/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:3225324/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:3350363/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:583835/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:633191/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:750204/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:95237/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:1339134/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:2463398/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:2328690/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:2134791/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:190741/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:1869409/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:1738100/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:1693868/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:1654249/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:1640692/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:1406653/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:1402590/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:1398562/1‑35 (MQ=255) cGGAACGTTTTCCAGCACACCGTCCATCATCCCaa > 1:134477/1‑35 (MQ=255) aaCGTTTTCCAGCACACCGTCCATCATACCAAGACCCGCACGCAGAATTGGCACAACGGTAATTTTCt < 1:769305/68‑1 (MQ=255) aCGTTTTCCAGCACACCGTCCATCATACCAAGACCCGCACGCAGAATTGGCACAACGGTAATTTTCt < 1:2142640/67‑1 (MQ=255) | CCGACAACGCTGATGCGCGCGCTCGGAACGTTTTCCAGCACACCGTCCATCATACCAAGACCCGCACGCAGAATTGGCACAACGGTAATTTTCT > minE/1653260‑1653353 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |