| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 228,181 | 0 | C | T | 63.2% | ‑6.6 / 15.0 | 20 | F191F (TTC→TTT) | yafD | endonuclease/exonuclease/phosphatase family protein |
| Reads supporting (aligned to +/- strand): ref base C (0/7); new base T (12/0); total (13/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.98e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-03 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTT > NZ_CP009273/228094‑228256 | gggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGGGGTTTTaa > 1:27177/1‑90 (MQ=255)gggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCCCACCACAGCGGCCCGGTCATTATGGC‑GGGAGATTTCaa < 2:768718/90‑1 (MQ=255) gTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGAGATTTCAATGCCTGGAGc < 2:586707/90‑1 (MQ=255) aGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGAGATTTCAATGCCTGGAGCCGTAGa < 2:1806/90‑1 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGGGGTTTTAATGCCCGGGGGCGGAGAAgg > 2:416959/1‑90 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGGGGGTTTAATTCCCGGGGGCGGAGAAgg > 2:568155/1‑90 (MQ=255) gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGGGAATTTAATTACTGGGGGCGGAGAAAgaagaa > 2:702279/1‑90 (MQ=255) ccTTCCTATTGGCGATCAGAAAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAAcgc < 2:617430/89‑1 (MQ=255) gATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGGGATTTTAATTCCTGGGGGCGGAGAAGGATGAACGCGGTTTATTGCttt > 2:27905/1‑90 (MQ=255) aTAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCgg < 1:702279/90‑1 (MQ=255) caccacAGCGGCCCGGTCATTATGGC‑GGGGGGTTTTTATTTCTGGGGGCCGAGAAGGGGGAAAGCGGTTTATTGCTTTTCGCGGGgaaag > 1:262527/1‑88 (MQ=255) caccacAGCGGCCCGGTCATTATGGC‑GGGAGATTTTAATGCCTGGGGCCCTAGAAGGAAGAAAGCGTTATATCGCTTTTTGCGGGgaaag > 2:766778/1‑88 (MQ=255) gCGGCCCGGTCATTATGGCGGGGGGGTTTAAAACCCTGGGGCC‑CAGAGGGAGGAAAGCCGTATATTGCTTTTTGGGGGGAATGTCGCTgg > 1:627571/1‑89 (MQ=255) gCGGCCCGGTCATTATGGC‑GGGAGGTTTTAATGCCTGGGGGCGGAGAAGGAGGAAAGCGTTATATCGCTTTTCGCGGGGAAAGGtgctgc > 2:450280/1‑90 (MQ=255) gCCCGGTCATTATGGC‑GGGGGGTTTTAATTCCCGGGGGCGGAGAAGGGGGAACGCGGTTTATTTCTTTTCGCGGGGAAAGGCGCTGCGcc > 1:351665/1‑90 (MQ=255) gCCCGGTCATTATGGC‑GGGGGGGTTTAATGCCCGGGGGCGGAGAAGGGGGAAAGCGGTTTATCGCTTTTTGCGGGGAAAGGCGCGGCGcc > 1:742450/1‑90 (MQ=255) gCCCGGTCATTATGGC‑GGGGGATTTTAATTCCCTGGGCCGCCGAAGGGTGAAACCGGTATATCCCTTTTCGCGGGGGAAGTCGCCGCGcc > 1:228866/1‑90 (MQ=255) gCCCGGTCATTATGGC‑GGGAGGTTTTAATTCCTGGGGCCGGAGAAGGGGGAAAGCGGTATATCGCTTTTTGCGGGGGAGGTCGCTGCGcc > 1:602450/1‑90 (MQ=255) aTTATGGC‑GGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGt < 1:755572/90‑1 (MQ=255) aTGGC‑GGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGtttt < 2:322974/90‑1 (MQ=255) | GGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTT > NZ_CP009273/228094‑228256 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |