Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 395,288 | 0 | A | T | 57.7% | 8.3 / 39.2 | 26 | F334I (TTC→ATC) | ampH | peptidoglycan DD‑carboxypeptidase/peptidoglycan DD‑endopeptidase |
Reads supporting (aligned to +/- strand): ref base A (7/4); new base T (15/0); total (22/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.21e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.47e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGAGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATAT > NC_000913/395227‑395349 | cGAGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATAcg > 1:2174461/1‑64 (MQ=255) cGAGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGCAATGAAAcg > 1:843538/1‑64 (MQ=255) 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) | CGAGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATAT > NC_000913/395227‑395349 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |