| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 395,291 | 0 | C | G | 57.7% | 11.6 / 36.9 | 26 | G333R (GGT→CGT) | ampH | peptidoglycan DD‑carboxypeptidase/peptidoglycan DD‑endopeptidase |
| Reads supporting (aligned to +/- strand): ref base C (6/5); new base G (15/0); total (21/5) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.02e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATAT > NC_000913/395232‑395349 | gACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGccac < 1:2606813/65‑1 (MQ=255) ccaccaCAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGcgacca > 1:3485809/1‑65 (MQ=255) ccaccaCAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGcgacca > 1:2934958/1‑65 (MQ=255) ccaccaCAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGcgacca > 1:2401454/1‑65 (MQ=255) ccaccaCAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGcgacca > 1:1819800/1‑65 (MQ=255) accacAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGaccacc > 1:1346865/1‑65 (MQ=255) accacAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGaccacc > 1:1577118/1‑65 (MQ=255) accacAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCaccacc > 1:736833/1‑65 (MQ=255) acaAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATGGGTAATGATACGGCGACCACCgag > 1:3990688/1‑63 (MQ=255) acaAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGcc > 1:3429598/1‑65 (MQ=255) acaAACGCCCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGcc < 1:2828285/65‑1 (MQ=255) cgcgCCAATGTTCTTGTGCGGAATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGtct < 1:1008072/65‑1 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCgcctgcct > 1:2352434/1‑65 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGCGTgcca > 1:458025/1‑62 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGCCTGtct > 1:3417383/1‑65 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGCCTGtct > 1:3088892/1‑65 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGACTgcca > 1:1367/1‑62 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACGGCGACCACCGACTGtct > 1:3586147/1‑65 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGtct < 1:3898207/65‑1 (MQ=255) cgcgCCAATATTCTTCTCCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGtct > 1:1794737/1‑65 (MQ=255) ttcttcTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGaataat < 1:3167190/65‑1 (MQ=255) cttctGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATcc < 1:2818084/65‑1 (MQ=255) tctGCGGGATCATCGCCATATAGGTAATGAACCCGCCACCACCGCCTGTCTTCTGAATAATcccc < 1:475985/65‑1 (MQ=255) tGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCgg > 1:3036509/1‑65 (MQ=255) gATCATCGCCATATAGGTAATGAACCCGCCGCCACCGCCTGTCTTCTGAATAATCCCCGgacgac < 1:2578187/65‑1 (MQ=255) gATCATCGCCATATAGGTAATGAAACCGCCCCCACCGCCTGTCTTCTGAATAATCCCCGgacgac < 1:3639538/65‑1 (MQ=255) atcatcGCCATATAGGTAATGATACGGCGACCACCGCGTGCTTTCTGACTCATCCCCGGACGAcc > 1:2091534/1‑65 (MQ=255) tataGGTAATGATACGGCGACCACCGCGTGCTTTATGAATACTTCCCGGACCACCCTCTTTCggg > 1:807648/1‑65 (MQ=255) tataGGTAATGAACCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCggg < 1:1135400/65‑1 (MQ=255) tataGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCggg > 1:1477819/1‑65 (MQ=255) ataGGTAATGATACGGCGACCACCGCCTGTCTTCTGAACAATCCCCGGACGACCCTCTTTCgggg > 1:1026338/1‑65 (MQ=255) aGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGcc > 1:3622805/1‑65 (MQ=255) aGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGcc > 1:2699906/1‑65 (MQ=255) ggTAATGATACGGCGACCACCGCCTTCTTTCCGAATAACCCCCGGACGACCCTCTTTCGGGGCCa > 1:3020578/1‑65 (MQ=255) taaaTAAAACCGCCGCCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:1576081/63‑1 (MQ=255) gTAATGAACCCGCCACCCCCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:1780960/65‑1 (MQ=255) gTAATGAACCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:1812559/65‑1 (MQ=255) ataataaCACCGCCACCCCCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:297854/64‑1 (MQ=255) tAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCata < 1:1098559/65‑1 (MQ=255) gaaTGAACCCCCCCCCCCCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCata < 1:2288260/64‑1 (MQ=255) gaaaaCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCatatat < 1:2704533/63‑1 (MQ=255) | GACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATAT > NC_000913/395232‑395349 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |