Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 395,294 | 0 | C | G | 30.8% | 43.4 / 20.0 | 26 | G332R (GGC→CGC) | ampH | peptidoglycan DD‑carboxypeptidase/peptidoglycan DD‑endopeptidase |
Reads supporting (aligned to +/- strand): ref base C (15/3); new base G (8/0); total (23/3) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.29e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAAC > NC_000913/395229‑395356 | aGTGCCCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGccac < 1:631781/68‑1 (MQ=255) tGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGCAATGATACGGcgacca > 1:1443134/1‑68 (MQ=255) ccaccaCAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCg > 1:22067/1‑68 (MQ=255) ccaccaCAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCg > 1:2648431/1‑68 (MQ=255) ccaccaCAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCg > 1:1075645/1‑68 (MQ=255) accacAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCgaa > 1:2217627/1‑66 (MQ=255) accacAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACGGCGACCACCGcc > 1:2628259/1‑68 (MQ=255) aaaCGCGCCAATATTCTTCTGCGGGATCATCGCCATATGGGTAATAAAACCGCCACCACCGCCTGtct < 1:834428/68‑1 (MQ=255) aaaCGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAACCCGCCACCACCGCCTGtct < 1:92292/68‑1 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGCCTGtcttgc > 1:195060/1‑66 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGCCTGtcttct > 1:2590782/1‑68 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGtcttct > 1:2465930/1‑68 (MQ=255) ttatTCTCCTGCGGGACCATCGCCATATAGGAAAACAACCCGCCACCACCGCCTGTCTTCTGaataat < 1:1870193/67‑1 (MQ=255) atTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGACTAATcc > 1:1534685/1‑68 (MQ=255) atTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATcc > 1:2357877/1‑68 (MQ=255) atTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATcc > 1:3394385/1‑68 (MQ=255) atTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATcc > 1:502131/1‑68 (MQ=255) atTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATcc > 1:2527771/1‑68 (MQ=255) atTCTTCTGCGGGATCATCGCCATATAGGCAATGATACGGCGACCACCGCCTGCCATCTGAATAATcc > 1:2350124/1‑68 (MQ=255) tcttctGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGCCTGCTTCCTGACTAATcccc > 1:1663763/1‑68 (MQ=255) tcttctGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGACAGTCATCCGACTCATtccc > 1:2350020/1‑68 (MQ=255) tctGCGGGATCATCGCCATATAGGTAATGAAACGGCGACCACCGCCTGTCTTCTGAATAATCCCCGga > 1:1876959/1‑68 (MQ=255) cGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGgacgac < 1:1921584/68‑1 (MQ=255) cGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGgacgac < 1:349838/68‑1 (MQ=255) cGGGATCATCGCCATATAGGTAAGAAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGgacgac < 1:110389/68‑1 (MQ=255) gggATCATCGCCATATAGNTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGAcc > 1:2573853/1‑68 (MQ=255) gCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCgg > 1:2442967/1‑68 (MQ=255) ccATATAGGTAATGAACCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCggg < 1:3332331/68‑1 (MQ=255) atatAGGAAATGAACCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGc < 1:1973965/68‑1 (MQ=255) tataGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGcc > 1:118048/1‑68 (MQ=255) ataGGTAGTGAACCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCa < 1:1302353/68‑1 (MQ=255) ataGGTAATGATACGGCGACCACCGCCTGCTTCCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCa > 1:1328409/1‑68 (MQ=255) ataGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCa > 1:701988/1‑68 (MQ=255) ataGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCa > 1:3449192/1‑68 (MQ=255) ataGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCa > 1:1262273/1‑68 (MQ=255) taGGTAATGAACCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:691374/68‑1 (MQ=255) taGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGNGGCCat < 1:2049158/68‑1 (MQ=255) agggTAATGAACCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:1190946/66‑1 (MQ=255) aaGGAAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:2874240/67‑1 (MQ=255) aGGTAATAAAGCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTGCGGGGCCata < 1:442218/68‑1 (MQ=255) ggTAATNATACGGCGACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCatat > 1:2041607/1‑68 (MQ=255) ggTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCatat > 1:3396719/1‑68 (MQ=255) ggTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCatat > 1:973751/1‑68 (MQ=255) tAATGAACCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCatatat < 1:2837012/68‑1 (MQ=255) aaTGAGACCGCCTCCCCCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCAtatata < 1:753600/68‑1 (MQ=255) aaTGAAACCGCCACCACCGCCTATCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCAtatata < 1:2433900/68‑1 (MQ=255) acccGCCCCCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCaa < 1:483066/66‑1 (MQ=255) cccGCCCCCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAAc < 1:1615219/67‑1 (MQ=255) aCCGCCCCCCCCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAAc < 1:452613/68‑1 (MQ=255) | AGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAAC > NC_000913/395229‑395356 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |