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% | 203.2 / 0.7 | 57 | 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 (28/26); total (31/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.94e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01 |
ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGC > REL606/3875550‑3875694 | aTAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAACACTACa > 2:573110/1‑76 (MQ=255) tAGCATTATTCACCNGACGCGCCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t > 1:1712085/1‑76 (MQ=255) tAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t < 2:1200123/76‑1 (MQ=255) tAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAAACGCTTCACCTTGAAAAAACTACA‑t < 1:1574798/76‑1 (MQ=255) aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGGACGAATTTAAACTGCTTCACCTNGAAAAAACTACA‑tt < 1:4847/76‑1 (MQ=255) attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tttt‑t > 1:1275160/1‑76 (MQ=255) tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tttt‑tt‑t < 2:126535/76‑1 (MQ=255) attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tttt‑tt‑tt > 1:1680277/1‑75 (MQ=255) tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tttt‑tt‑ttca < 1:611077/76‑4 (MQ=255) cATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGt > 1:978109/1‑76 (MQ=255) cATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGt > 1:774680/1‑76 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTg < 1:1381916/76‑1 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTg < 1:1318840/76‑1 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTg < 1:787253/76‑1 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTg > 1:537033/1‑76 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTAACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTg > 2:619976/1‑76 (MQ=255) aTACGCGTCCTGACTGAAACTTGANAACGAATTTAAACCGCTTCACATTGAAAAAACTACATTTTT‑TT‑TCATCGTg > 1:153990/1‑76 (MQ=255) cgcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTGACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGaaa < 1:1495633/76‑1 (MQ=255) cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATg > 1:944017/1‑76 (MQ=255) gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATgg < 1:1176222/76‑1 (MQ=255) ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGAc > 2:972981/1‑76 (MQ=255) cTGACTGTAATTTGAGAATGANTTNAAACCGCTTCACCTTGAAAAANCTACACTNTN‑TT‑TCATCGTGAAATGGACa > 2:255285/1‑76 (MQ=255) cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACa < 1:521876/76‑1 (MQ=255) cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACa < 1:645392/76‑1 (MQ=255) cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACa < 1:676616/76‑1 (MQ=255) gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACaga > 2:1015191/1‑76 (MQ=255) aCTGTAATTTGAGAACGAATTTAAACCGCTTGACNTTGAGAANACNACN‑TTTTGTT‑TGATCGTGAAATGGACagag > 1:1533910/1‑76 (MQ=255) cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTTAACAAACTACANTTTT‑TT‑TCATCGTGAAATGGACAGAgg < 1:1553740/76‑1 (MQ=255) cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTTCATTTTT‑TT‑TCATCGTGAAATGAACAGAgg > 1:26997/1‑76 (MQ=255) cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAgg < 2:1273065/76‑1 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGANAGAGGAt < 2:1417916/76‑1 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGAt < 1:1418829/76‑1 (MQ=255) aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATaa > 1:819745/1‑76 (MQ=255) aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATaa < 2:456995/76‑1 (MQ=255) tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAAt > 1:602791/1‑76 (MQ=255) ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAAtt < 1:1642311/76‑1 (MQ=255) agaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGtt < 1:1640478/76‑1 (MQ=255) aCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCaa > 2:1357435/1‑76 (MQ=255) cGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAAt > 1:1362445/1‑76 (MQ=255) aaTTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa > 2:341393/1‑76 (MQ=255) aaTTTAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa > 2:341378/1‑76 (MQ=255) ntttAAACCGCTTCACCTTGAAAAAACTACA‑TTTT‑TTGTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc > 1:1479712/2‑76 (MQ=255) tttAAACCGCTTCACCTNGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGAGAGAGGATAAACTGGTCAATAACg > 1:220763/1‑76 (MQ=255) tAAACCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa > 2:189413/1‑76 (MQ=255) aaaCCGCTACACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt > 1:1638013/1‑76 (MQ=255) aaCCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt < 2:1163067/76‑1 (MQ=255) aaCCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt < 1:1236860/76‑1 (MQ=255) aaCCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt > 2:1290337/1‑76 (MQ=255) aaCCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt < 2:516946/76‑1 (MQ=255) aaCCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt < 2:356570/76‑1 (MQ=255) aCCGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa < 2:610036/76‑1 (MQ=255) aCCGCTTCAACTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa > 2:694902/1‑76 (MQ=255) ccGCTTCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt < 1:1611046/76‑1 (MQ=255) ttCACCTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc < 2:127074/76‑1 (MQ=255) ccTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg > 2:1396907/1‑76 (MQ=255) cTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGANAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa > 1:1723510/1‑76 (MQ=255) cTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa > 2:1384579/1‑76 (MQ=255) cTTGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa > 2:1367392/1‑76 (MQ=255) ttGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc < 2:1671486/76‑1 (MQ=255) ttGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc < 1:352184/76‑1 (MQ=255) ttGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc < 1:384596/76‑1 (MQ=255) tGAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg > 1:785655/1‑76 (MQ=255) gAAAAAACTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGa < 1:1410645/76‑1 (MQ=255) aaaaaaCTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa < 2:1469752/76‑1 (MQ=255) aaaCTACATTTTT‑TT‑TCATCGTGAAATNGACAGAGGTTATATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc > 1:1557244/1‑76 (MQ=255) aaaCTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc > 1:1286014/1‑76 (MQ=255) aCTACATTTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcgc > 1:1164302/1‑76 (MQ=255) | ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGC > REL606/3875550‑3875694 |
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: ‑ |