Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,367,400 | 0 | A | C | 44.6% | 104.4 / 14.9 | 92 | G361G (GGT→GGG) | menE | o‑succinylbenzoate‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base A (28/23); new base C (0/41); total (28/64) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.69e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAG > NZ_CP009273/2367314‑2367485 | aaCTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATAccc < 1:2965392/90‑1 (MQ=255) tttttCGGCAACGGGGATCACAAACCCCTGCAAACCCGCCGGATTTGCAGAAATTACGCGCTCGCCTTCTTCCGGCTGAATCCCCTCTcc < 2:939181/86‑1 (MQ=255) cTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTcc < 2:40900/90‑1 (MQ=255) ttGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTcc < 1:790833/89‑1 (MQ=255) ttGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTcc > 2:790833/1‑89 (MQ=255) ttGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCg > 2:1525535/1‑90 (MQ=255) gTCGGCAAGGGGGACGTTAAACACCCGCAAAATCGCGGGATGTGCAGCAATTACCCGCTGGTCTTCTTCCGGCTGAATACCCTCTCCGcc < 2:2360287/90‑1 (MQ=255) gTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTCCAGCAATTACGCGCTCGCCTTCTTCCGGCTGAATCCCCTCTCCGcc < 2:169204/90‑1 (MQ=255) cGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCAc < 2:1584149/90‑1 (MQ=255) cGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCAc < 2:675518/90‑1 (MQ=255) cGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCAc < 2:1232771/90‑1 (MQ=255) cGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCAc > 2:3088905/1‑90 (MQ=255) cGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCAc > 1:2229205/1‑90 (MQ=255) gCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTg < 2:1043976/90‑1 (MQ=255) gCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTg < 2:3381784/90‑1 (MQ=255) gCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTg < 1:2050489/90‑1 (MQ=255) gCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTg > 2:694075/1‑90 (MQ=255) gCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTg < 2:2780535/90‑1 (MQ=255) gCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTg > 2:859254/1‑90 (MQ=255) gCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTg < 1:466299/90‑1 (MQ=255) gCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCCCTCCGCCACTg < 2:3306362/90‑1 (MQ=255) aaaCGGGGACGATAAACACCTGCAAAACCGCAGTATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTga < 1:2696908/89‑1 (MQ=255) aaCGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTgaa < 1:1700910/90‑1 (MQ=255) aaCGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTgaa > 2:3104145/1‑90 (MQ=255) aCGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGa < 1:3055015/90‑1 (MQ=255) cGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGAt < 1:363798/90‑1 (MQ=255) cGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCCCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGAt < 1:1123884/90‑1 (MQ=255) cGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGAt < 2:937763/90‑1 (MQ=255) cGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGAt < 1:1004922/90‑1 (MQ=255) cGATAAACACCTGCAAAACCGCAGGATGTACAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGAt < 1:1685511/90‑1 (MQ=255) aTAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTg > 2:1340081/1‑90 (MQ=255) ccTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAAc < 2:1492139/90‑1 (MQ=255) ctccAAACCCGCAGGATGTGCAGCAATTCCGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACg < 2:872537/87‑1 (MQ=255) cTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCAc < 1:2477470/68‑1 (MQ=255) cTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCAc > 2:2477470/1‑68 (MQ=255) cTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACg > 1:1081137/1‑90 (MQ=255) cTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACg > 1:2720304/1‑90 (MQ=255) tGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGt > 1:1841568/1‑90 (MQ=255) gCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTc > 1:284591/1‑90 (MQ=255) cAAACCCGCAGGATGTGCAGCAATTACGCGCTCGCCTTCTTCCGGCTAAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTcc < 2:2034395/90‑1 (MQ=255) cAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTcc < 1:2168631/90‑1 (MQ=255) cAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTcc < 2:2886922/90‑1 (MQ=255) cAAAACCGCAGGATGTGCAGCAATTACGCCCTCGCCTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTcc < 2:1755791/90‑1 (MQ=255) cAAAACCGAAGGATGTGCAGCAATTGCGCTCTGGTCTTTTTCCGGCTGAATCCCCTCCCCGCCACTGAAGAATAGATTGTCTAAACGTcc < 2:966025/90‑1 (MQ=255) cGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAAt < 1:803662/90‑1 (MQ=255) cGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAAt > 2:1459947/1‑90 (MQ=255) ggATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTc < 2:2226756/90‑1 (MQ=255) ggATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTc < 1:1402977/90‑1 (MQ=255) tgtgCAGCAATTCCGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGc < 1:371913/90‑1 (MQ=255) gtgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCt < 2:1237192/90‑1 (MQ=255) gtgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCt > 2:2097379/1‑90 (MQ=255) gtgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCt > 1:1858809/1‑90 (MQ=255) gtgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCt > 1:1465211/1‑90 (MQ=255) tgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCtt < 2:1623143/90‑1 (MQ=255) agcaATTCCGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGcc < 1:99761/90‑1 (MQ=255) caATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCAt < 2:864972/90‑1 (MQ=255) caATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCAt < 2:2686291/90‑1 (MQ=255) caATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCAc > 1:2969227/1‑45 (MQ=255) caATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCAc < 2:2969227/45‑1 (MQ=255) caATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCAt < 2:1081137/90‑1 (MQ=255) aaTTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTc < 1:2508425/79‑1 (MQ=255) aaTTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTc > 2:2508425/1‑79 (MQ=255) aTTACGCGCTCGATTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTa < 2:987411/90‑1 (MQ=255) aTTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTa > 2:701139/1‑90 (MQ=255) aCGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTg < 1:1229444/42‑1 (MQ=255) aCGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTg > 2:1229444/1‑42 (MQ=255) aCGCCCTCGCTTTTTTCCGGCTGAATCCCCTCCCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGc < 2:936795/90‑1 (MQ=255) cgcTCGATTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtc < 2:347730/90‑1 (MQ=255) cgcTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtc < 1:3278644/90‑1 (MQ=255) cgcTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtc < 2:258174/90‑1 (MQ=255) cgcTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtc < 1:2693403/90‑1 (MQ=255) cgcTCGACTTCTTCCGGCTGAATCCCCTCCCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtc < 2:1149208/90‑1 (MQ=255) gcTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtct > 1:1513221/1‑90 (MQ=255) gcTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtct > 2:2272782/1‑90 (MQ=255) gcTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtct > 2:488303/1‑90 (MQ=255) cGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCAc < 1:3088905/90‑1 (MQ=255) cGACTTCTTCCGGCTAAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCAc < 2:3001467/90‑1 (MQ=255) cttcttCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACcgc < 1:2675129/90‑1 (MQ=255) tcttcCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTAACCGCGa < 1:859254/90‑1 (MQ=255) tcttcCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTAACCGCGa < 1:694075/90‑1 (MQ=255) cTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgt < 2:3238020/90‑1 (MQ=255) cTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgt < 1:1340081/90‑1 (MQ=255) cTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgt < 1:1532886/90‑1 (MQ=255) gAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag < 2:52611/90‑1 (MQ=255) gAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag < 1:411469/90‑1 (MQ=255) gAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag > 2:2047825/1‑90 (MQ=255) gAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag > 2:2122355/1‑90 (MQ=255) gAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag < 1:1459947/90‑1 (MQ=255) gAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag < 2:3186087/90‑1 (MQ=255) gAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag > 1:3186087/1‑90 (MQ=255) gAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCACGATCGCgcgtag > 2:2622938/1‑90 (MQ=255) aaaTACCCTCCCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCTCCGCGATCGCgcgtag < 2:1932355/89‑1 (MQ=255) | AACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAG > NZ_CP009273/2367314‑2367485 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |