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 | 92.0% | 81.8 / 1.4 | 25 | 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 . (1/1); new base T (12/11); total (13/12) | |||||||||||
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.57e-01 |
AGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCAT > REL606/3875552‑3875697 | aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tt > 1:549432/1‑76 (MQ=255) attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt‑tt > 2:758166/1‑76 (MQ=255) tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt‑tttt < 2:734900/76‑1 (MQ=255) tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt‑tttt > 2:485022/1‑76 (MQ=255) attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt‑ttttt > 1:321825/1‑75 (MQ=255) attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt‑ttttt < 1:402955/76‑2 (MQ=255) attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt‑ttttt > 2:232841/1‑75 (MQ=255) tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt‑tttttca > 1:533403/1‑73 (MQ=255) aCATACGCGTCCTGACTGTNANTAGAGAANGAATTTAAACCGCTTCACCTTNAANAAGCTACATTTT‑TTTTCATCg < 2:120535/76‑1 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTg > 1:458389/1‑76 (MQ=255) aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGaa < 1:445616/76‑1 (MQ=255) gcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAAt < 2:370436/76‑1 (MQ=255) tCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGa > 2:365496/1‑76 (MQ=255) tCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGa > 2:365490/1‑76 (MQ=255) cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCGCCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACa < 2:158024/76‑1 (MQ=255) tGACTGTAATTTGAGAACGAATTTAAACCGCTGCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACag < 1:43557/76‑1 (MQ=255) gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACaga > 2:416436/1‑76 (MQ=255) cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAgg > 1:377320/1‑76 (MQ=255) cTGTAATTTGAGAACGAATGTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAgg > 1:715252/1‑76 (MQ=255) tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGa > 2:268978/1‑76 (MQ=255) gagaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGt < 2:553574/76‑1 (MQ=255) aaTTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa < 1:630388/76‑1 (MQ=255) aTTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc < 2:555616/76‑1 (MQ=255) tttAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg > 1:670570/1‑76 (MQ=255) ttAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGGTCAATAACgg < 1:64395/76‑1 (MQ=255) ttAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGACAAATTGTTCAATAACgg > 2:795701/1‑76 (MQ=255) tAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa > 2:40352/1‑76 (MQ=255) tAAACCGCTTCACCTTGAAAAAACTACA‑TTTCTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa > 2:376346/1‑76 (MQ=255) aaaCCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt < 2:466257/76‑1 (MQ=255) cGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATc > 2:556089/1‑76 (MQ=255) gCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcc > 1:20361/1‑76 (MQ=255) aaaaCTACATTTT‑TTTTCATCGTGAAATGGACAGAGGANAATTTGTTCAATAACGGATTATCCCCGCTGACGAAAg < 1:823542/76‑1 (MQ=255) tACA‑TTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCAt < 1:82995/76‑1 (MQ=255) | AGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCAT > REL606/3875552‑3875697 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |