Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 395,291 | 0 | C | G | 38.7% | 15.2 / 39.8 | 31 | G333R (GGT→CGT) | ampH | peptidoglycan DD‑carboxypeptidase/peptidoglycan DD‑endopeptidase |
Reads supporting (aligned to +/- strand): ref base C (15/4); new base G (12/0); total (27/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.39e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTC‑TGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAAC > NC_000913/395229‑395356 | aGTGACCACCACAAACGCGCCACTATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGccac > 1:2414567/1‑68 (MQ=255) aGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGcgac > 1:794467/1‑68 (MQ=255) aGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGcgac > 1:2626749/1‑68 (MQ=255) ccaccaCAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCg > 1:2549750/1‑68 (MQ=255) aaaCGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAACCCGCCACCACCGCCTGtct < 1:1822940/68‑1 (MQ=255) aaCGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGCAATGATACGGCGACCACCGCCTGtctt > 1:973272/1‑68 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACGGCCACCACCGCCTCtcttc‑t > 1:2845775/1‑68 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGtcttc‑t > 1:1360721/1‑68 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGtcttc‑t > 1:1774461/1‑68 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGCAATGATACGGCGACCACCCCCTGcctcc‑t > 1:3510654/1‑68 (MQ=255) atTCTTCTGCGGGATCATCGCCTTATAGGTAATGAAACCGCCACCACCGCCTGTCTTC‑TGAATAATcc > 1:2407378/1‑68 (MQ=255) atTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGCCTGTCTTC‑TGAATAATcc > 1:2913735/1‑68 (MQ=255) atTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGCCTGGCTTC‑TGAATAATcc > 1:21830/1‑68 (MQ=255) atTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTC‑TGAATAATcc > 1:1117427/1‑68 (MQ=255) tcttctGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGCGTGCCTTC‑TGAATAATcccc > 1:2836340/1‑68 (MQ=255) tcttctGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTC‑TGAATAATcccg > 1:2965326/1‑67 (MQ=255) tcttctGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTC‑TGAATAATcccc > 1:607763/1‑68 (MQ=255) tcttctGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTC‑TGAATAATcccc > 1:2886165/1‑68 (MQ=255) tcttctGCGGGATCATCGCCATATAGGTAATGAAACCCACACCACCGCCTGTCTTC‑TGAATAATcccc > 1:637125/1‑68 (MQ=255) ttctGCGGGATCATCGCCATATAGGAAAGGAAACCGCCACCACCGCCTGTCTTC‑TGAATAATCCCCgg < 1:2984531/68‑1 (MQ=255) tctGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTC‑TGAATAATCCCCGga > 1:2080567/1‑68 (MQ=255) ggggAACAGCGCCATAGAGGTAATGAACCCGCCACCACCGCCTGTCTTC‑TGAATAATCCCCGgacgac < 1:995501/67‑1 (MQ=255) cGGGATCATCGCCATATAGGTAATGAAACCGCCACCGCCGCCTGTCTTC‑TGAATAATCCCCGgacgac < 1:260508/68‑1 (MQ=255) cGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTC‑TGAATAATCCCCGgacgac < 1:3336212/68‑1 (MQ=255) atcatcGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTC‑TGAATAATCCCCGGACGACCctc > 1:2671101/1‑68 (MQ=255) ccATATAGGTAATGATACGGCGACCACCGACT‑TCTACACGCCTACTCCCCGGCCGACCCTCTTTCggg > 1:332788/1‑68 (MQ=255) ccATATAGGTAATGAAACCGCCACCACCGCCTGTCTTC‑TGAATAATCCCCGGACGACCCTCTTTCggg < 1:2056957/68‑1 (MQ=255) tgaGGAAATGAACCCGCCACCACCGGCTGTCTTC‑TGAATAATCCCCGGACGACCCTCTTTCGGGGCCa < 1:1100967/66‑1 (MQ=255) ataGGTAATGATACGGCGACCACCGCCTGTCTTC‑TGAATACTCCCCGGACGACCATCTTTCGGGGCCa > 1:2728866/1‑68 (MQ=255) taGGTAATGAACCCGCCACCACCGCCTGTCTTC‑TGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:543984/68‑1 (MQ=255) gaGGTAATGTACCCGCCACCACCGCCTGTCTTC‑TGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:1588308/67‑1 (MQ=255) agggTAATGAACCCGCCACCGCCGCCTGTCTTC‑TGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:2471094/66‑1 (MQ=255) aaGGAAATGAACCCCCCACCCCCGCCTGTCTTC‑TGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:884292/67‑1 (MQ=255) ggttaataaACCCGCCCCCACCGCCTGTCTTC‑TGAATAATCCCCGGACGACCCTCTTTCGGGGCCata < 1:1548403/65‑1 (MQ=255) atgTAATAAACCCGCCACCACCGCCTGTCTTC‑TGAATAATCCCCGGACGACCCTCTTTCGGGGCCata < 1:2175914/66‑1 (MQ=255) aGGTAATAAACCCGCCACCACCGCCTGTCTTC‑TGAATAATCCCCGGACGACCCTCTTTCGGGGCCata < 1:193625/68‑1 (MQ=255) ggTAATGATACGGCGACCACCGCGTTTTTTC‑CGAATCACCCCCGGACGACCCTCTTTCGGGGCCatat > 1:437751/1‑68 (MQ=255) ggTAATGATACGGCGACCACCGCCTGTCTTC‑TGAATAATCCCCGGACGACCCTCTTTCGGGGCCatat > 1:1514493/1‑68 (MQ=255) ggTAATGATACGGCGACCACCGCCTGTCTTC‑TGAATAATCCCCGGACGACCCTCTTTCGGGGCCatat > 1:2076275/1‑68 (MQ=255) ggTAATGATACGGCGACCACCGCCTGTCTTC‑TGAATAATCCCCGGACGACCCTCTTTCGGGGCCatat > 1:1799881/1‑68 (MQ=255) ggTAATGAAACCGCCACCACCGCCTGTCTTC‑TGAATAATCCCCGGACGACCCTCTTTCGGGGCCatat > 1:1477826/1‑68 (MQ=255) aaaCCGCCACCACCGCCTGTCTTC‑TGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCa > 1:499084/1‑68 (MQ=255) aaaCCGCCACCACCGCCTGTCTTC‑TGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCa > 1:1715207/1‑68 (MQ=255) acccGCCCCCTCCGCCTGTCTTC‑TGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCaa < 1:3039481/66‑1 (MQ=255) aCCGCCACCACCGCCTGTCTTC‑TGAATAATCCCCGGACGACCCTCTTTCGGGGCCGTATATACCCAAc < 1:438204/68‑1 (MQ=255) | AGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTC‑TGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAAC > NC_000913/395229‑395356 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |