Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,875,632 | (T)7→8 | intergenic (‑66/+287) | glmU ← / ← atpC | bifunctional N‑acetylglucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyltransferase/F0F1 ATP synthase subunit epsilon |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,875,625 | 1 | . | T | 97.6% | 162.8 / ‑2.8 | 42 | intergenic (‑59/+294) | glmU/atpC | bifunctional N‑acetylglucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyltransferase/F0F1 ATP synthase subunit epsilon |
Reads supporting (aligned to +/- strand): ref base . (0/1); new base T (21/20); total (21/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01 |
ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑T‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC > REL606/3875550‑3875695 | aTAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACa < 1:1841115/76‑1 (MQ=255) aTAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACa > 1:294831/1‑76 (MQ=255) tAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t > 1:11611/1‑76 (MQ=255) aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t < 2:275155/76‑1 (MQ=255) gCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑tt < 1:722919/76‑1 (MQ=255) cATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑ttt > 1:836623/1‑76 (MQ=255) tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑tttttt > 1:1145237/1‑76 (MQ=255) tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑tttttt > 1:219217/1‑76 (MQ=255) tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑tttttt > 1:1738144/1‑76 (MQ=255) cATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGt < 2:192063/76‑1 (MQ=255) cATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGt > 1:965289/1‑76 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTg > 1:1996503/1‑76 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTg < 1:306011/76‑1 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTg > 2:1533358/1‑76 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTg > 2:1710227/1‑76 (MQ=255) tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGa > 1:847097/1‑76 (MQ=255) gcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAAt < 2:8349/76‑1 (MQ=255) cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATg < 1:1339340/76‑1 (MQ=255) tCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGa < 1:227012/76‑1 (MQ=255) ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGAc < 1:1907112/76‑1 (MQ=255) ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGAc < 1:794677/76‑1 (MQ=255) ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGAc < 2:166621/76‑1 (MQ=255) cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACa < 1:1063998/76‑1 (MQ=255) tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACag > 1:813203/1‑76 (MQ=255) gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACaga < 1:496634/76‑1 (MQ=255) aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACagag < 1:1274171/76‑1 (MQ=255) cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAgg < 2:2069160/76‑1 (MQ=255) tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGa > 1:1331197/1‑76 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGAt > 2:126836/1‑76 (MQ=255) aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATaa > 1:1606027/1‑76 (MQ=255) aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATaaa > 1:916974/1‑76 (MQ=255) tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAAt > 1:1429306/1‑76 (MQ=255) gagaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGt > 1:1593587/1‑76 (MQ=255) gagaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGt < 1:882159/76‑1 (MQ=255) agaACGAATTTAAACCGCTTCGCCTTGAAAAAACTACA‑TGTTTGTTCATCGTGAAATGGACAGAGGATAAATTGtt < 2:1659235/76‑1 (MQ=255) gaACGAATTTAATCCGCTTCGCCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTc < 1:1608498/76‑1 (MQ=255) aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa > 2:1343226/1‑76 (MQ=255) aTTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc > 2:1833174/1‑76 (MQ=255) aTTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc > 2:1926076/1‑76 (MQ=255) aaaCCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt > 1:528498/1‑76 (MQ=255) aaCCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt < 1:205213/76‑1 (MQ=255) ccGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt > 1:1605688/1‑76 (MQ=255) cGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATc > 2:174306/1‑76 (MQ=255) cGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATc > 1:1860287/1‑76 (MQ=255) ttCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc < 1:1831413/76‑1 (MQ=255) cACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGc > 1:669473/1‑76 (MQ=255) ttGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc < 1:697666/76‑1 (MQ=255) ttGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc > 1:178953/1‑76 (MQ=255) aaaaaaCTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa < 1:206205/76‑1 (MQ=255) aaaaaCTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa < 1:915920/76‑1 (MQ=255) cTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc < 1:496334/76‑1 (MQ=255) | ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑T‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC > REL606/3875550‑3875695 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |