| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 236,101 | A→C | 55.6% | D12A (GAT→GCT) | dnaQ → | DNA polymerase III epsilon subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 236,101 | 0 | A | C | 55.6% | 68.3 / 62.5 | 63 | D12A (GAT→GCT) | dnaQ | DNA polymerase III epsilon subunit |
| Reads supporting (aligned to +/- strand): ref base A (18/10); new base C (13/22); total (31/32) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.37e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGC > W3110S.gb/236031‑236170 | aaCTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa < 1:2771/71‑1 (MQ=255) gTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAcc > 1:229342/1‑71 (MQ=255) cgtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAcc < 1:873740‑M2/62‑1 (MQ=255) tAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAcc > 1:937889/1‑70 (MQ=255) tacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACcgc < 1:610221‑M2/64‑2 (MQ=255) aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa < 1:535791/70‑1 (MQ=255) aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa > 1:572333/1‑70 (MQ=255) aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa > 1:1343632/1‑70 (MQ=255) aCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACAGAAaccac < 1:801053/71‑1 (MQ=255) aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAcc < 1:547749/61‑1 (MQ=255) tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa < 1:788515/71‑1 (MQ=255) cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt > 1:720620/1‑71 (MQ=255) tGACATAAATGACCGCTATGAGCACTGCAATTAGACGCCAGATCGTTCTCGATACCGAAaccac > 1:1415163/1‑64 (MQ=255) tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg < 1:1119568/67‑1 (MQ=255) aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt < 1:300616/68‑1 (MQ=255) aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa > 1:1318308/1‑71 (MQ=255) tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa < 1:817511/71‑1 (MQ=255) aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAc < 1:1096326/48‑1 (MQ=255) aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAc < 1:1559931/48‑1 (MQ=255) aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg > 1:1045822/1‑64 (MQ=255) aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt > 1:169967/1‑71 (MQ=255) aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg < 1:1126824/70‑1 (MQ=255) tATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgc > 1:1134971/1‑69 (MQ=255) aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt > 1:1320279/1‑60 (MQ=255) cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg < 1:1241498/71‑1 (MQ=255) cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa > 1:608437/1‑46 (MQ=39) cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa > 1:1182520/1‑46 (MQ=39) cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCAGTATGaa > 1:821606/1‑46 (MQ=255) aTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca > 1:115445/1‑71 (MQ=255) ttACACGCCAGATCGTTCTCTATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg < 1:634823/63‑1 (MQ=255) ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg < 1:260090/66‑1 (MQ=255) ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa < 1:988035/42‑1 (MQ=255) ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca > 1:1125859/1‑70 (MQ=255) tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt < 1:866224/61‑1 (MQ=255) tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac > 1:1243487/1‑70 (MQ=255) tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac > 1:1622425/1‑70 (MQ=255) acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTa > 1:571076/1‑37 (MQ=37) acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa > 1:104658/1‑46 (MQ=39) acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac > 1:1575206/1‑69 (MQ=255) acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa > 1:325930/1‑46 (MQ=255) acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg > 1:359616/1‑64 (MQ=255) acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac > 1:665423/1‑69 (MQ=255) acacGCCAGATCGTTCTCGATACCGAAACCACAGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac > 1:455722/1‑69 (MQ=255) acGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCGACAAGa < 1:227011/71‑1 (MQ=255) cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcaca < 1:498855/67‑1 (MQ=255) cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcaca < 1:512352/67‑1 (MQ=255) cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcaca < 1:1406534/67‑1 (MQ=255) cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcaca < 1:323610/67‑1 (MQ=255) gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg > 1:1592496/1‑61 (MQ=255) ccATATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcgc < 1:863819/50‑1 (MQ=255) ccAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcgc < 1:915677/50‑1 (MQ=255) gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAT‑AAGGCCACAAGATCATTga < 1:215379/71‑1 (MQ=255) gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg < 1:124974/53‑1 (MQ=255) aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt > 1:484483/1‑69 (MQ=255) aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga < 1:192226/71‑1 (MQ=255) aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg < 1:1503241/40‑1 (MQ=255) tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt > 1:1366007/1‑68 (MQ=255) tCGTTCTCGCTACCGCAACCACAGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga < 1:1157158/70‑1 (MQ=255) cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa > 1:32404/1‑65 (MQ=255) gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgaga < 1:483752/70‑1 (MQ=255) ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt < 1:189585/71‑1 (MQ=255) ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac < 1:1576465/56‑1 (MQ=255) ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg > 1:1386360/1‑71 (MQ=255) gATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGACACAAGATCATTGAGATTGGTGc > 1:451643/1‑71 (MQ=255) | AACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGC > W3110S.gb/236031‑236170 |
| 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: ‑ |