| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 895,856 | 0 | G | C | 43.9% | 44.2 / 36.0 | 41 | A26P (GCC→CCC) | ECB_00836 | putative Regulatory protein |
| Reads supporting (aligned to +/- strand): ref base G (22/1); new base C (0/18); total (22/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.40e-11 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGC‑TTAACAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAACAGCCAG > REL606/895783‑895930 | gACTCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGC‑TTAACAGCg‑cc < 1:5420992/76‑1 (MQ=255)gACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGC‑TTAACAGCg‑cc < 2:5523287/76‑1 (MQ=255) aCGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGC‑TTAACAGCg‑cca < 2:1316551/76‑1 (MQ=255) aCGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGC‑TTAACAGCg‑cca < 2:6118036/76‑1 (MQ=255) aCGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGC‑TTAACAGCg‑cca < 1:5888778/76‑1 (MQ=255) cgcgTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGC‑TTAACAGCg‑ccag < 2:4125897/76‑1 (MQ=255) cgcgTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGC‑TTAACAGCg‑ccag < 2:1748473/76‑1 (MQ=255) gcgTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGC‑TTAACAGCg‑ccagc < 1:6833539/76‑1 (MQ=255) cgTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGC‑TTAACAGCg‑ccagcc < 1:5328874/76‑1 (MQ=255) gTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGC‑TTAACAGCg‑ccagcca < 1:6857168/76‑1 (MQ=255) cgcgGGCAGTGTTGGTGGCGCTGGCACTCTTTGCGGTGGGCTGCCAGCGGTGGC‑GGC‑TAAACAGCG‑CCAGCCacac < 2:3210193/76‑1 (MQ=255) cgcTGGCAGTAGTGGTGGCGCTGGCCCTCGTTGCGCTGGGCTGGCAGTCGGGGC‑GGC‑TTAACAGCG‑CCAGCCacac < 2:3165773/76‑1 (MQ=255) gtGGCAGTGGGGGGGGGGCGGGCCCTCGTGGCGGGGGGCTGGCGGTGGTGGC‑GGC‑TAAACACCC‑CCAGCCACACCa < 2:2757616/75‑1 (MQ=255) tGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGG‑TTTACAACG‑CCCGCCccaccaa > 1:7170465/1‑75 (MQ=255) gCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGG‑GTTAAAGCG‑GCCGCCccaccaccc > 1:4707393/1‑73 (MQ=255) tggtggtgGCGCTGTCTATCGTTTCGCTGGGCGGGGGTTCGTGGC‑GTT‑TAACCACCC‑CCAGCCACACCATCGAAAc < 2:4875029/74‑1 (MQ=255) gggtggggGCTCTGGCGCTGGTGGCGCGGGGCTGGCGGGCGGGGC‑GTC‑TTAACAGCC‑CCAGCCACACCATCGAAAc < 2:447922/74‑1 (MQ=255) ggggtggCGCGCGCATTGGTTGCGGGGGGCGGGCAGGCGGGGC‑GGT‑TTAACGCCC‑CCAGCCACACCATCGAAACGc < 2:67725/74‑1 (MQ=255) gggtggCGCTGGCCCTGGTGGGGCGGGGCTGGCGGTGGGGGC‑GGC‑TAAACGCCC‑CCAGCCACACCATCGAAACGCa < 1:635765/75‑1 (MQ=255) ggtggCGCTGGCCCTCGTTGCGGTGGGCTGCCAGTGGTGGC‑GGC‑TAAACCGCC‑CCAGCCACACCATCGAAACGCAg < 1:5558879/76‑1 (MQ=255) ggtggCGCTGCCTCGCGTGCCGGGGGGCTGGCAGCCGGGCC‑GCC‑TTAACGTCC‑CCAGCCACACCATCGAAACGCAg < 2:6471819/76‑1 (MQ=255) ggtggCGCGGGCACTCGTGGGGCGGGGCTGGCAGGCGGGGC‑GCC‑TTAACAGCC‑CCAGCCACACCATCGAAACGCAg < 1:5179122/76‑1 (MQ=255) gggggCGCTGGCCCTCGTTGCGCTGGGCGGGCGGTGGGGGC‑GGC‑TTAACAGCC‑CCAGCCACACCATCGAAACGCAg < 1:4466008/73‑1 (MQ=255) gtggCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGG‑TTAAAAGCG‑CCCGCCACACCCTCGAAAAGCCgc > 2:1594273/1‑76 (MQ=255) tggCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGG‑GTTACAGCG‑CCAGCCACACCACCCAAAAGCAacg > 1:2291325/1‑76 (MQ=255) ggCGCTGGCCCTGGTTGCGGGGGGGGGGCGGGTGTGGG‑GGT‑TAAACAGCC‑CCAGCCACACCATCGAAACGCAgcgc < 2:1347091/76‑1 (MQ=255) cgcTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGGT‑TTAACACCG‑‑CAGCCCCACCCTCCAAAAGCCgcgcgc > 2:4958807/1‑76 (MQ=255) cgcTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGG‑TTTACAGCG‑CCAGCCCCACCCACAAAACCAAgcgcgc > 2:1615742/1‑76 (MQ=255) cgcTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGG‑TTTAAAGAG‑CCCGGCCCACCCTCGCAAAACAacgcgc > 2:2768167/1‑76 (MQ=255) ggCCCTCGTTGGGCGGGGGGGGCGGGCGGGGC‑GGC‑TAAACCGCC‑CCAGCCACACCATCGAAACGCAGCGCgcggcg < 1:1337250/76‑1 (MQ=255) tcACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGG‑TTTACCACG‑CCCGCCACACCATCGCAAAGCCACGCGCGggg‑c > 2:4575833/2‑76 (MQ=255) aCTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGG‑TTTACAGCG‑CCCGCCACACCAACCAAAAGCCGCGCGCCGCG‑C‑Tg > 1:995822/1‑76 (MQ=255) aCTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGG‑TTAACAGCG‑CCCGCCCCACCAACGCAACGCCGCGGCCGgcg‑cgg > 1:5568476/1‑74 (MQ=255) aCTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGCTTTAACA‑CG‑CCCGCCACACCATCGAAACGCCGCGCGCGggg‑ccg > 1:1883733/1‑74 (MQ=255) ctggTTGCGGTGGGGGGGCAGTCGTGGC‑GTC‑TAACCACCC‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGa < 1:6843782/73‑1 (MQ=255) cTCGTTGCGGTGGGCGGGCAGCCGGGGC‑GCC‑TAAACAGCC‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGa < 2:5259679/76‑1 (MQ=255) tCGTTGCGCGGGCGGGGCGGGCGTGGC‑GCC‑CCAACCCCC‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGaa < 1:5064621/76‑1 (MQ=255) gTTGCGCTGGGCTGGCAGTCGTGGC‑GGG‑TTTACAGCG‑CCAGCCACACCATCGCAACGCCGCGCGCGGCG‑C‑CGgaaa > 1:5435327/1‑76 (MQ=255) gTTGCGCTGGGCTGGCAGTCGTGGC‑GGG‑TTTAAAACG‑CCAGCCACACCATCCGCACGCAGCGCGCGGCG‑C‑GGcaaa > 2:2454813/1‑76 (MQ=255) gTTGCGCTGGGCTGGCAGTCGTGGC‑GGG‑GTTAAAGCG‑CCCGCCACACCAACGAAACCCACAGCGCGGCG‑C‑CGgaaa > 2:3373360/1‑76 (MQ=255) tGCGCTGGGCTGGCGGTCGGGGC‑GTC‑TAAACGGCC‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGt < 1:7112172/76‑1 (MQ=255) gcgcTGGGCTGGCAGTCGTGGC‑GGG‑TTTACAACG‑CCCGCCCCACCAACCAAACGACGCGCGCGGCG‑C‑CGgaaaaga > 1:1864515/1‑73 (MQ=255) gcgcTGGGCTGGCAGTCGTGGC‑GGG‑TTAACAGCG‑CCAGCCACACCCTCGAAACGCAGCGCGCGGGG‑C‑TGGAAAgaa > 1:6982822/1‑74 (MQ=255) gcgcGGGCCGGGCGGTGGGGGC‑GCT‑TTAACAGCC‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTa < 2:7117815/76‑1 (MQ=255) cgcTGGGCTGGCAGTCGTGGC‑GGG‑TTAACAGCG‑CCCGCCACACCACCCAAACGCCGCGCGCGGCG‑C‑CGAAAAGGaa > 1:4145490/1‑76 (MQ=255) cgcTGGGCTGGCAGTCGTGGC‑GGG‑TTAAAAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGGG‑C‑TGAAAAGaaa > 1:2076453/1‑76 (MQ=255) cgcTGGGCTGGCAGTCGTGGC‑GGC‑TTAACAGCG‑CCCGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAggca > 1:7094627/1‑73 (MQ=255) gtgggctggCGGTCGTGGC‑GCC‑TAACCAGCC‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAg < 2:6702776/75‑1 (MQ=255) ctgggctggCAGTCGTGGC‑GGG‑TTTAAAGCGCCCACCCCCCCCACCGAAACGCA‑CGCGCGGCG‑CGGgaaaagaaaa > 2:6171224/1‑76 (MQ=255) gggctggCAGTCGTGGCGGGC‑TTAAC‑GCG‑CCAACCCCACCCCCGAAAAGCAGCGCGCCGCG‑C‑Tgaaaagaaaaaag > 2:938375/1‑73 (MQ=255) gggctggCAGTCGTGGC‑GGG‑TTTACAGCG‑CCCGCCACCCCCCCCAAACGCCGCCCCCCGCG‑C‑TGGAaaagaaagag > 2:3698573/1‑74 (MQ=255) gggctggCAGTCGTGGC‑GGG‑TTTACAGCG‑CCCCCCCCCCCATCCAAAAGCAACGCCGCGCG‑C‑TGGAAAGaaaaaag > 2:2042943/1‑73 (MQ=255) gggctggCAGTCGTGGC‑GGG‑TTTACAGCG‑CCAGGCCCACCCACCAAAACCAGCGCGCCGCG‑C‑CGgaaaggaaagag > 2:4727709/1‑74 (MQ=255) gggctggCAGTCGTGGC‑GGG‑TTTAAAACG‑GCCGCCACACCCACCAAAAGCCGCGCGCGGCG‑C‑TGGAAAGaaaaaag > 2:5005524/1‑73 (MQ=255) gggctggCAGTCGTGGC‑GGG‑TTTAAAACG‑CCCACCACACCAACCAAACGCAGCGCGCGGCGCC‑Ggaaaagaaaacg > 2:2857771/1‑74 (MQ=255) gggctgcCGGTGGGGGC‑GGC‑TAAACCCCC‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAgcg < 2:2278168/76‑1 (MQ=255) tggCAGTCGTGGC‑GGG‑TTAAAAACG‑CCCGACCCACCCACGAAACGCAGCGCGCGGGG‑G‑TGCAAAGTAAAAAGCAgg > 2:1334242/1‑76 (MQ=255) gggCGGTCGGGGG‑GGC‑TGAACACCC‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAgg < 2:792444/75‑1 (MQ=255) gtcAGTTGTGGG‑GTC‑TTAACAGCC‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGa < 2:927856/74‑1 (MQ=255) gCAGTCGTGGC‑GGG‑TTAACAGCG‑CCCGCCCCACCCTCGAAAAGCAGCGCGCGGCG‑C‑TGGAAAGGAAAGAGCAGGaa > 2:950009/1‑76 (MQ=255) cGTGGC‑GGC‑TTAACAGCG‑CCAGCCCCACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGAGCAGGAACTGAc > 2:674456/1‑76 (MQ=255) gTGGC‑GGG‑TTTAAAGAG‑CCAGCCCCACCACCGAAACGCAGCGCGCGGGG‑C‑TGGAAAGTAAAACGCAAGAACTGaag > 2:248388/1‑74 (MQ=255) gTGGC‑GGC‑TTTAAAGCG‑CCCGCCACACCCACGAAAAGCAGCGCGCGGCG‑C‑TGGAAAGAAAAAAGCAGGAACTGACg > 2:3034851/1‑76 (MQ=255) gTGGC‑GGC‑TAAACCCCC‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACg < 1:3106930/76‑1 (MQ=255) gTGGC‑GCC‑TAAACACCC‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACg < 1:5100195/76‑1 (MQ=255) tGGC‑GGC‑TTAACAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACga > 2:2760927/1‑76 (MQ=255) tGGC‑GGC‑TTAACAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACga > 1:4046719/1‑76 (MQ=255) tGGC‑GCC‑TGACCAGCC‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACga < 2:3001131/76‑1 (MQ=255) ggc‑ggc‑TTAACAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACgaa > 1:2006109/1‑76 (MQ=255) ggc‑ggc‑TTAACAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACgaa > 1:3414306/1‑76 (MQ=255) ggc‑ggc‑TTAACAGCG‑CCAGCCACACCATCGAAACGCAGCACGCGGCG‑C‑TGACAAGTAAAGCGCAGGAACTGACgaa > 2:1782226/1‑76 (MQ=255) gc‑ggc‑TTAACAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACgaag > 2:1816301/1‑76 (MQ=255) c‑ggc‑TTACCCGCC‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACgaaga < 2:1001243/76‑1 (MQ=255) gc‑TTAACAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACGAAGaaa > 2:5926218/1‑76 (MQ=255) gc‑TTAACAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACGAAGaaa > 1:2343214/1‑76 (MQ=255) gc‑TTAACAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACGAAGaaa > 1:3257638/1‑76 (MQ=255) gc‑TTAACAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACGAAGaaa > 2:5498931/1‑76 (MQ=255) gc‑TTAACAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACGAAGaaa > 1:5468616/1‑76 (MQ=255) c‑TTAACAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACGAAGaaaa > 2:5910452/1‑76 (MQ=255) c‑TTAACAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACGAAGaaaa > 1:2567489/1‑76 (MQ=255) ttAACAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACGAAGaaaaa > 1:6295622/1‑76 (MQ=255) tAACAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAc > 2:726633/1‑76 (MQ=255) tAACAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAc > 2:1777221/1‑76 (MQ=255) aaCAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAca > 2:6590903/1‑76 (MQ=255) aaCAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAca > 2:1495508/1‑76 (MQ=255) aaCAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAca > 1:1936686/1‑76 (MQ=255) aCAGCC‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAcag < 1:1693546/76‑1 (MQ=255) cAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGAAGAAGAAAAAcagc > 2:6991308/1‑76 (MQ=255) gcg‑cCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAcagcca > 1:4354435/1‑76 (MQ=255) cg‑cCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAcagccag > 1:1873998/1‑76 (MQ=255) | GACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGC‑TTAACAGCG‑CCAGCCACACCATCGAAACGCAGCGCGCGGCG‑C‑TGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAACAGCCAG > REL606/895783‑895930 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |