Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 395,294 | 0 | C | G | 36.4% | 27.0 / 41.6 | 33 | G332R (GGC→CGC) | ampH | peptidoglycan DD‑carboxypeptidase/peptidoglycan DD‑endopeptidase |
Reads supporting (aligned to +/- strand): ref base C (15/6); new base G (12/0); total (27/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.48e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAAC > NC_000913/395232‑395356 | gACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGccac < 1:2673806/65‑1 (MQ=255) accacAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGaccacc > 1:772184/1‑65 (MQ=255) accacAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCaccacc > 1:1819468/1‑65 (MQ=255) accacAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCaccacc > 1:3489229/1‑65 (MQ=255) ccacAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCg > 1:1000720/1‑65 (MQ=255) acaAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGcc > 1:2938738/1‑65 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGCGTgcta > 1:737276/1‑62 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGCGTgcta > 1:2727145/1‑62 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGtct > 1:1373142/1‑65 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGtct > 1:2851399/1‑65 (MQ=255) cgcgCCAATATTCTCCGGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGtct < 1:526712/65‑1 (MQ=255) catatTCTTCGGCGGGATCTCCGCAATATAGGTAATAAACCCGCCCCCACCGCCTGTCT‑TCTGaa < 1:1663788/64‑1 (MQ=255) atatTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGaat < 1:2002476/65‑1 (MQ=255) atTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGCGTGCTT‑TCTGcatct > 1:75583/1‑63 (MQ=255) atTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGCCTGTCT‑TCTGaataa > 1:1233898/1‑65 (MQ=255) atTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGCAAGCCT‑TCTGcataa > 1:3096477/1‑65 (MQ=255) atTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGaatct > 1:135330/1‑63 (MQ=255) ttctgcTGAGAGACCATCACCAAATAAGCAAGGAAACCGCGACCACCGCCTGTCT‑TCTGaataat < 1:2312561/65‑1 (MQ=255) tcttctGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGCCTGCCC‑TACGAATAATc > 1:2544074/1‑65 (MQ=255) tcttctGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGACATTTA‑TATGAATAATc > 1:1474487/1‑65 (MQ=255) tcttctGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATc > 1:2986115/1‑65 (MQ=255) tctGCGGGATCATCGCCATATAGGTAATGATACGGCCACCACCGCCTGTCT‑TCTGAATAATcccc > 1:3304827/1‑65 (MQ=255) tctGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATcccc > 1:496432/1‑65 (MQ=255) tctGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATcccc > 1:2804034/1‑65 (MQ=255) tctGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATcccc > 1:1023435/1‑65 (MQ=255) tctGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATcccc > 1:1476955/1‑65 (MQ=255) gggATCATCGCCATATAGGTAAGGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGgacg < 1:1469263/65‑1 (MQ=255) ggtcatcGCCATAAAGGTAATGAACCCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGgacgac < 1:1126371/63‑1 (MQ=255) gATCATCGCCATATAGGTAATGAACCCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGgacgac < 1:1654284/65‑1 (MQ=255) gATCATCGCCAGATAGTTAATGAACCCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGgacgac < 1:1351786/65‑1 (MQ=255) ccATATAGGTAATGATACGGCGACCACCGAGT‑TCTACCTGAATCATCCCCGGACGACCCTCTTTc > 1:2091603/1‑65 (MQ=255) ccATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTc > 1:2328412/1‑65 (MQ=255) aGGTAATGATACGGCGACCACCGCCTGTCT‑TCTGAATACTCCCCGGACGACCCTCTTTCGGGGct > 1:2969091/1‑64 (MQ=255) ggTAATGATACGGCGACCACCGCCTTTCT‑TATGAATAATCCCCGGACGACCCTCTTTCGGGGCCa > 1:1966032/1‑65 (MQ=255) ggTAATGATACGGCGACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCa > 1:3449371/1‑65 (MQ=255) ggTAATGATACGGCGACCACCGCCTGCTT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCa > 1:3168179/1‑65 (MQ=255) ggTAATGAACCCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCa < 1:1692667/65‑1 (MQ=255) ggTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCa > 1:2455737/1‑65 (MQ=255) gTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCat < 1:3584818/65‑1 (MQ=255) tAAGGAACCCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCata < 1:464271/65‑1 (MQ=255) aaaTGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCata < 1:2573870/64‑1 (MQ=255) aaaTAAACCCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCata < 1:3083461/64‑1 (MQ=255) aTGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCatata > 1:1641086/1‑65 (MQ=255) tGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCatatat > 1:3696996/1‑65 (MQ=255) gCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAAc < 1:647843/65‑1 (MQ=255) gCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAAc < 1:2055170/65‑1 (MQ=255) | GACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAAC > NC_000913/395232‑395356 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |