| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 664,084 | C→T | 35.3% | G101D (GGC→GAC) | dacA ← | D‑alanyl‑D‑alanine carboxypeptidase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 664,084 | 0 | C | T | 35.3% | 32.2 / 30.3 | 34 | G101D (GGC→GAC) | dacA | D‑alanyl‑D‑alanine carboxypeptidase |
| Reads supporting (aligned to +/- strand): ref base C (6/16); new base T (6/6); total (12/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.65e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01 | |||||||||||
CCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCAT > W3110S.gb/664018‑664154 | cccGGTTTGAGGAACATCAGCGAAGAACCTTTAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGTCGa > 1:1731372/1‑70 (MQ=255)cccGGTTTGAGGAACATCAGCGAAGAACCTTTAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGTCGa > 1:496570/1‑70 (MQ=255) ggTTTGAGGAACATCAGCGAAGAACCTTTAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAg > 1:4073515/1‑70 (MQ=255) tCAGCGAAGAACCTTTAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGTCGa < 1:3108773/54‑1 (MQ=255) gaagaaCCTTTAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGTCGATa > 1:3851107/1‑51 (MQ=255) aagaaCCTTTAAACACCGGTTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTa < 1:3337854/71‑1 (MQ=255) aagaaCCTTTAAACACCGGTTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTa < 1:3333081/71‑1 (MQ=255) ccTTTAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGa > 1:205329/1‑43 (MQ=255) ccTTTAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAAttt < 1:2688727/71‑1 (MQ=255) ttAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGTCGATAGTGACTAAATCAGTTTCTTTAAATTTAcc > 1:1758579/1‑71 (MQ=255) acCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGtttt < 1:101253/71‑1 (MQ=255) acCCGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCttt < 1:3195864/71‑1 (MQ=255) ggTTACCGGTGGCCCATGCGTCGTTGTCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCAtt < 1:1610120/71‑1 (MQ=255) tACCGGTGGCCCATGCGTCGTTGTCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGcc > 1:538896/1‑71 (MQ=255) tACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCttt < 1:2597284/64‑1 (MQ=255) ggTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGc < 1:2173603/71‑1 (MQ=255) ccATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATaa < 1:3219245/71‑1 (MQ=255) ccATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATaa < 1:1054068/71‑1 (MQ=255) tGCGTCGTTGTCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGt < 1:391302/71‑1 (MQ=255) tGCGTCGTTGTCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGt < 1:3218279/71‑1 (MQ=255) tGCGTCGTTGTCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGt < 1:3563460/71‑1 (MQ=255) tGCGTCGTTGTCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGt < 1:3287158/71‑1 (MQ=255) cgtcgtTGCCGATAGTGACTAAATCAGTTTCTTTaaa > 1:3362276/1‑37 (MQ=255) cgtcgtTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCAtt > 1:1701209/1‑53 (MQ=255) cgtcgtTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCAtt > 1:419835/1‑53 (MQ=255) cgtcgtTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTa < 1:3262513/70‑1 (MQ=255) cgtcgtTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTa < 1:2272755/70‑1 (MQ=255) cgtcgtTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTa < 1:2189443/70‑1 (MQ=255) tcgtTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGt < 1:1046579/67‑1 (MQ=255) tcgtTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGt < 1:3158126/67‑1 (MQ=255) tcgtTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGt < 1:1563677/67‑1 (MQ=255) gtTGTCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCAt > 1:896737/1‑48 (MQ=255) gtTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGt < 1:3339886/65‑1 (MQ=255) gCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGa > 1:2363517/1‑56 (MQ=255) ccGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGtcat < 1:468424/71‑1 (MQ=255) | CCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAACACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCAT > W3110S.gb/664018‑664154 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |