Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,691,844 | 0 | T | G | 36.8% | 51.2 / 17.8 | 38 | N226T (AAC→ACC) | yiaY | predicted Fe‑containing alcohol dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (7/17); new base G (14/0); total (21/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.73e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATACGCT > REL606/3691772‑3691918 | gCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGgttt > 1:68139/1‑76 (MQ=255) gCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGgttt > 1:225387/1‑76 (MQ=255) cGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTc < 1:513343/76‑1 (MQ=255) cGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTc > 1:1496486/1‑76 (MQ=255) aGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTTgg > 2:1087371/1‑76 (MQ=255) aCTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAAt > 1:615796/1‑76 (MQ=255) aCTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAAt > 1:1587823/1‑76 (MQ=255) cTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATc > 2:435879/1‑76 (MQ=255) tGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCa > 1:543667/1‑76 (MQ=255) gggCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCAt < 1:1353659/76‑1 (MQ=255) aTAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTc > 1:713292/1‑76 (MQ=255) tAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCGATCATGGTCa > 1:1722348/1‑76 (MQ=255) aGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACg < 2:729298/76‑1 (MQ=255) gCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCCTGGTCCCGGCTTTct > 2:575984/1‑75 (MQ=255) cACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGc < 2:1546691/76‑1 (MQ=255) cACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTTAGTGc > 1:328442/1‑76 (MQ=255) aCGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGca < 1:236017/76‑1 (MQ=255) aCGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGca < 1:1627137/76‑1 (MQ=255) cgcgCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGCTTTTGGCAAACATGGGCACGGCTTTCAGTGccc > 1:724416/1‑74 (MQ=255) gcgcTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGcaca > 1:618146/1‑76 (MQ=255) gcTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGGCACGGGCTTCCGTGCCCaaa > 1:1597449/1‑75 (MQ=255) cTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGc > 2:1193042/1‑76 (MQ=255) tttCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCg > 1:274895/1‑76 (MQ=255) cGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCa < 2:981933/76‑1 (MQ=255) cATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCCTTCAGTGCACAAGCGTCAGt > 2:1489790/1‑76 (MQ=255) aTTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTg < 2:589832/76‑1 (MQ=255) ttACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATTGTCACGGCTTTTAGTGCACAAGCGGCAGTgg > 1:1294359/1‑75 (MQ=255) gCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGc > 1:172769/1‑76 (MQ=255) aTCTTCAACGGCTAACGGCAGGGTTTCGGCCATCATGGTCACGGCTTTCCGTGCCCAAGCGGCAGTGATCGGCGTg > 1:855710/1‑76 (MQ=255) cTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTggc < 1:1700114/76‑1 (MQ=255) ttCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGGCCGGGATCGGGCTggcg > 2:1737851/1‑76 (MQ=255) cAACGGCTAACGGCAGGGTTTTGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTggcggc > 2:1290445/1‑76 (MQ=255) cAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTggcggc > 2:698815/1‑76 (MQ=255) cGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAAt < 1:1428309/76‑1 (MQ=255) gCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAg < 1:452412/76‑1 (MQ=255) tAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGaa < 1:1426413/76‑1 (MQ=255) aaCGGCAGGTTTTCGGCCATCATGGTCACCGCTTTTTGTGCCCAAACGCCCGGGAACGGGGTGGGGGGAATAaaaa > 2:1211920/1‑76 (MQ=255) aaCGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGaaa < 2:178607/76‑1 (MQ=255) aCGGCAGGTTGTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAAc < 2:948368/76‑1 (MQ=255) cGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACa < 2:195267/76‑1 (MQ=255) gCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATa > 2:1193990/1‑76 (MQ=255) gCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATa < 1:976581/76‑1 (MQ=255) cAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATAc > 2:519472/1‑76 (MQ=255) aGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATagg < 1:717522/76‑3 (MQ=255) gTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATACGCt > 2:1658114/1‑76 (MQ=255) | GCGAGGAACTGGGCATAAGCCATTGCTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATACGCT > REL606/3691772‑3691918 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |