| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 553,240 | T→C | E234G (GAG→GGG) | lpxH ← | UDP‑2,3‑diacylglucosamine diphosphatase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 553,240 | 0 | T | C | 96.5% | 94.8 / NA | 29 | E234G (GAG→GGG) | lpxH | UDP‑2,3‑diacylglucosamine diphosphatase |
| Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (9/19); total (10/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.45e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAA > NC_000913/553167‑553307 | gagCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCCGCTCAACGt > 2:110774/1‑80 (MQ=255) aaGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATcc < 2:110734/81‑1 (MQ=255) aaGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATcc < 1:110733/81‑1 (MQ=255) aGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATccg < 1:110735/81‑1 (MQ=255) aGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATccg > 1:110779/1‑81 (MQ=255) ggAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATccgc > 1:110785/1‑81 (MQ=255) aaaaCGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATccgccg < 1:110736/81‑1 (MQ=255) aaCGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATCCGCCGTg < 2:110737/81‑1 (MQ=255) ttGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATCCGCCGTGACttt < 1:110738/81‑1 (MQ=255) gTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATCCGCCGTGACTTTCAcc < 1:110739/81‑1 (MQ=255) aaTCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATCCGCCGTGACTTTCACCATTGAACCtt < 2:110740/81‑1 (MQ=255) aGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCg > 2:110786/1‑80 (MQ=255) aGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCg > 1:110783/1‑80 (MQ=255) aGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGt < 2:110741/81‑1 (MQ=255) gCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTa < 1:110742/81‑1 (MQ=255) cAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTAt > 2:110791/1‑81 (MQ=255) aaaGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATg > 1:110801/1‑81 (MQ=255) aGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATCCTCCGTGACTTTCACCATTGAACCTTCCGTATGc > 1:110782/1‑80 (MQ=255) gTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCa < 2:110745/81‑1 (MQ=255) gTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCa < 1:110743/81‑1 (MQ=255) gTTACGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCa < 2:110744/81‑1 (MQ=255) ttGCGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAg > 2:110787/1‑81 (MQ=255) cGGGTTTTTTAAAACGGAAAATGAATCAGCCCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCa > 2:110795/1‑81 (MQ=255) ttttttAAAACGGAAAATGAATCAGCCCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCa < 2:110746/81‑1 (MQ=255) tAAAACGGAAAATGAATCAGCCCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTAc < 1:110747/80‑1 (MQ=255) aaaaCGGAAAATGAATCAGCCCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTAcc < 2:110749/80‑1 (MQ=255) aaaaCGGAAAATGAATCAGCCCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTAcc < 1:110748/80‑1 (MQ=255) cGGAAAATGAATCAGCCCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCAcgcg < 1:110750/81‑1 (MQ=255) aaaaTGAATCAGCCCAACGTCATCCGCCGTGACTTTCAACATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGaaa < 1:110751/81‑1 (MQ=255) | GAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTATGCCAGGCACCCAGTACCACGCGAAA > NC_000913/553167‑553307 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |