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% | 190.3 / 1.1 | 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 . (2/0); new base T (27/19); total (29/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.12e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CA‑TTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC > REL606/3875553‑3875695 | gCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CA‑ttt < 1:893146/76‑1 (MQ=255) cATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CA‑tttt > 2:482444/1‑76 (MQ=255) attattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CA‑ttttt < 2:388207/76‑1 (MQ=255) ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CA‑tttttt > 2:1366769/1‑76 (MQ=255) ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CA‑tttttt > 1:58065/1‑76 (MQ=255) tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CA‑ttttttt > 2:966160/1‑76 (MQ=255) tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CA‑ttttttt > 1:1758262/1‑76 (MQ=255) tattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CA‑ttttttt < 1:1958988/76‑1 (MQ=255) ttCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CA‑ttttttttc > 1:750742/1‑74 (MQ=255) tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CA‑ttttttttca > 1:1532667/1‑73 (MQ=255) aCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCg > 1:1112762/1‑76 (MQ=255) aCATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCg < 1:1232092/76‑1 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTg > 2:961872/1‑76 (MQ=255) tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGa > 2:758774/1‑76 (MQ=255) tACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGa > 2:665699/1‑76 (MQ=255) cgcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CA‑TTTTTTTCATCGTGAAAt > 1:1085497/1‑76 (MQ=255) gcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAAt < 2:1389155/76‑1 (MQ=255) tCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGa < 2:305124/76‑1 (MQ=255) cTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACa > 2:314625/1‑76 (MQ=255) tGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACag > 2:843282/1‑76 (MQ=255) aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACagag > 1:950339/1‑76 (MQ=255) tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑a < 2:1647420/76‑1 (MQ=255) gcnaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑At > 2:1498482/4‑76 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑At < 1:1051584/76‑1 (MQ=255) gTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑At > 2:627495/1‑76 (MQ=255) tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATa > 2:885372/1‑76 (MQ=255) aaTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATaa < 1:184322/76‑1 (MQ=255) aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATaaa < 1:140923/76‑1 (MQ=255) aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATaaa < 2:522919/76‑1 (MQ=255) aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATaaa > 1:1952013/1‑76 (MQ=255) tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAAt < 2:1635745/76‑1 (MQ=255) tttGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAAt > 2:1867952/1‑76 (MQ=255) tGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTg < 1:398463/76‑1 (MQ=255) tGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTg < 1:684712/76‑1 (MQ=255) agaACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGtt < 1:226894/76‑1 (MQ=255) aTTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAAc > 1:903138/1‑76 (MQ=255) aTTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAAc < 1:1380911/76‑1 (MQ=255) aTTTAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAAc > 1:807296/1‑76 (MQ=255) tttAAACCGCTTCACCTTGAAAAAACTACCT‑TTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACg > 2:1448845/1‑76 (MQ=255) ttAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACgg < 1:566937/76‑1 (MQ=255) ttAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACgg > 1:1107960/1‑76 (MQ=255) ttAAACAGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATGGGGAAAT‑GACCGAGGAATAATTTGTTCATTaagga > 2:346860/1‑73 (MQ=255) tAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACGGa > 1:690551/1‑76 (MQ=255) tAAACCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACGGa > 1:1326528/1‑76 (MQ=255) aaaCCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACGGAt < 1:1258034/76‑1 (MQ=255) aaCCGCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACGGAtt > 2:1656692/1‑76 (MQ=255) gCTTCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACGGATTATcc < 2:535266/76‑1 (MQ=255) tCACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACGGATTATCCCCg < 1:421340/76‑1 (MQ=255) cACCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACGGATTATCCCCGc > 1:134191/1‑76 (MQ=255) aCCTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACGGATTATCCCCGCt < 2:911766/76‑1 (MQ=255) cTTGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACGGATTATCCCCGCTGa > 1:425713/1‑76 (MQ=255) ttGAAAAAACTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACGGATTATCCCCGCTGAc > 2:2004964/1‑76 (MQ=255) aaaaaaCTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACGGATTATCCCCGCTGACGaa < 1:187986/76‑1 (MQ=255) aaaaaCTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa > 2:563116/1‑76 (MQ=255) aaaaaCTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACGGATTATCCCCGCTGACGaaa > 2:1147124/1‑76 (MQ=255) aaaCTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc > 2:1066563/1‑76 (MQ=255) aaCTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg > 2:366154/1‑76 (MQ=255) cTA‑CATTTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGcc > 1:892522/1‑76 (MQ=255) | GCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTA‑CA‑TTTTTTTCATCGTGAAATGGACAGAGG‑ATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCGCC > REL606/3875553‑3875695 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |