Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,556,841 | 0 | T | G | 57.9% | 3.7 / 16.6 | 19 | V490G (GTG→GGG) | glpD | glycerol‑3‑phosphate dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (4/4); new base G (11/0); total (15/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.81e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.36e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACA > NZ_CP009273/3556753‑3556930 | cgccgcgccgACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCtggtg < 2:510541/90‑1 (MQ=255) acgacgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgc < 1:612252/90‑1 (MQ=255) acgacgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACgc > 1:456562/1‑90 (MQ=255) cgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcagca > 1:369748/1‑90 (MQ=255) cgCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACgcagca < 2:247138/90‑1 (MQ=255) gtgGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAGGtt > 2:568744/1‑90 (MQ=255) tCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAGGTTATCGCt > 2:320088/1‑90 (MQ=255) tCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAGGGTATCGCt > 2:16202/1‑90 (MQ=255) gCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGGGGATACGCAGCAGAGGGTATCGCTgg > 2:497645/1‑90 (MQ=255) acaAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAGGGTCTCGCAGGcg > 2:112787/1‑90 (MQ=255) aaaaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCGGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGcgtc > 1:192325/1‑90 (MQ=255) aaCAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGcgtcgt < 2:360070/90‑1 (MQ=255) tgtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAAt < 1:38698/85‑1 (MQ=255) tgtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTAGAGAATACGGAGCAGAGGTTTTCGCTGGAGTCCTAAAATaaag > 1:446445/1‑88 (MQ=255) tgtgGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAAt > 2:38698/1‑85 (MQ=255) ggATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGGGATCAgg > 2:229927/1‑90 (MQ=255) ggATCAACAATCTCGTGTGAGTCAGTGGCTGGGGGAGGATACGCAGCAGAGGTTATCGCTGGGGGCGTAAATTAACGTAAGGTGGTCAgg > 2:569757/1‑90 (MQ=255) ctcGTGTGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCaa > 2:57802/1‑90 (MQ=255) gtgtGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGGGATCAGGTCAGATTTCAATCt > 2:144943/1‑90 (MQ=255) gtgtGAGTCAGTGGCTGGGGGAGTATACGCAGCAGAGGGTATCGCTGGGGGCGTAAATTAACGGAAGGTGATCAGGTCAGGTTTCAATCt > 2:403638/1‑90 (MQ=255) ggtggAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACtgatga > 1:367788/1‑90 (MQ=255) tggAGTATGCGCAGCAGAGGTTATCGCTGGCGTCGAAAATTAACGTAAGGTGATAAGGTCAGATTTCAATCTGGCCTGAGACTGATGACa < 2:370724/90‑1 (MQ=255) | CGCCGCGCCGACGACGCCCTGTGGCGTCGCACAAAACAAGGCATGTGGCTAAATGCGGATCAACAATCTCGTGTGAGTCAGTGGCTGGTGGAGTATACGCAGCAGAGGTTATCGCTGGCGTCGTAAATTAACGTAAGGTGATCAGGTCAGATTTCAATCTGGCCTGAGACTGATGACA > NZ_CP009273/3556753‑3556930 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |