Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,782,795 | 0 | T | C | 56.2% | 13.8 / 18.5 | 16 | V35A (GTC→GCC) | aroH | 3‑deoxy‑7‑phosphoheptulonate synthase AroH |
Reads supporting (aligned to +/- strand): ref base T (7/0); new base C (0/9); total (7/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.74e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAA‑CTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGA > NZ_CP009273/1782709‑1782882 | aCTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAA‑CTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCAc > 2:466140/1‑90 (MQ=255) ccggTAACCCCCCCCCAACCCCGCGCTAGGGTACCCCGAACCGCCTGGCGCCGCCCCCCATGCCCCCGACCCCCGCCGCAGAATTGaaaa < 2:744397/88‑1 (MQ=255) tGGTAACGCCCGCCGAA‑CTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGaaaaaa > 1:634866/1‑90 (MQ=255) ccccccgccGAA‑CTGGCGCTCCGGTACCCCATAACCCGCGGCGCCGCCCCCCAGGCCCCCGACTCCCGCCGCAGAATTGAAAAAATACTg < 1:453754/87‑1 (MQ=255) gcccgccGAA‑CTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGaa > 2:133224/1‑90 (MQ=255) ccgAC‑CCCGCGCTACGGTACCCCGTAACGCCTGGCGCCCCCCCCCATGCCCCCGCCTCCCGCCGCAGAATTGAAAAAATACTGAATGGTg < 1:140921/90‑1 (MQ=255) cgcTAGGGTACCCCGAACCGCCGGGCGCCCCCCCCCATGCCCCCGCCCCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGc < 1:192317/90‑1 (MQ=255) ggTACCCCGTAACGCCTGGCGCCGCCCCCCATGCCCCCCCCCCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGt < 1:481865/90‑1 (MQ=255) ccggTAACGCCGGGCGTCCCCCCCCATTCCCCCGCCTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTc < 2:714471/87‑1 (MQ=255) cccGTACCGCGGGGCGCCCCCCCCCATGCCCCCGCCCCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTc < 2:634866/90‑1 (MQ=255) cccGTACCGCCGGGCGCCCCCCCCCATGCCCCCCCCCCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTc < 1:219993/90‑1 (MQ=255) ccGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCa > 1:626115/1‑90 (MQ=255) aCGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCattatt > 2:49571/1‑90 (MQ=255) cTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGcc > 1:125260/1‑90 (MQ=255) cATGCCCCCGCCTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCAc < 2:316277/90‑1 (MQ=255) tGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGa > 2:526903/1‑90 (MQ=255) | ACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAA‑CTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGA > NZ_CP009273/1782709‑1782882 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |