Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 631,602 | A→G | 48.2% | *66* (TAA→TAG) | ybdD → | conserved hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 631,602 | 0 | A | G | 48.2% | ‑5.0 / 38.0 | 27 | *66* (TAA→TAG) | ybdD | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (8/6); new base G (6/7); total (14/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.06e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01 |
GAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGT > W3110S.gb/631537‑631671 | gAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGAt < 1:2865447/68‑1 (MQ=255) gAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTg < 1:2304092/71‑1 (MQ=255) agTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACg > 1:513230/1‑71 (MQ=255) aGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTaaa > 1:2277877/1‑71 (MQ=255) gCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGa < 1:3725709/70‑1 (MQ=255) gcgcTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTa < 1:2834630/44‑1 (MQ=255) cTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTg < 1:3218426/71‑1 (MQ=255) aCGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGa > 1:988417/1‑56 (MQ=255) cGGTGGAAAAGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGaa < 1:96560/70‑1 (MQ=255) gTGGAAAAGGCGGGGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGc > 1:1193605/1‑71 (MQ=255) gAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAAc < 1:3725025/48‑1 (MQ=255) gAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCt > 1:2207084/1‑69 (MQ=255) gAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTg > 1:4008560/1‑70 (MQ=255) aaaaGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGa > 1:3620717/1‑49 (MQ=255) aaGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACg > 1:1799114/1‑71 (MQ=255) aaGGCGGCGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACg > 1:4246167/1‑71 (MQ=255) aaGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGa < 1:2964106/47‑1 (MQ=255) ggcggcGCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCa < 1:4360994/71‑1 (MQ=255) gcggcgCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCaa > 1:5075362/1‑71 (MQ=255) cggcgCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAAc < 1:4211185/71‑1 (MQ=255) ggcgCGCGCTGCTGCTAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACg < 1:3407895/71‑1 (MQ=255) gcgcgcgcTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGt < 1:1677446/71‑1 (MQ=255) gctgctgctAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCt > 1:682914/1‑54 (MQ=255) tgctAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTg < 1:4004061/43‑1 (MQ=255) gctAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTa > 1:814548/1‑60 (MQ=255) gctAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTa > 1:820503/1‑60 (MQ=255) ctAGATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGaa > 1:4982321/1‑43 (MQ=255) aaaTTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGc > 1:327695/1‑69 (MQ=255) aaTTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGt < 1:1905221/70‑1 (MQ=255) | GAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGT > W3110S.gb/631537‑631671 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |