| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 195,303 | A→G | 13.5% | N134S (AAT→AGT) | ispU → | undecaprenyl pyrophosphate synthase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 195,303 | 0 | A | G | 13.5% | 86.9 / 7.2 | 37 | N134S (AAT→AGT) | ispU | undecaprenyl pyrophosphate synthase |
| Reads supporting (aligned to +/- strand): ref base A (19/13); new base G (2/3); total (21/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.34e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01 | |||||||||||
GCGATACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAA > W3110S.gb/195234‑195373 | gCGATACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGaa < 1:4321901/70‑1 (MQ=255) cGATACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAAt < 1:3790712/70‑1 (MQ=255) tACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAGTACCg > 1:3618449/1‑71 (MQ=255) aCCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCgg > 1:3962204/1‑71 (MQ=255) ccAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGt < 1:4135436/71‑1 (MQ=255) aGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCt < 1:1410309/71‑1 (MQ=255) gTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCt > 1:2122158/1‑70 (MQ=255) ttAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCt > 1:4389315/1‑70 (MQ=255) aaCTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAGTACCGGTCTGACGCTGaa > 1:2239276/1‑71 (MQ=255) tCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCt < 1:4635504/59‑1 (MQ=255) cgcgTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAat > 1:4251686/1‑68 (MQ=255) gcgTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTg < 1:3261197/71‑1 (MQ=255) gcgTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTg > 1:3856400/1‑71 (MQ=255) cgTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACCGCCGGGAATACCGGTCTGACGCTGAATATTgc < 1:3334847/71‑1 (MQ=255) tGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGa > 1:462748/1‑54 (MQ=255) aaGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAAc < 1:110797/71‑1 (MQ=255) aGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGaa > 1:4102643/1‑69 (MQ=255) aGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGaa > 1:4097219/1‑69 (MQ=255) gAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCt < 1:2614733/47‑1 (MQ=255) gAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGa > 1:1613049/1‑49 (MQ=255) aaCGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAGTa < 1:2453870/39‑1 (MQ=255) aCGTATTCGTAAATCTGATGCGCTAACAGCCGGGAGTa < 1:2402838/38‑1 (MQ=38) cGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCt < 1:1313683/44‑1 (MQ=255) aTTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAata > 1:4529510/1‑52 (MQ=255) aTTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAata > 1:1828276/1‑52 (MQ=255) ttCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTg < 1:3083531/47‑1 (MQ=255) cGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGa > 1:1267360/1‑40 (MQ=255) cGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGa > 1:3083556/1‑40 (MQ=255) cGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGATTATTGCGGCGAACTACGGTGGa > 1:2959511/1‑69 (MQ=255) cGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACgg > 1:2276729/1‑65 (MQ=255) gTAAATCTGAAGCGCTAACAGCCGGGAATACCGGt < 1:2627209/35‑1 (MQ=255) aaaTCTGAAGCGCTAACAGCCGGGAATACCGGTCTGa > 1:275808/1‑37 (MQ=255) aTCTGAAGCGCTAACAGCCGGGAGTACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTgg < 1:2610355/70‑1 (MQ=255) tCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGa > 1:3941883/1‑71 (MQ=255) gcgcTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGt < 1:2989890/71‑1 (MQ=255) aaCAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAAc > 1:4147965/1‑43 (MQ=255) ggAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGa > 1:3048741/1‑44 (MQ=255) gAATACCGGTCTGACGCTGAATATTGCGGCGAACTAGGGTGGACGTTGGGATATAGTcc < 1:2874342/59‑1 (MQ=255) aaTACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAgg < 1:3177983/61‑1 (MQ=255) aTACCGGTCTGACGCTGAATATTGCGGCGAACTACgg > 1:913094/1‑37 (MQ=255) aTACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAggatc > 1:3642653/1‑68 (MQ=255) aTACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAggatc > 1:1276452/1‑68 (MQ=255) | GCGATACCAGTCGCTTTAACTCGCGTTTGCAAGAACGTATTCGTAAATCTGAAGCGCTAACAGCCGGGAATACCGGTCTGACGCTGAATATTGCGGCGAACTACGGTGGACGTTGGGATATAGTCCAGGGAGTCAGGCAA > W3110S.gb/195234‑195373 |
| 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: ‑ |