Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_012967 | 2,016,091 | 0 | T | A | 46.9% | 12.2 / 26.9 | 32 | intergenic (‑123/+56) | yeeE/plaP | thiosulfate utilization transporter TsuA/YeeE/putrescine/proton symporter PlaP |
Reads supporting (aligned to +/- strand): ref base T (6/11); new base A (0/15); total (6/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.92e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AATGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTT > NC_012967/2016020‑2016160 | aaTGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAAt > 1:3178035/1‑72 (MQ=255) gCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATaaa > 1:6752106/1‑72 (MQ=255) cAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATaaaa > 2:247510/1‑72 (MQ=255) cAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATaaaa > 1:7716824/1‑72 (MQ=255) aGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATaaaaa > 1:5473463/1‑72 (MQ=255) aaTAACCAAATTCTCTAAATTAGAACAAATTGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCg < 2:9476223/72‑1 (MQ=255) aaaTTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCGGAGCCaa < 2:7792694/72‑1 (MQ=255) ttCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCGGAGCCAAata < 1:2408627/72‑1 (MQ=255) tctAAATTAGAAAAAATGTTTATTAATGAGGATTTCGAAAAAAAAGAATAAAAAAA‑CCGGAGCCAAATattc < 1:5315050/72‑1 (MQ=255) tcaaaaTAAGAAAAAATGGTAATTAATGAGGATTTCGATAAAAAAAAATAAAAAAA‑CCGGAGCCAAATattc < 2:330811/69‑1 (MQ=255) ctaaaaTTAAAAAAAATGGTTTTTAAAGGGGATTTCAAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATatt < 1:6156457/69‑1 (MQ=255) taaaaTTAAAAAAAAGTGTTATTAAAGGGGATTTAAAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATattc < 1:3447722/70‑1 (MQ=255) aaaTTAGAAAAAATGGTTTTTAATGAGGTTTTCGAAAAAAAAGAAAAAAAAAA‑CCGGAGCCAAATattcatt < 2:7390903/72‑1 (MQ=255) aTTAGAAAAAGTGGTTATTAATGAGGATTACGATAAAAAAGAAAAAAAAAG‑CCGGAGCCAAATATTCATTCg < 2:8157229/72‑1 (MQ=255) tAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCGGAGCCAAATATTCATTCGAc < 2:7769801/72‑1 (MQ=255) tAGAAAAAATGGTTATTAATGAGGATTTCAATAAAAAAAAAAAAAAAAA‑CCGGAGCCAAATATTCATTCGAc < 1:5664925/72‑1 (MQ=255) aaaaaaGGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCGGAGCCAAATATTCATTCGACTcc < 2:7832947/69‑1 (MQ=255) aaaaaTGGTTTATAAAGGGGGTTTCGAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATATTCATTCTACTcc < 1:6779932/70‑1 (MQ=255) aaaaTGGTTTTTAAGGGGGTTTTGAAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATATTCATTCGACTCCg < 2:3208190/71‑1 (MQ=255) aaaTGGTTTATAAAGGGGATTTCAAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATATTCATTCGACTCCgg < 1:8488087/72‑1 (MQ=255) aaaTGGTTATTAATGGGAATTTCGATAAAAAAGAAAAAAAAAG‑CCGGAGCCAAATATTCATTCGACTCCGGt < 2:7190886/72‑1 (MQ=255) aaaTGGTTATTAAAGAGGATTTAAAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATATTCATTCGACTCCgg < 1:8281487/72‑1 (MQ=255) aaTGGTTATTAATAGTGATTTCGATAAAAAAGAATAAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCGGtt < 2:9845289/72‑1 (MQ=255) aTGGTTATTAATGAGGATTTCAATAAAAAAAAAAAAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCGGTTg < 1:9193961/72‑1 (MQ=255) gtttataaAGGAGGATTTCGATAAAAAAAAATAAAAAAC‑CCGGAGCCAAATATTCATTCGACTCCGGTTGat < 2:280828/70‑1 (MQ=255) gtttataaAGGAGGATTCGGAAAAAAAAAAAAAAAAAAG‑CCGGAGCCAAATATTCATTCGACTCAGGTTGat < 1:5501469/70‑1 (MQ=255) ggTTATGAATGAGGATTTCGATAAAAAAGAAAAAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCGGTTGat < 1:1961538/72‑1 (MQ=255) ttattaATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCGGTTGatta < 1:4616695/72‑1 (MQ=255) ttattaATGAGGATTTCAATAAAAAAGAATAAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCGGTTGatta < 2:8575834/72‑1 (MQ=255) taaaTGA‑GATTTCCAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATATTCATTCGACTCCGGTTGAttatta < 1:8467722/70‑1 (MQ=255) taATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCCGATCCAAAAAATTATTTGACTCCCGTTGGTTATTAg > 1:2132961/1‑72 (MQ=255) aaTGAGGATTTCGATAAAAAAAAAAAAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGa < 2:7628389/72‑1 (MQ=255) gagtttCGAAAAAAAAAAAAAAAAAAACCCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAg < 2:4234711/69‑1 (MQ=255) gATTTAGAAAAAAAAGAAAAAAAAAACCCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAga < 2:3384463/72‑1 (MQ=255) ttCGATAAAAAAGAATAAAAAAA‑CCGGAGACAAATATTCATTCGACTCCGGTTGATTATTAGAAAAgatgat > 2:9371540/1‑72 (MQ=255) gaaaaaaaaaaaaaaaaaaC‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTAt < 1:9424381/69‑1 (MQ=255) aaaaaaaaaGAAAAAAAAAC‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTAt < 1:2006600/69‑1 (MQ=255) aaaaaaaaaaaaaaaaaaC‑CCGGGGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATg < 2:3613189/70‑1 (MQ=255) aaaaaaaaaaaaaaaaaaC‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATg < 2:2580260/70‑1 (MQ=255) gaaaaaaaaaaaaaaaaC‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGc < 2:3462807/71‑1 (MQ=255) aaaaaaaaaaaaaaA‑CCGGAGCCAAATATTCATTCGANTCCGGTTGATTATTAGAAAAGATGATTATGCTAc < 1:529526/72‑1 (MQ=255) aaTAAAAAAC‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCtt < 2:8908456/72‑1 (MQ=255) | AATGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTT > NC_012967/2016020‑2016160 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |