| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 228,175 | 0 | A | G | 55.6% | 15.2 / 21.0 | 18 | G189G (GGA→GGG) | yafD | endonuclease/exonuclease/phosphatase family protein |
| Reads supporting (aligned to +/- strand): ref base A (1/7); new base G (10/0); total (11/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.51e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGA > NZ_CP009273/228094‑228261 | gggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCaa < 1:151200/90‑1 (MQ=255)gggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCaa < 2:484656/90‑1 (MQ=255) ggATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGAATTTAATGcccg > 2:70502/1‑88 (MQ=255) gTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGc < 1:247543/90‑1 (MQ=255) gTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGGTTTAATGCCCGGGGGCGGAGaa > 2:274621/1‑90 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAATGCCCGGAGGCGTAGAAgg > 2:220347/1‑90 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGAGTTAAAT‑CCTGGGGCCGTAGAAggg > 2:284283/1‑89 (MQ=255) cAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAATGCCTGGAGCCGGAGAAAGATg > 2:171897/1‑90 (MQ=255) gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGaa < 2:171593/90‑1 (MQ=255) aCTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCTGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAAcgc < 2:243955/90‑1 (MQ=255) tCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAATTCCCGGGGCCGGAGAAAGGTGAAAGCGGTATATTGCttttc > 2:384163/1‑88 (MQ=255) aTAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTgcgctg < 1:123693/90‑3 (MQ=255) tAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCTGGAGGCGGAGAAGGATGAAAGCGGTATATCGCTTTTGGCggg > 2:7392/1‑90 (MQ=255) caccacAGCGGCCCGGTCATTATGGCGGGGGGTTTTAATGCCTGGGGCCGGAGAAGGAAGAAAGCGGTATATTGCTTTTCGCGGGgaaag > 2:225469/1‑88 (MQ=255) gCGGCCCGGTCATTATGGCGGGGGGTTTTAAAACCTGGGGCCGGAGAAGGGAGAAAGCGGTATATTGCTTTTTGCGGGGAAAATcgccgc > 1:137093/1‑90 (MQ=255) gCGGCCCGGTCATTATGGCGGGGGATTTTAATGCCCGGGGCCGGAGAAGGAGGAAAGCGGTGTATTGCTTTTCGCGGGGAAAGTCGCTgc > 1:14478/1‑90 (MQ=255) ccGGTCATTATGGCGGGGGATTTTAATGCCTGGGGGCGGAGAAGGGTGAACGCGGTATATTGCTTTTCGCGGGGAATGTCGCTGCGcccg > 2:333431/1‑88 (MQ=255) gggAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCga < 1:70502/90‑1 (MQ=255) | GGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGA > NZ_CP009273/228094‑228261 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |