Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 2,855,732 | G→A | 23.3% | intergenic (‑270/‑17) | 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,732 | 0 | G | A | 23.3% | 55.6 / 14.2 | 30 | intergenic (‑270/‑17) | ygbA/mutS | hypothetical protein/methyl‑directed mismatch repair protein |
Reads supporting (aligned to +/- strand): ref base G (12/11); new base A (4/3); total (16/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01 |
CTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATG > W3110S.gb/2855662‑2855790 | cTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAggg < 1:140837/71‑1 (MQ=255) tATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCgg < 1:2781456/71‑1 (MQ=255) tATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCg > 1:3890201/1‑70 (MQ=255) aTGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACa > 1:1950529/1‑70 (MQ=255) aTGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACa > 1:4139530/1‑70 (MQ=255) aTGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACAt < 1:2017205/71‑1 (MQ=255) tgCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATaa < 1:488559/70‑1 (MQ=255) ttATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCAt > 1:2141176/1‑70 (MQ=255) tgtgATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGa < 1:4083592/56‑1 (MQ=255) ttACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATaa > 1:59231/1‑56 (MQ=255) ttACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa > 1:998240/1‑70 (MQ=255) cGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGaaaa > 1:2132749/1‑71 (MQ=255) gAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAAt > 1:1297512/1‑71 (MQ=255) aTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCAt < 1:732696/51‑1 (MQ=255) aTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTc > 1:314423/1‑70 (MQ=255) aaaaCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGac < 1:1433666/70‑1 (MQ=255) aaCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATg < 1:823573/47‑1 (MQ=255) ccATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCCacgccca < 1:1666740/71‑1 (MQ=255) tCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcccatac > 1:4036708/1‑71 (MQ=255) cacaCCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcccatacg > 1:1058097/1‑71 (MQ=255) cacCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcccatacgcc < 1:45715/71‑1 (MQ=255) cccATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcccatacgcccat > 1:1134987/1‑71 (MQ=255) taattTAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgggcagggggc > 1:779844‑M1/3‑31 (MQ=255) aatTTAATATCAGTAAGCCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgggcag > 1:842587‑M1/2‑30 (MQ=255) aatTTAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgggcag > 1:2185561‑M1/2‑30 (MQ=255) aTTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgccca < 1:3265051/59‑1 (MQ=255) aTTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgc < 1:2156569/56‑1 (MQ=255) tttAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgggcagggggcaagt < 1:2498356‑M1/71‑44 (MQ=255) tttAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaa > 1:2976159‑M1/1‑28 (MQ=255) tttAATATCAGTAAACCGGACATAACCCatgaagaaaaat < 1:3009564‑M1/40‑13 (MQ=255) tatCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgg < 1:2862287‑M1/53‑31 (MQ=255) | CTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATG > W3110S.gb/2855662‑2855790 |
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: ‑ |