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.4% | 96.3 / ‑2.7 | 28 | 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 (13/14); total (13/15) | |||||||||||
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.97e-01 |
ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTTT‑T‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGC > REL606/3875550‑3875692 | aTAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACa > 1:289268/1‑76 (MQ=255) tAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t > 2:626815/1‑76 (MQ=255) aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tt > 2:54557/1‑76 (MQ=255) attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tttt‑t‑ttt < 2:369494/76‑2 (MQ=255) attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tttt‑t‑ttt < 1:798245/76‑2 (MQ=255) tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tttt‑t‑tttca < 2:599586/76‑4 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTg < 1:89836/76‑1 (MQ=255) cgTCCTGACTGTAATTTGAGAACGAATTTTAACCGCTTCACCTTGAAAAAATTACA‑TTTT‑TNTTCATCGTGAAATg < 1:6842/76‑1 (MQ=255) gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑T‑TTCCTCGTGAAATgg < 1:30874/76‑1 (MQ=255) cccgACTGTCATTTGATAACGCCTCGAAACCGCTTCCCCCTGAAAACACTACA‑TTTTCT‑TTCATCGTGAAATGCAc > 2:1076209/4‑76 (MQ=255) cTGACTGTAATTTGAGAACGCATTTAAACCGCTACACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACa < 1:1199521/76‑1 (MQ=255) tGACTGTAATTTGAGAACGAATTTAAANCGCTTCACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACag > 1:1059328/1‑76 (MQ=255) tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTNGAGAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACag < 1:1177897/76‑1 (MQ=255) gTAATTTNCGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACAGAGGAt > 2:977917/1‑76 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAANAACTACATTTTT‑T‑TTCATCGTGAAATGGACAGAGGAt < 1:1175287/76‑1 (MQ=255) aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACAGAGGATaa > 1:499991/1‑76 (MQ=255) aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACAGAGGATaa > 2:148500/1‑76 (MQ=255) tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACAGAGGATAAAt < 1:125103/76‑1 (MQ=255) ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACAGAGGATAAAtt > 2:1140397/1‑76 (MQ=255) aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCa < 2:123224/76‑1 (MQ=255) aCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCaa < 1:582460/76‑1 (MQ=255) gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa > 2:125238/1‑76 (MQ=255) aaTTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa < 1:834511/76‑1 (MQ=255) aTTTAAACCGCTTCGCCTTCAAAAAACTATATTTTT‑T‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc < 2:24534/76‑1 (MQ=255) tAAACCGCTTCACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa > 2:522340/1‑76 (MQ=255) tAAACCGCTTCACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa > 1:207657/1‑76 (MQ=255) aagccGCTTCACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt > 1:366883/4‑76 (MQ=255) aaaCCGCTTCACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt > 1:366858/1‑76 (MQ=255) aaaCCGCTTCACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt < 1:231566/76‑1 (MQ=255) aCCGCTTCACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTGAGATGGACAGAGGANAAATTGTCGANTAATGGATTa < 2:1194421/76‑1 (MQ=255) gCTTCACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcc > 1:421289/1‑76 (MQ=255) ttCACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc < 2:1078021/76‑1 (MQ=255) tCACCTTGAAAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg > 2:868835/1‑76 (MQ=255) tGAAAAAACTACATTTTT‑T‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg < 1:773112/76‑1 (MQ=255) aaaCTACATTTTT‑T‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc > 1:165009/1‑76 (MQ=255) | ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTTT‑T‑TTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGC > REL606/3875550‑3875692 |
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: ‑ |