Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 489,981 | 0 | A | G | 53.8% | 7.8 / 13.8 | 13 | M41V (ATG→GTG) | recR | recombination mediator RecR |
Reads supporting (aligned to +/- strand): ref base A (1/5); new base G (7/0); total (8/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.66e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATCGGCCACTGCGCCGATTGCCGCACTTTCACCGAACAGGAAG > NZ_CP009273/489893‑490062 | ggAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGAt < 1:199654/90‑1 (MQ=255) gTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGGGGCGTCTGGCGCCgg > 1:925949/1‑90 (MQ=255) tGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGgcgc < 2:12609/90‑1 (MQ=255) tGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGGTGCGGCTGGGGCAGGGGCCCCCCCcggg > 2:715714/1‑89 (MQ=255) gCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGa < 1:298557/90‑1 (MQ=255) tGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGGGGGGGCTGGGGCCGGCGCGCCCCCCGGGGGGGTCGGGAAACCGCCCCTgcc > 1:786772/1‑89 (MQ=255) ttCACGCTGCTTCAGCGCGATCGTAGCGGCGGGGTGGGGCTGGCGCCGGCGCTCCCCCCGGCGGGGTCGGGAAAAGGCCCCCGCGcccct > 1:153141/1‑87 (MQ=255) ttCACGCTGCTTCAGCGCGATCGTAGCGGCGGGGTGCGTCTGGCGCCGGGGCTCCCCCCGGCGGGGGCGGGAAACGGCCCCCTCGCCggt > 1:152843/1‑88 (MQ=255) ttCACGCTGCTTCAGCGCGATCGTAGCGGCGGGGTGCGGCTGGCGCCGGGGCCCCCCCCGGGGGGGTCGGGAAACGGGCCCCGCGCCggt > 1:344882/1‑88 (MQ=255) gATCGTAGCGGCGGGGTGCGTCTGGCGCAGGCGCTCACCCGGGCGGTGTCGGAAAACGGCCCCCGCGCCGATTGCCCCACTTTCACCGaa > 2:522504/1‑90 (MQ=255) gATCGTAGCGGCGGGATGCGTCTGGCGCCGGCGCTCACCCCGGCGGTGTCGGAAAACGGCCCCTGCGCCGATTGCCGCCCTTTTACCCaa > 2:368026/1‑90 (MQ=255) tAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATCGGCCACTGCGCCGATTGCCGCACTTTCACCGAACAGGa < 2:52286/90‑1 (MQ=255) gcggcgGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATCGGCCACTGCGCCGATTGCCGCACTTTCACCGAACAGGAAg < 2:745303/90‑1 (MQ=255) | GGAAGCACTGCGCTGTCTGCCGGGCGTTGGCCCGAAGTCGGCGCAGCGTATGGCGTTCACGCTGCTTCAGCGCGATCGTAGCGGCGGGATGCGTCTGGCGCAGGCGCTCACCCGGGCGATGTCGGAAATCGGCCACTGCGCCGATTGCCGCACTTTCACCGAACAGGAAG > NZ_CP009273/489893‑490062 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |