Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 2,855,731 | G→T | 17.8% | intergenic (‑269/‑18) | ygbA ← / → mutS | hypothetical protein/methyl‑directed mismatch repair protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 2,855,731 | 0 | G | T | 17.8% | 85.4 / 5.0 | 28 | intergenic (‑269/‑18) | ygbA/mutS | hypothetical protein/methyl‑directed mismatch repair protein |
Reads supporting (aligned to +/- strand): ref base G (10/13); new base T (2/3); total (12/16) | |||||||||||
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.85e-01 |
CATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCT > W3110S.gb/2855671‑2855801 | tctcgatctggtgcagaagaaaggcggaagttgcctctggcgagcataaTTTAATATCAGTAAACCGGACa > 1:791875‑M2/50‑71 (MQ=255) tgtgCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATa < 1:3135383/71‑1 (MQ=255) gtgCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATaa > 1:2319088/1‑71 (MQ=255) cgcCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCC > 1:1059870‑M1/1‑71 (MQ=255) ccTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGAc < 1:57209/62‑1 (MQ=255) cTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCAt < 1:2964988/71‑1 (MQ=255) gtgcagaagaaaggcggaagttgcctctggcgagcgtaaTTTAATATCAGTAAACCGGACATAAcc < 1:1471675‑M2/27‑1 (MQ=255) tATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCC < 1:2904546‑M1/66‑1 (MQ=255) agaagaaaggcggaagttgcctctggcgagcataaTTTAATATCAGTAAACCGGACATAAc < 1:957376‑M2/26‑1 (MQ=255) aTTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa < 1:1140891/71‑1 (MQ=255) tACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa > 1:2200929/1‑69 (MQ=255) tACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa > 1:1553681/1‑69 (MQ=255) aaCGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAAc < 1:3071333/43‑1 (MQ=255) gAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGaaa > 1:50747/1‑69 (MQ=255) aaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACAt < 1:953858/44‑1 (MQ=255) aTAAAAACCATCACACCCCATTTAATATCAGGGAAc < 1:1369009/36‑1 (MQ=255) aTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCAt < 1:1329413/51‑1 (MQ=255) acacCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcccatacgc < 1:1161956/71‑1 (MQ=255) ccATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcccatacgcccat > 1:2280436/1‑70 (MQ=255) ccATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcccatacgcccat > 1:238663/1‑70 (MQ=255) cATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCatga < 1:315788/71‑1 (MQ=255) tAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgggcagggggcaagta < 1:1014154‑M1/70‑45 (MQ=255) tAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCatgatg > 1:499699/1‑69 (MQ=255) tatCAGTAAACCGGACATAACCCatgaagaaaaatc > 1:159597‑M1/1‑23 (MQ=255) tatCAGGGAACCGGACATAACCCCATGAGTGCAATAGaa < 1:1000298/39‑1 (MQ=255) tCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCatgatg > 1:1649716/1‑64 (MQ=255) aGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTc < 1:356251/41‑1 (MQ=255) aGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGa > 1:1738932/1‑43 (MQ=255) ggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTAtct < 1:2713093/71‑1 (MQ=255) | CATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCT > W3110S.gb/2855671‑2855801 |
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: ‑ |