Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 5,159 | 0 | G | A | 57.1% | 1.5 / 28.2 | 28 | intergenic (+139/‑75) | thrC/yaaX | threonine synthase/DUF2502 domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base G (12/0); new base A (0/16); total (12/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.29e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-06 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAA > NZ_CP009273/5070‑5248 | cTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGgag < 1:612443/90‑1 (MQ=255) aTGAAAAAATTAGGAAAAAAAATAAAAAGAAAAAAGGAAAAATCCCAAAAAAATGGGGTAACTTAGAGATTGGGATTGGGGAAAAAaaca < 2:186210/90‑1 (MQ=255) aaaaaaaTGACAGGAAAAAAGGAAAAATTCTCAATAAAGGCGGTAACTTAGAGATTAGGATTGCGGAAAATAACAACCGCCGTTCTCATc < 2:1229805/88‑1 (MQ=255) aaaaaaaGGAAAGGAAAAAAGAAAAATTCCAAAAAAAAGGGGGTAATTAAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATc < 2:225741/88‑1 (MQ=255) gACAGGAAAAAAGGAAAAATTCTCAATAAATGGGGTAACTTAGAGATTAGGATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAAtc < 1:443782/90‑1 (MQ=255) gggAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCg > 2:129873/1‑90 (MQ=255) aaGGAAAATTTCTAAAAAAAGGGGGTAATTTAGAGTTTGGGTTGGGGGAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATc < 2:364714/90‑1 (MQ=255) aaGGAAAAATTCAAAAAAAAGGGGGAAATTAAGAAATTAGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc < 2:485703/90‑1 (MQ=255) agaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGgag < 2:1110314/46‑1 (MQ=255) agaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGgag > 1:1110314/1‑46 (MQ=255) agaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc > 1:1284726/1‑90 (MQ=255) agaAATTCCCAATAAATGCGGTAACTTAGAGATTGGGTTTGCGGAAAAAACCACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc < 2:24785/90‑1 (MQ=255) aaaaTGGGGTAATTTAGAGATTGGGTTGGGGGAAAAAAAAACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa < 2:165782/89‑1 (MQ=255) aaaaTGCGTTAACTAAGAGATTAGGATGGGGAAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa < 2:1080826/89‑1 (MQ=255) aaaaGGCGTTAATTTAGAGATAAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa < 1:129873/89‑1 (MQ=255) aaaaGGCGGTAATTTAGAGATAAGGTTGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa < 2:728166/89‑1 (MQ=255) aaaaGGCGGTAATTAAGAGATAAGGATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa < 1:79709/89‑1 (MQ=255) gCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATa > 1:12933/1‑90 (MQ=255) ggaaaCTGAGAGATTGGGATCGCGAGAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa < 1:1145849/87‑1 (MQ=255) ggTAACTTAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa < 1:81562/90‑1 (MQ=255) tAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTa > 1:749372/1‑90 (MQ=255) ttgggATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgt > 2:945785/4‑90 (MQ=255) aGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc > 2:102354/1‑56 (MQ=255) aGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa > 2:1066572/1‑90 (MQ=255) aGGATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc < 1:102354/56‑1 (MQ=255) ttGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaaa > 2:714318/1‑90 (MQ=255) gCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa > 2:432068/1‑62 (MQ=255) gCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGa > 2:975037/1‑90 (MQ=255) gCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGa > 1:73267/1‑90 (MQ=255) gCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGa > 1:103734/1‑90 (MQ=255) gCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa < 1:432068/62‑1 (MQ=255) gaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTa > 2:588571/1‑71 (MQ=255) gaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCaa > 2:349943/1‑90 (MQ=255) | CTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAA > 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 < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |