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.3% | 133.7 / ‑2.4 | 37 | 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); new base T (17/19); total (18/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.86e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01 |
AGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑T‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC > REL606/3875552‑3875695 | aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGNATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t < 2:1227079/76‑1 (MQ=255) attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t‑ttt > 2:261123/1‑76 (MQ=255) ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t‑tttt < 2:1263110/76‑1 (MQ=255) ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t‑tttt > 1:1065110/1‑76 (MQ=255) tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t‑ttttt > 1:695903/1‑76 (MQ=255) attCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t‑tttttt < 2:270723/76‑2 (MQ=255) ttCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t‑ttttttc < 2:157174/76‑3 (MQ=255) tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t‑ttttttca > 1:556960/1‑73 (MQ=255) aCATACGGGTCCTGACTGTAATTTGAGCACGAATTTAAACCGCTTCACCTTGTAAAAACTACATT‑T‑TTTTT‑CATCg < 1:72602/76‑1 (MQ=255) aCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCg > 2:1152083/1‑76 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGACCGAATTTAAACCGCTTCACCTTGAAAACACGACATT‑T‑TTTTT‑CATCGTg > 1:1364069/1‑76 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTg > 2:1383890/1‑76 (MQ=255) aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGaa < 1:646719/76‑1 (MQ=255) gcgTCCTGACTGTACTTTGAGAACTAATGTANACCGCTTCTCCTTGNAAAAANTACATT‑T‑TTTTT‑CGTCGTGAAAt < 2:185897/76‑1 (MQ=255) tCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGa < 2:696644/76‑1 (MQ=255) ccTGACTGTAATTTGAGAANGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGAc > 1:1339745/1‑76 (MQ=255) ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGAc < 1:285329/76‑1 (MQ=255) ccTGACTGTAATTTGAGAACGAATTGAANCCGCTTCACTGTCAAAAAGCTACATT‑T‑TTTTT‑CATCGTGAAATGGAc < 1:133924/76‑1 (MQ=255) cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACa > 1:714912/1‑76 (MQ=255) cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACa > 1:714937/1‑76 (MQ=255) aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACagag > 2:604627/1‑76 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTATA‑TGA‑TATTT‑CATCGTGAAATGGACAGAGCAt > 1:1306483/1‑76 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGAt > 2:749597/1‑76 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑T‑TGTATTT‑CATAGTGCAATGGACAGAGGAt > 1:1306480/1‑76 (MQ=255) aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATaa < 2:470503/76‑1 (MQ=255) aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATaaa < 2:330071/76‑1 (MQ=255) tttGAGAACGAATTTAAAACGATTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATAAAt < 1:1354195/76‑1 (MQ=255) cGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAAt < 1:1279058/76‑1 (MQ=255) gAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATa < 2:176362/76‑1 (MQ=255) gAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATa < 2:929420/76‑1 (MQ=255) gAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATa < 1:370198/76‑1 (MQ=255) aaTTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATaa > 2:556828/1‑76 (MQ=255) ttAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACgg < 1:763514/76‑1 (MQ=255) tAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa < 1:806866/76‑1 (MQ=255) aCCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa < 1:758786/76‑1 (MQ=255) ccGCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CNTCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTAt > 2:1216847/1‑76 (MQ=255) gCTTCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcc < 2:948004/76‑1 (MQ=255) gCTTCACCATGAAAAAACTACATT‑T‑TTTTT‑CATCGTGATATGGAAATACGCTAAATTATCCCATAACAGATTATcc > 1:1308933/1‑76 (MQ=255) ttCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc > 2:1038384/1‑76 (MQ=255) ttCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc > 1:611039/1‑76 (MQ=255) tCACCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg > 2:394051/1‑76 (MQ=255) aCCTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCt > 1:587016/1‑76 (MQ=255) ccTTGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg > 2:618506/1‑76 (MQ=255) ttGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc < 2:18208/76‑1 (MQ=255) tGAAAAAACTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATAAATTGGTCAATAACGGATTATCCCCGCTGACg < 1:1360176/76‑1 (MQ=255) aCTACATT‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcgc > 1:959521/1‑76 (MQ=255) cTACA‑T‑T‑TTTTTNCATCGTGTAGTGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc > 1:20115/1‑76 (MQ=255) | AGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑T‑T‑TTTTT‑CATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC > REL606/3875552‑3875695 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |