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.8% | 184.5 / 0.6 | 48 | 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/1); new base T (28/18); total (29/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01 |
ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑T‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCG > REL606/3875550‑3875693 | aTAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACa > 1:857707/1‑76 (MQ=255) tAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t < 1:680762/76‑1 (MQ=255) aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t < 1:650347/76‑1 (MQ=255) aGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t < 1:1275355/76‑1 (MQ=255) gCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t‑t < 2:749808/76‑1 (MQ=255) cATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t‑tt > 1:1494914/1‑76 (MQ=255) attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t‑ttt < 1:1759145/76‑1 (MQ=255) attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑t‑t‑ttt < 1:1044226/76‑1 (MQ=255) cATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGt > 2:1521626/1‑76 (MQ=255) cATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGt > 1:222149/1‑76 (MQ=255) atccgcgTCCTGACTGAAATTTGAGAACGAATTTAAGCCGCCTCACCGTGAAATAGCTACA‑T‑TCTGTTTCATCGTg < 2:807129/73‑1 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTg > 2:799430/1‑76 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTg < 2:807117/76‑1 (MQ=255) aCGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGaa > 1:149528/1‑76 (MQ=255) gcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAAt < 1:538122/76‑1 (MQ=255) cgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATg > 2:223038/1‑76 (MQ=255) gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATgg < 1:657163/76‑1 (MQ=255) gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATgg > 1:1292216/1‑76 (MQ=255) ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGAc < 1:1154799/76‑1 (MQ=255) ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGAc < 2:1346378/76‑1 (MQ=255) tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACag < 1:1085560/76‑1 (MQ=255) tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACag < 1:1085568/76‑1 (MQ=255) tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACag < 1:1454788/76‑1 (MQ=255) gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACaga > 2:1518534/1‑76 (MQ=255) aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATaaa > 1:239723/1‑76 (MQ=255) ttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAAtt > 1:146212/1‑76 (MQ=255) tGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTg < 1:717160/76‑1 (MQ=255) gagaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGt > 2:1516070/1‑76 (MQ=255) agaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGtt < 2:1096477/76‑1 (MQ=255) gaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTc < 2:190134/76‑1 (MQ=255) gaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTc < 1:66089/76‑1 (MQ=255) aCGAATTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCaa < 1:1096813/76‑1 (MQ=255) aTTTAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc > 1:1556444/1‑76 (MQ=255) tAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa > 2:628574/1‑76 (MQ=255) tAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa > 2:1203998/1‑76 (MQ=255) tAAACCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa > 2:675721/1‑76 (MQ=255) aaaCCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt > 1:1312288/1‑76 (MQ=255) aaaCCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAt > 1:146594/1‑76 (MQ=255) aaCCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt < 1:1500253/76‑1 (MQ=255) aCCGCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa < 2:460716/76‑1 (MQ=255) gCTTCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcc < 1:22740/76‑1 (MQ=255) ttCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc > 1:1269421/1‑76 (MQ=255) tCACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCg > 2:1524656/1‑76 (MQ=255) cACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGc > 1:1749256/1‑76 (MQ=255) cACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGc > 2:161661/1‑76 (MQ=255) cACCTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGc > 1:408238/1‑76 (MQ=255) ccTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg < 2:1268995/76‑1 (MQ=255) ccTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTg > 1:1533439/1‑76 (MQ=255) cTTGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa < 2:1310449/76‑1 (MQ=255) ttGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGAc > 1:197005/1‑76 (MQ=255) tGAAAAAACTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACg > 1:1305621/1‑76 (MQ=255) aaaaCTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAg > 1:806483/1‑76 (MQ=255) aaaCTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc > 2:124807/1‑76 (MQ=255) aaaCTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc > 2:457542/1‑76 (MQ=255) aaaCTACATT‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc > 2:1581421/1‑76 (MQ=255) aaCTACA‑TCT‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg > 1:15739/1‑76 (MQ=255) | ATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑T‑T‑TTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCG > REL606/3875550‑3875693 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |