| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_012967 | 2,016,098 | 1 | . | C | 24.3% | 83.1 / 28.5 | 37 | intergenic (‑130/+49) | yeeE/plaP | thiosulfate utilization transporter TsuA/YeeE/putrescine/proton symporter PlaP |
| Reads supporting (aligned to +/- strand): ref base . (2/26); new base C (0/9); total (2/35) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTTCGTACTGCG > NC_012967/2016030‑2016169 | 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) aaaaaaaaac‑CCGGAGCCAAATNTTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCtt < 1:4349351/69‑1 (MQ=255) aaTAAAAAAC‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCtt < 2:8908456/72‑1 (MQ=255) aaac‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTTCGTACt < 1:4444176/72‑1 (MQ=255) aaaa‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTTCGTACt < 1:7904473/72‑1 (MQ=255) aaaa‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTTCGTACt > 1:5711624/1‑72 (MQ=255) aaa‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTTCGTACTg > 2:2082428/1‑72 (MQ=255) aaa‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTTCGTACTg < 2:574055/72‑1 (MQ=255) aaa‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTTCGTACTg < 1:4761457/72‑1 (MQ=255) aaa‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTTCGTACTg > 1:4111941/1‑72 (MQ=255) aaa‑CCGCAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTTCGTACTg > 1:5941333/1‑72 (MQ=255) aa‑CCGGAGCCATATATTCATTCGACTCCGGATGATTATTAGAAAAGATGATTATGCTACGTCTTCGTACTGc > 2:3809448/1‑72 (MQ=255) aa‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTTCGTACTGc > 2:3217212/1‑72 (MQ=255) aa‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTTCGTACTGc > 1:698611/1‑72 (MQ=255) aa‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTTCGTACTGc > 2:3215623/1‑72 (MQ=255) a‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTTCGTACTGCg > 2:4407852/1‑72 (MQ=255) | AATAACCAAATTCTCTAAATTAGAACAAATGGTTATTAATGAGGATTTCGATAAAAAAGAATAAAAAAA‑CCGGAGCCAAATATTCATTCGACTCCGGTTGATTATTAGAAAAGATGATTATGCTACGTCTTCGTACTGCG > NC_012967/2016030‑2016169 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |