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 | 88.5% | 84.4 / 6.7 | 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 . (1/2); new base T (12/11); total (13/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01 |
GCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC > REL606/3875553‑3875695 | gCATTATTCAGCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt < 2:8371/76‑1 (MQ=255) ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tttttt > 2:433350/1‑76 (MQ=255) attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tttttt‑tt > 1:3937/1‑75 (MQ=255) ttCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tttttt‑ttc < 1:629987/76‑3 (MQ=255) cAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTTTTT‑TCATc < 2:402827/76‑1 (MQ=255) aCATACGCGTCCCGTCTCTAATTGNAGAACGAATTTAAACCGCTTCACCTTGAAACAANTACATTTTTTT‑TCATCg < 2:655227/76‑1 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTg > 1:308010/1‑76 (MQ=255) cgcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGaaa < 1:520278/76‑1 (MQ=255) gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATgg < 1:143576/76‑1 (MQ=255) cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACa < 1:182962/76‑1 (MQ=255) gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAACAACTACATTTTTTT‑TCATCGTGAAATGGACaga > 2:644559/1‑76 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGAt > 2:595772/1‑76 (MQ=255) tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATa < 1:13486/76‑1 (MQ=255) tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAAt < 2:353580/76‑1 (MQ=255) gagaACGAATTTAAACCGCTTCACNTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGt > 2:134522/1‑76 (MQ=255) agaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGtt > 2:530467/1‑76 (MQ=255) agaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGtt > 1:225827/1‑76 (MQ=255) cGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAAt > 2:339211/1‑76 (MQ=255) gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAAATGTCCAATa < 2:87008/76‑1 (MQ=255) aTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc > 1:412346/1‑76 (MQ=255) aTTTAAACCGCTTCACCTTGAAAAAACTACA‑TTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg > 2:208188/1‑76 (MQ=255) ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg > 1:646766/1‑76 (MQ=255) tAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa < 2:539509/76‑1 (MQ=255) aaaCCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt < 1:311084/76‑1 (MQ=255) tCACCTAGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg > 1:121236/1‑76 (MQ=255) ccTTGAGAAAACTACA‑TTTTTTNTCATCGTGAANTGGACAGANGATAAATTGTTCAATAACTGATNATCCCCGCTg < 2:686044/76‑1 (MQ=255) aaaaaCTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCACTGACGaaa < 2:397560/76‑1 (MQ=255) aaaCTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc > 1:150936/1‑76 (MQ=255) aaCTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg > 2:21511/1‑76 (MQ=255) cTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc < 1:45330/76‑1 (MQ=255) | GCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC > REL606/3875553‑3875695 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |