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 | 94.1% | 182.2 / 4.1 | 51 | 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 . (3/0); new base T (21/27); total (24/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.72e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01 |
ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑T‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCA > REL606/3875550‑3875696 | aTAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACa > 1:1608652/1‑76 (MQ=255) aTAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACa < 1:1854733/76‑1 (MQ=255) aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t < 1:1488086/76‑1 (MQ=255) aGCATTATTCAACANACGCGTCCTTACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t > 1:11526/1‑76 (MQ=255) gCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t‑t > 1:1149784/1‑76 (MQ=255) cATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t‑tt < 2:866785/76‑1 (MQ=255) attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t‑ttt < 1:1818754/76‑1 (MQ=255) ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t‑ttt‑t < 1:738137/76‑1 (MQ=255) aCATAAGCGTCCTGACTGTAATTTGAAAACGAATTTAAACCGCTTCCCCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCg < 1:131522/76‑1 (MQ=255) aTACGCTTCCTGTCTTTAATTTGAGACCGAACCCAAACCGCTTCACCTTTTACCCACTCCATT‑T‑TTT‑TT‑CATCGtc < 2:1394768/76‑2 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTg < 1:871192/76‑1 (MQ=255) tACGCGTCCTGACTGTAATTTGAGAACNAATTTAANNCGCTTGACCTTGAACACACTCCA‑T‑T‑TTTCTT‑AATCGTGa > 1:2005225/1‑76 (MQ=255) tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGa < 1:1391436/76‑1 (MQ=255) cgcgTCCTGGCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGaaa > 1:153531/1‑76 (MQ=255) cgcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGaaa < 1:505080/76‑1 (MQ=255) gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATgg < 2:1801739/76‑1 (MQ=255) gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATgg < 2:1513553/76‑1 (MQ=255) ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGAc < 1:1038043/76‑1 (MQ=255) ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGAc > 1:1846195/1‑76 (MQ=255) ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGAc < 1:548071/76‑1 (MQ=255) ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑T‑T‑TTT‑TT‑CATCGTGAAATGGACa > 2:1325431/1‑76 (MQ=255) cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACa > 1:1406034/1‑76 (MQ=255) cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACa < 1:718824/76‑1 (MQ=255) tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACag < 2:1831079/76‑1 (MQ=255) tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACag < 1:961292/76‑1 (MQ=255) tGACTGGATTTTGAGAACTAATTAAACCCGCTTCACCTTGAAAAACCTACATT‑T‑TTT‑TT‑CATCGGGAAATGGACag < 2:749862/76‑1 (MQ=255) gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACaga > 2:1191510/1‑76 (MQ=255) cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAgg < 1:1728299/76‑1 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGAt < 1:693492/76‑1 (MQ=255) aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCATGAAATGGACAGAGGATaa < 2:828117/76‑1 (MQ=255) ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAAtt > 1:998800/1‑76 (MQ=255) ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAAtt > 1:1400266/1‑76 (MQ=255) tGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTg > 2:1944062/1‑76 (MQ=255) tGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTg < 1:566796/76‑1 (MQ=255) gaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTc < 1:1548505/76‑1 (MQ=255) gaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTc < 1:1548513/76‑1 (MQ=255) aCGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCaa < 1:1010579/76‑1 (MQ=255) cGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAAt < 1:647828/76‑1 (MQ=255) aaTTTAAACCGCTTCACCTTGAAACAACNACA‑T‑T‑TTT‑TTGCATCGTGACANGGACAGAGGATAAATTGTTCAATaa > 1:292642/1‑76 (MQ=255) tttAAACCGCTTCACCTTGAAAAACCTNCA‑T‑T‑TTTNNT‑CATCCTGAAATGGACAGAGCATAAATTGTTCAATAACg > 1:2008070/1‑76 (MQ=255) ttAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg < 1:1635011/76‑1 (MQ=255) tAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAAT‑GACAGAGGATAAATTGTTCAATAACGGAt > 1:486571/1‑76 (MQ=255) aaaCCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt > 2:1785935/1‑76 (MQ=255) aaaCCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt > 2:145035/1‑76 (MQ=255) aCCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa > 2:1773448/1‑76 (MQ=255) aCCGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa < 2:1563240/76‑1 (MQ=255) ccGCTTCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt > 2:589623/1‑76 (MQ=255) ttCACCTTGAAAAAACTAGA‑TGT‑TCT‑GT‑CATCGTGTAATGGGCAGAGGATAAATAGTTCAGTGACGGATTCTcccc > 1:2041323/1‑76 (MQ=255) ttCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc < 1:510243/76‑1 (MQ=255) ttCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc > 1:1297832/1‑76 (MQ=255) tCACCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTATTCAATAACGGATTATCCCCg > 2:189406/1‑76 (MQ=255) aCCTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCt < 2:593241/76‑1 (MQ=255) cTTGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa < 1:1470636/76‑1 (MQ=255) ttGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc > 2:965736/1‑76 (MQ=255) ttGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc > 1:1929827/1‑76 (MQ=255) tGAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg > 1:807002/1‑76 (MQ=255) gAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACgn < 1:575597/76‑2 (MQ=255) gAAAAAACTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGa > 1:1240206/1‑76 (MQ=255) aaaaaCTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGCTAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa < 1:1699573/76‑1 (MQ=255) aaaaaCTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa > 1:928946/1‑76 (MQ=255) aaaaCTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTCTCCCCGCTGACGAAAg > 2:115815/1‑76 (MQ=255) aaaaCTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg > 1:514370/1‑76 (MQ=255) aaaCTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc < 1:990185/76‑1 (MQ=255) cTACATT‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc < 1:763360/76‑1 (MQ=255) tACA‑T‑TGTTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGTCGAAAGCGCCa < 1:63841/76‑1 (MQ=255) | ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑T‑T‑TTT‑TT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCCA > REL606/3875550‑3875696 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |