Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_012967 | 2,016,091 | 0 | T | A | 31.9% | 66.3 / 26.7 | 47 | intergenic (‑123/+56) | yeeE/plaP | thiosulfate utilization transporter TsuA/YeeE/putrescine/proton symporter PlaP |
Reads supporting (aligned to +/- strand): ref base T (13/19); new base A (0/15); total (13/34) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.86e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AATGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCGGAGCCAAATATTCATTCGACT‑CCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTTCG > NC_012967/2016020‑2016162 | aaTGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAAt < 2:6837437/72‑1 (MQ=255) aaTGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAAt > 2:5233231/1‑72 (MQ=255) aaTGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAAt > 1:8953/1‑72 (MQ=255) aaTGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAAt < 1:7968204/72‑1 (MQ=255) aTGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATa > 1:6864402/1‑72 (MQ=255) aTGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATa > 2:101288/1‑72 (MQ=255) cAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGNATTTCGATAAAAAAGAATaaaa > 1:2452633/1‑72 (MQ=255) cAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATaaaa > 1:5523546/1‑72 (MQ=255) ccAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATaaaaaaaa > 1:5076201/1‑71 (MQ=255) ccAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑c > 2:7953560/1‑72 (MQ=255) ccAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑c < 1:4729758/72‑1 (MQ=255) ccAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑c > 1:4575235/1‑72 (MQ=255) cAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑cc > 1:7010629/1‑72 (MQ=255) cAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑cc < 2:2184841/72‑1 (MQ=255) cAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑cc < 2:4491204/72‑1 (MQ=255) aaTAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCg < 1:6617236/72‑1 (MQ=255) aaTAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAAAATAAAAAAA‑CCg < 1:1305039/72‑1 (MQ=255) aaaaCAAAATTCTCTAAATTAGAAAAAATGGTTATTAATGAGGATTTCGATAAAAAAAAAAAAAAAAA‑CCgg < 2:70292/70‑1 (MQ=255) aTAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑cccg > 1:6666322/1‑70 (MQ=255) aTAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCgg < 1:8605308/72‑1 (MQ=255) aaCCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCGGAg < 1:1699197/72‑1 (MQ=255) aaCAAAATACTTAAAATAAGAAAAAATGTTAATTAATGAGGATTCCGAAAAAAAAGAAAAAAAAAA‑CCGGAg < 2:4444243/72‑1 (MQ=255) aCCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAAAAAAAAAA‑CCGGAGc < 1:101430/72‑1 (MQ=255) aCCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGAAAAAAAAGAAAAAAAAAA‑CCGGAGc < 1:4245438/72‑1 (MQ=255) aaTTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCGGAGCCaaa < 1:1414631/72‑1 (MQ=255) ttCTCTAAATTAGAACAAATGTTTATAAATGAGGATTTCGATAAAAAAGAAAAAAAAAA‑CCGGAGCCAAata < 2:8164758/72‑1 (MQ=255) ttCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCGGAGCCAAata < 2:6181009/72‑1 (MQ=255) ctctAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCGGAGCCAAATatt < 2:8108163/72‑1 (MQ=255) ctctAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAAAAAAAAAA‑CCGGAGCCAAATatt < 1:5934363/72‑1 (MQ=255) tAAATAAGAAAAAATGTTTATTAATGAGGATTTGGAAAAAAAAAAAAAAAAAAA‑CCGGAGCCAAATattcat < 1:213986/72‑1 (MQ=255) aaaTTAAAAAAAATTGTTTATAAAGGGGATTTTAAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATattcat < 2:32468/72‑1 (MQ=255) aaTTAGAACAAATGTTTATTAATAAGGATTTCGATAAAAAAGAAAAAAAAAA‑CCGGAGCCAAATattcattc < 2:7008616/72‑1 (MQ=255) ttAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCCGAGCCAAAAATTCATTCGa > 2:8219125/1‑72 (MQ=255) tAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAAAATAAAAAAA‑CCGGAGCCAAATATTCATTCGAc < 2:5100076/72‑1 (MQ=255) aaaaaaTGGTTATAAATGAGGTTTAAAAAAAAAAAAAAAAAAAAAC‑CCGGAGCCAAATATTCATTCGACT‑cc < 1:5066656/69‑1 (MQ=255) atgggTTATAAAGAAGGATTTGAAAAAAAAAAAAAAAAAAAC‑CCGGAGCCAAATATTCATTCGACT‑CCGGtt < 2:4253634/69‑1 (MQ=255) aTGGTTAATAAAGGAGATTTCGAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATATTCATTCGACT‑CCGGtt < 2:7174096/72‑1 (MQ=255) ggTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAC‑CCGGAGCCAAATATTCATTCGACT‑CCGGTTGat < 2:8204672/72‑1 (MQ=255) tagtaAGGAGGATTTCAAAAAAAAAGAAAAAAAAAC‑CCGGAGCCAAATATTCATTCGACT‑CCGGTTGattat < 2:4205316/69‑1 (MQ=255) attaATGAGGATTTCGATAAAAAAGAATAAAAAAC‑CCGGAGCCAAATATTCATTCTACT‑CCGGTTGattatt < 2:519293/72‑1 (MQ=255) attaAGGAGGATTTCGATAAAAAAAAAAAAAAAAA‑CCGGAGCCAAATATTCATTCGACT‑CCGGTTGattatt < 1:4555130/72‑1 (MQ=255) ataaaTGGGGATTCGAATAAAAAAGAATAAAAAAA‑CCGGAGCCAAATATTCATTCGACT‑CCGGTTGattatt < 1:6572098/69‑1 (MQ=255) taATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCCGAGCCCAAAAATTATTTTACTCCCGGTTG‑ttattat > 2:8568075/1‑71 (MQ=255) taATGAGGATTTCGATAAAAAAGAAAAAAAAAA‑CCGGAGCCAAATATTCATTCGACT‑CCGGTTGATTATTAg < 1:8765876/72‑1 (MQ=255) aaTGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCGGAGCCAAATATTCATTCGACT‑CCGGTTGATTATTAGa < 2:8707608/72‑1 (MQ=255) aaTGAGAATTTCGAAAAAAAAAAAAAAAAAAA‑CCGGAGCCAAATATTCATTCGACT‑CCGGTTGATTATTAGa < 1:7758110/72‑1 (MQ=255) aGGATTTCAAAAAAAAAAAAAAAAAAAA‑CCGGAGCCAAATATTCATTCGACT‑CCGGTTGATTATTAGAAAAg < 1:548292/72‑1 (MQ=255) ggATTTCGATAAAAAAAAATAAAAAAA‑CCGGAGCCAAATATTCATTCGACT‑CCGGTTGATTATTAGAAAAga < 1:8492574/72‑1 (MQ=255) gATTTCGATAAAAAAGAATAAAAAAA‑CCGGAGCCAAATATTCATTCGACT‑CCGGTTGATTATTAGAAAAgat < 2:6398714/72‑1 (MQ=255) tCGATAAAAAAGAATAAAAAAA‑CCGGAGCCAAATATTCATTCGACT‑CCGGTTGATTATTAGAAAAGATGAtt < 1:4007395/72‑1 (MQ=255) gATAAAAAAGAATAAAAAAA‑CCCGGGCCAAAAAATTATTTTACT‑CCCGTTGATTATTTGAAAAGAAGATTAt > 2:8562559/1‑72 (MQ=255) aaaaaacaAAAAAAAAAAC‑CCGGAGCCAAATATTCATTCGACT‑CCGGTTGATTATTAGAAAAGATGATTATg < 2:617585/70‑1 (MQ=255) aaaaaaaaaaaaaaaaaC‑CCGGAGCCAAATATTCATTCGACT‑CCGGTTGATTATTAGAAAAGATGATTATGc < 2:2703433/71‑1 (MQ=255) aaaaaaaaaaaaaaaa‑CCGGAGCCAAATATTCATTCGACT‑CCGGTTGATTATTAGAAAAGATGATTATGCTa < 1:7974186/72‑1 (MQ=255) aaGAATAAAAAAA‑CCCGGGGCAAAAAATCATTCTACT‑CCCGGTGGTTTTTTGAAAAAAAGAATATGCTACGt > 1:8585459/1‑72 (MQ=255) aaTAAAAAAA‑CCGGAGCCAAATATTCATTCGACT‑CCGGTTGATTATTAGAAAAGATGATTATGCTACGTCtt < 2:8461884/72‑1 (MQ=255) tAAAAAAA‑CAGGAGCCAAATATTCATTCGACT‑CCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTTCg < 1:7279582/72‑1 (MQ=255) | AATGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCGGAGCCAAATATTCATTCGACT‑CCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTTCG > NC_012967/2016020‑2016162 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |