Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,225,948 | G→C | 21.6% | V129V (GTG→GTC) | dusB → | tRNA‑dihydrouridine synthase B |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,225,948 | 0 | G | C | 21.6% | 177.2 / 46.3 | 88 | V129V (GTG→GTC) | dusB | tRNA‑dihydrouridine synthase B |
Reads supporting (aligned to +/- strand): ref base G (28/41); new base C (13/6); total (41/47) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.94e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01 |
ATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAA > minE/2225878‑2226017 | aTCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTg < 1:3783823/71‑1 (MQ=255) tCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTg > 1:371905/1‑70 (MQ=255) aaGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGc < 1:1197369/71‑1 (MQ=255) aaGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGc < 1:824064/71‑1 (MQ=255) aaGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGc < 1:2968828/71‑1 (MQ=255) gCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGc < 1:3736900/69‑1 (MQ=255) tCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTg > 1:249402/1‑71 (MQ=255) gCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGa < 1:846568/71‑1 (MQ=255) aGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTg < 1:149032/67‑1 (MQ=255) aGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTg < 1:3351605/67‑1 (MQ=255) ggTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGt < 1:1773993/71‑1 (MQ=255) tCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTcc > 1:474043/1‑71 (MQ=255) cAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCt > 1:1836973/1‑71 (MQ=255) aGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCt < 1:1100646/70‑1 (MQ=255) cccTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGaa > 1:3746082/1‑52 (MQ=255) tctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGAAGTTc < 1:1390836/63‑1 (MQ=255) ttGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAg > 1:844470/1‑52 (MQ=255) ttGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGa > 1:109602/1‑56 (MQ=255) ttGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCaa > 1:737736/1‑47 (MQ=255) ttGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCaa > 1:4067667/1‑47 (MQ=255) tGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTg < 1:3999125/44‑1 (MQ=255) tGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGt < 1:3767477/58‑1 (MQ=255) tGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCt > 1:3194082/1‑70 (MQ=255) gCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTg > 1:2407788/1‑70 (MQ=255) gCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTg < 1:3061687/70‑1 (MQ=255) gCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTg > 1:4089697/1‑70 (MQ=255) aGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGc < 1:302763/46‑1 (MQ=255) aGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTg < 1:4049192/68‑1 (MQ=255) aGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTg < 1:3317275/68‑1 (MQ=255) tACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGa < 1:3367987/50‑1 (MQ=255) tACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTa < 1:3619858/61‑1 (MQ=255) tACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATg < 1:563923/43‑1 (MQ=255) cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGATTGCAGTGGACGTTCCTGTTACCCTGAAGAtt < 1:1251302/70‑1 (MQ=255) cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTa < 1:603562/59‑1 (MQ=255) cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTa < 1:904980/59‑1 (MQ=255) cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTc > 1:552476/1‑71 (MQ=255) cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTc > 1:3070399/1‑71 (MQ=255) cccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGt > 1:1257279/1‑57 (MQ=255) ccGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCg > 1:961502/1‑71 (MQ=255) cGGATGTCGTTAAATCGATCCTTACCGAGGTCGTg < 1:2527596/35‑1 (MQ=255) cGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGa < 1:978580/66‑1 (MQ=255) cGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGc < 1:480107/71‑1 (MQ=255) ggATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGc > 1:1019297/1‑70 (MQ=255) gATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGc > 1:2336690/1‑38 (MQ=255) gATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGaa < 1:862912/62‑1 (MQ=255) gATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGaa < 1:1559226/62‑1 (MQ=255) aTGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCGTGAAGATTCGCAcc > 1:3568028/1‑71 (MQ=255) tGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTAcc > 1:2867657/1‑55 (MQ=255) tGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCg > 1:872906/1‑71 (MQ=255) tCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGtt < 1:3296814/44‑1 (MQ=255) tCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGtt < 1:4118384/44‑1 (MQ=255) tCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCg < 1:2935189/69‑1 (MQ=255) tCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCg < 1:638173/69‑1 (MQ=255) tCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGc < 1:2194316/71‑1 (MQ=255) tCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGc > 1:1892846/1‑71 (MQ=255) cGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCt > 1:4072572/1‑71 (MQ=255) gTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGa < 1:26493/55‑1 (MQ=255) gTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTg < 1:371425/71‑1 (MQ=255) gTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTg < 1:2442608/71‑1 (MQ=255) gTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGATg > 1:1704525/1‑71 (MQ=255) gTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGtt > 1:834632/1‑48 (MQ=255) tAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTg > 1:2925710/1‑69 (MQ=255) aaaTCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTggg > 1:902497/1‑70 (MQ=255) aaaTCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTggg < 1:948761/70‑1 (MQ=255) aaaTCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTggg > 1:1040573/1‑70 (MQ=255) aaTCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGc < 1:4113760/70‑1 (MQ=255) cGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTAcc > 1:1707126/1‑44 (MQ=255) cGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGa < 1:861811/51‑1 (MQ=255) cGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCg > 1:2639476/1‑71 (MQ=255) cGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCg > 1:2734110/1‑71 (MQ=255) gATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGa > 1:3753712/1‑50 (MQ=255) gATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCgg < 1:4043621/71‑1 (MQ=255) gATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCgg < 1:894718/71‑1 (MQ=255) gATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCgg < 1:3449094/71‑1 (MQ=255) gATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCg > 1:3700855/1‑70 (MQ=255) aTCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGt < 1:706093/38‑1 (MQ=255) tCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGt > 1:2208856/1‑37 (MQ=255) tCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCAc > 1:2570900/1‑66 (MQ=255) cTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGt < 1:1812174/35‑1 (MQ=255) tACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAacac > 1:3741284/1‑70 (MQ=255) tACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAaca > 1:3100097/1‑69 (MQ=255) tACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACAcc > 1:759367/1‑71 (MQ=255) aCCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCg > 1:3348407/1‑71 (MQ=255) ccGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGt < 1:1447235/71‑1 (MQ=255) ccGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGt < 1:3423814/71‑1 (MQ=255) ccGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGt < 1:3193126/71‑1 (MQ=255) cGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTa < 1:2980963/71‑1 (MQ=255) tCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGAt < 1:3172347/37‑1 (MQ=255) cGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGc > 1:2270427/1‑46 (MQ=255) cGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGc > 1:766059/1‑70 (MQ=255) gTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCg > 1:1238930/1‑70 (MQ=255) tGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGaa < 1:3317308/71‑1 (MQ=255) | ATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAA > minE/2225878‑2226017 |
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: ‑ |