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.9% | 126.9 / 9.2 | 36 | 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/2); new base T (18/14); total (20/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑T‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC > REL606/3875560‑3875695 | tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑ttt‑tt‑ttca < 2:126946/76‑4 (MQ=255) aCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCg > 1:403520/1‑76 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTg < 1:804449/76‑1 (MQ=255) tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGa < 2:424327/76‑1 (MQ=255) aCGCGTCCTGACTGTAATCTGAGAACGNATNTCNACTGCTTNACATTGGAAAAACTACATT‑TTT‑TT‑TCATCGTGaa < 1:1098435/76‑1 (MQ=255) cgcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGaaa > 1:292232/1‑76 (MQ=255) gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATgg < 1:374055/76‑1 (MQ=255) gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACaga > 2:764626/1‑76 (MQ=255) aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACagag > 1:698634/1‑76 (MQ=255) cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAgg > 2:676946/1‑76 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGAt > 1:999709/1‑76 (MQ=255) aaTTTGAGAACGAAGTCAAACCGCTTCACCGTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATaa < 1:21481/76‑1 (MQ=255) tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAAt < 1:895961/76‑1 (MQ=255) tttGAGAAAGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑T‑TTTCTT‑TCATCGTGAAATCGACAGAGGATAAAt > 2:29250/1‑76 (MQ=255) ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAAtt > 2:1188371/1‑76 (MQ=255) agaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGtt < 1:416402/76‑1 (MQ=255) gaACGAATTTAANCCNCTTCACCTTGAAAAAACTNCA‑T‑TTTATT‑TCATCGTGAAATGGACAGAGGATAAATTGTTc < 1:1037010/76‑1 (MQ=255) aCGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCaa < 1:991582/76‑1 (MQ=255) cGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAAt > 2:258051/1‑76 (MQ=255) gAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATa > 2:902744/1‑76 (MQ=255) gAATTTAAACAGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATa > 1:104448/1‑76 (MQ=255) aaTTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa < 1:728932/76‑1 (MQ=255) aTTTAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc > 1:351083/1‑76 (MQ=255) tttAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACg > 1:1024422/1‑76 (MQ=255) tAAACCGCTTCANCTNGAAAAAACTACA‑T‑TTT‑TTNTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa < 1:1032019/76‑1 (MQ=255) tAAACCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa < 1:1051784/76‑1 (MQ=255) aaaCCGCTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt < 2:595213/76‑1 (MQ=255) cTTCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc > 2:489980/1‑76 (MQ=255) tCACCTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg > 1:818280/1‑76 (MQ=255) ccTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg < 1:914826/76‑1 (MQ=255) cTTGAAAAAACTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa > 1:771043/1‑76 (MQ=255) aaaaaaCTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa > 2:1129379/1‑76 (MQ=255) aaaaaaCTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa < 1:587794/76‑1 (MQ=255) aaaaCTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg > 2:149862/1‑76 (MQ=255) aaCTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg < 1:489058/76‑1 (MQ=255) aaCTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg < 2:842889/76‑1 (MQ=255) aaCTACA‑TCTTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCACGCTGACGAAAgcg > 2:27233/1‑76 (MQ=255) cTACATT‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc > 2:822552/1‑76 (MQ=255) | TCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑T‑TTT‑TT‑TCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC > REL606/3875560‑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: ‑ |