Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_012967 | 2,016,091 | 0 | T | A | 42.4% | 32.1 / 28.5 | 33 | intergenic (‑123/+56) | yeeE/plaP | thiosulfate utilization transporter TsuA/YeeE/putrescine/proton symporter PlaP |
Reads supporting (aligned to +/- strand): ref base T (6/13); new base A (0/14); total (6/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.72e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AATGCAGACCAATAACCAAA‑TTCTCTAAATTAGAACAAATGG‑TTATTAATGAGGATTTCGATAAAAAAGAAT‑AAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGT > NC_012967/2016020‑2016157 | aaTGCAGACCAATAACCAAA‑TTCTCTAAATTAGAACAAATGG‑TTATTAATGAGGATTTCGATAAAAAAGAAt < 1:5273357/72‑1 (MQ=255) aTGCAGACCAATAACCAAA‑TTCTCTAAATTAGAACAAATGG‑TTATTAATGAGGATTTCGATAAAAAAGAAT‑a < 2:4936936/72‑1 (MQ=255) aTGCAGACCAATAACCAAA‑TTCTCTAAATTAGAACAAATGG‑TTATTAATGAGGATTTCGATAAAAAAGAAT‑a > 2:6240835/1‑72 (MQ=255) gCAGACCAATAACCAAA‑TTCTCTAAATTAGAACAAATGG‑TTATTAATGAGGATTTCGATAAAAAAGAAT‑aaa < 2:6424666/72‑1 (MQ=255) gCAGACCAATAACCAAA‑TTCTCTAAATTAGAACAAATGG‑TTATTAATGAGGATTTCGATAAAAAAGAAT‑aaa > 2:3322217/1‑72 (MQ=255) gCAGACCAATAACCAAA‑TTCTCTAAATTAGAACAAATGG‑TTATTAATGAGGATTTCGATAAAAAAGAAT‑aaa < 1:5446632/72‑1 (MQ=255) aCCAATAACCAAA‑TTCTCTAAATTAGAACAAATGG‑TTATTAATGAGGATTTCGATAAAAAAGAAT‑aaaaaaa < 1:6249074/72‑1 (MQ=255) ccAATAACCAAC‑TTCTCTAAATTAGAACAAATGG‑TTATTAATGAGGATTTCGATAAAAAAGAAT‑AAAAAAA‑c > 2:3592569/1‑72 (MQ=255) ccAATAACCAAA‑TTCTCTAAATTAGAACAAATGG‑TTATTAATGAGGATTTCGATAAAAAAGAAT‑AAAAAAA‑c < 2:2620151/72‑1 (MQ=255) aTAACCAAA‑TTCTCTAAATTAGAACAAATGG‑TTATTAATGAGGATTTCGATAAAAAAGAAT‑AAAAAAA‑CCgg < 1:5548338/72‑1 (MQ=255) aTAACCAAA‑TTCTCTAAATTAGAACAAATGG‑TTATTAATGAGGATTTCGATAAAAAAGAAT‑AAAAAAA‑CCgg < 1:597942/72‑1 (MQ=255) aaCCAAA‑TTCTCTAAATTAGAACAAATGG‑TTATTAATGAGGATTTCGATAAAAAAGAAT‑AAAAAAA‑CCGGAg > 2:3148487/1‑72 (MQ=255) aCCAAATTTCTCAAAATTAAAAAAAAGGG‑TTATTAATGAGAATTCCAAAAAAAAAAAAA‑AAAAAAA‑CCGGAg < 2:5616824/72‑1 (MQ=255) ccAAA‑TTCTCTAAATTAGAACAAATGG‑TTATTAATGAGGATTTCGATAAAAAATAAA‑AAAAAAA‑CCGGAGcc < 2:7075097/72‑1 (MQ=255) cAAA‑TTCTCAAAATTAGAACAAATGG‑TTATTAATGAGGATTTCGAAAAAAAAGAAA‑AAAAAAA‑CCGGAGCCa < 2:3326696/72‑1 (MQ=255) aatttCTCAAAATTAGAAAAAATGT‑TTATTAAGGAGGATTTCGAAAAAAAAAAAA‑AAAAAAA‑CCGGAGCCaa < 1:6963499/69‑1 (MQ=255) a‑TTCTATAAATAAAAAAAAATGG‑TTATAAATGAGGATTTAAAAAAAAAAAAAA‑AAAAAAA‑CCGGAGCCAAat < 1:6284688/72‑1 (MQ=255) ttCTCTAAATTAGAACAAATGG‑TTATTAATGAGGATTTCGATAAAAAAGAAT‑AAAAAAA‑CCGGAGCCAAata < 2:5440139/72‑1 (MQ=255) tctcTAAATGAGAACAAGTGT‑TTATAAATGAGGATTTCAAAAAAAAAAAAA‑AAAAAAA‑CCGGAGCCAAatat < 2:6538303/72‑1 (MQ=255) taaaaTTAAAAAAAAAGGTTTATAAAGGAGATTTTAAAAAAAAAAAAAA‑AAAAAAC‑CCGGAGCCAAATattc < 1:2757421/70‑1 (MQ=255) ctAAATTAGAACAAATGG‑TTATTAATGAGGATTTCGAAAAAAAAAAAA‑AAAAAAC‑CCGGAGCCAAATattca < 2:6311955/72‑1 (MQ=255) aaaaTTAGAAAAAATGT‑TTATTAATGAGGTTTTCAAAAAAAAAAAAA‑AAAAAAA‑CCGGAGCCAAATattcat < 1:2425705/71‑1 (MQ=255) aTTAGAACAAATGG‑TTATTAATGAGGATTTCGATAAAAAAAAAT‑AAAAAAC‑CCGGAGCCAAATATTCATTCg < 1:4899508/72‑1 (MQ=255) aGAACAAATGT‑TTATTAATGAGGATTTCGATAAAAAAGAAA‑AAAAAAA‑CCGGAGCCAAATATTCATTCGACt < 2:5980824/72‑1 (MQ=255) aCAAATGT‑TAATTAATGAGGATTTCAATAAAAAAAAAA‑AAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCg < 2:6992517/72‑1 (MQ=255) aaaaTGG‑TTAATAAAGGAGATTTAAAAAAAAAAAAAA‑AAAAAAACCCGGAGCCAAATATTCATTCGACTCCg < 1:5046183/71‑1 (MQ=255) aaggg‑TAATTAAGGAGAATTCCGAGAAAAAAGAAA‑AAAAAAC‑CCGGAGCCAAATATTCATTCGACTCCGGtt < 2:7044853/69‑1 (MQ=255) aaTGG‑TTATTAATGAGGATTTCGATAAAAAAGAAT‑AAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCGGtt < 1:7429088/72‑1 (MQ=255) taATGAGGATTTCGATAAAAAAGAAA‑AAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAg < 1:6033109/72‑1 (MQ=255) taAAGAGGATTTCGAAAAAAAAAAAA‑AAAAAAACCCGGAGCCAAATATTCATTCGACTCCGGTTGAttatta < 2:2984816/72‑1 (MQ=255) aTGAGGATTTCGATAAAAAAGAATAAAAAAAA‑CCGGGGCCAAAAATTCATTAACCACCGGTTGtttatttaa > 2:3616108/1‑69 (MQ=255) ggATTTCGATAAAAAAGAAT‑AAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAga < 2:5778340/72‑1 (MQ=255) aTTTCGATAAAAAAGAAT‑AAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAgatg < 2:5663926/72‑1 (MQ=255) ttCGATAAAAAAGAAT‑AAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAgatgat < 2:7238140/72‑1 (MQ=255) gATAAAAAAGAAT‑AAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTAt < 1:415537/72‑1 (MQ=255) aaaaaaaaaaaaa‑aaaaACC‑CCGGACCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTAt < 2:1180671/69‑1 (MQ=255) tAAAAAAGAAT‑AAAAAAA‑CCGGAGCCAAATATTTATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGc > 1:5699713/1‑72 (MQ=255) aaaaaaaaaa‑aaaaACC‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCt < 1:4358966/72‑1 (MQ=255) aaaaaaaa‑aaaaaaC‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTAc < 1:1711800/72‑1 (MQ=255) aaGAAT‑AAAAAGC‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGt < 2:2969796/72‑1 (MQ=255) | AATGCAGACCAATAACCAAA‑TTCTCTAAATTAGAACAAATGG‑TTATTAATGAGGATTTCGATAAAAAAGAAT‑AAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGT > NC_012967/2016020‑2016157 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |