Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,103,887 | 4 | . | C | 50.0% | ‑1.5 / 39.6 | 14 | coding (154/216 nt) | ECB_01992 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base . (7/0); new base C (7/0); total (14/0) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAAACATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑CAG‑‑‑‑‑CCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAAAT > REL606/2103812‑2103947 | aaaaCATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTca < 1:5890612/76‑1 (MQ=255) aaaaCATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTca < 1:5747872/76‑1 (MQ=255) aaaaCATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTca > 1:4758898/1‑76 (MQ=255) aaaaCATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTca < 1:3954854/76‑1 (MQ=255) aaaaCATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTca > 2:2837404/1‑76 (MQ=255) aaCATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTca‑‑‑‑ca < 1:59631/76‑1 (MQ=255) aaCATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTACATTTca‑‑‑‑ca > 1:2313694/1‑76 (MQ=255) aCATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag < 2:6012195/76‑1 (MQ=255) cATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑c < 2:1925421/76‑1 (MQ=255) cATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑c > 1:7772345/1‑76 (MQ=255) aTGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑cc < 1:2050482/76‑1 (MQ=255) aTGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑cc < 1:5809380/76‑1 (MQ=255) aTGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑cc < 1:2639174/76‑1 (MQ=255) tGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑cca < 1:4199790/76‑1 (MQ=255) tGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑cca < 2:573511/76‑1 (MQ=255) tGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑cca > 1:561784/1‑76 (MQ=255) tGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑cca > 1:7781897/1‑76 (MQ=255) tGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑cca < 2:2571866/76‑1 (MQ=255) ggCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccag < 2:5308889/76‑1 (MQ=255) ggCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccag < 1:870918/76‑1 (MQ=255) ggCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccag < 1:5405752/76‑1 (MQ=255) gCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagc < 1:7828993/76‑1 (MQ=255) gCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagc > 2:3026587/1‑76 (MQ=255) gCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagc > 2:1499639/1‑76 (MQ=255) cTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagcc > 1:1466159/1‑76 (MQ=255) cTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagcc < 1:6881815/76‑1 (MQ=255) cTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagcc > 2:2390280/1‑76 (MQ=255) cTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagcc > 2:47443/1‑76 (MQ=255) tataCAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagcca < 2:5892223/76‑1 (MQ=255) tataCAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagcca < 2:910397/76‑1 (MQ=255) tataCAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATGTCA‑‑‑‑cag‑‑‑‑‑ccagcca < 2:6127315/76‑1 (MQ=255) ataCAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccag < 1:6428056/76‑1 (MQ=255) ataCAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccag < 2:5005874/76‑1 (MQ=255) aCAATGAAATCTGCTTCTTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGTTTTCC‑‑‑‑cag‑‑‑‑‑ccagccagcc < 2:3619108/76‑1 (MQ=255) aCAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagcc > 2:2286287/1‑76 (MQ=255) aCAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagcc > 2:6014110/1‑76 (MQ=255) aCAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagac > 1:1747232/1‑74 (MQ=255) aCAATGAAAACTGATTCGTTATTGCCTTCGAATATAAAGGGGTGCACTTTGCG‑TTAGATTTAA‑‑‑‑cag‑‑‑‑‑ccagcgagca > 2:2286288/1‑75 (MQ=255) cAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagcca > 1:5955464/1‑76 (MQ=255) cAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagcca > 1:2910600/1‑76 (MQ=255) cAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagcca < 2:2879970/76‑1 (MQ=255) cAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagcca > 1:4306365/1‑76 (MQ=255) aaTGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccag < 2:3338317/76‑1 (MQ=255) aaTGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccag < 2:1847587/76‑1 (MQ=255) ataaaacCTGCTTCGTTATTGACTTCGAATATAAGGGGGTGAAATTTGCA‑TTAGTTTTCC‑‑‑‑cag‑‑‑‑‑ccagccagccagc < 2:4648037/73‑1 (MQ=255) aTGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccag > 1:1750677/1‑71 (MQ=255) aTGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccag < 2:1750677/71‑1 (MQ=255) tGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagcc > 2:474259/1‑76 (MQ=255) gAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagcca > 1:6795143/1‑76 (MQ=255) aaaaCTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccag < 2:5285303/76‑1 (MQ=255) aaaaCTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccag > 2:7459005/1‑76 (MQ=255) aaaaCTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccag < 1:5794161/76‑1 (MQ=255) aaaaCTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccag < 1:6166908/76‑1 (MQ=255) aaCTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccagcc > 1:4614620/1‑76 (MQ=255) aaCTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccagcc > 2:6215596/1‑76 (MQ=255) aaCTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccagcc > 2:1423600/1‑76 (MQ=255) aaCTGCTTCGTTATTGACTTCGAATATAAAGGGGCGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccagcc > 2:7240896/1‑76 (MQ=255) aaCTGCTTCGTTATTGACTTCAAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccagcc < 1:7580280/76‑1 (MQ=255) aaCTGCTCCTTTATGGACTTCGAATATAAAGGGGTGAATTTTGCT‑TTAGTTTTCA‑‑‑‑cgg‑‑‑‑‑ccagccagccagccagcc < 2:6151823/76‑1 (MQ=255) aCTGCTTGTTTATTGATTTCGAATATAAGGGGGTAAAATTGGCA‑TTAGTTTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccagcca < 2:4029336/76‑1 (MQ=255) aCTGCTTCTTTTTTGATTTGGAATAAAAAGGGGGAAATTTTGCA‑TTAGTTTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccagcca < 2:1527879/76‑1 (MQ=255) cTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccagccag > 1:5107521/1‑76 (MQ=255) tGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccagccagc > 1:2958515/1‑76 (MQ=255) ttttCGTTTTTGACTTCGAATATAAAGGGGGAAATTTTGAT‑TTAGTTTGCC‑‑‑‑cag‑‑‑‑‑ccagccagccagccagccagcc < 2:1020676/74‑1 (MQ=255) gCTTCGTTATTGATTTCGAATAAAAAGGGGTAAAATTTGCA‑TTAAATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccagccagcc < 1:1270144/76‑1 (MQ=255) gCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccagccagcc < 1:1624366/76‑1 (MQ=255) gCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccagccagcc > 1:1381109/1‑76 (MQ=255) gCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccagccagcc > 2:7857884/1‑76 (MQ=255) gCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccagccagcc > 2:1607870/1‑76 (MQ=255) gCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccagccagcc > 2:4227789/1‑76 (MQ=255) cTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCT‑TTAGATTTCC‑‑‑‑cag‑‑‑‑‑ccagccagccagccagccagcca < 2:3932536/76‑1 (MQ=255) ttCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccagccagccag > 2:763463/1‑76 (MQ=255) ttCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑cag‑‑‑‑‑ccagccagccagccagccagccag > 1:1308988/1‑76 (MQ=255) tttATTGACTTCGAATATAAGGGGGTAAATTTTGCA‑TTAAATTCCA‑‑‑‑CAg‑‑‑‑‑ccagccagccagccagccagccagcca < 1:1255352/75‑1 (MQ=255) gctATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑CAg‑‑‑‑‑ccagccagccagccagccagccagcca > 1:3338605/3‑76 (MQ=255) gTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑CAg‑‑‑‑‑ccagccagccagccagccagccagcca > 2:7414632/1‑76 (MQ=255) tttttGACTTCAAATAAAAAGGGG‑GAAATTTTCATTTAGTTTCCc‑‑‑‑cag‑‑‑‑‑ccagccagccagccagccagccagccag < 1:6538088/73‑1 (MQ=16) ttATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑CAG‑‑‑‑‑ccagccagccagccagccagccagccag > 1:684043/1‑76 (MQ=255) ttATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑CAG‑‑‑‑‑ccagccagccagccagccagccagccag > 1:1768988/1‑76 (MQ=255) tATTGACTTTGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑CAG‑‑‑‑‑ccagccagccagccagccagccagccagc > 1:3318601/1‑75 (MQ=255) tATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑CAG‑‑‑‑‑ccagccagccagccagccagccagccagc > 1:298705/1‑75 (MQ=255) aaGGGGTGAAATTTGCA‑TTAGATTTCACAGCCAG‑‑‑‑‑CCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGAt > 2:1521059/1‑76 (MQ=255) aaGGGGTGAAATTTGCA‑TTAGATTTCACAGCCAG‑‑‑‑‑CCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGAt > 1:2510230/1‑76 (MQ=255) tGAAATTTGCA‑TTAGATTTCACAGCCAG‑‑‑‑‑CCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGTATGTGGATTTTGta > 2:2841146/1‑76 (MQ=25) taaaTTTGCA‑TTAGTTTTCA‑‑‑‑CAGCCCGCCCACCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGta < 2:3608518/75‑1 (MQ=19) gAAATTTACA‑TTAGATTTCATAGCCAG‑‑‑‑‑CCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGtac > 2:4558389/1‑76 (MQ=25) aaaTTTGCA‑TTAGATTTCACAGCCAG‑‑‑‑‑CCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGtact > 2:6240042/1‑76 (MQ=37) aTTTGCA‑TTAGATTTCACAGCCAG‑‑‑‑‑CCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGtactac > 2:2898388/1‑76 (MQ=37) tttGCA‑TTAGATTTCACAGCCAG‑‑‑‑‑CCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACg > 2:4332680/1‑76 (MQ=37) tttGCA‑TTAGATTTCA‑‑‑‑CAG‑‑‑‑‑CCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAaat > 2:7611092/1‑76 (MQ=255) | AAAACATGGCTATACAATGAAAACTGCTTCGTTATTGACTTCGAATATAAAGGGGTGAAATTTGCA‑TTAGATTTCA‑‑‑‑CAG‑‑‑‑‑CCAGCCAGCCAGCCAGCCAGCCAGCCAGGTGTGAATGTGGATTTTGTACTACGAAAT > REL606/2103812‑2103947 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |