| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,964,452 | 0 | T | G | 56.2% | 10.9 / 17.9 | 16 | G320G (GGT→GGG) | wecB | UDP‑N‑acetylglucosamine 2‑epimerase (non‑hydrolyzing) |
| Reads supporting (aligned to +/- strand): ref base T (2/5); new base G (9/0); total (11/5) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.81e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCA > NZ_CP009273/3964368‑3964532 | cATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGGAcagg > 1:256240/1‑87 (MQ=255) aTTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACgggg > 2:249132/1‑88 (MQ=255) aTTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACgggg > 2:421703/1‑88 (MQ=255) aTTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTg < 2:153739/90‑1 (MQ=255) aGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGTc < 1:249132/90‑1 (MQ=255) aagCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGGACGGGGCGGCTggggg > 2:72450/1‑87 (MQ=255) cgcCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGGACGGGGCGTCTGGGAggga > 2:219191/1‑88 (MQ=255) aaCCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGGAGGGGGGGGCTGGGTGGGACGGGTAAGACGCGaa > 1:492354/1‑90 (MQ=255) gtgCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGGGCGGGGGGGCTGGGGGGGGCGGGTAAGCAGCGGAATGTc > 1:485498/1‑90 (MQ=255) gtgCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGGACGGGGGGTGTGGGAGGGGGGGAGAAGCAGGGAATTGtt > 2:311897/1‑89 (MQ=255) gCGTCCGGAAGCGGTGACGGCGGGGACGGGGCGGCTGGGAGGGACGGATAAACAGCGAAATTTCGGGGGAGGGGGGCGGCTTTTtaaaaa > 2:281575/1‑88 (MQ=255) gTCCGGAAGCGGTGACGGCGGGGACGGTGCGGCTGGGAGGGACGGAGAAGCAGGGGAATGTCGGGGGAGGGGGGCGGCTTTTTAAAGaag > 1:380588/1‑88 (MQ=255) gTGACGGCGGGGACAGGGCGTCTGGTAGGGACGGGGAAGCAGCGAATTGTTGGGGGAGGGACGCGGCTTTTTAAAAAAGAAAAAGAaaat > 1:34766/1‑90 (MQ=255) gACGGCGGGTACGGTGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGTATATCa < 1:421703/90‑1 (MQ=255) gACGGCGGGTACGGTGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCa < 1:281575/90‑1 (MQ=255) gACGGCGGGTACGGTGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCa < 1:515702/90‑1 (MQ=255) | CATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCA > NZ_CP009273/3964368‑3964532 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |