Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,964,466 | 0 | T | G | 53.8% | 11.0 / 16.2 | 13 | V325G (GTA→GGA) | wecB | UDP‑N‑acetylglucosamine 2‑epimerase (non‑hydrolyzing) |
Reads supporting (aligned to +/- strand): ref base T (2/4); new base G (7/0); total (9/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.10e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAACC > NZ_CP009273/3964378‑3964555 | gaagCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGTCTGGTa < 2:299684/90‑1 (MQ=255) gtgCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGGACGGGGGGGCTGGGAGGGACGGGGAAGCAGCGGAATGTc > 2:119511/1‑90 (MQ=255) gATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGGGGGTCTGGGAGGCACGGATAAGCAGCGAATTTTCGAGGGAGAGAcgcgg > 1:210585/1‑89 (MQ=255) aCTGAGCGTCCGGAAGCGGTGACGGCGGGGACGGGGCGGCTGGGGGGGACGGAAAAGCAGGGGAATGTGGGGGGAGGGAAGCGGCttttt > 1:519649/1‑89 (MQ=255) gAGCGTCCGGAAGCGGTGACGGCGGGGTGGGGGTGTCTGGGGGGCCAGGGTTACCAGGGAAATGTGGAGGGAGGGGGGCGGCTTTTtaaa > 2:511142/1‑90 (MQ=255) gCGTCCGGAAGCGGTGACGGCGGGGACGGGGCGTCTGGGAGGGACGGGTAAAAAGCGAAATTTCGAGGGAGGGAAGCGTCTTTTtaaaaa > 2:114463/1‑88 (MQ=255) ccGGAAGCGGTGACGGCGGGGACGGGGCGGCTGGGAGGGACGGGTAAGCAGCGGATTGTCGGGGGAGGGGCGCGGGTTTTTAAAAAAGaa > 2:9118/1‑90 (MQ=255) gCGGTGACGGCGGGGACGGGGCGGCTGGGAGGGACGGGTAAACAGCGAAATTTCGGGGGAGGGGAGCGGCTTTTTAAAaaagaaaaagaa > 1:512068/1‑90 (MQ=255) gTGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCcgc < 1:279302/90‑1 (MQ=255) gTGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCcgc < 1:391815/90‑1 (MQ=255) cGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGcc < 1:98139/90‑1 (MQ=255) gTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGccc > 1:372491/1‑90 (MQ=255) gTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGccc > 2:447738/1‑90 (MQ=255) tAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAAcc > 2:178602/1‑90 (MQ=255) | GAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAACC > NZ_CP009273/3964378‑3964555 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |