Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 228,175 | 0 | A | G | 57.1% | 10.0 / 14.9 | 14 | G189G (GGA→GGG) | yafD | endonuclease/exonuclease/phosphatase family protein |
Reads supporting (aligned to +/- strand): ref base A (1/5); new base G (8/0); total (9/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.00e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGA > NZ_CP009273/228094‑228261 | gggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCaa < 1:310480/90‑1 (MQ=255) gTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGc < 1:468791/90‑1 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTAATGCCTCGGGCCGCAGAAgg > 2:453461/1‑90 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGGTTTAATTCCTGGGGGCGGAGAAgg > 1:178015/1‑90 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTaa > 1:247590/1‑71 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCaa < 2:247590/71‑1 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGAATTTAAATCCTGGCGGCCTAGAAgg > 2:589864/1‑90 (MQ=255) cAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAATGCCTGGGGGCGGAGAAGGaag > 2:221618/1‑88 (MQ=255) gCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTAATTCCCGGAGCCGGAGAAGGGTGAACGCGTTATATCGCt > 1:328828/1‑90 (MQ=255) gATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTAATACCCGGGGGCGCAGAAGGGGGAAAGCGGTTTATTGCttt > 1:605135/1‑90 (MQ=255) caccacAGCGGCCCGGTCATTATGGCGGGGGATTTTAATGCCTGGGGGCGGAGAAGGGGGAAAGCGGTATATTGCTTTGCGCGGGgaaag > 2:38808/1‑88 (MQ=255) cGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAgg < 2:257950/90‑1 (MQ=255) ggTCATTATGGCGGGGGATTTTAATGCCTGGAGCCGTAGAAGGATGAACGCGGTATATCGCTTTGCGCGGGGAATGTCGCTGCGCCAggg > 2:152192/1‑89 (MQ=255) gggAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCga < 2:335363/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 < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |