| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 5,159 | 0 | G | A | 45.5% | 11.6 / 10.2 | 22 | intergenic (+139/‑75) | thrC/yaaX | threonine synthase/DUF2502 domain‑containing protein |
| Reads supporting (aligned to +/- strand): ref base G (9/3); new base A (0/10); total (9/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.62e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGG‑ATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGG > NZ_CP009273/5027‑5291 | ttCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGG‑ATTGCGGAGAATaa > 2:215488/1‑138 (MQ=255) tCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGG‑ATTGCGGAGAATaac > 1:5356/1‑138 (MQ=255) tCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGG‑ATTGCGGAGAATAACAACCGCCGTtct > 2:354637/1‑138 (MQ=255) ggTTTGTTTTTAAGCAGCCGGGTTTTTTTATGAAAAAATGAGGAAAAAAAATGAGAGGAAAAAAGGAAAAATTCTC‑AATAAATGCGGTAACTGAGAGATTAGG‑ATGGCGGAAAATAAAAACCGCCGTTCTCATCGAGTa < 2:169900/139‑1 (MQ=255) ggTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGG‑ATTGCGGAGAATAACAACCGCCGTTCTCATCGAGt > 2:32531/1‑138 (MQ=255) gaaGCCCGCTTTTTTTAGGAA‑AAAATACGGGAAAAAACGGCCGGGAAAAAAGGAAAATTCCCAAACAAAGGCATTAAATTAGAAATCGGA‑TTTGCGGAAAACAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc < 1:44651/137‑1 (MQ=255) gCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGG‑ATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAc > 2:44651/1‑139 (MQ=255) cccGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGG‑ATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc > 1:52098/1‑139 (MQ=255) tttttATGAAAAAATTATGGAAAAAAAGGAAAGGAAAAAAGGAAAAATTCTA‑AAAAAAGGCGTAAACTTAGAGATTAGG‑ATTGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc < 1:469348/138‑1 (MQ=255) tttttATAAAAAAATTATGGAAAAAAATGAGAGGAAAAAAGGAAAAATTCAA‑AAAAAAGGCGGTAATTTAGAGATTAGG‑ATGGGGGAGAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc < 1:398365/138‑1 (MQ=255) ataaaaaaaTAATGGAAAAAAAGGAAAGAAAAAAAGGAAAAATTCAA‑AAAAAAGGTGGTAATTAAGAAATTAGG‑AGGGCGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa < 2:611794/132‑1 (MQ=255) aggaaaaaATTAGGGAAAAAAATGAAAGGAAAAAGGGAAAAATTCCA‑AAAAAGGGGGGAAACTAAGA‑ATTGGGTATGGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa < 1:462423/136‑1 (MQ=255) aaaaaaaaaTTAGGAAAAAAAAAGACGGGAAAAAAGGAGAAACTCTA‑AAAAATGGCTGAAACATAGAAATTGGG‑GTTGGGGAGAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa < 1:600859/132‑1 (MQ=255) ggaaaaaaGGAAAAATTCTA‑AAAAAATGCTGAAACTAAGAGTTTAGG‑ATGGGGGAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAg < 2:52098/136‑1 (MQ=255) gaAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGG‑ATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt > 1:290809/1‑139 (MQ=255) aaaaTTCTA‑AAAAAAGGCGGTAACTAAGAGATTAGG‑ATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt < 1:354637/138‑1 (MQ=255) aTGCGGTAACTTAGAGATTAGG‑ATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACttt > 2:108572/1‑139 (MQ=255) aataTTAGG‑ATTGCGTAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGt < 1:215488/136‑1 (MQ=255) agaTTAGG‑ATGGGGGAGAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc < 1:632121/139‑1 (MQ=255) agaTAAGG‑TTTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc < 1:253238/139‑1 (MQ=255) aaaTTAGG‑ATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc < 2:5356/137‑1 (MQ=255) cGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATgg > 2:194235/1‑137 (MQ=255) | TTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGG‑ATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGG > NZ_CP009273/5027‑5291 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |