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 | 95.0% | 227.7 / 4.5 | 60 | 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/1); new base T (25/32); total (27/33) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.83e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01 |
ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC > REL606/3875550‑3875695 | aTAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACa > 1:728538/1‑76 (MQ=255) tAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t > 1:1002485/1‑76 (MQ=255) tAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t < 2:1824018/76‑1 (MQ=255) ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt‑ttt > 2:811302/1‑76 (MQ=255) attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt‑ttttt > 2:105826/1‑75 (MQ=255) ttCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt‑tttttc < 1:387553/76‑3 (MQ=255) aCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCg < 1:1192170/76‑1 (MQ=255) aCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTT‑TTTTCATCGt > 2:43304/1‑76 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTg > 2:1899778/1‑76 (MQ=255) tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGa < 2:1012412/76‑1 (MQ=255) cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATg > 2:17967/1‑76 (MQ=255) gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATgg < 1:1262908/76‑1 (MQ=255) tCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGa > 1:1408272/1‑76 (MQ=255) tCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGa > 1:149309/1‑76 (MQ=255) tCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGa > 1:1408257/1‑76 (MQ=255) tCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGa < 1:123656/76‑1 (MQ=255) ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGAc < 1:1754049/76‑1 (MQ=255) tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACag > 2:231160/1‑76 (MQ=255) gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACaga > 1:1643577/1‑76 (MQ=255) cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAgg < 1:815348/76‑1 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGAt < 1:1970925/76‑1 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGAt < 1:272809/76‑1 (MQ=255) tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATa > 2:1967294/1‑76 (MQ=255) tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATa > 1:1314890/1‑76 (MQ=255) tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATa > 1:243767/1‑76 (MQ=255) tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAAt < 1:920534/76‑1 (MQ=255) tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAAt > 2:96093/1‑76 (MQ=255) tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAAt > 1:1069198/1‑76 (MQ=255) tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAAt < 1:1261499/76‑1 (MQ=255) gagaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGt < 1:1154338/76‑1 (MQ=255) agaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGtt < 1:535776/76‑1 (MQ=255) aaCGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCa < 1:641539/76‑1 (MQ=255) cGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAAt < 2:1502813/76‑1 (MQ=255) cGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAAt < 1:942865/76‑1 (MQ=255) gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa < 2:112906/76‑1 (MQ=255) gAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa > 2:923160/1‑76 (MQ=255) aaTTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa < 1:169089/76‑1 (MQ=255) aaTTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa > 1:1710156/1‑76 (MQ=255) aaTTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa > 2:1031769/1‑76 (MQ=255) ttAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg > 2:1947227/1‑76 (MQ=255) tAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa > 1:346644/1‑76 (MQ=255) aaCCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt > 2:249571/1‑76 (MQ=255) aCCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa > 1:566304/1‑76 (MQ=255) aCCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa < 2:1858978/76‑1 (MQ=255) aCCGCTTCACCTTGAAAAAACTACA‑TTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt < 2:461300/76‑1 (MQ=255) ccGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt < 2:1798808/76‑1 (MQ=255) ccGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt < 2:284755/76‑1 (MQ=255) gCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcc < 1:542897/76‑1 (MQ=255) cTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc > 1:1527841/1‑76 (MQ=255) ttCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc < 1:1649595/76‑1 (MQ=255) ttCACCTTGAAAAAACTACA‑TTTATTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc > 2:128166/1‑76 (MQ=255) tCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg < 1:217564/76‑1 (MQ=255) cACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGc < 2:100222/76‑1 (MQ=255) ccTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg > 2:1714190/1‑76 (MQ=255) ccTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg < 1:1369820/76‑1 (MQ=255) cTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa > 1:224916/1‑76 (MQ=255) tGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg < 1:315828/76‑1 (MQ=255) tGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg < 1:1309007/76‑1 (MQ=255) tGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg < 1:195069/76‑1 (MQ=255) tGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg < 1:84495/76‑1 (MQ=255) aaaaaCTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa < 1:984195/76‑1 (MQ=255) aaaaCTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg < 2:334517/76‑1 (MQ=255) aaaCTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc > 2:605663/1‑76 (MQ=255) aaCTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg > 1:1988611/1‑76 (MQ=255) aCTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcgc < 2:85542/76‑1 (MQ=255) cTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc > 1:1093215/1‑76 (MQ=255) | ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC > REL606/3875550‑3875695 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |