Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 895,856 | 0 | G | C | 36.0% | 34.7 / 19.5 | 25 | A26P (GCC→CCC) | ECB_00836 | putative Regulatory protein |
Reads supporting (aligned to +/- strand): ref base G (14/2); new base C (0/9); total (14/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.69e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGCTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAACAGC > REL606/895781‑895927 | aTGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGCTTAACAgcg < 1:4301117/76‑1 (MQ=255) aTGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGCTTAACAgcg < 1:570242/76‑1 (MQ=255) aTGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGCTTAACAgcg < 1:1118428/76‑1 (MQ=255) tGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGCTTAACAgcgc < 2:6878055/76‑1 (MQ=255) gACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGCTTAACAGCgcc < 1:2444987/76‑1 (MQ=255) aCGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGCTTAACAGCgcca < 1:2618797/76‑1 (MQ=255) cgTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGCTTAACAGCgccagcc < 1:4192337/76‑1 (MQ=255) gTGCGCTGGCGGTAGGGGTGGCCCGGGCACTCGTTGCGGTGGGCTGGCAGTGGTGGC‑GGCTTAACAGCgccagcca < 2:2798849/76‑1 (MQ=255) gTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCGGTCGTGGC‑GGCTTAACAGCgccagcca < 1:1525177/76‑1 (MQ=255) gTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGCTTAACAGCgccagcca < 1:7053990/76‑1 (MQ=255) gTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGCTTAACAGCgccagcca < 1:4692034/76‑1 (MQ=255) cgcTGGCAGTAGTGGTGGCGCTGGCCCTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGCTTAACAGCGCCAGCCacac < 1:2239360/76‑1 (MQ=255) cgcTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGCTTAACAGCGCCAGCCacac < 1:7120021/76‑1 (MQ=255) cTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGG‑GGGGTGAAAGCGCCAGACCCAccc > 2:2281928/1‑75 (MQ=255) ggtggggGCGCGGGCCCTGGTGGGGTGGGGCGGGGGGTGGGGGC‑GGCTAACCAGCCCCAGCCACACCATCGAAACg < 1:6991037/75‑1 (MQ=255) agtGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGGGTTACAACGCCCGCCCCACCCTCGCAAcc > 1:5498246/1‑75 (MQ=255) tggtggCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGGGTTAAAACCCCCGCCCCACCCTCGAAAAGCa > 1:4663695/1‑76 (MQ=255) gtgggCGCGGGCCCTGGTTGCGCGGGGCTGCCAGTGGTGGC‑GGCTAACCACCCCCAGCCACACCATCGAAACGCAg < 2:4329536/73‑1 (MQ=255) ggttgCGCGGGCCCTCGTTGCGCGGGGCGGGCAGTCGGGGC‑GGCTTAACGGCCCCAGCCACACCATCGAAACGCAg < 1:7814245/76‑1 (MQ=255) ggggCGCGGCCCCTGGTGGCGCGGGGGGGGCGGTGGGGGC‑GGCTAACCGGCCCCAGCCACACCATCGAAACGCAgc < 1:4635135/74‑1 (MQ=255) ggCGCGGGCCCTCGTGGCGGGGGGCGGGGGGTGGGGGC‑GTCTAAACGCCCCCAGCCACACCATCGAAACGCAgcgc < 1:6190093/76‑1 (MQ=255) ggtGGCCCTGGTTGCGGGGGGCTGGGGGGGGGGGC‑GGTTAAACAGCCCCAGCCACACCATCGAAACGCAGcgtgcg < 1:1606230/74‑1 (MQ=255) gcTGGCCCTAGTTGGGCGGGGCGGACAGTCGTGGG‑CGATAAACAGCCCCAGCCACACCATCGAAACGCAGcgcgcg < 1:1777875/76‑1 (MQ=255) gcTGGCCCGGGTTGCGGGGGGCTGGCTGTCGTGGC‑GCTTAAACGCCCCCAGCCACACCATCGAAACGCAGcgcgcg < 1:321635/76‑1 (MQ=255) gcTGCCACTCGTCGCGCGGGGGGGGCGGTGGGGGC‑GGCTTACCCCCCCCAGCCACACCATCGAAACGCAGcgcgcg < 1:5264043/76‑1 (MQ=255) cTGGCACTGGTTGGGCGGGGCGGGCGGGGGGGGC‑GGCTAACCCCCCCCAGCCACACCATCGAAACGCAGCGCgcgg < 2:743918/76‑1 (MQ=255) tGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGGTTAACAGCGGCAGCCACACCATCCAAACGCAGCGCgcgg‑c > 1:7756555/1‑76 (MQ=255) gCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGGTTAACAGCGCCCGCCACACCATCGAAACGCAGAGCGCGg‑cg‑c > 1:1454038/1‑76 (MQ=255) cACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGGGTTACCAG‑GCCCGCCCCACCCTCGCAACGCAGCGCGCGg‑gg‑cg > 2:6412490/1‑75 (MQ=255) cACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGGTTAACAACGCCCGCCACACCATCCAAACGCCGCGCCCGG‑GG‑Ct > 1:2720577/1‑76 (MQ=255) aCTCGTTGCGCTGGGCTGGCAGTCGTGGG‑GGGCTTACACCGCCAGCCCCACCACCGAAACGCAGCGCGCGg‑cg‑ccg > 1:6718542/1‑74 (MQ=255) cTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGGTTTAAAGCGCCAGCCACAACCTCCAAAAGCAGCGCGCGG‑CG‑CTgg > 2:4658897/1‑75 (MQ=255) tCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGGTTTACAGCGCCAGCCCCCCCATCGAAACGCAGCGCGCGG‑GG‑CTgca > 1:2018252/1‑74 (MQ=255) tCGGTGCGGGGGGGTGGCAGGGGGGGC‑GGCTAAACCCCCCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGaa < 1:3707732/76‑1 (MQ=255) cGTTGCGCTGGGCTGGCAGTCGTGGCGGGTTTAAAAGGCCCACCACCACCACCCAAA‑GCAGCGCGCGG‑CG‑CTGgaa > 2:6264791/1‑76 (MQ=255) gTTGCGCTGGGCTGGCAGTCGTGGC‑GGCTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGG‑GG‑CTGgaaa > 1:3423658/1‑76 (MQ=255) ttGCGCTGGGCTGGCAGTCGTGGC‑GGCTTTAAAACGCCCACCACACCCTTCAAAAGCACCGGGCGG‑GG‑CCGgaaaa > 1:1502060/1‑75 (MQ=255) gtGCGCTGGGCGGGCAGGGGGGGC‑GGCTAAACACCCCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAg < 1:5970545/75‑1 (MQ=255) tGCGCTGGGCTGGCAGTCGTGGC‑GGCTTAACAGCGCCCGCCACACCCTCGAAACCCAGGGCGCGG‑CG‑CCGCAAAgg > 1:3249693/1‑75 (MQ=255) tGCGCTGGGCTGGCAGTCGTGGC‑GGCTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGGAAAGt > 1:5921732/1‑76 (MQ=255) tGCGCTGGGCTGGCAGTCGTGGC‑GGCTTAACACCGCCCGCCACACCAACTAAACGCAGCGCGCGGCCG‑CGGAAAAg > 1:6231680/1‑76 (MQ=255) gcgcTGGGCTGGCAGTCGTGGC‑GGGTTAACAGCCCCCGCCCCACCACCGCAACGCAGCGCGCGG‑CG‑CGGACAAgaa > 1:121938/1‑74 (MQ=255) ctgggctggCAGTCGTGGC‑GGGTTTAAAACGCCAGCCCCCCCCCCCAAAAGCAGCGCGCCG‑GG‑CCGCAAAGaaaaa > 2:5070884/1‑75 (MQ=255) tgggctggCAGTCGTGGC‑GGGTTTACAAAGCCAGCCCCACCCCAGAAACGCAGCGCGCGG‑GG‑CCGGTAAGTAAAgc > 2:5778958/1‑76 (MQ=255) tgggctggCAGTCGTGGC‑GGGCTTACAGGGCCCGCCACCCCAACGACAAGCCACGCGCGG‑CGCCGgaaaagaaaac > 2:5501881/1‑75 (MQ=255) ggggctggCTGTGGGGGG‑GGCTAACCCCCCCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAgc < 1:2919305/75‑1 (MQ=255) gggctggCAGTCGTGGC‑GGCTTAACAGCGCCAGCCACACCCTCCAAAAGCAGCGCGCGG‑CG‑CGGGAAAGGAAAacg > 1:3019022/1‑76 (MQ=255) gggcgggcAGTGGGGGC‑GGCTAAACCCCCCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAgcg < 1:6961813/76‑1 (MQ=255) gctggCGGTGGGGGC‑GGCTTAACCCCCCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCa < 2:7366380/76‑1 (MQ=255) cgggCAGGCGGGGC‑GGCTAAACCCCCCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAg < 1:5734875/74‑1 (MQ=255) gggCGGTCGGGGC‑GCCTTACCCGCCCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAgg < 1:1598272/75‑1 (MQ=255) gCAGTCGTGGC‑GGGTTAAAAACCCCAGCCCCACCCACGGAAAGCAGCGCGCGG‑CG‑CCGGAAAATAAAAAGCCGGaa > 2:5042339/1‑76 (MQ=255) cAGTCGTGGC‑GGGTTTAAAGCGCCCGCCACACCACCCGAAAGCAGCGCGCCG‑CG‑CCGGAAAAAAAAACGCAGGaaa > 2:41124/1‑75 (MQ=255) cAGTCGTGGC‑GGCTTAACAGGGCCCCCCCCACCAACGCAAAGCAGCGCGCGG‑GC‑CTGTAAAGTAAAACGCCGGAAc > 2:7977004/1‑76 (MQ=255) cAGTCGTGGC‑GGCTTAACAGCGCCAGCCCCACCACCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGAACAAAGCAGGaga > 2:4302842/1‑74 (MQ=255) aGTCGTGGC‑GGGTTTAAAGCGCCCGCCCCACCCCCGAAACGCAGCGCGCGG‑CG‑CTGAAAACAAAAGAGCAGGaaat > 2:2031722/1‑74 (MQ=255) gTCGTTGC‑GGCTTTTTAGCGCCAGCCACACCATCGAAAAGCAGCGCGCGG‑CG‑CTGGAAAGTAAAACGCAGGGAccg > 2:173394/1‑74 (MQ=255) gTCGTGGC‑GGGTTTACAGCGCCCGCCCCCCCATCCAAAAGCAGCGCGCGG‑CG‑CTGAAAAGAAAAACGCAGGAAcag > 2:259980/1‑74 (MQ=255) gTCGTGGC‑GGGTTTACAACGCCAGCCCCACCATCGAAACACAGCGCGCGG‑GG‑CCGGAAAGGAAAAAGCAGGAAcag > 2:4054601/1‑74 (MQ=255) gTCGTGGC‑GGGTTAACAGAGCCAGCCCCACCCTCCAAACGCAGCGCGCGG‑CG‑CTGGAAAGTAAAGCGGAGGGACTg > 2:3369061/1‑76 (MQ=255) cGTGGC‑GGGTTAACCGCGCCAGCCACACCACCGAAAAGCAGCGCGCGG‑CG‑CTGAAAAGCAAAGAGCAGGAACTggc > 2:2894511/1‑74 (MQ=255) cGTGGC‑GGCTTTACAGCGCCCGCCCCACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGAc > 2:7354989/1‑76 (MQ=255) cGTGGC‑GGCTTAACAGCGCCCGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGGAAAGGAAAACGCAGGAAATGAc > 2:3922755/1‑76 (MQ=255) ggggc‑ggcTAAACCCCCCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACg < 1:3961676/74‑1 (MQ=255) gTGGC‑GGGTTTACAGCGCCCGCCACACCATCGAAAACCAGCGCGCCG‑CG‑GTGGAAAAAAAAACGCAGGGACTGACg > 2:3458548/1‑76 (MQ=255) tGGC‑GGCTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACga > 1:907277/1‑76 (MQ=255) tGGC‑GGCTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACga > 2:2017357/1‑76 (MQ=255) tGGC‑GGCTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCCG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACga > 2:1979081/1‑76 (MQ=255) ggc‑ggcTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACgaa > 2:3929213/1‑76 (MQ=255) gc‑ggcTAACCACCCCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACgaag < 1:7066826/76‑1 (MQ=255) c‑ggcTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACgaaga > 1:3873238/1‑76 (MQ=255) ccTTACCACCGCAAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACGAAGaaa < 1:474720/75‑1 (MQ=255) ttaaaCAGCCCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaa < 2:5348196/73‑1 (MQ=255) cTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAACAGAAAAGCGCAGGAACTGACGAAGaaaa > 2:1939545/1‑76 (MQ=255) cTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaa > 1:7613698/1‑76 (MQ=255) ttAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaaa > 1:75942/1‑76 (MQ=255) ttAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaaa > 2:4743263/1‑76 (MQ=255) ttAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaaa > 2:6639940/1‑76 (MQ=255) tAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAc > 1:1591238/1‑76 (MQ=255) aaCAGCGCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAca > 2:2713753/1‑76 (MQ=255) aaCAGCGCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAca > 1:1943529/1‑76 (MQ=255) aaCACCCCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAca < 2:2851941/76‑1 (MQ=255) aaCACCCCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAca < 2:2918345/76‑1 (MQ=255) gcccgcCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAcagc < 1:3183742/73‑1 (MQ=255) | ATGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGC‑GGCTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGG‑CG‑CTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAACAGC > REL606/895781‑895927 |
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: ‑ |