Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,983,700 | 0 | T | G | 31.6% | 71.6 / 10.7 | 39 | T397P (ACC→CCC) | ECB_03713 | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (6/20); new base G (12/0); total (19/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.86e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTGC > REL606/3983625‑3983775 | gAAGTTTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGt < 1:152771/76‑1 (MQ=255) cGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACaaa < 2:405801/76‑1 (MQ=255) cGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCCCCCCCaaa > 2:936556/1‑76 (MQ=255) ggggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAAt < 2:1155178/76‑1 (MQ=255) ggggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCCCAAt > 2:1893393/1‑76 (MQ=255) ggggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCCCCCCCaaaa > 2:705559/1‑75 (MQ=255) ggggTTACCACCGAGGGAGATTGCTTGCATCCTACGGCTGTTATTGGTGGCGTCGGCAGGGGGGGGTCCCACaaaa > 2:158670/1‑75 (MQ=255) gggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAAtt < 1:1238645/76‑1 (MQ=255) ttACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGa < 2:789539/76‑1 (MQ=255) tACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCACAAATTTGAc > 2:528127/1‑76 (MQ=255) accaccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCCAAAATTTACg > 1:1622474/1‑76 (MQ=255) ccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCAAAATTTACGATCg > 1:1024662/1‑76 (MQ=255) ccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCACAAATTTTACGATCg > 1:1787610/1‑76 (MQ=255) cangggTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCCCCCCCAAAATTTACGATCGc > 2:1582285/4‑76 (MQ=255) gTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGgag < 2:1432441/76‑1 (MQ=255) ttgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGa < 2:1558382/76‑1 (MQ=255) ttgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGa < 2:772151/76‑1 (MQ=255) tgcttgcATCCTACGGCTGTTATTGGCGGCATCGGCAGGGGGGGCCCCCCCAAAATTTACGAACGCCGGGGGTGAc > 2:1324677/1‑76 (MQ=255) cttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCAAAATTTACCATCCCCGGAGAATACGc > 1:1964495/1‑76 (MQ=255) tgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCa < 1:1708129/76‑1 (MQ=255) cATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATc < 2:1838266/76‑1 (MQ=255) tCCTGCGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATcgc < 1:187130/76‑1 (MQ=255) cTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGa < 1:699828/76‑1 (MQ=255) cTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGa < 1:2058333/76‑1 (MQ=255) cTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCAAATTTGACGATCGCTGGAGATGACGCCATCGCGa > 2:1681951/1‑76 (MQ=255) cGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCCCAAAATTTACCACCGCTGGAGATGCAGCCACCGCGAAGc > 1:1386632/1‑76 (MQ=255) cTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCACAAATTTTACGATCGCTGGAGATGACGCCATCGCGGAGccac > 1:406815/1‑76 (MQ=255) tGTTATTGGTGGCATCGGCAGGGGGGCCTCCCACAAATTTTACGATCGCTGGGGATGAAGCCAACGCGGAGccacc > 2:989037/1‑76 (MQ=255) gTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCAAATTTGACGAACGCTGGAGATGACGCCATCGCGGAGccaccc > 1:1080411/1‑76 (MQ=255) ttATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccaccca < 1:1829563/76‑1 (MQ=255) ttATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccaccca < 1:2204483/76‑1 (MQ=255) tATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccacccac < 2:441588/76‑1 (MQ=255) ttGGTGGCATCGGCAGGGGGGCCTCCCACAAATTTGAAGATCGCTGGGGATGACGCCATCGCGGAGCCAACCacac > 1:1879400/1‑76 (MQ=255) gtggGATGGGCGGGGGTGTCGTCCCCAAAGTGGGCGATCGCTGGAGATGACGCCATCGCGAAGCCACCCCCAccgg < 2:324927/76‑3 (MQ=255) tggCATCGGCAGGGGAGCCTCCCACAAATTTGACGAACGCTTGAGATGACCCCATCGCGAAGCCACCCACCCCTGc > 2:2033379/1‑76 (MQ=255) gCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCCGCgg > 2:88631/1‑76 (MQ=255) cATCGGTAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGt > 2:816308/1‑76 (MQ=255) tCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGttt < 1:1775909/76‑1 (MQ=255) cGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTc < 2:1051480/76‑1 (MQ=255) ggCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTcg > 1:1007125/1‑75 (MQ=255) acnggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTa > 2:1592450/4‑76 (MQ=255) aGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTa > 1:509059/1‑76 (MQ=255) ggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAAtt < 2:1879400/76‑1 (MQ=255) ggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAAtt < 1:2250265/76‑1 (MQ=255) ggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAAtt < 1:1977803/76‑1 (MQ=255) tGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTgc > 2:1788838/1‑76 (MQ=255) tGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTgc > 1:1677134/1‑76 (MQ=255) | GAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTGC > REL606/3983625‑3983775 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |