Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,631,750 | G→T | 26.4% | A14A (GCG→GCT) | livM → | leucine/isoleucine/valine transporter subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,631,750 | 0 | G | T | 26.4% | 46.8 / 13.5 | 34 | A14A (GCG→GCT) | livM | leucine/isoleucine/valine transporter subunit |
Reads supporting (aligned to +/- strand): ref base G (11/14); new base T (7/2); total (18/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.25e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01 |
CTGGGTCGCCCGGAGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTGGATGGCACCAAACTGGTG > minE/2631683‑2631816 | cTGGGTCGCCCGGAGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGa > 1:2599055/1‑69 (MQ=255) cTGGGTCGCCCGGAGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGa > 1:2676293/1‑69 (MQ=255) cTGGGTCGCCCGGAGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGa > 1:2294475/1‑69 (MQ=255) cTGGGTCGCCCGGAGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGa > 1:2693919/1‑69 (MQ=255) cTGGGTCGCCCGGAGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGa > 1:3006672/1‑69 (MQ=255) cTGGGTCGCCCGGAGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGa > 1:3257243/1‑69 (MQ=255) cTGGGTCGCCCGGAGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGa > 1:1647182/1‑69 (MQ=255) gAGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCTATGTTCTTTGTGc < 1:889389/69‑1 (MQ=255) aGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCTATGTTCTTTGTGc < 1:3078406/68‑1 (MQ=255) gagaAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCTATGTTCt > 1:634277/1‑57 (MQ=255) gagaAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCTATGTTCt > 1:938691/1‑57 (MQ=255) gagaAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCTATGTTCt > 1:2083148/1‑57 (MQ=255) gagaAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCTATGTTCt > 1:1335351/1‑57 (MQ=255) gagaAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCTATGTTCt > 1:1221124/1‑57 (MQ=255) gagaAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCTATGTTCt > 1:2345044/1‑57 (MQ=255) gagaAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCTATGTTCt > 1:2500035/1‑57 (MQ=255) gaAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTggcgggc < 1:1161865/69‑1 (MQ=255) tATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCtt < 1:270570/69‑1 (MQ=255) tATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCtt < 1:2440938/69‑1 (MQ=255) aaCCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTggcgggc < 1:1191421/59‑1 (MQ=255) gATGCATATTGCAATGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTggcgggc < 1:1184634/55‑1 (MQ=255) ttGGAATGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCtt < 1:888944/52‑1 (MQ=255) ttGCAATGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCtt < 1:2095395/52‑1 (MQ=255) ttGCAATGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCtt < 1:808427/52‑1 (MQ=255) aTGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTg < 1:1917672/69‑1 (MQ=255) gcgcTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTg < 1:2413408/66‑1 (MQ=255) gcgcTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTg < 1:2938006/66‑1 (MQ=255) gcgcTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTg < 1:1784076/66‑1 (MQ=255) gcgcTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTg < 1:1746936/66‑1 (MQ=255) gcgcTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTg < 1:1382218/66‑1 (MQ=255) gcgcTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTGGAt > 1:81463/1‑69 (MQ=255) ctctGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTGGATGGCACCaa > 1:724202/1‑69 (MQ=255) ctctGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTGGATGGCACCaa > 1:1373647/1‑69 (MQ=255) ctctGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTGGATGGCACCa > 1:1538133/1‑68 (MQ=255) gcgATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTGGATGGCACCAAACtggtg < 1:1957312/69‑1 (MQ=255) | CTGGGTCGCCCGGAGGTAGAGAAAGTATGAAACCGATGCATATTGCAATGGCGCTGCTCTCTGCCGCGATGTTCTTTGTGCTGGCGGGCGTCTTTATGGGCGTGCAACTGGAGCTGGATGGCACCAAACTGGTG > minE/2631683‑2631816 |
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: ‑ |