Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 5,162 | 0 | T | A | 45.5% | 19.9 / 10.0 | 23 | intergenic (+142/‑72) | thrC/yaaX | threonine synthase/DUF2502 domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base T (10/2); new base A (0/10); total (10/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.04e-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‑ATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCA > NZ_CP009273/5027‑5296 | 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) aaTAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagca > 2:501435/1‑137 (MQ=255) | TTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGG‑ATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCA > NZ_CP009273/5027‑5296 |
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: ‑ |