| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 5,159 | 0 | G | A | 44.8% | 30.9 / 14.3 | 29 | intergenic (+139/‑75) | thrC/yaaX | threonine synthase/DUF2502 domain‑containing protein |
| Reads supporting (aligned to +/- strand): ref base G (16/0); new base A (0/13); total (16/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.47e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAA > NZ_CP009273/5070‑5248 | cTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGgag < 1:719871/90‑1 (MQ=255) ggAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCa > 1:94166/1‑90 (MQ=255) aaaaaaaTGACAGGAAAAAAGGAAAAATTCTAAATAAATGCGGTAACTTAGAGATTAGGATTGGGGAAAAAAACAA‑CCGCCGTTCTCATc < 1:837772/88‑1 (MQ=255) aaaaaaaGGACAGAAAAAAAGAAAAAATCCAAAAAAAAGGGGGTAATTAAGAGATTAGGATGGGGAAAAAAAACAA‑CCGCCGTTCTCATc < 2:812742/88‑1 (MQ=255) aaaaaaaGGAAAGGAAAAAAGGAAAAATTCTAAAAAAATGCGGTAACTTAGAGATTAGGATTGGGGAAAATAACAA‑CCGCCGTTCTCATc < 1:85723/88‑1 (MQ=255) aaaaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAg > 1:256834/1‑90 (MQ=255) aaaaGACAGGAAAAAAGGAAAAATTCTCAAAAAATGCGGTAACTTAGAGATTAGGATTGCGGAAAATAACAA‑CCGCCGTTCTCATCGAGt < 1:400079/90‑1 (MQ=255) aaaaGAAAGGAAAAAAGGAAAAATTCAAAAAAAAGGGGTTAACTAAGAAATTAGGATGGGGGAAAAAAACAA‑CCGCCGTTCTCATCGAGt < 2:1017638/90‑1 (MQ=255) aaaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGt > 1:472737/1‑90 (MQ=255) aGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTcc > 1:812742/1‑90 (MQ=255) aGGAAAAAAGAAAAAATTCAAAAAAAAGGGGGAAATTAAGAGATTAGG‑TTGGGGAAAAAAAAAACCCGCCGTTCTCATCGAGTAATCTcc < 2:1308903/90‑1 (MQ=255) ggaaaaaaGAAAAAATTCAAAAAAAAGGCGTTAATTAAGAGATTAGGTTTGGGAAAAAAAACAC‑CCGCCGTTCTCATCGAGTAATCTCCg < 2:754881/87‑1 (MQ=255) aaGGAAAAATTCAAAATAAAGGCGGAAATTAAGAAATTAGGATGGGGAAAAAAAACAC‑CCGCCGTTCTCATCGAGTAATCTCCGGATATc < 2:1549491/90‑1 (MQ=255) ggAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGa > 2:648055/1‑90 (MQ=255) aaaaaaTGGGGTAATTTAGAAATTGGGATGGGGGAAAAAAACAC‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc < 1:995295/87‑1 (MQ=255) aaaaaaTGCGGTAATTTAGAGATTAGGATGGGGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc < 1:989250/87‑1 (MQ=255) aaaaGGCGGTAACTTAGAGATTAGGATGGCGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa < 1:14709/89‑1 (MQ=255) ggTAATTAAGAGATTAGAATGGGGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa < 1:1293273/90‑1 (MQ=255) ggTAATTAAGAAATTAGGATTGCGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa < 1:548961/90‑1 (MQ=255) ggTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa > 2:317094/1‑90 (MQ=255) ggTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa > 1:1071813/1‑90 (MQ=255) ttAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgt > 2:1185412/1‑90 (MQ=255) ttAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgt > 2:304190/1‑90 (MQ=255) ttAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgt > 1:156136/1‑90 (MQ=255) aGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa > 2:1375250/1‑90 (MQ=255) aTTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa > 2:836367/1‑90 (MQ=255) gCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGa > 2:1050889/1‑90 (MQ=255) gCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGa > 1:652384/1‑90 (MQ=255) gCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGa > 2:433179/1‑90 (MQ=255) gCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGa > 1:1036804/1‑90 (MQ=255) gaATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCaa > 1:721694/1‑90 (MQ=255) gaATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCaa > 2:430517/1‑90 (MQ=255) gaATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCaa > 1:428744/1‑90 (MQ=255) | CTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAA > NZ_CP009273/5070‑5248 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |