Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 3,473,895 | A→G | 71.4% | E162G (GAG→GGG) | trmA → | tRNA (uracil‑5‑)‑methyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 3,473,895 | 0 | A | G | 71.4% | 62.7 / 29.2 | 42 | E162G (GAG→GGG) | trmA | tRNA (uracil‑5‑)‑methyltransferase |
Reads supporting (aligned to +/- strand): ref base A (4/8); new base G (17/13); total (21/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.06e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01 |
CGCGATGCACTGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAG > W3110S.gb/3473828‑3473963 | cgcgATGCACTGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCt > 1:2709679/1‑71 (MQ=255) cACTGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAg > 1:553833/1‑71 (MQ=255) cACTGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAg < 1:273863/71‑1 (MQ=255) cACTGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAg < 1:1954498/71‑1 (MQ=255) tGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATc < 1:824313/66‑1 (MQ=255) cgcgcgcAGAATCTGAATGTGCATCTGATTGGTCGGGCAACTAAAACCTAAATCGAGCTGGATCAGGATTa < 1:2683055/71‑1 (MQ=255) cgcgcgcAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGAtt > 1:2746583/1‑70 (MQ=255) cgcgcgcAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTa > 1:1220863/1‑71 (MQ=255) cgcgcgcAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTa < 1:261753/71‑1 (MQ=255) cgcgcgcAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTa > 1:2596303/1‑71 (MQ=255) cgcgcgcAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTa < 1:1614877/71‑1 (MQ=255) gAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATg < 1:78377/71‑1 (MQ=255) tCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAAc > 1:1585510/1‑71 (MQ=255) cTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGa > 1:1237045/1‑56 (MQ=255) cTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGa > 1:1971792/1‑56 (MQ=255) tGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGt > 1:96775/1‑71 (MQ=255) gtgCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGaa > 1:166411/1‑63 (MQ=255) tgCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTgg < 1:57469/42‑1 (MQ=255) tgCATCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGa > 1:617624/1‑58 (MQ=255) tgCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACg < 1:72337/64‑1 (MQ=255) tgCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACg < 1:1399398/64‑1 (MQ=255) tCTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGAt < 1:2432326/46‑1 (MQ=255) cTGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGt > 1:1799806/1‑60 (MQ=255) tGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCg > 1:569648/1‑69 (MQ=255) tGATTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCg > 1:566358/1‑69 (MQ=255) tGATTGGACGGGCAACGAAAACCAAAATAGGTCTGGAta > 1:1698628/1‑38 (MQ=21) aaTTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGaa > 1:623519/2‑55 (MQ=255) aTTGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAgg > 1:1066419/1‑71 (MQ=255) tGGTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTa > 1:333794/1‑42 (MQ=255) ggTCGGGCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGa < 1:2453458/70‑1 (MQ=255) gggCAACGAAAACCAAAATCGGGCTGGATCATGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAaga < 1:2018838/70‑1 (MQ=255) gggCAACGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAaga < 1:1016961/70‑1 (MQ=255) aCGAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGa < 1:529890/61‑1 (MQ=255) gAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATc > 1:1811202/1‑70 (MQ=255) gAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATc > 1:2515356/1‑70 (MQ=255) gAAAACCAAAATCGGGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCt < 1:2449333/71‑1 (MQ=255) aCCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTAc < 1:1204936/69‑1 (MQ=255) aaaaTCGGGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAaga < 1:758166/56‑1 (MQ=255) aaaTCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAaga < 1:519022/55‑1 (MQ=255) aaaTCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAaga < 1:809293/55‑1 (MQ=255) aaTCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAgatga < 1:46842/59‑1 (MQ=255) gAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAg > 1:2791733/1‑70 (MQ=255) | CGCGATGCACTGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAG > W3110S.gb/3473828‑3473963 |
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: ‑ |