Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 895,856 | 0 | G | C | 39.5% | 51.9 / 31.1 | 38 | A26P (GCC→CCC) | ECB_00836 | putative Regulatory protein |
Reads supporting (aligned to +/- strand): ref base G (21/2); new base C (0/15); total (21/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.79e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑G‑G‑CTTAA‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAACAGCCA > REL606/895781‑895929 | aTGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑G‑G‑CTTAA‑CAgcg < 1:6427796/76‑1 (MQ=255) gACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑G‑G‑CTTAA‑CAGCgcc < 2:3255078/76‑1 (MQ=255) gACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑G‑G‑CTTAA‑CAGCgcc < 1:5108071/76‑1 (MQ=255) aCGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑G‑G‑CTTAA‑CAGCgcca < 2:5811306/76‑1 (MQ=255) aCGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑G‑G‑CTTAA‑CAGCgcca < 2:2042993/76‑1 (MQ=255) aCGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑G‑G‑CTTAA‑CAGCgcca < 2:4648721/76‑1 (MQ=255) gTGCGCTGGCAGTAGTGCTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑G‑G‑CTTAA‑CAGCgccag‑cca < 1:684027/76‑1 (MQ=255) cgcTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑G‑G‑CTTAA‑CAGCGCCAG‑CCacac < 1:1692167/76‑1 (MQ=255) cTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGG‑G‑G‑GGTTT‑AAGCGCCAC‑CCCCAccc > 1:1264288/1‑75 (MQ=255) tGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑G‑G‑GGTGA‑CAAAGCCCG‑CCccaccac > 1:2303303/1‑75 (MQ=255) tGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGGCGTGGCGG‑G‑GTTAA‑CACGCCCAG‑CCCCAccc > 2:3221757/1‑75 (MQ=255) gggCTGTAGTGGTGGCGCTGACGCGCGTTGCGTGGGGCCGCCGGTCGTGGC‑G‑C‑CTAAA‑CCGCCCCAG‑CCACACCAt < 2:5521310/75‑1 (MQ=255) agtGGTGGCGCTGACCCTCGTTGGGTGTGGCGGCCCGTCGGGGG‑G‑G‑TTAAG‑CACCGCCAG‑CCACACCATCGAAACg < 2:2388854/76‑1 (MQ=255) ggtggCGCGGGCCCTGTTTGCGCGGGGCGGGAGGGGGTGGG‑G‑T‑CTAAA‑CACCCCCAG‑CCACACCATCGAAACGCAg < 1:272311/76‑1 (MQ=255) gggggccctggccctGGTTGCGCGGGGCTGGCGGTGGGGGC‑G‑G‑CTAAA‑CAGCCCCAG‑CCACACCATCGAAACGCAg < 1:1903708/73‑1 (MQ=255) tggCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑G‑G‑GTTTA‑CAGAGCCGGCCCCCACCA‑CGAAAAGCAggg > 1:1435526/1‑74 (MQ=255) gcgcTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑G‑G‑GTTTA‑CAGCGCCCG‑CCCCCCCCACCAAAAGCAgcccg > 1:2043896/1‑76 (MQ=255) gcgcTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑G‑G‑CGTAA‑AACGGCCAG‑GCCCACCAACCAAACGCcgcgcg > 1:2860476/1‑76 (MQ=255) ggtGGCCCTCGTGGCGGTGGGCGGGCAGGCGGGGC‑G‑G‑CTAAA‑CAGCCCCAG‑CCACACCATCGAAACGCAGcgcgcg < 2:6425317/74‑1 (MQ=255) cTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑G‑G‑GTTAA‑AAGGGCCCA‑CCACAACACACAAACCCAACGCgcgg > 2:5020522/1‑76 (MQ=255) tGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑G‑G‑CTTAA‑CAGCGCCCG‑CCACACCCTCGCAACGCCGCGCgcggc > 1:1537323/1‑76 (MQ=255) ggCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑G‑G‑GTTAA‑CAGCGCCAG‑CCACACCACCGAAACGCAGGGCgcggcg > 1:6888571/1‑76 (MQ=255) ccTCGGTGCGCGGGGCGGGCCGTCGTGGC‑G‑T‑CTAAA‑CGCCCCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTg < 2:1077335/75‑1 (MQ=255) aCTCGTTGCGCTGGGCTGGCAGTCGTGGC‑G‑G‑GTTTA‑CAGCGCCCG‑CCACACCATCGAAACGCCGCGCGCGGCGCTg > 1:6736550/1‑76 (MQ=255) tCGTTGCGCTGGGCTGCAGGTCG‑GGC‑GTA‑GTTAA‑CGGCCCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGaa < 2:1046116/76‑1 (MQ=255) gcGTGGCGCGGGGCGGGCGGTCGGGGG‑G‑G‑CTTAA‑CCCCCCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGaa < 2:863729/75‑1 (MQ=255) gtgtGCGCGGGGCTGGCGGTGGGGGC‑G‑G‑CTAAA‑CGCCCCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGaaa < 2:1521049/73‑1 (MQ=255) gtGCGCGGGGGGGGCTGTGGGGGC‑G‑G‑TTAAC‑CAGCCCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAg < 2:5896771/75‑1 (MQ=255) ggcggTGGGCGGGGGGTGGGGGC‑G‑G‑CTAAA‑CCCCCCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGt < 1:909987/75‑1 (MQ=255) ggcgcGGGGCGGGCGGGCGGGGC‑G‑G‑CTAAC‑CACCCCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGt < 1:3213542/75‑1 (MQ=255) gcgcTGGGCTGGCAGTCGTGGC‑G‑G‑CTTAA‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGGAAAGTa > 1:3340256/1‑76 (MQ=255) ggtggggtggCGGTCGCGGC‑G‑T‑CAAAA‑CGGCCCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTaaa < 2:6162326/74‑1 (MQ=255) tgggctggCAGTCGTGGC‑G‑G‑GTTAA‑CAGCGCCCG‑CCACACCATCCAAAAGCCGCCGGCGGGGCCGGAAAGGaacag > 1:6065845/1‑73 (MQ=255) tgggctggCAGTCGTGGC‑G‑G‑CTTAA‑CAAGGCCAC‑CCCCACCCCCCAAACGCAGCGCGCGGCGCTGGAAAGACAAgc > 2:6838014/1‑76 (MQ=255) gggctgtGAGTCGTGGC‑G‑G‑CTAAA‑CAGCTTTAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAgcg < 2:1005491/76‑1 (MQ=255) gggctggGAGT‑GGGGC‑G‑GCCTAAA‑CAGCCCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAgcg < 2:5622197/76‑1 (MQ=255) gggctggCGGTGGGGGC‑G‑C‑CTAAA‑CAGCCCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAgcg < 2:1031020/76‑1 (MQ=255) gggctggCAGTCGTGGC‑G‑G‑GTTTA‑AAACGCCAG‑CCACACCAACGAAAACCCGCGCGCGGCGGCGGAAAATaaaaag > 2:2407361/1‑73 (MQ=255) ggggCAGTGGGGGC‑G‑G‑CTTAA‑CGGCCCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAg < 2:1671813/74‑1 (MQ=255) ctggCAGTCGTGGC‑G‑G‑GTTAA‑CAGCGCCCG‑CCCCCCCCACCAAAAGGAGGGCGCGGCGCTGGAAAGGAAAGAGCAg > 2:5821447/1‑76 (MQ=255) ctggCAGTCGTGGC‑G‑G‑CTTTA‑AAACGCCCG‑CCCCACCATCGAAAACCCGCGCGGCGCGCTGCGAAATAAAGCGCAg > 2:5201028/1‑76 (MQ=255) tggCAGTCGTGGC‑G‑G‑GCTAA‑CAGCGGCCC‑ACCCCCCATCGAAACGCACCGCGCGGCGGCGAAAAGTAAAACGCAgg > 2:2167045/1‑76 (MQ=255) tggCAGGGGTGGC‑G‑G‑CTAAC‑CAGCCCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAgg < 2:6170418/76‑1 (MQ=255) ggCGGTGGTGGC‑G‑T‑CTTAC‑CAGCCCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGa < 2:826508/76‑1 (MQ=255) cAGTCGTGGC‑G‑G‑CTTAA‑CAGCGCCAG‑CTACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAAc < 1:4587318/76‑1 (MQ=255) cAGTCGTGGC‑C‑T‑CTAAA‑CCCCCCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAAc < 2:5011720/76‑1 (MQ=255) aGTCGTGGC‑G‑G‑CTTTA‑AAGCGCCAG‑CCACCCCCTCGAAAAGCAGCGCGCGGCGCGGAAAAGAAAAGAGCCGGAAcc > 2:4932830/1‑75 (MQ=255) gcgggc‑g‑g‑cCTAACCGGCCCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGa < 2:1278065/73‑1 (MQ=255) cGTGGC‑G‑G‑CTTAA‑CAGCGCCCG‑CCCCCCCAACGAAACGCAGCGCGCGGCGCTGCAAAGGAAAGGCCAGGGACttgc > 2:5041527/1‑73 (MQ=255) cGTGGC‑G‑G‑CTTAA‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGAc > 1:373040/1‑76 (MQ=255) ggc‑g‑g‑cTTAA‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACgaa > 2:1384842/1‑76 (MQ=255) ggc‑g‑g‑cTTAA‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACgaa > 2:3514381/1‑76 (MQ=255) ggc‑g‑g‑cTTAA‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACgaa > 1:1336332/1‑76 (MQ=255) ggc‑g‑g‑cTTAA‑CAGCGACAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACgaa > 1:2631580/1‑76 (MQ=255) gc‑g‑g‑cTTAA‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACgaag > 2:568296/1‑76 (MQ=255) g‑cTTAA‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGaaa > 2:852152/1‑76 (MQ=255) g‑cTTAA‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGaaa > 1:1282941/1‑76 (MQ=255) cTTAA‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaa > 1:1672410/1‑76 (MQ=255) cTTAA‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGATGAAAAGTAAAGCGCAGGAACTGACGAAGaaaa > 2:6850218/1‑76 (MQ=255) ttAA‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaaa > 2:4569027/1‑76 (MQ=255) ttAA‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaaa > 1:2476318/1‑76 (MQ=255) taaa‑CCGCCCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaaa < 1:6838014/74‑1 (MQ=255) tAA‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAc > 2:1604276/1‑76 (MQ=255) tAA‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAc > 1:1409370/1‑76 (MQ=255) caa‑CACCCCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAc < 2:1291541/75‑1 (MQ=255) aa‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAca > 2:1138499/1‑76 (MQ=255) aa‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAca > 2:5470578/1‑76 (MQ=255) aa‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAca > 1:5092399/1‑76 (MQ=255) aa‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAca > 1:4611489/1‑76 (MQ=255) a‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAcag < 2:5841124/76‑1 (MQ=255) cAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAcagc > 2:3945441/1‑76 (MQ=255) aGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAcagcc > 2:1096711/1‑76 (MQ=255) aGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAcagcc > 2:6629565/1‑76 (MQ=255) gcgcCAG‑ACACACCATCGAAACGCAGCGCGCGGCGCTGAAAAATAAAGCGCAGGAACTGACGAAGAAAAAcagcca < 1:4934270/76‑1 (MQ=255) | ATGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑G‑G‑CTTAA‑CAGCGCCAG‑CCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAACAGCCA > REL606/895781‑895929 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |