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.6% | 160.7 / ‑3.5 | 42 | 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 . (0/1); new base T (18/23); total (18/24) | |||||||||||
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.95e-01 |
TAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCG > REL606/3875551‑3875693 | tAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACa < 1:247999/75‑1 (MQ=255) tAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACa < 1:1308110/75‑1 (MQ=255) tAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACa > 1:1530061/1‑75 (MQ=255) gCATTATTTAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tt > 1:1592606/1‑75 (MQ=255) gCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑tt > 1:1592618/1‑75 (MQ=255) cATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt > 1:2071962/1‑75 (MQ=255) ttattCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt‑tt > 1:2067381/1‑75 (MQ=255) ttCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt‑ttttt > 1:990556/1‑74 (MQ=255) ttCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt‑ttttt < 1:434809/75‑2 (MQ=255) ttCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt‑ttttt > 1:613805/1‑74 (MQ=255) tCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑ttt‑tttttc < 2:1868850/75‑3 (MQ=255) aTACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGt > 1:278829/1‑75 (MQ=255) gcgTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGaaa > 2:1192452/1‑75 (MQ=255) gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATg > 2:978149/1‑75 (MQ=255) gTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATg < 2:1828996/75‑1 (MQ=255) gTCCTGACTGTAATTCGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATg < 1:1240683/75‑1 (MQ=255) tCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATgg < 1:1211587/75‑1 (MQ=255) ccTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGa < 1:1678376/75‑1 (MQ=255) gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACag < 2:1889319/75‑1 (MQ=255) gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACag < 2:600160/75‑1 (MQ=255) gACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACag > 1:965135/1‑75 (MQ=255) aCTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACaga > 2:1464609/1‑75 (MQ=255) cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACagag < 2:1399217/75‑1 (MQ=255) cTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACagag < 2:1399233/75‑1 (MQ=255) cTGTAAGTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACagag < 2:807381/75‑1 (MQ=255) tGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAgg < 1:1397268/75‑1 (MQ=255) tAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGAt < 1:307935/75‑1 (MQ=255) aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATaa > 1:236881/1‑75 (MQ=255) aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATaa > 1:236876/1‑75 (MQ=255) aTTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATaa > 1:207223/1‑75 (MQ=255) gagaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTg > 1:1860153/1‑75 (MQ=255) gaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGtt > 1:259249/1‑75 (MQ=255) gaACGAATTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGtt > 2:731574/1‑75 (MQ=255) aaTTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATa < 2:1530061/75‑1 (MQ=255) aTTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa < 1:460296/75‑1 (MQ=255) aTTTAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATaa < 1:1017849/75‑1 (MQ=255) tttAAACCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAAc < 2:751931/75‑1 (MQ=255) aaaCCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGa < 2:1898526/75‑1 (MQ=255) aCCGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGAtt > 2:687787/1‑75 (MQ=255) ccGCTTCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTa < 2:2150659/75‑1 (MQ=255) ttCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATccc > 2:1826089/1‑75 (MQ=255) tCACCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATcccc < 2:888178/75‑1 (MQ=255) aCCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGATGATAAATTGTTCAATAACGGATTATCCCCGc > 1:1786222/1‑75 (MQ=255) aCCTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGc < 2:743071/75‑1 (MQ=255) aCCTTGAAAAAACTACA‑TTTCTTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGc < 1:108106/75‑1 (MQ=255) ccTTGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCt < 1:558836/75‑1 (MQ=255) ttGAAAAAACTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGa > 2:1490768/1‑75 (MQ=255) aaaaaaCTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGa < 1:1648207/75‑1 (MQ=255) aaaaaaCTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGa < 1:305825/75‑1 (MQ=255) aaaaaaCTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGa < 2:1667292/75‑1 (MQ=255) aaCTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgc > 2:825003/1‑75 (MQ=255) aCTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg > 1:1545343/1‑75 (MQ=255) aCTACATTTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAgcg > 1:1276423/1‑75 (MQ=255) | TAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTAAACCGCTTCACCTTGAAAAAACTACA‑TTT‑TTTTCATCGTGAAATGGACAGAGGATAAATTGTTCAATAACGGATTATCCCCGCTGACGAAAGCG > REL606/3875551‑3875693 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |