Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 571,979 | 0 | C | G | 42.9% | 13.4 / 10.6 | 21 | R20R (CGG→CGC) A504P (GCA→CCA) | ECB_00529 ECB_00530 | hypothetical protein conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (7/5); new base G (9/0); total (16/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.51e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGACGGCCCGTCAGGTCGCTCACCGCATGATTCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGCCGTATCCATGTTCTTCTGCGCGTGGATGTAGACCTGTTCTCTCCCTGTCGCATCATCAAACTTCAGCTCATTAAA > REL606/571905‑572054 | ggACGGCCCGTCAGGTCGCTCACCGCATGATTCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGcc > 2:3726695/1‑76 (MQ=255) gACGGCCCGTCAGGTCGCTCACCGCATGATTCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGcgg > 2:300508/1‑74 (MQ=255) ggCCCGTCAGGTCGCTCACCGCATGATTCGTCCGCGACACCCACTGCGGGGCCACGGAAGGATGCAGATGCCgcag > 2:922825/1‑73 (MQ=255) gCCCGTCAGGTCGCTCACCGCATGATTCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGGCGGATc > 1:2507135/1‑76 (MQ=255) ccGTCAGGTCGCTCACCGCATGATTCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGCCGGATcac > 1:1387806/1‑74 (MQ=255) tCAGGTCGCTCACCGCATGATTCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGGCGGGTCCCTGt > 2:755402/1‑76 (MQ=255) aGGGCGCTCACCGCATGATTCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGGCGTATCCCTGttc > 1:1968996/1‑76 (MQ=255) ggTCGCTCACCGCATGATTCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGGGGGATCCCTGgtct > 2:2797203/1‑76 (MQ=255) gCTCACCGCATGATTCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGGCGGATCCCTGTtcttctt > 2:1070551/1‑75 (MQ=255) gCTCACCGCATGATTCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGCCGTATCCATGTtcttctt > 2:1314620/1‑75 (MQ=255) cTCACCGCATGATTCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGGCGTATCCCTTTTCTTTTgt > 2:1328071/1‑75 (MQ=255) cACCGCATGATTCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGGCGGGTCCATCTTCTTTTtcgc > 2:4258633/1‑76 (MQ=255) cACCGCATGATTCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGGCGGATCCCTGTTCTTCTgcgc > 2:544599/1‑76 (MQ=255) cACCGCATGATTCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGGCGCATCCCTGTTCTTTTgtgc > 2:1645606/1‑76 (MQ=255) cACCGCATGATTCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGGCCGGTCCACGTTCTTCTccgc > 2:1306134/1‑76 (MQ=255) cACCGCATGATTCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGGCCGATTCATTTTCTTTTgcgc > 2:587142/1‑76 (MQ=255) cATGATTCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGGCGGATCCATGTTCTTTTTCGCGCggg > 2:557423/1‑75 (MQ=255) aTGATTCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGGCCTATTCATGTTTTTCTGCTCGTgggc > 1:4051357/1‑74 (MQ=255) aTGATTCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGGCCTATCCCTGTTCTTCTGCGCgggggg > 1:1426842/1‑74 (MQ=255) tCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGGCGTATCCATGTTTTTTTGCGCGTGGGTGTAGa > 1:31638/1‑76 (MQ=255) cGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGGCCCATCCATGTTATTCTTCTCGCGGATCTAAAc > 2:2245413/1‑76 (MQ=255) gcccGCGTCCCCCCCGGGGGGGCCAGGAACGGATGCGCCTGCCGTATCCATGTTCTTCTGCGCGTGGATGTAGAcc < 1:380867/74‑1 (MQ=21) gTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGGCGGATCCATGTTTTTCTGCGCGTGGATGGAGAcc > 1:1278662/1‑76 (MQ=38) gcgTCACCCACTGCGGGGCCACGGACGGATGCGGCTGCCGTATCCATGTTCTTCTGCGCGGGGATGTAGACCCGtt > 1:3291341/1‑76 (MQ=255) cgTCACCCACTGCGGGGCCACGGACGGATGCGGCTGGCGGATCCCTGTTCTTCTTCGCGCGGGGGTAGACCTGttt > 2:3351710/1‑75 (MQ=255) cgTCACCCACTGCGGGGCCACGGACGGATGCGGCTGCCGGATCCCCATTCGTCTGCGCGGGGGTGTAGACCTGtgc > 2:1077689/1‑74 (MQ=255) cACCCACTGCGGGGCCACGGACGGATGCGGCTGGCGGATCCATGTTCTTCTGCGCGCGGATGTAGACCTGTtctcc > 1:1854629/1‑75 (MQ=37) cTGCGGGGCCACGGACGGATGCGGCTGGCGGATCCCTGTTTTTTTGGGCGGGGATGTAGAACTGGTCTCTCCCCGt > 2:2194406/1‑76 (MQ=255) cTGCGGGGCCACGGACGGATGCGGCTGGCGGATCCATGTTCTTTTGCGCGTGGATGTAGACCTGTTCTCTCCCTGt > 2:1267660/1‑76 (MQ=25) cTGCGGGGCCACGGACGGATGCGGCTGCCGTATCCATGTTCTTCTGCGCGTGGGTGTAGACCTGTTTTCTCCCTGt > 1:3158350/1‑76 (MQ=255) tGCGGGGCCACGGACGGATGCGGCTGGCGTTTCCCTGTTCTTCTGCCGGCGGATGTAGACCTGTTTCCTCCCttgc > 1:2186741/1‑73 (MQ=25) tGCGGGGCCACGGACGGATGCGGCTGCCGTATCCCTGTTTTTCTGCGCGTGGGTGTAGACCTGTTTTCTCCCTGTc > 1:4313775/1‑76 (MQ=255) tGCGGGGCCACGGACGGATGCGGCTGCCGGAACCGTCATATTTTGCGCGGCGAGGGCGGCGTGTTCTCTCCCCGTc > 1:2099589/1‑76 (MQ=255) cGGGGCCACGGACGGATGCGGCTGCCGTATCCATGTTCTTCTGCGCGTGGGTGTAGACCTGTTCTCTCCCTGTCGc > 1:3433752/1‑76 (MQ=255) ggggCCACGGACGGATGCGGCTGGCCCAACCATGTCCTTCTGCGCGTGGGTGGGGAACTGTTCTTTTCCCGTCGca > 2:725322/1‑76 (MQ=12) gcccACGGAGTGCTCCGCCTGCCGTATCCATGTTCTTCTGCGCGTGGATGTAGACCTGTTCTCTCCCTGTCGcatc < 1:1706661/74‑1 (MQ=17) acggacggATGCGGCTGGCGGATCCCTGTTTTTTTTCGCGTGGGTGTAGACCCGTTCTCTCCCCGTCGCCtcatca > 2:1374076/1‑76 (MQ=255) gggaGGCGCCTGCCGTATCCATGTTCTTCTGCGCGTGGATGTAGACCTGTTCTCTCCCTGTCGCATCATCAAACtt < 1:3885028/75‑1 (MQ=255) ggaTGCGGCTGGCGTATCCATGTTTTTCTGGGCGGGGATGTGGACCCGTTTTCTCCCTGTCGCCTCATCAAAATTc > 2:4877827/1‑76 (MQ=255) tGCGGCTGGCGTATCCATGTTCTTCTGCGCGTGGATGTAGACCTGTTCTCTCCCTGTCGCATCATCAAACTTCAGc > 2:231409/1‑76 (MQ=255) cGGCTGCCGTATCCATGTTCTTCTGCGCGTGGATGTAGACCTGTTCTCTCCCTGTCGCATCATCAAACTTCAGCTc > 2:2603124/1‑76 (MQ=255) cGGCTGCCGTATCCATGTTCTTCTGCGCGTGGATGTAGACCTGTTCTCTCCCTGTCGCATCATCAAACTTCAGCTc > 1:1433795/1‑76 (MQ=255) gCTGCCGTATCCATGTTCTTCTGCGCGTGGATGTAGACCTGTTCTCTCCCTGTCGCATCATCAAACTTCAGCTCAt > 1:3174437/1‑76 (MQ=255) ccTGCCGTATCCATGTTCTTCTGCGCGTGGATGTAGACCTGTTCTCTCCCTGTCGCATCATCAAACTTCAGCTCAt < 2:116840/75‑1 (MQ=255) tGCCGTATCCATGTTCTTCTGCGCGTGGATGTAGACCTGTTCTCTCCCTGTCGCATCATCAAACTTCAGCTCATTa < 2:2277438/76‑1 (MQ=255) tGCCGTATCCATGTTCTTCTGCGCGTGGATGTAGACCTGTTCTCTCCCTGTCGCATCATCAAACTTCAGCTCATTa < 1:5005068/76‑1 (MQ=255) ccGTATCCATGTTCTTCTGCGCGTGGATGTAGACCTGTTCTCTCCCTGTCGCATCATCAAACTTCAGCTCATTaaa < 2:3835091/76‑1 (MQ=255) ccGTATCCATGTTCTTCTGCGCGTGGATGTAGACCTGTTCTCTCCCTGTCGCATCATCAAACTTCAGCTCATTaaa > 2:1982071/1‑76 (MQ=255) | GGACGGCCCGTCAGGTCGCTCACCGCATGATTCGTCCGCGTCACCCACTGCGGGGCCACGGACGGATGCGGCTGCCGTATCCATGTTCTTCTGCGCGTGGATGTAGACCTGTTCTCTCCCTGTCGCATCATCAAACTTCAGCTCATTAAA > REL606/571905‑572054 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |