| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 552,981 | C→T | 15.1% | A61A (GCG→GCA) | lpxH ← | UDP‑2,3‑diacylglucosamine pyrophosphatase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 552,981 | 0 | C | T | 15.1% | 75.7 / 6.0 | 33 | A61A (GCG→GCA) | lpxH | UDP‑2,3‑diacylglucosamine pyrophosphatase |
| Reads supporting (aligned to +/- strand): ref base C (16/12); new base T (3/2); total (19/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01 | |||||||||||
CAGAAAATCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGC > W3110S.gb/552915‑553045 | cAGAAAATCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCTCTTTGATcgccgcc < 1:3086948/70‑1 (MQ=255)cAGAAAATCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATcgccgcc < 1:1894532/70‑1 (MQ=255)cAGAAAATCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATcgccgac < 1:1324609/70‑3 (MQ=255) aTCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCt < 1:3260661/71‑1 (MQ=255) aTCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCt < 1:433104/71‑1 (MQ=255) gCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATg > 1:2577939/1‑69 (MQ=255) aatgaaATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTg > 1:3194369/1‑58 (MQ=255) atgaaATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGt < 1:3693105/70‑1 (MQ=255) tgaaATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATgg < 1:341538/62‑1 (MQ=255) aaaTAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTgg > 1:3661787/1‑65 (MQ=255) aaaTAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTggg < 1:1605388/71‑1 (MQ=255) tAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATgg < 1:302185/57‑1 (MQ=255) aaCAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGAtc < 1:2358767/70‑1 (MQ=255) aGGGAACGCCGGAATCGGACACCGCTTTGATCgccgccgcc > 1:3596802/1‑41 (MQ=255) ccGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAAt > 1:2104609/1‑71 (MQ=255) cGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGAtcgtc > 1:2124701/1‑61 (MQ=255) aTCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAg > 1:2476390/1‑42 (MQ=255) aTCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAg > 1:1517047/1‑42 (MQ=255) aTCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATg > 1:2250739/1‑71 (MQ=255) aTCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATg > 1:743152/1‑71 (MQ=255) tCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTggg < 1:3728935/50‑1 (MQ=255) ggACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGAt > 1:906814/1‑50 (MQ=255) ggACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGAt > 1:514977/1‑50 (MQ=255) ggACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGAt > 1:3370827/1‑50 (MQ=255) acacCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCa < 1:3441602/71‑1 (MQ=255) tacCGCTTTGATTGCTGCGGCCATCTGGTGATGAAGTGGGTTGGGATCGTCGTCGCCGATCCAGGCTTCa > 1:2131296/2‑70 (MQ=255) tacCGCTTTGATTGCTGCGGCCATCTGGTGATGAAGTGGGTTGGGATCGTCGTCGCCGATCCAGGCTTCa > 1:380728/2‑70 (MQ=255) tacCGCTTTGATTGCTGCGGCCATCTGGTGATGAAGTGGGTTGGGATCGTCGTCGCCGATCCAGGCTTCa > 1:700069/2‑70 (MQ=255) cacCGCTTTGAGCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCa < 1:1702986/70‑1 (MQ=255) ccGCTTTGATTGCTGCGGCCATCTGGTGATGAAGTGGGTTGGGATCGTCGTCGCCGATCCAGGCTTCAAAc < 1:1991868/71‑1 (MQ=255) ccGCTTTGATTGCTGCGGCCATCTGGTGATGAAGTGGGTTGGGATCGTCGTCGCCGATCCAGGCTTCAAAc < 1:141476/71‑1 (MQ=255) ccGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCaaa > 1:1629352/1‑70 (MQ=255) cGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACa < 1:235084/71‑1 (MQ=255) cGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACa < 1:1679859/71‑1 (MQ=255) gCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCaaa > 1:1808479/1‑68 (MQ=255) tGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCg > 1:16782/1‑71 (MQ=255) gATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCAATGCTACAAACAGATCGc > 1:3305373/1‑71 (MQ=255) | CAGAAAATCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGC > W3110S.gb/552915‑553045 |
| 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: ‑ |