Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 395,291 | 0 | C | G | 47.2% | ‑0.8 / 64.8 | 36 | G333R (GGT→CGT) | ampH | peptidoglycan DD‑carboxypeptidase/peptidoglycan DD‑endopeptidase |
Reads supporting (aligned to +/- strand): ref base C (14/5); new base G (17/0); total (31/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.73e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.79e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATAC > NC_000913/395228‑395351 | gAGTGCCCACCCCAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCg < 1:963725/65‑1 (MQ=255) gAGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCg < 1:1795491/65‑1 (MQ=255) tGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGcga > 1:1124017/1‑63 (MQ=255) gACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGcgac > 1:2348582/1‑65 (MQ=255) ccaccaCAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGcgacca > 1:72739/1‑65 (MQ=255) ccaccaCAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGccacca > 1:929451/1‑65 (MQ=255) ccaccaCAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGccacca > 1:612536/1‑65 (MQ=255) accacAAACTCTCCAATATTCTTCTTCTTGATCATCGCCATATAGGTAATGATACGGCGaccacc > 1:2917928/1‑65 (MQ=255) accacAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGaccacc > 1:973947/1‑65 (MQ=255) accacAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGaccacc > 1:840983/1‑65 (MQ=255) accacAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGaccacc > 1:2099465/1‑65 (MQ=255) accacAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCaccacc > 1:380301/1‑65 (MQ=255) gcaACCGCGCCAATATGCGTCTGCGGGATCATCGCCATATAGGAAATGAAACCGCCACCACCGcc < 1:2346973/64‑1 (MQ=255) acaAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACGGCGACCACCGcc > 1:438077/1‑65 (MQ=255) aaaCGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTg < 1:58176/65‑1 (MQ=255) aCGCGCGAATATTCTCCTGCGGGATCATCGCCATATAGGTAATAAAACCGCCCCCACCGCCTGtc < 1:553785/65‑1 (MQ=255) cgcgCCAATTTTCTTCTGCGGGATCATCACCATATAGGTAATAAAACCGCCACCACCGCCTGtct < 1:2885826/65‑1 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGACTgaca > 1:521539/1‑62 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACAACCGCCTGtct > 1:1831534/1‑65 (MQ=255) gcCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGCCTGtcg‑tct > 1:2488706/1‑65 (MQ=255) atTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGCGTGCCA‑TCTGaataa > 1:2509574/1‑65 (MQ=255) atTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGAGTTTCA‑TATGaatac > 1:1038002/1‑64 (MQ=255) atTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGaatct > 1:1643290/1‑63 (MQ=255) atTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGaataa > 1:2071739/1‑65 (MQ=255) ttctccTGCGGGATCATCGCCATATAGGTAATGAACCCGCCACCACCGCCTGTCT‑TCTGaataat < 1:2009926/65‑1 (MQ=255) tcttctGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATc > 1:2016510/1‑65 (MQ=255) tctGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGCCCGCAT‑TCCGAATAATcccc > 1:626656/1‑65 (MQ=255) tctGCGGGATCATCGCCATATAGGTAATGAAACCGCGACCACCGCCTGTCT‑TCTGAATAATcccc > 1:363350/1‑65 (MQ=255) tctGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATcccc > 1:388423/1‑65 (MQ=255) tcggCGGAATCATCGCCATAAAGGTAATGCAACCGCCACCACCGCCTGTCT‑TCTGAATAATcccc < 1:2104704/62‑1 (MQ=255) gATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGgacgac < 1:151716/65‑1 (MQ=255) gATCATCCCCATATAGGTAATAAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGgacgac < 1:2767276/65‑1 (MQ=255) aatcatcGCCATATAGGAAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGgacgac < 1:2302572/64‑1 (MQ=255) atcatcGCCATATAGGTAATGATACGGCGACCACCGCCTGTCT‑TCTGAATAATCCCCGGCCGAcc > 1:566167/1‑65 (MQ=255) cGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCtt > 1:2926060/1‑65 (MQ=255) gCCATATAGGTAATGATACGGCGACCACCG‑AGATCTATCTGAATACCCCCCGGACGACCCTCttt > 1:2837422/1‑65 (MQ=255) ccATATAGGTAATGATACGGCGACCACCGCGATTTT‑CATGAATAACCCCCGGACGACACTCTTTc > 1:997129/1‑65 (MQ=255) aataGCTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCggg < 1:2735316/64‑1 (MQ=255) ataGGTAATGATACGGCGACCACCGCCTGTCT‑TCCGAATAATCCCCGGACGACCCGCTTTAgggg > 1:2446874/1‑65 (MQ=255) ataGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCgggg > 1:1317297/1‑65 (MQ=255) ataGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCgggg > 1:128928/1‑65 (MQ=255) taGGTAATGAACCCGCCCCCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGc < 1:2278494/65‑1 (MQ=255) ggTAATGATACGGCGACCACCGCGTGTCT‑TCTGAATAATCCCCGGACGACACTCTTTCGAGGCCa > 1:1681029/1‑65 (MQ=255) ggTAATGATACGGCGACCACCGCCTGCCT‑TCCGAATAATCCCCGGACGCCCCTCTTTCGGGGCCa > 1:1458128/1‑65 (MQ=255) ggTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCa > 1:1388078/1‑65 (MQ=255) ggTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCa > 1:1466848/1‑65 (MQ=255) gaaaTGAACCCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:232805/63‑1 (MQ=255) gaaaTAAACCCGCCACCCCCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:1711743/63‑1 (MQ=255) gaaaTAAACCCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:1529276/63‑1 (MQ=255) gTAATGAACCCCCCCCCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:2034211/65‑1 (MQ=255) gTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:2920010/65‑1 (MQ=255) gTAATAAACCCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:1622440/65‑1 (MQ=255) gTAATAAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:1874440/65‑1 (MQ=255) tAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCata < 1:2805153/65‑1 (MQ=255) aactGAACCCGCCACCCCCGCCTGTCT‑TCTCAATAATCCCCGGGCGACCCTCTTTCTGGGCCata < 1:2744076/62‑1 (MQ=255) aaaTGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCata < 1:1827262/64‑1 (MQ=255) gAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCAtatata < 1:18773/65‑1 (MQ=255) aaaCCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCCTATATAc > 1:1804400/1‑65 (MQ=255) aaaCCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATAc > 1:2019708/1‑65 (MQ=255) | GAGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATAC > NC_000913/395228‑395351 |
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: ‑ |