Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 333,793 | 0 | A | G | 57.1% | 11.0 / 16.4 | 14 | G435G (GGA→GGG) | fdrA | acyl‑CoA synthetase FdrA |
Reads supporting (aligned to +/- strand): ref base A (3/3); new base G (8/0); total (11/3) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.49e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGAT > NZ_CP009273/333714‑333882 | ccaccaACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGGGATCGGGGTGGGTTTTTTGGTCGGGTTTGGATCgg > 1:204417/1‑89 (MQ=255) caccaACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGGGATACGGGTGGGTTTTTGGGTCGGGTTTGGGTCGCa > 1:308419/1‑90 (MQ=255) accaACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGGTTTTTTGCTCGGGTTTGGGTCGCAt > 2:7555/1‑90 (MQ=255) aTCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGGTTTTGTGCTCGGGTTTGGATCGCATGAAGATCCg > 2:353265/1‑90 (MQ=255) cAGAAGTGGCGGGGATCGTGATGGATTTTTTGCTCGGGTTTTGGTCGCCTGAAGATCCGGGCGGCCCCCCCATCGAGACGAGCaaaaaaa > 1:382474/1‑89 (MQ=255) ggCGGTGGTCGTGGTGGGTTTTGTGCTTGGGATTGGAGCGCATGAAGGTCCGGGGGGGGCCACCCTCGTGACGATCAAAGaaacgaaagc > 1:280911/1‑90 (MQ=255) ggCGGTGGTCGGGATGGATTTTGTGCTCGGGTTTGGGTCGCATGAAGAGCCGGGGGGCTCCACCCTCGAGAGGATCAAAGaagggaaagg > 1:10306/1‑89 (MQ=255) ggCGGTGATCGTGATGGGTTTTTTGCTCGGGTTTTGGTCGCATGAAGATCCGGGCGGCTCCACCATCGGGACGATCAAAGaaacggaaac > 1:206454/1‑88 (MQ=255) aTCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATcgcc < 2:308419/90‑1 (MQ=255) gTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGc < 1:24860/48‑1 (MQ=255) gTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGc > 2:24860/1‑48 (MQ=255) tGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTg < 1:7555/90‑1 (MQ=255) ttGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGAc > 1:31802/1‑90 (MQ=255) cTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAg > 1:15444/1‑90 (MQ=255) aTTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGAt < 2:168029/90‑1 (MQ=255) | CCACCAACCGCATCAGTCGCTTGATCGAAGAGGCGCGCGATCCAGAAGTGGCGGTGATCGTGATGGATTTTGTGCTCGGATTTGGATCGCATGAAGATCCGGTCGGCTCCACCATCGAGACGATCAAAGAAGCGAAAGCGATCGCCGCTGCCGAAGGACGCGAGTTGAT > NZ_CP009273/333714‑333882 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |