Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,964,452 | 0 | T | G | 52.9% | 20.4 / 24.1 | 17 | G320G (GGT→GGG) | wecB | UDP‑N‑acetylglucosamine 2‑epimerase (non‑hydrolyzing) |
Reads supporting (aligned to +/- strand): ref base T (1/7); new base G (9/0); total (10/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.11e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGA > NZ_CP009273/3964372‑3964520 | cAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGt < 2:394592/90‑1 (MQ=255) aGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGTc < 1:329488/90‑1 (MQ=255) aagaagCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGTCTgg < 2:415123/90‑1 (MQ=255) gaagCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGTCTGGTa < 2:158573/90‑1 (MQ=255) aagCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGGACGGGGGGGCTGGTAg > 2:81428/1‑90 (MQ=255) ccTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGTCTGGTAGGCACGGATAAGCAGCGAAtt < 1:511949/90‑1 (MQ=255) gtgCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTc < 1:154272/90‑1 (MQ=255) gtgCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGGACGGGGCGTCTGGGAGGGACGGATAAGCAGCGAATTGtg > 2:144042/1‑89 (MQ=255) gCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGGACGGGGCGTCTGGTAGGGACGGATAAGCAGCGAAATGTCGa > 2:71404/1‑90 (MQ=255) ggTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGGACGGGGCGGCTGGGAGGGACGGATAAGCAGCGGAATGTCGAggg > 2:394535/1‑89 (MQ=255) tGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGGACGGGGGGGCTGGGGGGGACGGGTAAACAGGGGAATGTTGAGGGAgggg > 2:243417/1‑87 (MQ=255) tGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGGACGGGGCGTCTGGGAGGGCCGGGTAAGGAGCGGAATGTGGAgggaggga > 1:523051/1‑90 (MQ=255) gATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGTCTGGTAGGCACGGATAAGCAGCGAAtt < 1:111743/72‑1 (MQ=255) gATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGTCTGGTAGGCACGGATAAGCAGCGAAtt > 2:111743/1‑72 (MQ=255) gAGCGTCCGGAAGCGGTGACGGCGGGGACGGGGGGTCTGGGAGGGACGGGTAAACAGCGAAATTTTGGGGGAGGGACGCGGCTTTTtaaa > 1:115042/1‑90 (MQ=255) gTCCGGAAGCGGTGACGGCGGGGACGGGGGGTCTGGGAGGGACGGGTAAACAGCGAATTGTCGGGGGAGGGGCGCGGCTTTTTAAAGaag > 1:199430/1‑88 (MQ=255) tCCGGAAGCGGTGACGGCGGGGACGGGGCGTCTGGGAGGCACGGATAAGCAGGGAAATGTCGGGGGAGGGACGCGGCTTTTTAAAAACGa > 2:528326/1‑90 (MQ=255) | CAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGA > NZ_CP009273/3964372‑3964520 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |