| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | CP009273 | 4,185,493 | C→G | 24.9% | E179Q (GAG→CAG) | thiC ← | phosphomethylpyrimidine synthase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP009273 | 4,185,493 | 0 | C | G | 24.9% | 28.0 / 4.3 | 16 | E179Q (GAG→CAG) | thiC | phosphomethylpyrimidine synthase |
| Reads supporting (aligned to +/- strand): ref base C (7/5); new base G (2/2); total (9/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GGCTCCGATTCCGGATGATTAATGTTGGCCGGGATAATCGCACGTCCGGCAGCAACTTCATCACGGACAAATTCCGCAGTGATATTTTCCGGCAGATGTGCGCCAAAGCTCATTCCCGGATGCTGGTGGCGTAAAACCTCGCTACGGATGCGCTCGCGGCCCATATTCTCGCGGATGGCGATGAATTCCATTTCCGGCGTGATGATGCCCTGGCGGGCGTAGTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGC > CP009273/4185354‑4185638 | ggCTCCGATTCCGGATGATTAATGTTGGCCGGGATAATCGCACGACCGGCAGCAACTTCATCACGGACAAATTCCGCAGTGATATTTTCCGGCAGATGTGCGCCAAAGCTCATTCCCGGATGCTGGTGGCGTAAAACCTCGCTACGGAt < 2:226222/149‑1 (MQ=255) tatggtTAATGTTGGCCGGGATAATCGCACGTCCGGCAGCAACTTCAACACGGACAAATTCCGCAGTGCTATTTTTCGGCAGATGTGCGCCAAAGCTCATTCCCGGATGCTGGTGGCGTAAACCCTCGCTCCGCATACGCTAGCggggc < 1:207571/148‑4 (MQ=255) gatgatTAATGTTGGCCGGGATAATCGCACGTCCGGCAGCAACTTCATCACGGACAAATTCCGCAGTGATATTTTCCGGCAGATGTGCGCCAAAGCTCATTCCCGGATGCTGGTGGCGTAAAACCTCGCTACGGATGCGCTCGCGGccc < 2:76636/149‑1 (MQ=255) tCGCACGTCCGGCAGCAACTTCATCACGGACAAATTCCGCAGTGATATTTTCCGGCAGATGTGCGCCAAAGCTCATTCCCGGATGCTGGTGGCGTAAAACCTGGCTACGGATGCGCTCGCGGCCCATATTCTCGCGGATGGCGATGAAt > 1:281116/1‑149 (MQ=255) cACGTCCGGCAGCAACTTCATCACGGACAAATTCCGCAGTGATATCCTCCCGCACATCTGAGCAACCCCTCCTCCCCGGCTGTCAGCGGTCACAAACCAGCATACGGACGCGCTCGCGGGTCATATTCGCGCGGTTTGCCATGAAtctc > 1:27492/1‑146 (MQ=255) gattgcatcaCGCACAGATTCCGCAGTGATACTTCCCGGCAGATGTGCGCCAAAGCTCATTCCCGGATGCTGGTGGCGTAAAACCTCGCTACGGAGGCGCTCGCGGCTCATATTCCCGCGGATGGCGATGAATTCCATTTCCGGCGtga < 2:32118/144‑1 (MQ=255) catcaCGGACAAATTCCGCAGTGATATTTTCCGGCAGATGTACGCCAAAGCTCATTCCCGGATGCTGGTGGCGTAAAACCTCGCTACGGATGCGCTCGCGGCCCATATTCTCGCGGATGGCGATGAATTCCATTTCCGGCGGGatgatg > 2:103153/1‑149 (MQ=255) ccGCAGGGATATTTTCCGGCAGATGTGCGCCAAAGCTCATTCCCGGATGCTGGTGGCGTAAAACCTCGCTACGGATGCGCTCGCGGCCCATATTCTCGCGGATGGCGATGAATTCCATTTCCGGCGTGATGATGCCCTGGCGGGCGTAg > 1:295419/1‑149 (MQ=255) atTTTCCGGCAGATGTGCGCCAAAGCTCATTCCCGGATGCTGGTGGCGTAAAACCTCGCTACGGATGCGCTCGCGGCCCATATTCTCGCGGATGGCGATGAATTCCATTTCCGGCGTGATGATGCCCTGGCGGGAGTAGTGCGGTgggg > 1:237950/1‑149 (MQ=255) ttttCCGGCAGAGGTCCGCCAAAGCTCATTCCCGGATGCTGGTCGCGTACAACCTCGCTACGGATGCGCTCGCGGCCCATATTCTCGCGGATGGCGATGAATTCCATTTCCGGCGTGATGATGCCCTGGCGGGCGTAGTGCAGTTGGGt < 1:157630/149‑1 (MQ=255) aTTCCCGGATGCTGGTGGCGTAAAACCTCGCTACGGATGCGCTCGCGGCCCATATTCTCGCGGATGGCGATGAATTCCATTTCCGGCGTGATGATGCCCTGGCGGGCGTAGTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGttt > 1:303353/1‑149 (MQ=255) gCTGGTGGCGTAAAACCTCGCTACGGATGCGCTCGCGGCCCATATTCTCGCGGATGGCGATGAATTCCATTTCCGGCGTGATGATGCCCTGGCGGGCGTAGTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGGGTTCGt > 1:158936/1‑149 (MQ=255) gtggCGTAAAACCTCGCTACGGATGCGCTCGCGGCCCATATTCTCGCGGATGGCGATGAATTCCATTTCCGGCGTGATGATGCCCCGTCGCGCGGAGTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTAcgc > 2:73041/1‑149 (MQ=255) gCGTAAAGCCTGGCTACGGATGCGCTCGCTGCCCATATTCTCGCGGATGGCGATGAATTCCATTTCCGGCGTGATGATGCCCTGGCGGGCGTAGTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTAcgccgc < 2:2408/149‑1 (MQ=255) gCGTAAAACCTCGCTACGGATGCGCTCGCGGCCCATATTCTCGCGGATGGCGATGAATTCCATTTCCGGCGTGATGATGCCCTGGCGGGCGTAGTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTAcggagc > 2:66996/1‑145 (MQ=255) ccTGGCTACGGATGCGCTCGCGGCCCATATTCTCGCGGATGGCGATGAATTCCATTTCCGGCGTGATGATGCCCTGGCGGGCGTAGTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGc < 2:281116/149‑1 (MQ=255) | GGCTCCGATTCCGGATGATTAATGTTGGCCGGGATAATCGCACGTCCGGCAGCAACTTCATCACGGACAAATTCCGCAGTGATATTTTCCGGCAGATGTGCGCCAAAGCTCATTCCCGGATGCTGGTGGCGTAAAACCTCGCTACGGATGCGCTCGCGGCCCATATTCTCGCGGATGGCGATGAATTCCATTTCCGGCGTGATGATGCCCTGGCGGGCGTAGTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGC > CP009273/4185354‑4185638 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |