Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 87,506 | G→A | 10.9% | L178L (TTG→TTA) | ddlB → | D‑Alanine:D‑alanine ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 87,506 | 0 | G | A | 10.9% | 202.4 / 12.6 | 73 | L178L (TTG→TTA) | ddlB | D‑Alanine:D‑alanine ligase |
Reads supporting (aligned to +/- strand): ref base G (36/29); new base A (6/2); total (42/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.54e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01 |
AAAAGTAGTAGCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTC > minE/87437‑87571 | aaaaGTAGTAGCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTattga > 1:1910941/1‑71 (MQ=255) aaaaGTAGTAGCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTattg < 1:3926107/70‑1 (MQ=255) aaGTAGTAGCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTattgatt > 1:561298/1‑71 (MQ=255) agtagtAGCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTattgatt > 1:2736590/1‑70 (MQ=255) tagtagCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGaa < 1:2202706/71‑1 (MQ=255) tagtagCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGaa < 1:3522590/71‑1 (MQ=255) gtagCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTAATTGaaa > 1:1901573/1‑70 (MQ=255) gCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAAt < 1:2796325/70‑1 (MQ=255) gCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAAt < 1:3783991/70‑1 (MQ=255) gCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAAt < 1:3653775/70‑1 (MQ=255) gCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATg > 1:1022684/1‑71 (MQ=255) gCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATg > 1:49429/1‑71 (MQ=255) gCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATg > 1:3663242/1‑71 (MQ=255) gAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGc > 1:496292/1‑70 (MQ=255) aaaaTGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCt > 1:3328293/1‑70 (MQ=255) aaaaTGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTa < 1:523555/71‑1 (MQ=255) aaaTGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATg > 1:1423266/1‑66 (MQ=255) aaaTGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTaa < 1:688994/71‑1 (MQ=255) aaaTGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTa < 1:2203262/70‑1 (MQ=255) aaaTGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTa < 1:399325/70‑1 (MQ=255) aaTGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAg > 1:2776234/1‑71 (MQ=255) ctctACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTggg < 1:714186/71‑1 (MQ=255) tACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGc < 1:1853581/69‑1 (MQ=255) tACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCg > 1:952001/1‑71 (MQ=255) aCAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGaaaaa < 1:1668945/55‑1 (MQ=255) aCAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCg > 1:2252793/1‑70 (MQ=255) cAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCg < 1:850058/69‑1 (MQ=255) cAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCg > 1:874656/1‑69 (MQ=255) cAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGa < 1:1086445/71‑1 (MQ=255) gATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGtt > 1:2212623/1‑71 (MQ=255) gATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGtt > 1:2697714/1‑71 (MQ=255) gATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGtt < 1:2061174/71‑1 (MQ=255) gATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGtt > 1:873590/1‑71 (MQ=255) tGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGcc > 1:1283155/1‑63 (MQ=255) gCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCAc > 1:3363820/1‑71 (MQ=255) cATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCAc < 1:1503440/70‑1 (MQ=255) ttAAGATTGGCATTTCAGCACGATGAAGAAGTAttgattga < 1:3033176/41‑1 (MQ=255) tAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGcc > 1:717445/1‑58 (MQ=255) tAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGtt > 1:3404997/1‑71 (MQ=255) tAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGtt > 1:2015352/1‑71 (MQ=255) tAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGtt > 1:416369/1‑71 (MQ=255) tAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGt < 1:724142/70‑1 (MQ=255) aaGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACgg > 1:1952423/1‑68 (MQ=255) aaGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGtt > 1:755904/1‑70 (MQ=255) aaGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTg > 1:1783941/1‑71 (MQ=255) aaGATTGGCATTTCAGCACGATGAAGAAGTATTAATTGaa > 1:1389254/1‑40 (MQ=255) aaGATTGGCATTTCAGCACGATGAAGAAGTATTAATTGaa > 1:1066798/1‑40 (MQ=255) gATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTa > 1:169781/1‑47 (MQ=255) gATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCg < 1:1394584/71‑1 (MQ=255) gATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCg < 1:880899/71‑1 (MQ=255) gATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCg < 1:2446568/71‑1 (MQ=255) aTTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGa < 1:3248154/71‑1 (MQ=255) aTTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGa < 1:760291/71‑1 (MQ=255) tGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTgg > 1:70090/1‑49 (MQ=255) ggCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATAc > 1:1439692/1‑71 (MQ=255) aTTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCg > 1:3496096/1‑71 (MQ=255) tttCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAATTGGGCCGGAGTTCACGGTTGCGATACTCg > 1:3614355/1‑70 (MQ=255) ttCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGt > 1:1017626/1‑71 (MQ=255) ttCAGCACGATGAAGAAGTATTAATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGt > 1:2909913/1‑71 (MQ=255) ttCAGCACGATGAAGAAGTATTAATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGt > 1:1226640/1‑71 (MQ=255) tCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGt < 1:2991558/70‑1 (MQ=255) cAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGa > 1:2146666/1‑47 (MQ=255) cAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGa > 1:3112522/1‑47 (MQ=255) aGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTga < 1:3895381/70‑1 (MQ=255) gCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCgg < 1:1948444/66‑1 (MQ=255) cACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGtt < 1:2065807/47‑1 (MQ=255) cACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTgaag > 1:926755/1‑70 (MQ=255) aCGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGtt > 1:2352595/1‑46 (MQ=255) gATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAAt > 1:2343965/1‑71 (MQ=255) tGAAGAAGTATTAATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAAttt > 1:3328773/1‑71 (MQ=255) agaagTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTgaagaa < 1:1334574/64‑1 (MQ=255) gaagTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTAcc > 1:2772154/1‑71 (MQ=255) aagTATTAATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCg < 1:1257842/71‑1 (MQ=255) agTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGa > 1:513025/1‑47 (MQ=255) gTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGt > 1:3549608/1‑70 (MQ=255) gTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTc < 1:1529224/71‑1 (MQ=255) gTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTc < 1:466497/71‑1 (MQ=255) gTATTAATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTc < 1:952263/71‑1 (MQ=255) gattgaAAAATGGCTAAGTGGGCCGGAGTTCACgg > 1:1096686/1‑35 (MQ=255) | AAAAGTAGTAGCAGAAAATGCTCTACAAGATGCATTAAGATTGGCATTTCAGCACGATGAAGAAGTATTGATTGAAAAATGGCTAAGTGGGCCGGAGTTCACGGTTGCGATACTCGGTGAAGAAATTTTACCGTC > minE/87437‑87571 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |