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 | 97.6% | 155.6 / ‑5.6 | 43 | 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/0); major base T (29/12); minor base C (1/0); total (31/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.95e-01 |
AGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGC > REL606/3875552‑3875694 | aGCCTTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tt > 1:1247215/1‑76 (MQ=255) ttattCCACATTCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tttttt > 1:1223714/1‑76 (MQ=255) ttCAANACACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAANCTACA‑tttttt‑ttc > 2:1341544/1‑74 (MQ=255) ttCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tttttt‑ttc > 2:601287/1‑74 (MQ=255) cATACGCGTCCTNACNGTAATTTGAGAACGAAGTNAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGt < 1:1086346/76‑1 (MQ=255) cATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTATATTTTTTT‑TCATCGt > 1:1102036/1‑76 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTg < 1:103819/76‑1 (MQ=255) tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGa < 2:1166165/76‑1 (MQ=255) aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGaa > 2:1148752/1‑76 (MQ=255) cgcgTCCTGACTGTACTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGaaa > 1:122798/1‑76 (MQ=255) cgcgTCCTGACTGTAATTTGACAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGaaa < 1:647802/76‑1 (MQ=255) gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATgg > 2:831647/1‑76 (MQ=255) ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGAc > 1:914957/1‑76 (MQ=255) aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACagag < 1:981964/76‑1 (MQ=255) tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGa > 1:301710/1‑76 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGAt > 1:1038978/1‑76 (MQ=255) aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATaa > 1:407002/1‑76 (MQ=255) aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATaaa > 1:584594/1‑76 (MQ=255) aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATaaa > 1:842074/1‑76 (MQ=255) ttGAGAACGAATTTACACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAAtt > 1:127038/1‑76 (MQ=255) ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAAtt > 1:628861/1‑76 (MQ=255) gagaACGAATTTAAACCGCTTCACCTTNAAAAAACTACACTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGt > 1:1167924/1‑76 (MQ=255) gagaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGt > 2:690242/1‑76 (MQ=255) tttAAACCGCTTCACCTTGAAAAAACTNCA‑TTTTTTGTCTTCGTGAAATGGACAGAGGATAAATTGTTCAATAACg > 1:1181535/1‑76 (MQ=255) ttAAACCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg > 1:157552/1‑76 (MQ=255) aaaCCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt < 2:834700/76‑1 (MQ=255) aCCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa > 2:407249/1‑76 (MQ=255) aCCGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa < 1:941432/76‑1 (MQ=255) ccGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt > 1:904652/1‑76 (MQ=255) cGCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATc < 1:298206/76‑1 (MQ=255) gCTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcc > 2:1043063/1‑76 (MQ=255) cTTCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc > 1:783345/1‑76 (MQ=255) ttCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc < 2:1106814/76‑1 (MQ=255) ttCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc > 2:863420/1‑76 (MQ=255) tCACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg > 2:550346/1‑76 (MQ=255) cACCTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGc > 1:837016/1‑76 (MQ=255) cTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa > 1:823936/1‑76 (MQ=255) cTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa > 2:174017/1‑76 (MQ=255) cTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa > 1:422150/1‑76 (MQ=255) cTTGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa > 1:1109303/1‑76 (MQ=255) ttGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGGTAAATTGTTCAATAACGGATTGTCCCCGCTGAc < 1:1228834/76‑1 (MQ=255) ttGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc > 1:652706/1‑76 (MQ=255) tGAAAAAACTACATTTTTTT‑TCATCGTGAAATGGACAGAGGCTAAATTGTTCAATAACGGATTATCCCCGCTGACg < 2:92634/76‑1 (MQ=255) aaaaCTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg > 1:815078/1‑76 (MQ=255) aaaNTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc > 1:1208230/1‑76 (MQ=255) aaCTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg < 2:933049/76‑1 (MQ=255) aCTACATTTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcgc > 1:500578/1‑76 (MQ=255) | AGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTTTTT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGC > REL606/3875552‑3875694 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |