| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,530,398 | G→A | 23.8% | A26T (GCA→ACA) | mutM → | formamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,530,398 | 0 | G | A | 23.8% | 121.7 / 37.1 | 63 | A26T (GCA→ACA) | mutM | formamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase |
| Reads supporting (aligned to +/- strand): ref base G (26/22); new base A (7/8); total (33/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.69e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01 | |||||||||||
AATTACCCGAAGTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGACCAAC > minE/2530330‑2530464 | aaTTACCCGAAGTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATg < 1:2144665/69‑1 (MQ=255) tACCCGAAGTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATg > 1:1485971/1‑66 (MQ=255) tACCCGAAGTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATACAgtg > 1:727380/1‑71 (MQ=255) tACCCGAAGTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATACAgtg > 1:2527213/1‑71 (MQ=255) aCCCGAAGTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCa > 1:2618841/1‑67 (MQ=255) aaGTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTgcgc < 1:3740012/71‑1 (MQ=255) gAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCa > 1:2926562/1‑57 (MQ=255) aCCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATACAGTGGTGCGCAACGGACg > 1:3102137/1‑71 (MQ=255) ccAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGc > 1:1103848/1‑71 (MQ=255) cAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATACAgt < 1:1001801/54‑1 (MQ=255) cAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATACAgt < 1:2823865/54‑1 (MQ=255) cAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATACAgt < 1:1380848/54‑1 (MQ=255) cAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATACAgt < 1:1296541/54‑1 (MQ=255) cAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATACAgt < 1:447697/54‑1 (MQ=255) aGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATACAGTGGTGCGCAACGGACGCtt > 1:2064581/1‑71 (MQ=255) gccgccGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATg < 1:3433241/48‑1 (MQ=255) ccgccgCGGCATAGAACCTCATCTCGTTGGTGCAACCATTCTTCATg < 1:2346359/47‑1 (MQ=255) cATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATACAgt < 1:1340505/42‑1 (MQ=255) aCCGCATCTCGGTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCag < 1:2759591/71‑1 (MQ=255) cGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCagaa < 1:3620065/71‑1 (MQ=255) cGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCag > 1:3105212/1‑69 (MQ=255) gCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTgcg > 1:3853736/1‑54 (MQ=255) gCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCag > 1:3416963/1‑68 (MQ=255) gCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCag > 1:1652250/1‑68 (MQ=255) aTCTCGTTGGTGCAACCATTCTTCATACAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAagag < 1:2982327/71‑1 (MQ=255) tctcGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCt < 1:99970/47‑1 (MQ=255) tctcGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCt < 1:1918508/47‑1 (MQ=255) tctcGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCagaa < 1:3116760/67‑1 (MQ=255) tctcGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCagaa < 1:1066919/67‑1 (MQ=255) tctcGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCagaa < 1:1989796/67‑1 (MQ=255) tctcGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCagaa < 1:173081/67‑1 (MQ=255) tctcGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCagaa < 1:1766491/67‑1 (MQ=255) tcGTTGGTGCAACCATTCTTCATACAGTGGTGCGCaa > 1:2834623/1‑37 (MQ=255) tcGTTGGTGCAACCATTCTTCATACAGTGGTGCGCaa > 1:766729/1‑37 (MQ=255) cAACCATTCTTCATGCAGTGGTGCGCAACGGACGCt > 1:1087059/1‑36 (MQ=255) cAACCATTCTTCATGCAGTGGTGCGCAACGGACGCt > 1:825576/1‑36 (MQ=255) cAACCATTCTTCATGCAGTGGTGCGCAACGGACGCt > 1:1242384/1‑36 (MQ=255) cAACCATTCTTCATGCAGTGGTGCGCAACGGACGCt > 1:1577113/1‑36 (MQ=255) cAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGttt > 1:1743514/1‑70 (MQ=255) cAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGttt > 1:2179586/1‑70 (MQ=255) aCCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTAc > 1:1180372/1‑63 (MQ=255) aCCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTAc > 1:3107193/1‑63 (MQ=255) aCCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTAc > 1:1988374/1‑63 (MQ=255) ccATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCag < 1:3440108/51‑1 (MQ=255) ccATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCag < 1:1232888/51‑1 (MQ=255) ccATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCag < 1:3370631/51‑1 (MQ=255) ccATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCag < 1:3199546/51‑1 (MQ=255) ccATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCag < 1:304886/51‑1 (MQ=255) cATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGc > 1:1080884/1‑70 (MQ=255) ttcttcATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGAtg > 1:3172528/1‑55 (MQ=255) ttcttcATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATc > 1:3793447/1‑56 (MQ=255) ttcttcATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATc > 1:3931505/1‑56 (MQ=255) ttcttcATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATc > 1:245939/1‑56 (MQ=255) ttcttcATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATc > 1:1890768/1‑56 (MQ=255) ttcttcATACAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCg < 1:1109372/61‑1 (MQ=255) tcttcATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGac > 1:1757309/1‑70 (MQ=255) cttcATACAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGtt > 1:2709450/1‑61 (MQ=255) tcATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGaccaac > 1:3593368/1‑71 (MQ=255) tcATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGaccaac > 1:260340/1‑71 (MQ=255) tcATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGaccaac > 1:697491/1‑71 (MQ=255) cATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGaccaac < 1:3675804/70‑1 (MQ=255) cATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGaccaa < 1:203123/69‑1 (MQ=255) cATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGaccaa < 1:542499/69‑1 (MQ=255) cATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGaccaa < 1:1119055/69‑1 (MQ=255) cATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGaccaa < 1:1251667/69‑1 (MQ=255) tGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAagag < 1:1873923/46‑1 (MQ=255) tGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAagag < 1:1876533/46‑1 (MQ=255) | AATTACCCGAAGTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGACCAAC > minE/2530330‑2530464 |
| 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: ‑ |