Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 395,288 | 0 | A | T | 40.6% | 26.1 / 44.7 | 32 | F334I (TTC→ATC) | ampH | peptidoglycan DD‑carboxypeptidase/peptidoglycan DD‑endopeptidase |
Reads supporting (aligned to +/- strand): ref base A (15/4); new base T (13/0); total (28/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.28e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATAT > NC_000913/395228‑395349 | gAGTGCCCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCg < 1:1945992/65‑1 (MQ=255) aGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGc > 1:2888230/1‑65 (MQ=255) 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:2727145/1‑62 (MQ=255) cgcgCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGATACGGCGACCACCGCGTgcta > 1:737276/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:1023435/1‑65 (MQ=255) tctGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATcccc > 1:496432/1‑65 (MQ=255) tctGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATcccc > 1:1476955/1‑65 (MQ=255) tctGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATcccc > 1:2804034/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) | GAGTGACCACCACAAACGCGCCAATATTCTTCTGCGGGATCATCGCCATATAGGTAATGAAACCGCCACCACCGCCTGTCT‑TCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATAT > NC_000913/395228‑395349 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |