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.6% | 83.2 / 0.4 | 27 | 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 . (0/2); new base T (19/6); total (19/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.98e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
ATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGC > REL606/3875555‑3875694 | attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttttt > 1:156749/1‑76 (MQ=255) cATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGt > 1:503634/1‑76 (MQ=255) tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGa > 1:105625/1‑76 (MQ=255) tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGa > 1:304172/1‑76 (MQ=255) tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGa > 1:406170/1‑76 (MQ=255) tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGa > 1:105619/1‑76 (MQ=255) tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGa > 1:469611/1‑76 (MQ=255) aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaa < 1:95406/76‑1 (MQ=255) cgcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGaaa > 2:362556/1‑76 (MQ=255) tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATa < 2:347648/76‑1 (MQ=255) tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATa < 2:511975/76‑1 (MQ=255) aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa > 2:213592/1‑76 (MQ=255) aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa > 1:411114/1‑76 (MQ=255) aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATaa > 1:38996/1‑76 (MQ=255) gagaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGt > 2:145353/1‑76 (MQ=255) gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa > 1:107005/1‑76 (MQ=255) aaTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa > 2:69083/1‑76 (MQ=255) aaTTTAAACCGCCTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa < 1:176366/76‑1 (MQ=255) aTTTAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc > 1:611754/1‑76 (MQ=255) tAAACCGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa > 1:411321/1‑76 (MQ=255) ccGCTTCACCTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt < 1:196603/76‑1 (MQ=255) cTTGAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa > 1:424693/1‑76 (MQ=255) cTTGAAAAAACTACA‑TTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc < 1:492350/76‑1 (MQ=255) cTTGAAAAAACTACA‑TTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc < 1:492353/76‑1 (MQ=255) gAAAAAACTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGa < 1:302795/76‑1 (MQ=255) aaaaaCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa > 2:738700/1‑76 (MQ=255) aCTACATTTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcgc > 1:91125/1‑76 (MQ=255) aCTACATTTTTTTTCATCGTGAAATGGACAGAGGACAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcgc > 2:184857/1‑76 (MQ=255) | ATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTTTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGC > REL606/3875555‑3875694 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |