Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,983,700 | 0 | T | G | 25.0% | 90.8 / 12.6 | 40 | T397P (ACC→CCC) | ECB_03713 | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (4/26); new base G (10/0); total (14/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.18e-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 | gAAGGGGAAATGCGGGGTTGCCGCCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGt < 2:1688377/76‑1 (MQ=255) aaGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCAACGGCAGGGGTg > 2:153535/1‑76 (MQ=255) gtAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGtctc > 1:1989281/1‑76 (MQ=255) tAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGgccgcc > 2:853063/1‑76 (MQ=255) aTGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCcac < 1:14953/76‑1 (MQ=255) gCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACaa < 2:314162/76‑1 (MQ=255) gCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGCGGCCCCCCCaa > 1:1782448/1‑76 (MQ=255) cGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCCCCCaaa > 1:146988/1‑76 (MQ=255) cGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCCCCCCCaaa > 1:490950/1‑76 (MQ=255) ggggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTTCCACCAAt > 2:127375/1‑76 (MQ=255) gggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAAtt < 2:120936/76‑1 (MQ=255) ggTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAAttt < 2:1894633/76‑1 (MQ=255) gTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTg < 1:228471/76‑1 (MQ=255) gTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTg < 1:672341/76‑1 (MQ=255) gtACCACCGAGAGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGa < 2:1143638/75‑1 (MQ=255) tACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGAc < 1:1584528/76‑1 (MQ=255) ccaccGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCACAAATTTTACGa > 2:938138/1‑76 (MQ=255) accGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATc < 1:574020/76‑1 (MQ=255) cGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCCCCCACAAAATTTACGATTGc > 2:1637163/1‑76 (MQ=255) gggTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTgg < 2:760586/76‑1 (MQ=255) ggTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGga < 1:1711134/76‑1 (MQ=255) gtgcttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGa < 1:330844/75‑1 (MQ=255) cttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGGGCCCCCCCCAAAATTTACGATCGCTGGGGAAGaacg > 1:766902/1‑73 (MQ=255) ttgcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGcc < 1:143411/76‑1 (MQ=255) gcATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCAt < 1:151434/76‑1 (MQ=255) tCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATcgc < 1:1993820/76‑1 (MQ=255) ccTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATcgcg < 1:37859/76‑1 (MQ=255) tACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGaa < 1:1405345/76‑1 (MQ=255) aCGGCTGTTATTGGTGGCATCGGCAGGGGGGGCTCTCCCAAAATTGACGACCGCTGGGGATGACGCGCTAGCGAAg > 2:381596/1‑76 (MQ=255) aCGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCACAAAATTTACGATCCCTGGAGATGAAGCCACCGCGCAg > 1:434006/1‑76 (MQ=255) cGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGc < 2:1883530/76‑1 (MQ=255) cGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGc < 2:468501/76‑1 (MQ=255) cGGCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCCCAAATTTGACGATCGCCGGAGATGACGCCATCCCGAAGc > 2:2012509/1‑76 (MQ=255) gCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCACAAATTTTACGATCCCTGGAGATGACGCCATCCCGAAGcca > 2:1347463/1‑76 (MQ=255) gCTGTTATTGGTGGCATCGGCAGGGGGGCCTCCCACAAATTTGACGATCGCTGGGGATGACGCCATCGCGAAGccc > 1:1218804/1‑75 (MQ=255) gCTATTATTGGTGGCATTGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGcca < 2:101542/76‑1 (MQ=255) cTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccac < 2:652222/76‑1 (MQ=255) ttATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGccaccca < 1:1784852/76‑1 (MQ=255) tATTGGTGGCATCGGCAGGGGTGCCTCCCACAAAATTGACGATCGCTGGGGATGACGCCATCCCGAAGccacccac > 1:1193192/1‑76 (MQ=255) tATTGGTGGCATCGGCAGGGGGGCCTCCCCCCAATTTGACGATCGCTGGAGAAGACGCCATCGCGAAGccacccac > 2:1481293/1‑76 (MQ=255) aTTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCcaca < 1:384222/76‑1 (MQ=255) tggtggCATCGGCAGGGGGGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCAcccc > 2:811037/1‑76 (MQ=255) gtggCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTg < 2:284696/76‑1 (MQ=255) ggCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCa > 2:2103094/1‑76 (MQ=255) ggCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCa > 2:2055525/1‑76 (MQ=255) ggCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCa < 2:1791567/76‑1 (MQ=255) cAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGt < 1:1626933/76‑1 (MQ=255) aGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTa > 2:551773/1‑76 (MQ=255) tGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTgc > 1:95542/1‑76 (MQ=255) tGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTgc > 2:1470107/1‑76 (MQ=255) | GAAGTGTAAATGCGGGGTTACCACCGAGGGTGATTGCTTGCATCCTACGGCTGTTATTGGTGGCATCGGCAGGGGTGCCTCCCACAAATTTGACGATCGCTGGAGATGACGCCATCGCGAAGCCACCCACACCTGCGGTTTCAGTAATTGC > REL606/3983625‑3983775 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |