Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,983,700 | 0 | T | G | 27.8% | 69.6 / 15.4 | 36 | T397P (ACC→CCC) | ECB_03713 | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (5/21); new base G (10/0); total (15/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.18e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATT > REL606/3983625‑3983773 | gAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGt < 1:974602/76‑1 (MQ=255) tgtAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCt < 1:805190/76‑1 (MQ=255) aTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCcac > 2:133413/1‑76 (MQ=255) aTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTccccc > 1:1209491/1‑74 (MQ=255) gggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAAta < 2:862089/76‑2 (MQ=255) gTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTg < 1:16241/76‑1 (MQ=255) gTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCCCCCCCAAATTTg > 1:886905/1‑76 (MQ=255) tACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCCAAATTTTAc > 2:787809/1‑76 (MQ=255) accaccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACg < 1:1030750/76‑1 (MQ=255) accaccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCACAAAATTTACg > 2:1138656/1‑76 (MQ=255) accaccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCAGCAGGGGGGGCCCCCCCAAAATTTACg > 2:205907/1‑76 (MQ=255) caccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGAt < 2:612960/76‑1 (MQ=255) caccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGAt < 2:15086/76‑1 (MQ=255) gggTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTgg < 1:553190/76‑1 (MQ=255) gggTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTgg < 1:617460/76‑1 (MQ=255) tGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGgaga < 2:197065/76‑1 (MQ=255) tGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCCAAAATTTACCATCGCCggggg > 2:282525/1‑73 (MQ=255) gATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGAt < 2:823756/76‑1 (MQ=255) gATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGAt < 1:836863/76‑1 (MQ=255) ttgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGa < 1:875081/76‑1 (MQ=255) aTCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATcg < 1:432548/76‑1 (MQ=255) cTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCCAAAATTTACGATCGCTGGGGATGACGCCCAcgcgg > 1:258614/1‑75 (MQ=255) cTACGGCTGTTATTGGTGGCAACGGCAGGGGGGGCTCCCACAAAATTGAAGAAAGCTGGGGAAGAAGACAACGCGa > 1:123456/1‑76 (MQ=255) tGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccacc > 2:246039/1‑76 (MQ=255) tGTTATTGGTGGCATCGGCAGGGGGGCCTCCCACAAAATTGACGATCGCTGGAGATGAAGCCATCGCGGAGccccc > 2:842770/1‑76 (MQ=255) gTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCAAAATTGACGATCGCTGGAGATGACGCCCTCGCGCAGccaccc > 2:269901/1‑76 (MQ=255) ttATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccaccca < 2:626310/76‑1 (MQ=255) ttATTGGTGGCATCGGCAGGGGGGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccaccca > 1:492532/1‑76 (MQ=255) ttATTGGTGGCATCGGCAGGGGGGCCTCCCACAAAATTGACGATCGCTGGAGAAGACGCCATCGCGAAGccacccc > 2:1032188/1‑75 (MQ=255) tATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCCAAGccaccccc > 2:503886/1‑74 (MQ=255) tATTGGTGGCATCGGCAGGGGGGCCTCCCACAAAATTGACGATCGCTGGAGGTGACGCCCTCGCGAAGacacccac > 2:740796/1‑76 (MQ=255) aTTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCcaca < 2:1210377/76‑1 (MQ=255) ttGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCacac < 1:1176721/76‑1 (MQ=255) ttGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCacac < 2:268228/76‑1 (MQ=255) ttGGTGGCATCGGCAGGGGGGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCacac > 1:405236/1‑76 (MQ=255) ggtggCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCt < 2:452357/76‑1 (MQ=255) gtggCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTg < 2:123366/76‑1 (MQ=255) tCGGCA‑GGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTc < 1:47550/76‑1 (MQ=255) cGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTc > 2:808465/1‑76 (MQ=255) cGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGGTTc > 2:808475/1‑76 (MQ=255) cAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGt > 2:45554/1‑76 (MQ=255) ggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAAtt < 1:643145/76‑1 (MQ=255) | GAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATT > REL606/3983625‑3983773 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |