Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 5,159 | 0 | G | A | 42.9% | 26.5 / 11.1 | 28 | intergenic (+139/‑75) | thrC/yaaX | threonine synthase/DUF2502 domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base G (15/1); new base A (0/12); total (15/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.27e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCA > NZ_CP009273/5022‑5296 | atttATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGTTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTAAATAAATGCGGTAACTTAGAGATTAGGATTGCGgaga < 2:64172/136‑1 (MQ=255) cAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGtt > 1:56589/1‑139 (MQ=255) aaTCAGGCCGGGTTGTTTTTAATCCAGCCCGGTTTTTTTATAAAAAAATTAGGGAAAAAAATGACGGGAAAAAAGGAGAAATTCACAAAAATGGCGTTAACTTAGAGATTAGGATTGGGGAAAAAAAAAACCGCCGTtc < 1:268755/139‑1 (MQ=255) gCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATc > 1:137002/1‑138 (MQ=255) gCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATc > 1:395986/1‑138 (MQ=255) gCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATc > 1:45112/1‑138 (MQ=255) ggTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGt > 1:96335/1‑138 (MQ=255) tttttATGCAGCCCGGTTTTTTTAAGAAAAAATTAGGAAAAAAAAGGACGGGAAAAAAGGAAAAATTCTAAAAAAAGGGGGTAATTTAGAGATTAGGATTGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTcc < 1:263465/138‑1 (MQ=255) aTGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGata > 2:521286/1‑139 (MQ=255) aTGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGata > 1:240140/1‑139 (MQ=255) gCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAc > 1:566757/1‑139 (MQ=255) ataaagaaATTATGGAAAAAAATGAAAGGAAAAAAGGAAAAATTCTAAATAAATGGGGTAATTGAGAAATTAGGATTCCGGAGAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa < 1:1181/135‑1 (MQ=255) ataaaaaaaTTAGGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAATAAATGCGGTAACTTAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa < 1:606750/132‑1 (MQ=255) aaaaaaaaaTTAGGAAAAAAAAGGAAAGGAAAAAAGAAAAAATTCTAAAAAAAGGCGGTAATTGAAAAATTAGGATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa < 2:96335/132‑1 (MQ=255) aTGAAAAAATGAGGAAAAAAAAGGAAAGGAAAAAAGGAAAAATTCAAAAAAAATGCGGTAACTAAGAGATTAGGATGGCGAAAAATAAAACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa < 2:137002/138‑1 (MQ=255) gaaaaaaaaTGAAAGGAAAAAAGGAAAATTCCTCAATAAATGCGGTATTTAAAAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt < 1:458510/135‑1 (MQ=255) gaAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAc > 2:365527/1‑139 (MQ=255) aaaaaaGGACAGGAAAAAAGGAAAAATTCTAAATAAATGCGTTAATTAAGAGTTTAGGTTGGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAc < 2:590252/138‑1 (MQ=255) aaaaaTGACAGGAAAAAAGGAAAAATTCTAAAAAAATGCGGTAACTTAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc < 1:9486/139‑1 (MQ=255) aaaaaGAAAAAATTCCAAAAAAGTGCTGTAACTTAGAGATAAGGATTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg < 2:45112/139‑1 (MQ=255) gaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt > 1:222840/1‑139 (MQ=255) gTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTcc > 2:324426/1‑136 (MQ=255) gTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg > 1:487311/1‑139 (MQ=255) gTAACTTAGAGATTAGGATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTcc < 1:324426/136‑1 (MQ=255) agagTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc < 2:363067/139‑1 (MQ=255) agaTTAGGAGAGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc < 2:222840/139‑1 (MQ=255) aGGATTGCGGAGTATAACAACCGCCGTTCTCATAGAGTAATCTCCGGATATCGACCCCTAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc > 2:424577/1‑139 (MQ=255) cGGAGAATGACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATgg > 2:112987/1‑137 (MQ=255) cGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATgg > 1:207744/1‑137 (MQ=255) agaAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagca < 1:112987/139‑1 (MQ=255) | ATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCA > NZ_CP009273/5022‑5296 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |