| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_012967 | 2,105,084 | 0 | T | C | 23.3% | 68.0 / 11.8 | 30 | intergenic (+45/+167) | ECB_RS10540/ECB_RS10545 | glycosyltransferase family 2 protein/phage tail protein |
| Reads supporting (aligned to +/- strand): ref base T (1/22); new base C (7/0); total (8/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.93e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TATTCTTCTTTTTGCGAAAAAATATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑CCC‑TAAATTAAGGGTGAGGGGCAAAAAAAATAATTATTTGATGTTATATTTATATG > NC_012967/2105013‑2105142 | tattCTTCTTTTTGCGAAAAAATATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCt > 2:7251335/1‑72 (MQ=255) attCTTCTTTTTGCGAAAAAATATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTc < 2:54071/72‑1 (MQ=255) attCTTCTTTTTGCGAAAAAATATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTc < 1:9218933/72‑1 (MQ=255) attCTTCTTTTTGCGAAAAAATATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CACCTc < 2:4069704/72‑1 (MQ=255) ttcttcTTTTTGCGAAAAAATATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCa < 1:2600086/72‑1 (MQ=255) ttcttcTTTTTGCGAAAAAATATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCa < 2:1626548/72‑1 (MQ=255) ttcttcTTTTTGCGAAAAAATATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCa < 1:7763961/72‑1 (MQ=255) ttcttcTTTTTGCGAAAAAATATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCa < 1:7238511/72‑1 (MQ=255) tcttctTTTTGCGAAAAAATATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑c > 2:8916671/1‑72 (MQ=255) cttcttTTTGCGAAAAAATATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑cc < 2:7866532/72‑1 (MQ=255) ttcttTTTGCGAAAAAATATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑ccc < 2:3967993/72‑1 (MQ=255) tcttTTTGCGAAAAAATATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑CCC‑t < 2:3603823/72‑1 (MQ=255) tcttTTTGCGAAAAAATATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCCCC‑CCC‑t > 2:9075237/1‑72 (MQ=255) tttttGCGAAAAAATATANAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑CCC‑Taa < 1:3563654/72‑1 (MQ=255) tttttGCGAAAAAATATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑CCC‑Taa < 1:3896068/72‑1 (MQ=255) ttttGCGAAAAAATATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑CCC‑Taaa > 2:9029105/1‑72 (MQ=255) tttGCGAAAAAATATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑CCC‑TAAAt < 1:9823/72‑1 (MQ=255) tGCGAAAAAATATAAAAANACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑CCC‑TAAATTa < 1:3548724/72‑1 (MQ=255) tGCGAAAAAATATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑CCC‑TAAATTa < 2:5654480/72‑1 (MQ=255) cGAAAAAATATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑CCC‑TAAATTAAg < 2:3776320/72‑1 (MQ=255) aaaaaaTATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑CCC‑TCAAATAAagg > 2:4236255/1‑72 (MQ=255) aaaaaaTATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑CCC‑TAAATTAAggg > 1:3289749/1‑72 (MQ=255) aaaaaTATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑CCC‑TAAATTAAGGGt < 1:3415644/72‑1 (MQ=255) aaaaaTATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑CCC‑TAAATTAAGGGt < 1:343421/72‑1 (MQ=255) aaaaTATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑CCG‑TAAATTAAGGGTg < 1:3696790/72‑1 (MQ=255) aaaTATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑CCC‑TAAATTAAGGGTGa < 1:4522317/72‑1 (MQ=255) aaaTATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑CCC‑TAAATTAAGGGTGa < 2:3059848/72‑1 (MQ=255) aaaTATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCCCACCCC‑CAAATAAAgggga > 1:3957889/1‑70 (MQ=255) tataAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑CCC‑TAAATTAAGGTTGAggg < 2:8059197/72‑1 (MQ=255) cataAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCCCA‑CCC‑CCAATTAAAGGTGAggg > 1:36794/2‑72 (MQ=255) cataAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑ATCCCCA‑CCCTTAACTTAAGGGTgggg > 1:36784/2‑72 (MQ=255) aaaaaCACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCCCA‑CCC‑CCCAATAAAAGTGAgggggg > 2:4719867/1‑70 (MQ=255) aaaCACTTAACTACAGATACTCCTAACTTT‑GTATCTGACGC‑CCCCCCC‑CCC‑CAAATTAAGGTGAAGGGGCaa > 1:4208992/1‑72 (MQ=255) aaacaacAGATACT‑CTAATTTTGGTATCGGACC‑‑CCCCCCC‑CCA‑AAAATAAGGGGGGGGGGGCAAAAAAaata < 2:4904487/71‑1 (MQ=255) aCTACGGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCCCC‑CCC‑CCAACTCAATGAGAGGGGGGaaaaaaaaaaa > 1:8884678/1‑71 (MQ=255) cTACAGATACTCCTAACTTT‑GTATCTGACGC‑CCCCCCC‑CCC‑CAAATAAAAGGGGGGGGGAAAAAAAaaaaat > 2:6472831/1‑72 (MQ=255) tACAGATACTCCTAACTTT‑GTATCTGACGC‑CCCCCCC‑CCC‑TAAATAAAGGG‑GAGGGGGGAAAAAAANAaata > 2:403792/1‑72 (MQ=255) aCAGATACTCCTAACTTT‑GTATCTGACGCCCCCCCCC‑CCC‑CAAAATAAGGGCGGGGGGaaaaaaaaaaaaac > 2:7429373/1‑70 (MQ=255) aCAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCCCC‑CCC‑CCAATTAAAGGTGAGGGGGCAAAAAAATAaatat > 2:4790134/1‑72 (MQ=255) aTACTCCTAACTTT‑GTATCGGCCC‑‑CCCCCCC‑CCA‑AAAATAAGGGGGGGGGGGCAAAAAAAATAATTATTTGa < 1:7206267/72‑1 (MQ=255) ccTAACTTT‑GTATCTGACG‑‑CCCCCCC‑CCC‑CCAATTAAAGGTGAGGGGGAAAAAAAAAAAATATTTGATGttt > 1:2366947/1‑71 (MQ=255) cTAATTTT‑GTATCGGACG‑‑CCCCTCA‑CCC‑AAAATAAGGGGGGGGGGGCAAAAAAAATAATTATTTGATGttat < 2:9053764/72‑1 (MQ=255) t‑GTATCTGACG‑‑CCCCTCA‑CCC‑AAAATTAAGGGTGAGGGGCAAAAAAAATAATTATTTGATGttatatttata < 1:9085436/72‑1 (MQ=255) tATCTGACGC‑CCCCCCC‑CCC‑CAAATTAAAGGGGGGGGGGAAAAAAAA‑AAATATTTTTTGttttatttttatt > 2:7416522/1‑71 (MQ=255) | TATTCTTCTTTTTGCGAAAAAATATAAAAACACTTAACTACAGATACTCCTAACTTT‑GTATCTGACG‑‑CCCCTCA‑CCC‑TAAATTAAGGGTGAGGGGCAAAAAAAATAATTATTTGATGTTATATTTATATG > NC_012967/2105013‑2105142 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |