Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 966,025 | 0 | A | G | 60.0% | 6.6 / 38.6 | 30 | Q394R (CAA→CGA) | ycaQ | winged helix DNA‑binding protein YcaQ |
Reads supporting (aligned to +/- strand): ref base A (8/4); new base G (10/8); total (18/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.09e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCTCTCTGTGGTTACAGGAAGGTATTAAACCAACGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTG > NZ_CP009273/965901‑966156 | aTCTCTCTGTGGTTACAGGAAGGTATTAAACCAACGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACt > 2:435400/1‑138 (MQ=255) tctctcTGTGGTTACAGGAAGGTATTAAACCAACGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCGAGGCCTCTTTACTGa > 2:243099/1‑139 (MQ=255) tgtgGTTACAGGAAGGTATTAAACCAACGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCGAGGCCTCTTTACTGATTGCCg > 1:167392/1‑139 (MQ=255) cAGGAAGGTATTAAACCAACGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCGAGGCCTCTTTACTGATTGCCGCACCGGCt > 1:378165/1‑139 (MQ=255) ggAAGGTATTAAACCAACGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTgg > 2:104884/1‑139 (MQ=255) aTTAAACCAACGACAACGCTGCAAAGAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGaaa > 2:534038/1‑136 (MQ=255) aCGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATaaa < 1:104884/139‑1 (MQ=255) cGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCGAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAAt > 2:374851/1‑139 (MQ=255) cGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAAt > 1:138280/1‑139 (MQ=255) cAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCGAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGAtt < 1:274127/138‑1 (MQ=255) cAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCGAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGAt < 2:167392/137‑1 (MQ=255) aaaGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCGAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGa < 1:402299/139‑1 (MQ=255) aaGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGATATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGAt < 2:175955/139‑1 (MQ=255) tACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCGAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGAtt > 2:199646/1‑139 (MQ=255) tACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCGAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCAT‑AATGAATTTTATAAAAttt > 1:240004/1‑139 (MQ=255) ggCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCAt > 1:43350/1‑139 (MQ=255) ttACTGATTTCGCTGACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCAt > 1:69243/1‑138 (MQ=255) ttACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCGAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCAt > 2:553530/1‑138 (MQ=255) tttCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCGAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGAt < 2:378165/139‑1 (MQ=255) tttCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCGAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGAt < 1:173082/139‑1 (MQ=255) cTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCGAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGAccc > 1:354642/1‑88 (MQ=255) cTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCGAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGAccc < 2:354642/88‑1 (MQ=255) cTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCGAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTgga < 1:81769/139‑1 (MQ=255) cTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCGAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTg > 2:81769/1‑137 (MQ=255) agcagGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAAc < 1:435400/139‑1 (MQ=255) agcagGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAAc > 1:178665/1‑139 (MQ=255) agGCAACGCGGGTGACATTAGGATGCTGCCCGCGAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTAt < 1:243099/139‑1 (MQ=255) ctgcCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTggg > 1:321528/1‑41 (MQ=255) ctgcCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTggg < 2:321528/41‑1 (MQ=255) tgcCCGCGAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATtg > 1:357398/1‑139 (MQ=255) | ATCTCTCTGTGGTTACAGGAAGGTATTAAACCAACGACAACGCTGCAAAAAGGGCTACGTCAGGCGATTACTGATTTCGCTAACTGGCAGCAGGCAACGCGGGTGACATTAGGATGCTGCCCGCAAGGCCTCTTTACTGATTGCCGCACCGGCTGGGAAATAGACCCCGTCGCATAAATGATTATGATAAGATTTATGCATTCATAACCCGATCCATCTGGAGGAACTATGGATCATCGTCTGCTTGAAATCATTG > NZ_CP009273/965901‑966156 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |