| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 1,489,255 | G→A | 12.1% | A103T (GCG→ACG) | ydcF → | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 1,489,255 | 0 | G | A | 12.1% | 90.4 / 3.6 | 33 | A103T (GCG→ACG) | ydcF | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base G (18/11); new base A (2/2); total (20/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01 | |||||||||||
TCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTCAACAAA > W3110S.gb/1489193‑1489322 | tCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGa < 1:1860960/67‑1 (MQ=255) gCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGc > 1:2329129/1‑71 (MQ=255) gCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCtcatc < 1:2456140/71‑1 (MQ=255) gCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCtcat > 1:2473706/1‑70 (MQ=255) cGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTc > 1:1652289/1‑71 (MQ=255) aCTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTgg < 1:2567930/71‑1 (MQ=255) caacaCTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGcaca > 1:1690895/1‑71 (MQ=255) caacaCTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGcaca < 1:2671654/71‑1 (MQ=255) caacaCTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGACGGATATCGCACATCAGTTCTGGcaca > 1:1095799/1‑71 (MQ=255) aacaCTATCCGCACCACTGGCAGAGCAGAAGCGAACATCCTGGCGGATATCGCTCATCAGTTCTGGCACAt > 1:2704403/1‑71 (MQ=255) tATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCg > 1:1558501/1‑70 (MQ=255) tCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCa > 1:2346358/1‑70 (MQ=255) tCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCa > 1:2470282/1‑70 (MQ=255) tCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCa > 1:2473502/1‑70 (MQ=255) tCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCAt > 1:2565878/1‑71 (MQ=255) cGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATg < 1:2771454/70‑1 (MQ=255) gcagagcagaAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTg < 1:2121350/70‑1 (MQ=255) gcagagcagaAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTg < 1:2447550/70‑1 (MQ=255) agcataAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTg < 1:601688/71‑1 (MQ=255) gcagaAGCGACCATCCTGACGGATATCGCTCAGCAGTTCTGGCACATTCCGCATGAAAAAATCTGGAGTg > 1:673160/1‑70 (MQ=255) cagaAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTg > 1:1833634/1‑69 (MQ=255) aaGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGAc < 1:82935/71‑1 (MQ=255) gCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAAt > 1:1123342/1‑56 (MQ=255) gCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAAt > 1:1802891/1‑56 (MQ=255) gCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGAc > 1:1783121/1‑69 (MQ=255) gCGACCATCCTGACGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCa < 1:751481/71‑1 (MQ=255) gCGACCATCCTGACGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCa < 1:1466059/71‑1 (MQ=255) ccATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTg < 1:2197598/60‑1 (MQ=255) tCCTGGCGGATATCGCTCATCAGTTCTTGCACATTCCGCATGAAAAAATCTGGATTGAAGAc < 1:1515297/62‑1 (MQ=255) cTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTcaacaa > 1:655746/1‑70 (MQ=255) cTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTCAACaaa > 1:2253260/1‑71 (MQ=255) cTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTCAACaaa > 1:59222/1‑71 (MQ=255) cTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTCAACaaa > 1:962958/1‑71 (MQ=255) | TCGCACAGCATCCGCACTACAACACTATCCGCACCACTGGCAGAGCAGAAGCGACCATCCTGGCGGATATCGCTCATCAGTTCTGGCACATTCCGCATGAAAAAATCTGGATTGAAGACCAGTCAACAAA > W3110S.gb/1489193‑1489322 |
| 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: ‑ |