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.9% | 150.9 / 3.7 | 42 | 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/1); new base T (22/17); total (24/18) | |||||||||||
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.96e-01 |
ATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑T‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC > REL606/3875555‑3875695 | attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑ttt‑t > 2:695731/1‑76 (MQ=255) tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑ttt‑tt‑t > 1:665300/1‑76 (MQ=255) tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAAATACA‑t‑ttt‑tt‑t > 1:665316/1‑76 (MQ=255) attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑ttt‑tt‑tt > 2:1066384/1‑75 (MQ=255) attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑ttt‑tt‑tt < 2:1140513/76‑2 (MQ=255) attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑ttt‑tt‑tt < 2:1304928/76‑2 (MQ=255) ttCAACATACGCGTCCTGACTGTAAGTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑ttt‑tt‑ttc < 1:596233/76‑3 (MQ=255) aCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCg < 1:1499835/76‑1 (MQ=255) cATACGCGTCCTGCCTGTAATTTGAGAACGAATNTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGt > 2:1373734/1‑76 (MQ=255) gTCCTGACTGTAATTTGAGNACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATgg > 1:1513293/1‑76 (MQ=255) gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATgg < 1:1212249/76‑1 (MQ=255) cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCTTGAAATGGACa < 2:321300/76‑1 (MQ=255) cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAgg > 2:1431606/1‑76 (MQ=255) tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGa < 2:1090871/76‑1 (MQ=255) tGTAATTTGAGAACGAATTTAAACCGCNTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGa > 1:1613193/1‑76 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAANTGGACAGAGGAt < 1:13958/76‑1 (MQ=255) aaTTTGAGAACGAATTTAAACCGCTTCACCTTGATATAACTACA‑T‑TTN‑TTCTCATCGTGAAATGGACAGAGGATaa > 1:133617/1‑76 (MQ=255) aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATaaa > 2:1037550/1‑76 (MQ=255) aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATaaa > 2:1037559/1‑76 (MQ=255) ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAAtt > 2:38740/1‑76 (MQ=255) ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAAtt < 2:735432/76‑1 (MQ=255) agaACGAATTTAAACCGATTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGACGATACATTNtt > 2:1526374/1‑76 (MQ=255) gaACGAATTTAAACCGCGTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATNGACAGAGGATAAATTGTTc < 1:1512494/76‑1 (MQ=255) aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCa < 2:1411575/76‑1 (MQ=255) aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCa > 2:402945/1‑76 (MQ=255) aCGAATTTAAACCGCTTCACCTTGAAAAATCTACA‑TCTTT‑TT‑TCATCGTGAAATGGANAGAGGATAAATTGTTCaa < 1:1541487/76‑1 (MQ=255) gAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATa < 1:864369/76‑1 (MQ=255) aaTTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa < 1:1084663/76‑1 (MQ=255) tttAAACCGCTTCACCTTGAAAAAACTATA‑T‑TTTCTA‑CCATCGTGAANTGCACAGAGGATAAATTGTTCAATAACg > 1:138429/1‑76 (MQ=255) tttAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg < 1:964918/76‑1 (MQ=255) ttAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGAGAGAGGATAAANTGTTCAATAACgc > 1:207181/1‑75 (MQ=255) tAAACCGCGTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTAAATAACGGa > 1:54665/1‑76 (MQ=255) aaaCCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAACTGGACAGAGGATAAATTGTTCAATAACGGAt > 2:263263/1‑76 (MQ=255) aaaCCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt > 2:263291/1‑76 (MQ=255) aaaCCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt < 1:628275/76‑1 (MQ=255) ccGCTTCACCTTGANAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTan < 2:217421/76‑2 (MQ=255) gCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcc < 1:923074/76‑1 (MQ=255) tCACATTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg > 1:97109/1‑76 (MQ=255) aCCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGNGGATAAATTGTTCAATAACGGATTATCCCCGCt < 2:1365268/76‑1 (MQ=255) cTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa > 1:435408/1‑76 (MQ=255) cTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa < 2:936681/76‑1 (MQ=255) ttGAAAAAACTACATT‑TTT‑TT‑TCATGGTGAAATGGACAGANNATAAATTGTTNNACAACGGATTATCCCCGCTGAc > 1:1449999/1‑76 (MQ=255) ttGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc < 1:503811/76‑1 (MQ=255) ttGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc > 1:1496761/1‑76 (MQ=255) tGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATNGATAGAGGNTAAATAGTTCAATAACGGANTATCCCCGCTGACg > 2:1334088/1‑76 (MQ=255) aaaaaaCTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa > 2:91587/1‑76 (MQ=255) aaaaaCTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa > 1:122871/1‑76 (MQ=255) aaaCTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAGGAAAgc > 1:1470669/1‑76 (MQ=255) cTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc > 2:1444877/1‑76 (MQ=255) | ATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑T‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC > REL606/3875555‑3875695 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |