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 | 96.2% | 192.6 / 1.1 | 52 | 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 (23/27); total (25/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.26e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01 |
TAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑A‑TTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC > REL606/3875551‑3875695 | tAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑A‑t > 1:313997/1‑76 (MQ=255) gCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑A‑ttt > 1:1082904/1‑76 (MQ=255) attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑A‑ttttt > 1:294946/1‑76 (MQ=255) tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑A‑ttttttttca > 1:1694670/1‑73 (MQ=255) cATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGt > 1:53799/1‑76 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTg < 1:1344368/76‑1 (MQ=255) tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGa < 2:958277/76‑1 (MQ=255) aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGaa < 1:1840672/76‑1 (MQ=255) aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGaa > 1:1962511/1‑76 (MQ=255) cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATg < 2:747977/76‑1 (MQ=255) cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATg < 1:806722/76‑1 (MQ=255) cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCAGCGTGAAATg < 2:1608607/76‑1 (MQ=255) gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATgg < 2:425536/76‑1 (MQ=255) tCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGa > 1:755491/1‑76 (MQ=255) cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACa < 2:1933988/76‑1 (MQ=255) aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACagag > 1:1162596/1‑76 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGAt < 1:1442374/76‑1 (MQ=255) aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATaa > 2:73023/1‑76 (MQ=255) aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATaaa < 2:116731/76‑1 (MQ=255) tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAAt < 1:740329/76‑1 (MQ=255) ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt < 2:278812/76‑1 (MQ=255) ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAAtt < 2:987316/76‑1 (MQ=255) tGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTg < 2:430579/76‑1 (MQ=255) tGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTg < 2:1912940/76‑1 (MQ=255) gaACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTc < 1:1110324/76‑1 (MQ=255) aaTTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa < 2:827329/76‑1 (MQ=255) aaTTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa < 2:616894/76‑1 (MQ=255) aTTTAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc > 1:1314597/1‑76 (MQ=255) tttAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg < 1:1539523/76‑1 (MQ=255) tttAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg > 1:1899080/1‑76 (MQ=255) tAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa < 2:438674/76‑1 (MQ=255) tAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa > 2:227717/1‑76 (MQ=255) tAAACCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa > 1:243434/1‑76 (MQ=255) tAAACCGCTTCACCTCGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa > 2:701153/1‑76 (MQ=255) aaaCCGCTTCACCTTGAAAACACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt > 1:1898190/1‑76 (MQ=255) aaaCCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt > 1:697548/1‑76 (MQ=255) aaCCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt < 1:645445/76‑1 (MQ=255) aaCCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt > 2:42854/1‑76 (MQ=255) aaCCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt < 1:921963/76‑1 (MQ=255) aCCGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa < 1:1680294/76‑1 (MQ=255) ccGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt > 1:1492940/1‑76 (MQ=255) ccGCTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt < 1:1589801/76‑1 (MQ=255) cTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc < 2:1855510/76‑1 (MQ=255) cTTCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc > 1:873539/1‑76 (MQ=255) ttCACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc < 2:1195474/76‑1 (MQ=255) cACCTTGAAAAAACTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGc > 1:447695/1‑76 (MQ=255) aaaaaaCTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa > 1:1672696/1‑76 (MQ=255) aaaaaCTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa > 2:1644609/1‑76 (MQ=255) aaaaaCTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa > 2:1127040/1‑76 (MQ=255) aaaaaCTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCAGCTGACGaaa > 2:1688107/1‑76 (MQ=255) aaaaaCTAC‑A‑TTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg > 2:1368165/1‑76 (MQ=255) aaaaCTACAA‑TTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg > 1:330301/1‑76 (MQ=255) aaaCTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc > 1:1842572/1‑76 (MQ=255) aaaCTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc < 2:1345947/76‑1 (MQ=255) cTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc > 2:1062487/1‑76 (MQ=255) cTAC‑ATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc > 1:963001/1‑76 (MQ=255) | TAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTAC‑A‑TTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC > REL606/3875551‑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: ‑ |