| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 2,855,732 | G→A | 13.1% | 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 | 13.1% | 145.3 / 11.2 | 54 | intergenic (‑270/‑17) | ygbA/mutS | hypothetical protein/methyl‑directed mismatch repair protein |
| Reads supporting (aligned to +/- strand): ref base G (17/29); new base A (4/3); total (21/33) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.15e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01 | |||||||||||
TATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTC > W3110S.gb/2855668‑2855802 | tATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGTAACCgg < 1:9062/71‑1 (MQ=255)tATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCgg < 1:2189242/71‑1 (MQ=255) tCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGa < 1:1337240/70‑1 (MQ=255) tCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGa < 1:3889494/70‑1 (MQ=255) tCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGa < 1:1935422/70‑1 (MQ=255) tCATGTGCGCCTTATGTGATTACAACCAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGa < 1:4637267/70‑1 (MQ=255) aTGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGAc > 1:1746979/1‑69 (MQ=255) aTGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACa > 1:5106685/1‑70 (MQ=255) aTGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACa > 1:3919250/1‑70 (MQ=255) aTGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACAt < 1:1072738/71‑1 (MQ=255) gcCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAAcccc < 1:2470756/71‑1 (MQ=255) cTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCAt < 1:95789/71‑1 (MQ=255) ttATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGa < 1:5076217/53‑1 (MQ=255) aTGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGaa > 1:2851972/1‑52 (MQ=255) aTGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATg < 1:3557281/69‑1 (MQ=255) aTGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAg > 1:1965695/1‑71 (MQ=255) aCAACTAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa > 1:4070476/1‑68 (MQ=255) aCAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa > 1:4079594/1‑68 (MQ=255) aCAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa > 1:1512268/1‑68 (MQ=255) cAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAg < 1:4323522/70‑1 (MQ=255) cAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGa < 1:1109101/71‑1 (MQ=255) aaCGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGAc < 1:3140541/48‑1 (MQ=255) aaCGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGaa < 1:257582/71‑1 (MQ=255) aaCGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGaa < 1:1948917/71‑1 (MQ=255) aaaaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGC‑ATAGAAAAttt > 1:584102/1‑71 (MQ=255) aaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTc < 1:3945470/71‑1 (MQ=255) aTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATaa > 1:5094849/1‑45 (MQ=255) aaaCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgc < 1:4336931/71‑1 (MQ=255) aCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGa > 1:1392546/1‑66 (MQ=255) ccATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgccca < 1:4721304/71‑1 (MQ=255) tCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa < 1:2908267/49‑1 (MQ=255) tCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa < 1:382216/49‑1 (MQ=255) tCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcc < 1:1587774/66‑1 (MQ=255) acacCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcccata < 1:5018441/68‑1 (MQ=255) cacCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAAt > 1:888996/1‑54 (MQ=255) cacCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcccatacgcc < 1:1635658/71‑1 (MQ=255) cacCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcccat > 1:3563832/1‑66 (MQ=255) cacCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgccca > 1:1798739/1‑65 (MQ=255) taattTAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgggca < 1:85384‑M1/62‑34 (MQ=255) aTTTAATATCAGGGAACCGGACATAACCCCATGAGTGc < 1:3489432/38‑1 (MQ=255) tttAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgggcagggggcaagt < 1:2030762‑M1/71‑44 (MQ=255) tttAATATCAGTAAACCGGACATAACCCatgaagaaaaa > 1:3642286‑M1/1‑28 (MQ=255) tttAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTc < 1:3207574/50‑1 (MQ=255) tttAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTc < 1:3213705/50‑1 (MQ=255) tttAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcccatacgcc < 1:4001933/64‑1 (MQ=255) tttAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCatgat < 1:1546536/70‑1 (MQ=255) ttAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgggcagggggcaagta > 1:1119572‑M1/1‑27 (MQ=255) tAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgggcagggggcaagtat > 1:4048855‑M1/1‑26 (MQ=255) tAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgggcagggggcaagtat > 1:2371928‑M1/1‑26 (MQ=255) tAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCatga > 1:1636736/1‑67 (MQ=255) aaTATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAAt > 1:4122962/1‑44 (MQ=255) aaTATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATgca < 1:4417024/70‑1 (MQ=255) atatCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgggcagggggcaagtatc < 1:1480321‑M1/70‑47 (MQ=255) cAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTa > 1:3843244/1‑71 (MQ=255) gAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcccatacgccca > 1:1146162/1‑54 (MQ=255) gAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTAtctc < 1:880781/71‑1 (MQ=255) | TATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTC > W3110S.gb/2855668‑2855802 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |