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 | 93.9% | 118.6 / 0.8 | 33 | 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 . (2/0); new base T (17/14); total (19/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.96e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01 |
AGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑T‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGC > REL606/3875552‑3875692 | aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t > 1:865179/1‑76 (MQ=255) aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t < 1:52584/76‑1 (MQ=255) aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t < 1:448144/76‑1 (MQ=255) gCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑tt < 2:959329/76‑1 (MQ=255) gCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCACTTCACCTTGAAAAAACTACA‑t‑tt > 1:329270/1‑76 (MQ=255) tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑tttttttca > 1:1200608/1‑73 (MQ=255) cATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGt > 1:499301/1‑76 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTg < 1:315718/76‑1 (MQ=255) tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGa > 2:2740/1‑76 (MQ=255) tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGa > 1:548989/1‑76 (MQ=255) cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATg > 1:683159/1‑76 (MQ=255) ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGAc < 2:72772/76‑1 (MQ=255) tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACag < 1:596280/76‑1 (MQ=255) tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TATTTTTTCATAGTGAAATGGACag > 1:1113449/1‑76 (MQ=255) gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACaga < 1:846173/76‑1 (MQ=255) cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAgg < 1:827871/76‑1 (MQ=255) cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAgg < 1:855729/76‑1 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGANAAAACTAGA‑TGTTTTATCATCGTGAAATGGACAGNGGAt > 1:151561/1‑76 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGAt > 2:938170/1‑76 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGAt > 2:56214/1‑76 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGAt > 2:938161/1‑76 (MQ=255) aaTTTGAGACCGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATaa > 2:1115588/1‑76 (MQ=255) ttGAGAATGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAAtt > 1:1267052/1‑76 (MQ=255) agaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTCAAATGGACAGAGGATAAATTGtt > 1:1122226/1‑76 (MQ=255) aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa < 1:773545/76‑1 (MQ=255) aTTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc < 1:1059853/76‑1 (MQ=255) ttAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg > 1:782797/1‑76 (MQ=255) ttAAACCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg > 2:1266119/1‑76 (MQ=255) aaCCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCTATAACGGAtt < 1:1198979/76‑1 (MQ=255) aaCCGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt < 2:529601/76‑1 (MQ=255) ccGCTTCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt < 1:283300/76‑1 (MQ=255) tCACCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg > 1:519249/1‑76 (MQ=255) aCCTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCt < 1:423557/76‑1 (MQ=255) ccTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg < 1:1253795/76‑1 (MQ=255) ccTTGAAAAAACTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg < 1:379554/76‑1 (MQ=255) aaaaaCTACATT‑TTTTTTCATCGTGAANTGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa > 1:1297893/1‑76 (MQ=255) aaaaaCTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa > 2:274706/1‑76 (MQ=255) aaaCTACATT‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc > 2:280857/1‑76 (MQ=255) aaaCTACATT‑TTTTTTCATCGTCAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc > 2:145901/1‑76 (MQ=255) | AGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑T‑TTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGC > REL606/3875552‑3875692 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |