Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,499,798 | T→G | T336P (ACC→CCC) | glgC ← | glucose‑1‑phosphate adenylyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,499,798 | 0 | T | G | 98.0% | 197.3 / ‑4.9 | 50 | T336P (ACC→CCC) | glgC | glucose‑1‑phosphate adenylyltransferase |
Reads supporting (aligned to +/- strand): ref base T (0/1); new base G (20/29); total (20/30) | |||||||||||
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.93e-01 |
CGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTG > REL606/3499725‑3499869 | cgcgAGAACAGACCGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCa < 2:1020770/76‑1 (MQ=255) cgcgAAAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCa < 1:1339077/76‑1 (MQ=255) gcgAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCAt > 1:376373/1‑76 (MQ=255) gcgAGAACAGAACGGACTGCACCACCACCGAACCAGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCAt > 2:283743/1‑76 (MQ=255) cgAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATc > 1:875334/1‑76 (MQ=255) agaacagaacGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATccc > 1:777970/1‑76 (MQ=255) acagaacGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTg < 1:244865/76‑1 (MQ=255) cagaacGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTgg > 2:1209154/1‑76 (MQ=255) aacGGACTGCCCCACCCCCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTa < 1:200752/76‑1 (MQ=255) aacGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTa < 2:459250/76‑1 (MQ=255) aacGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTa < 1:1671177/76‑1 (MQ=255) ggACTGCCCCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCg < 2:803014/76‑1 (MQ=255) cTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAg < 1:154513/76‑1 (MQ=255) cTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAg > 1:895468/1‑76 (MQ=255) cTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAg > 2:620236/1‑76 (MQ=255) tGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGc > 1:745534/1‑76 (MQ=255) gCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCg > 1:1894111/1‑76 (MQ=255) caccaccacCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGa > 1:970105/1‑76 (MQ=255) accaccaccGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGAt < 2:1249354/76‑1 (MQ=255) cccaccGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCt < 2:1504237/75‑1 (MQ=255) accaccGAATCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGCGGCTACCGGAGCGATccc > 1:1677857/1‑75 (MQ=255) accaccGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCt < 2:1663146/76‑1 (MQ=255) cccGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCa < 1:1415767/75‑1 (MQ=255) cGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACg > 1:1751925/1‑76 (MQ=255) cGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGGTCCTGCACGAAttt > 2:108341/1‑76 (MQ=255) cGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAAttt < 2:1296295/76‑1 (MQ=255) ggAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTc < 2:1340186/76‑1 (MQ=255) gagaTCACACAACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCg < 1:1539907/76‑1 (MQ=255) aTCACACAACCTCCCGCAACCAGTGAGTGAAAGGGCATCCCGCGGCTATCATTGAGATTCTGCACGAATTTCggcg > 2:5668/1‑73 (MQ=255) cacacaACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGc < 2:513073/76‑1 (MQ=255) acaACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGt > 1:1209684/1‑76 (MQ=255) acaACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGt < 1:155122/76‑1 (MQ=255) caACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTa > 1:1523666/1‑76 (MQ=255) caACCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTa < 1:1546688/76‑1 (MQ=255) accccgccgGAAACCAGTGAGTGAAGGGGCACCTCGTCGCTACCGGAGCGATCCTGCGCGAAGTTCCCTGGCGGTa < 2:227277/73‑1 (MQ=255) aaCCGCCGGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTaa > 1:260297/1‑76 (MQ=255) gccgGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGAt < 2:1417386/76‑1 (MQ=255) cgGAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGGAATGATTc > 2:443019/1‑76 (MQ=255) gAAACCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCAt > 1:1448402/1‑76 (MQ=255) aaaCCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCAtt > 2:915066/1‑76 (MQ=255) aCCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGt < 1:699442/76‑1 (MQ=255) aCCAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGt < 1:699436/76‑1 (MQ=255) cAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAg < 1:1028668/76‑1 (MQ=255) cAGTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAg < 1:859590/76‑1 (MQ=255) gTGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGt < 2:929936/76‑1 (MQ=255) tGAGTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTg < 1:130281/76‑1 (MQ=255) gTTAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTGcca > 1:1021458/1‑74 (MQ=255) ttAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGaa < 1:1209154/76‑1 (MQ=255) tAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAAt < 2:421177/76‑1 (MQ=255) tAAGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAAt < 2:971931/76‑1 (MQ=255) aGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTg < 1:1508944/76‑1 (MQ=255) aGGGGCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTg < 2:590673/76‑1 (MQ=255) | CGCGAGAACAGAACGGACTGCACCACCACCGAACCGGAGATCACACAACCGCCGGAAACCAGTGAGTTAAGGGTCATCCCGTGGCTACCGGAGCGATCCTGCACGAATTTCGCTGGCGGTAATGATTCATTGTAGGTGCGAATTG > REL606/3499725‑3499869 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |