Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,983,700 | 0 | T | G | 32.8% | 108.0 / 18.5 | 58 | T397P (ACC→CCC) | ECB_03713 | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (7/32); new base G (19/0); total (26/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.94e-10 | |||||||||||
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 | gAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGt < 1:1615339/76‑1 (MQ=255) gAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGt < 2:767222/76‑1 (MQ=255) aaGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTg < 2:783085/76‑1 (MQ=255) aaGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTg < 2:715270/76‑1 (MQ=255) aGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGc < 2:771882/76‑1 (MQ=255) aaaTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTccc > 1:1815356/1‑76 (MQ=255) tGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCcaca < 2:843498/76‑1 (MQ=255) gCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACaa < 2:945419/76‑1 (MQ=255) cGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCCCCCACaaa > 2:1656899/1‑76 (MQ=255) cGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCaaa > 1:1180578/1‑76 (MQ=255) ggggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCCaaaa > 2:1614422/1‑75 (MQ=255) ggggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCCCCCCCaaaa > 2:1477448/1‑75 (MQ=255) ggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCACAAAatt > 2:247852/1‑76 (MQ=255) gTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTg < 1:1180643/76‑1 (MQ=255) gTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTg < 1:1461751/76‑1 (MQ=255) tACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGAc < 1:138175/76‑1 (MQ=255) tACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCCAAATTTGAc > 2:526973/1‑76 (MQ=255) caccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGAt < 2:600793/76‑1 (MQ=255) caccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGAt < 2:575794/76‑1 (MQ=255) cGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCCAAAATTTACGATCGc > 2:1872205/1‑76 (MQ=255) gAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCAAATTTGACGATCGCt > 2:443544/1‑76 (MQ=255) taaTTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGgaga < 2:1306963/74‑1 (MQ=255) tGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCCCCCCCAAAATTGACGATAGCTGgagg > 1:2042289/1‑75 (MQ=255) gATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGAt < 2:214021/76‑1 (MQ=255) gATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGAt < 2:1234970/76‑1 (MQ=255) aTTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCAAAATTTACGATCGCTGgggagg > 2:757752/1‑74 (MQ=255) ttgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCAAATTTGACGATCGCTGGAGATGa > 1:1025579/1‑76 (MQ=255) ttgcttgcATCCTACGGCTGGTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGa < 2:2039675/76‑1 (MQ=255) tgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGAc < 2:304041/76‑1 (MQ=255) tgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCCCCCACAAAATTTACGATCGCTGGGGATGAc > 2:638540/1‑76 (MQ=255) gcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCCCCCCCAAAATTTACGAACGCTGGGGGTGACg > 1:1567549/1‑76 (MQ=255) gcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCCCCCCCAAAATTTACGAACGCTGGAGATGGCg > 1:983091/1‑76 (MQ=255) cttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGGTCCCACAAATTTGACGATGGCTGGCGGTGACGc > 1:979463/1‑76 (MQ=255) cttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCCCCCCCAAAATTGACGACCGCTGGAGATGgcgc > 2:1422399/1‑76 (MQ=255) ctCCTGCGGCTGGTATTGGTGGCATCGGCAGTGGTGCCTCCCACAAATTTGGCGATCGCTGGAGATGACGCCATcg < 1:316574/75‑1 (MQ=255) aTCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCCCCCCCCAAATTGACGATCGCCGGGGATGACGCCATcc > 1:1387944/1‑75 (MQ=255) ccTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATcgcg < 1:1957634/76‑1 (MQ=255) tACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCCCCCACAAAATTTACAATCGCTGGAGAGGACGCCCTCGCGaa > 2:1874589/1‑76 (MQ=255) cGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGTCGAGCGCTGGAGATGTCGCCATCGCGAAGc < 2:1953189/76‑1 (MQ=255) cGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGc < 2:808343/76‑1 (MQ=255) ggCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGcc < 1:78982/76‑1 (MQ=255) gCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGcca < 2:1105446/76‑1 (MQ=255) cTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccac > 2:1198581/1‑76 (MQ=255) gTTATTGGTGGCATCGGCAGGGGGGCCTCCCACAAATTTGACGGTCGCTGGAGATGACGCCCTCGCGGAGccaccc > 1:1154564/1‑76 (MQ=255) ttATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccaccca < 1:1782223/76‑1 (MQ=255) ttATTGGTGGCATCGGCAGGGGGGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccacccc > 2:1420512/1‑75 (MQ=255) ttATTGGTGGCATCGGCAGGGGGGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccacccc > 2:336648/1‑75 (MQ=255) aTTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCcaca < 2:1815356/76‑1 (MQ=255) ttGGTGGCCTCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCacac < 2:1009186/76‑1 (MQ=255) ttGGTGGCATCGGCAGGGGTGCCTCCCTCAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCacac < 1:1138549/76‑1 (MQ=255) tggtggCATCGGCAGGGGGGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACAcc > 2:1359656/1‑76 (MQ=255) ggtggCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCt > 1:762238/1‑76 (MQ=255) ggtggCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCt < 1:1024495/76‑1 (MQ=255) cATCGGCAGGGGGGCCTCCCACAAAATTGACGATCGCTTGGGATGACGCCATCGCGGAGCCACCCCCACCTCCggg > 2:1011679/1‑75 (MQ=255) cGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTc < 1:182915/76‑1 (MQ=255) cGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTc < 1:1976154/76‑1 (MQ=255) cGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTc < 1:85041/76‑1 (MQ=255) cGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTc < 2:85088/76‑1 (MQ=255) ggCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCa < 1:904480/76‑1 (MQ=255) gCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAg > 1:1387903/1‑76 (MQ=255) cAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGt > 1:1756219/1‑76 (MQ=255) aGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTa > 1:990826/1‑76 (MQ=255) aGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTa > 2:1635373/1‑76 (MQ=255) aGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTa > 2:894663/1‑76 (MQ=255) ggTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAAtt < 2:1161248/76‑1 (MQ=255) gTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTg > 1:1153618/1‑76 (MQ=255) tGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTgc < 2:738243/76‑1 (MQ=255) | GAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTGC > REL606/3983625‑3983775 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |