Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 85,795 | C→A | 36.8% | A97A (GCC→GCA) | murC → | UDP‑N‑acetylmuramate:L‑alanine ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 85,795 | 0 | C | A | 36.8% | 69.2 / 50.6 | 73 | A97A (GCC→GCA) | murC | UDP‑N‑acetylmuramate:L‑alanine ligase |
Reads supporting (aligned to +/- strand): ref base C (25/21); new base A (20/7); total (45/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.35e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01 |
TGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGT > minE/85735‑85861 | tGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATg < 1:1475059/68‑1 (MQ=255) ccAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAgcg > 1:2657213/1‑68 (MQ=255) ccAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAgcg > 1:418154/1‑68 (MQ=255) ccAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAgcg > 1:1895214/1‑68 (MQ=255) ccAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAgcg > 1:2461887/1‑68 (MQ=255) ccAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAgcg > 1:679783/1‑68 (MQ=255) ccAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAgcg > 1:2465818/1‑68 (MQ=255) ccAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAgcg > 1:1139901/1‑68 (MQ=255) ccAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAgc > 1:716253/1‑67 (MQ=255) aGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGcgcg < 1:914139/68‑1 (MQ=255) tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc > 1:2113430/1‑51 (MQ=255) tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc > 1:3279938/1‑51 (MQ=255) tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc > 1:2063991/1‑51 (MQ=255) tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc > 1:1950353/1‑51 (MQ=255) tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc > 1:1912017/1‑51 (MQ=255) tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc > 1:2980834/1‑51 (MQ=255) tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc > 1:2205583/1‑51 (MQ=255) tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc > 1:1893919/1‑51 (MQ=255) tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc > 1:1849214/1‑51 (MQ=255) tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc > 1:2966344/1‑51 (MQ=255) tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc > 1:3324169/1‑51 (MQ=255) tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc > 1:2851824/1‑51 (MQ=255) tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc > 1:2833015/1‑51 (MQ=255) tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc > 1:652277/1‑51 (MQ=255) tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc > 1:668186/1‑51 (MQ=255) tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc > 1:1503255/1‑51 (MQ=255) tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc > 1:2744764/1‑51 (MQ=255) tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc > 1:926644/1‑51 (MQ=255) tCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTc > 1:270735/1‑51 (MQ=255) tCGTTGTTTCCAGCGCGATTTCTGCCCATAACCCGGAAATTGACGCAGCTc > 1:1114910/1‑51 (MQ=38) cGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGt > 1:339875/1‑69 (MQ=255) cGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGt > 1:582912/1‑69 (MQ=255) gttgttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTg > 1:1047129/1‑69 (MQ=255) gttgttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTg > 1:1731210/1‑69 (MQ=255) gttgttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTg > 1:1550106/1‑69 (MQ=255) tgttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGa < 1:829224/68‑1 (MQ=255) tgttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGa < 1:966003/68‑1 (MQ=255) tgttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGa < 1:1045018/68‑1 (MQ=255) gttTCCAGCGCTATTTCTGCCGATAACCCGGAAATTGTCGCAGCTCATGAAGcgcg < 1:2251874/56‑1 (MQ=255) gttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTCATGAAGcgcg < 1:47038/56‑1 (MQ=255) gttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTCATGAAGcgcg < 1:97616/56‑1 (MQ=255) gttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTCATGAAGcgcg < 1:2635868/56‑1 (MQ=255) gttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTCATGAAGcgcg < 1:1252370/56‑1 (MQ=255) gttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTCATGAAGcgcg < 1:1828749/56‑1 (MQ=255) gttTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCAGCTCATGAAGcgcg < 1:2065945/56‑1 (MQ=255) tCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgt > 1:1341358/1‑69 (MQ=255) tCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgt > 1:2470052/1‑69 (MQ=255) tCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgt > 1:2193091/1‑69 (MQ=255) tCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcg > 1:1099762/1‑68 (MQ=255) tCCAGCGCGATTTCTGCCGATAACCCGGAAATTGCCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgt > 1:2812007/1‑69 (MQ=255) aaCCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGCTCCGTCGTGCCGAACTGGTGGCTGa > 1:62575/1‑68 (MQ=255) aaCCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGa > 1:2453500/1‑68 (MQ=255) aaCCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGa > 1:2601393/1‑68 (MQ=255) aaCCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGa > 1:1315032/1‑68 (MQ=255) aaCCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGa > 1:2735061/1‑68 (MQ=255) aaCCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGa > 1:3245296/1‑68 (MQ=255) aaCCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGa > 1:2239244/1‑68 (MQ=255) ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt < 1:2464190/68‑1 (MQ=255) ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt < 1:2881622/68‑1 (MQ=255) ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt < 1:728799/68‑1 (MQ=255) ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt < 1:2700594/68‑1 (MQ=255) ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt < 1:2696828/68‑1 (MQ=255) ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt < 1:1358890/68‑1 (MQ=255) ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt < 1:1766462/68‑1 (MQ=255) ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt < 1:1764246/68‑1 (MQ=255) ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt < 1:517736/68‑1 (MQ=255) ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt < 1:1518461/68‑1 (MQ=255) ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt < 1:2882793/68‑1 (MQ=255) ccGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt < 1:1520824/68‑1 (MQ=255) aaTTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt < 1:607783/63‑1 (MQ=255) aaTTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt < 1:1177481/63‑1 (MQ=255) aaTTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt < 1:2883544/63‑1 (MQ=255) ttGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGtt < 1:226387/69‑1 (MQ=255) gccgcTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGt < 1:1661053/69‑1 (MQ=255) gccgcTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGt < 1:182589/69‑1 (MQ=255) gccgcTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGt < 1:3353200/69‑1 (MQ=255) gccgcTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGt < 1:2887549/69‑1 (MQ=255) gccgcTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGt < 1:2593624/69‑1 (MQ=255) gccgcTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGt < 1:3068612/69‑1 (MQ=255) gccgcTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGt < 1:1151842/69‑1 (MQ=255) | TGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGT > minE/85735‑85861 |
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: ‑ |