| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_012967 | 2,016,088 | 0 | G | A | 38.2% | 34.9 / 14.0 | 34 | intergenic (‑120/+59) | yeeE/plaP | thiosulfate utilization transporter TsuA/YeeE/putrescine/proton symporter PlaP |
| Reads supporting (aligned to +/- strand): ref base G (9/12); new base A (0/13); total (9/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.13e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGGAATGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTAC > NC_012967/2016017‑2016155 | tggAATGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAg > 2:8107373/1‑72 (MQ=255) ggAATGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAATAAGa > 2:4586583/1‑72 (MQ=255) ggAATGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGa > 1:9468919/1‑72 (MQ=255) 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) | TGGAATGCAGACCAATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTAC > NC_012967/2016017‑2016155 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |