| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 395,291 | 0 | C | G | 44.4% | 10.4 / 35.1 | 27 | G333R (GGT→CGT) | ampH | peptidoglycan DD‑carboxypeptidase/peptidoglycan DD‑endopeptidase |
| Reads supporting (aligned to +/- strand): ref base C (10/5); new base G (12/0); total (22/5) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.70e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AACGAGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCA > NC_000913/395225‑395354 | aaCGAGTGCCCACCACAACCGCGCCAGTATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCg < 1:2975177/68‑1 (MQ=255) aGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGccac > 1:3253341/1‑68 (MQ=255) aGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGccac < 1:1231947/68‑1 (MQ=255) tGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGcgacca > 1:1917223/1‑68 (MQ=255) ccaccaCAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCg > 1:1055064/1‑68 (MQ=255) ccaccaCAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCg > 1:2549293/1‑68 (MQ=255) ccaccaCAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCg > 1:2341083/1‑68 (MQ=255) ccaccaCAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCg > 1:868137/1‑68 (MQ=255) accacAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGcg > 1:3006656/1‑67 (MQ=255) accacAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGcc > 1:1421073/1‑68 (MQ=255) cgcCAATATTGTTCGGCGGGATCATCGCCATATGGGTAATGAAACCGCCACCACCGCCTGTCTTCTGa < 1:3789297/68‑1 (MQ=255) cgcCAATATTCTTCTGGGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGa < 1:179137/68‑1 (MQ=255) cgcCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGa < 1:3173210/68‑1 (MQ=255) atTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATcc > 1:3031139/1‑68 (MQ=255) ttctGCGGGGTCAGCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCgg < 1:1481193/68‑1 (MQ=255) ttctGCGGGATTATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCgg < 1:112047/68‑1 (MQ=255) gggATCATCGCCATATAGGTAATGATACGGCGACCACCGCCAGTTTCACGAATAACATCCGgacgatc > 1:3073755/1‑66 (MQ=255) atcatcGCAATATAGGTAATGAAACCGCCGCCACCGCCTGTCTTCTGAATAATCCCCGGACGACCctc < 1:872840/68‑1 (MQ=255) ccATATAGGTAATGATACGGCGACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCggg > 1:1536821/1‑68 (MQ=255) ccATATAGGTAATGAAACGGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCggg > 1:99357/1‑68 (MQ=255) ccATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCggg > 1:1912560/1‑68 (MQ=255) ataGGTATTGAACCCGCCCCCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCa < 1:3652203/68‑1 (MQ=255) ataGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGctg > 1:1446110/1‑66 (MQ=255) ataGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGCGCCa > 1:3182536/1‑68 (MQ=255) tgggTAATGAACCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:2130849/66‑1 (MQ=255) taGGTAATAAACCCGCCACCCCCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:3378978/68‑1 (MQ=255) ggggTAATAAAGCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:1930003/66‑1 (MQ=255) aaGGTAATGAACCCGCCACCCCCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:1004304/67‑1 (MQ=255) aGGTAATGATACGGCGACCACCGCCGGCTATACGAATACTCCCCGGACGACCCTCCTTCGGGGCCaaa > 1:3080467/1‑66 (MQ=255) aGGAAATGAACCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCata < 1:412433/68‑1 (MQ=255) ggTAATGATACGGCGACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCatat > 1:3214570/1‑68 (MQ=255) ggTAATGATACGGCGACCACCGAGATTCACCTGACTAACTCCCGGACGACCCTCTTTCGGGGACatat > 1:1179298/1‑68 (MQ=255) ggTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCatat > 1:3051743/1‑68 (MQ=255) gTAATGATACGGCGACCACCGCCTGCCTCCTGCATCATCCCCGGACGACCCTCTTTCGGGTCCatata > 1:3283546/1‑68 (MQ=255) aaTGAAACCGGCTCCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCAtatata < 1:3592229/68‑1 (MQ=255) aaTGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCatatatc > 1:2156154/1‑67 (MQ=255) aaaCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATCTATAcctg > 1:3475512/1‑66 (MQ=255) | AACGAGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCA > NC_000913/395225‑395354 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |