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% | 95.1 / ‑3.0 | 26 | 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 (17/8); total (17/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.46e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCA > REL606/3875551‑3875696 | tAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t > 2:288384/1‑76 (MQ=255) cATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt‑t > 2:356849/1‑76 (MQ=255) cATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt‑t < 2:832161/76‑1 (MQ=255) cATTATTCAACATACGCGTCCTGACCGTAATGTGAGAACANNTTTAAACCGCNNCGCCTTGAAAAAACTACA‑ttt‑t > 1:1088671/1‑76 (MQ=255) tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACGGCTTCACCTTGAAAAAACTACA‑ttt‑tttt < 2:197174/76‑1 (MQ=255) tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt‑tttt > 1:195386/1‑76 (MQ=255) tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTGAAACCGCTTCACCTTGAAAAAACTACA‑ttt‑tttt < 1:582884/76‑1 (MQ=255) tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt‑tttttca < 1:1171382/76‑4 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTg > 1:1016464/1‑76 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTg > 1:1016450/1‑76 (MQ=255) tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGa > 1:200004/1‑76 (MQ=255) aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGaa > 2:628511/1‑76 (MQ=255) gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATgg < 1:573429/76‑1 (MQ=255) tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACag > 1:894193/1‑76 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGAt < 1:481512/76‑1 (MQ=255) aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATaa > 1:1110578/1‑76 (MQ=255) gagaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGt > 2:982031/1‑76 (MQ=255) agaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGtt > 2:543936/1‑76 (MQ=255) aCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCaa > 2:136236/1‑76 (MQ=255) ttAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg < 2:405011/76‑1 (MQ=255) ttAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTAAATAACgg < 1:467069/76‑1 (MQ=255) tAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa > 1:954342/1‑76 (MQ=255) aaaCCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt > 2:392965/1‑76 (MQ=255) aaaCCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt > 2:392977/1‑76 (MQ=255) aaCCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt < 2:490559/76‑1 (MQ=255) cGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATc > 1:353548/1‑76 (MQ=255) cTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc > 2:54084/1‑76 (MQ=255) ttGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc > 2:361274/1‑76 (MQ=255) tGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg < 2:693069/76‑1 (MQ=255) tGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg < 1:262949/76‑1 (MQ=255) gAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGa > 1:452316/1‑76 (MQ=255) aaCTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg < 1:263920/76‑1 (MQ=255) cTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc > 2:786530/1‑76 (MQ=255) tACA‑TTTGTGTTCATCGTGAAATGGACAGAGGATAAATTGTTCAGTAACGGATTATCCCCGCTGACGAAAGCGCCa < 1:864622/76‑1 (MQ=255) | TAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCA > REL606/3875551‑3875696 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |