| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 1,557,824 | C→T | 22.8% | Q57Q (CAG→CAA) | bdm ← | biofilm‑dependent modulation protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 1,557,824 | 0 | C | T | 22.8% | 65.5 / 17.0 | 35 | Q57Q (CAG→CAA) | bdm | biofilm‑dependent modulation protein |
| Reads supporting (aligned to +/- strand): ref base C (10/17); new base T (4/4); total (14/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.85e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01 | |||||||||||
TAGGTTGCCGATCAAGCATAGCACCTTAAAGCGTAGGGTGCTGGCCACTGACCACATAATTGATCGTTTGCTGGTAGATATCACTGAGGATGTCGTTATCAGAAGCTTCAACCCATTTGGTCAGCTCCATGAGAATGT > W3110S.gb/1557754‑1557891 | tAGGTTGCCGATCAAGCATAGCACCTTAAAGCGTAGGGTGCTGGCCACTGACCACATAATTGATCGTTTGc > 1:2300717/1‑71 (MQ=255) gCCGATCAAGCATAGCACCTTAAAGCGTAGGGTGCTGGCCACTGACCACATAATTGATCGTTTGCTGGTa < 1:64540/70‑1 (MQ=255) ccGATCAAGCATAGCACCTTAAAGCGTAGGGTGCTGGCCACTGACCACATAATTGATCGTTTGCTGGTAGa < 1:2627483/71‑1 (MQ=255) aTCAAGCATAGCACCTTAAAGCGTAGGGTGCTGGCCACTGACCACATAATTGATCGTTTGCTGGTAGatat < 1:2810862/71‑1 (MQ=255) cAAGCATAGCACCTTAAAGCGTAGGGTGCTGGCCACTGACCACATAATTGATCGTTTGCTGGTAGATATCa > 1:887538/1‑71 (MQ=255) cAAGCATAGCACCTTAAAGCGTAGGGTGCTGGCCACTGACCACATAATTGATCGTTTGCTGGTAGATATCa > 1:23775/1‑71 (MQ=255) cATAGCACCTTAAAGCGTAGGGTGCTGGCCACTGACCACATAATTGATCGTTTGCTGGTAGATATCACt > 1:390661/1‑69 (MQ=255) aTAGCACCTTAAAGCGTAGGGTGCTGGCCACTGACCACATAATTGATCGTTTGCTGGTAGATATCACTGAg < 1:50090/71‑1 (MQ=255) aTAGCACCTTAAAGCGTAGGGTGCTGGCCACTGACCACATAATTGATCGTTTGCTGGTAGATATCACTGAg < 1:711814/71‑1 (MQ=255) gCACCTTAAAGCGTAGGGTGCTGGCCACTGACCACATAATTGATCGTTTGTTGGTAGATATCACTGAGGAt > 1:795897/1‑71 (MQ=255) gCACCTTAAAGCGTAGGGTGCTGGCCACTGACCACATAATTGATCGTTTGTTGGTAGATATCACTGAGGAt < 1:1130968/71‑1 (MQ=255) gCACCTTAAAGCGTAGGGTGCTGGCCACTGACCACATAATTGATCGTTTGCTGGTAGATATCACTGAGGAt < 1:536226/71‑1 (MQ=255) tAAAGCGTAGGGTGCTGGCCACTGACCACATAATTGATCGTTTGCTGGTAGATATCACTGAGGa > 1:1257657/1‑64 (MQ=255) aaaGCGTAGGGTGCTGGCCACTGACCACATAATTGATCGTTTGCTGGTAGATATCACTGAGGATGTCGTTa > 1:652591/1‑71 (MQ=255) aGCGTAGGGTGCTGGCCACTGACCACATAATTGATCGTTTGCTgg < 1:377884/45‑1 (MQ=255) gTGCTGGCCACTGACCACATAATTGATCGTTTGCTGGTAGATATCACTGAGGATGTCGTTATCAGAAGCtt < 1:1620527/71‑1 (MQ=255) cTGGCCACTGACCACATAATTGATCGTTTGTTGGTAGATATCACTGAGGATGTCGTTATCAGAAGCTTCaa > 1:1332451/1‑71 (MQ=255) gCCACTGACCACATAATTGATCGTTTGCTGGTAGATATCACTGAGGATGTCGTTATCAGAAGCTTCAAccc < 1:985115/71‑1 (MQ=255) cACTGACCACATAATTGATCGTTTGCTGGTAGatat < 1:1856464/36‑1 (MQ=255) aCTGACCACATAATTGATCGTTTGCTGGTAGatat < 1:593850/35‑1 (MQ=255) cacaTAATTGATCGTTTGTTGGTAGATATCACTGAg < 1:2809724/36‑1 (MQ=255) cacaTAATTGATCGTTTGTTGGTAGATATCACTGAg < 1:1964678/36‑1 (MQ=255) cacaTAATTGATCGTTTGTTGGTAGATATCACTGAg < 1:1467363/36‑1 (MQ=255) caTAATTGATCGTTTGTTGGTAGATATCACTGAGGa > 1:634407/1‑36 (MQ=255) caTAATTGATAGTTTGTTGGTAGATATCACTGAGGa > 1:550309/1‑36 (MQ=38) tAATTGATCGTTTGCTGGTAGATATCACTGAGGATg < 1:658296/36‑1 (MQ=255) tAATTGATCGTTTGCTGGTAGATATCACTGAGGATGTCGTTATCAGAAGCTTCAACCCATTTGGTCAGCTc > 1:965693/1‑71 (MQ=255) aTTGATCGTTTGCTGGTAGATATCACTGAGGATGTCGTTATCAGAAGCTTCAAccc > 1:636290/1‑56 (MQ=255) tGATCGTTTGCTGGTAGATATCACTGAGGATGTCGTTATCAGAAGCTTCAACCCATTTGGTCAGCTCCAt > 1:2277024/1‑70 (MQ=255) aTCGTTTGCTGGTAGATATCACTGAGGATGTCGTTATCAGAAGCTTCAACCCa < 1:2474742/53‑1 (MQ=255) aTCGTTTGCTGGTAGATATCACTGAGGATGTCGTTATCAGAAGCTTCAACCCa < 1:1593794/53‑1 (MQ=255) aTCGTTTGCTGGTAGATATCACTGAGGATGTCGTTATCAGAAGCTTCAACCCATTTGGTCAGCTCCATgag < 1:293422/71‑1 (MQ=255) aTCGTTTGCTGGTAGATATCACTGAGGATGTCGTTATCAGAAGCTTCAACCCATTTGGACAGCTCCATgag < 1:2394005/71‑1 (MQ=255) tCGTTTGCTGGTAGATATCACTGAGGATGTCGTTATCAGAAGCTTCAACCCATTTGGTCAGCTCCATgag > 1:1241722/1‑70 (MQ=255) cGTTTGCTGGTAGATATCACTGAGGATGTCGTTATCAGAAGCTT‑AACCCATTTGGTCAGCTCCATGAGAAt > 1:1713047/1‑71 (MQ=255) tGCTGGTAGATATCACTGAGGATGTCGTTATCAGAAGCTTCAACCCATTTGGTCAGCTCCATGAGAATGt < 1:505361/70‑1 (MQ=255) | TAGGTTGCCGATCAAGCATAGCACCTTAAAGCGTAGGGTGCTGGCCACTGACCACATAATTGATCGTTTGCTGGTAGATATCACTGAGGATGTCGTTATCAGAAGCTTCAACCCATTTGGTCAGCTCCATGAGAATGT > W3110S.gb/1557754‑1557891 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |