| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 5,159 | 0 | G | A | 53.3% | 1.4 / 14.1 | 15 | intergenic (+139/‑75) | thrC/yaaX | threonine synthase/DUF2502 domain‑containing protein |
| Reads supporting (aligned to +/- strand): ref base G (6/1); new base A (0/8); total (6/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.40e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-03 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGC‑GGTAACTTA‑GAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGA > NZ_CP009273/5074‑5243 | tttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGC‑GGTAACTTA‑GAGATTAGGATTGCGGAGAATa > 1:112460/1‑90 (MQ=255) tATGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAATAAATGG‑GGTAACTTA‑GAGATTAGGATTGGGGAAAAAAACAACCGCCGTtc < 2:459723/90‑1 (MQ=255) tATGGAAAAAAATGAAAGGAAAAAAGGAAAAATTCTCAAAAAATGC‑GGTAACTTA‑GAGATTAGGATTGCGGAAAATAACAACCGCCGTtc < 2:144222/90‑1 (MQ=255) ataaaaaaGAAAGGGAAAAAGGA‑AAA‑TCTCAAAATATACGGGAAATTTAGGAGATTAGGGTGGGGAAGAATAACACCCGCCGTTCTCATc < 2:520689/88‑1 (MQ=255) aaaaaaaTGAAAGGAAAAAAGGAAAAATTCTCAAAAAATGC‑GGTAATTTA‑GAGATTAGGATGGGGGAAAAAAAAAACCGCCGTTCTCATc < 2:272493/88‑1 (MQ=255) agaaaaaaaGGAAAAATACAAAAAAAATGG‑GGTAACTAA‑GAGATGGGGATGGAGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTcc < 1:799363/86‑1 (MQ=255) ggaaaaaaGAAAAAATTCTAAAAAAATGC‑GTTAACTTA‑GAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCg < 1:286908/87‑1 (MQ=255) aaGGAAAAATTCTAAATAAATGG‑GGTAACTTA‑GAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc < 2:588641/90‑1 (MQ=255) aaGGAAAAATTCTAAAAAAATGG‑GGTAATTTA‑GAAATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc < 2:888116/90‑1 (MQ=255) gaAATTCTCAATAAATGC‑GGTAACTTA‑GAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcct > 2:799363/1‑89 (MQ=255) cAATGAATGC‑GGTAACTTA‑GAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACggg > 1:865237/1‑90 (MQ=255) aaaaaaTGG‑GGTAATTTA‑GAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc < 2:918976/87‑1 (MQ=255) gagaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAc > 1:892835/1‑90 (MQ=255) gaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCt > 1:685538/1‑90 (MQ=255) gCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGa > 1:911993/1‑90 (MQ=255) | TTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGC‑GGTAACTTA‑GAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGA > NZ_CP009273/5074‑5243 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |