| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 4,574,252 | 0 | T | G | 30.6% | 65.5 / 19.7 | 36 | S5S (TCA→TCC) | hpaX | hypothetical 4‑hydroxyphenylacetate permease |
| Reads supporting (aligned to +/- strand): ref base T (2/23); new base G (11/0); total (13/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.30e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTTAATAACCGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATGCGGCGTAA > REL606/4574177‑4574327 | cttAATAACCGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCaggt < 2:928613/76‑1 (MQ=255) ataaCCGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCaggggagg > 1:1701661/1‑76 (MQ=255) taaCCGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCaggggtggt > 1:559101/1‑76 (MQ=255) cGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCg < 2:1745761/76‑1 (MQ=255) cGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCg < 1:1764274/76‑1 (MQ=255) tgttgtCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGGGGCGCTcc > 2:1674943/1‑75 (MQ=255) ttgtCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTGGCTCGtt > 1:72196/1‑76 (MQ=255) ttgtCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCAtt < 1:361537/76‑1 (MQ=255) tCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGGGTGGCTCAttttt > 2:1572093/1‑76 (MQ=255) ggCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTttctt < 2:1660944/76‑1 (MQ=255) tCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCctct < 1:685745/76‑1 (MQ=255) cAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTa < 1:1031535/76‑1 (MQ=255) aGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTaa < 1:488024/76‑1 (MQ=255) cgcTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATg < 1:798930/76‑1 (MQ=255) tttATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGa < 2:1094902/76‑1 (MQ=255) tttATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGa < 1:1514451/76‑1 (MQ=255) tttATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGAGGGGTCGCTCATTTTTTTTCCTCTTATgggg > 1:763408/1‑75 (MQ=255) tGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGGGTCGCTCATTTTTTTTTCTCTTATGGGGtttt > 1:83207/1‑76 (MQ=255) cTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGAGGGGTCGCTCCTTTTTTTTTCTCTTATGGGGtttttt > 2:1652544/1‑75 (MQ=255) gATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGCGGTGGTCCGTTTTTTTTCTCTTATGGGGTTTTTGtt > 2:31238/1‑75 (MQ=255) tCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAggg < 2:21459/76‑1 (MQ=255) cGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCNTCTAATGGGATTTTTGTAGGGt < 2:261266/76‑1 (MQ=255) cGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGGGTCGCTCCTTTTTTTTTCTCTTATGGGGTTTTTTTTgggg > 1:1700772/1‑75 (MQ=255) gCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTa < 2:1347223/76‑1 (MQ=255) gCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTa < 2:1404417/76‑1 (MQ=255) ccGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTAc < 2:1296144/76‑1 (MQ=255) ccGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTAc < 2:1462465/76‑1 (MQ=255) cGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACg < 1:693051/76‑1 (MQ=255) cGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACg < 2:10901/76‑1 (MQ=255) aTCGATACTCTCCGGTATGGCAGGGGGGGGGGCGCTCCTTTTTTTTTCTCTTATGGGGTTTTTTTTGGGGACGGCa > 1:1133764/1‑76 (MQ=255) aTCGATACTCTCCGGTATGGCAGGGGAGGGGTCGCTCATTTTTTTTCCTCTAATGGGGTTTTTTTAGGGGACGaaa > 2:1489359/1‑74 (MQ=255) cGATACTCTCCGGTATGGCAGGGGAGGGGTGGCTCATTTTTCTTTCTCTTATTGGATTTTTTTTGGGTACGACAAg > 1:1300386/1‑76 (MQ=255) aCTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGttt < 2:5930/76‑1 (MQ=255) tctcCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTAt < 2:83207/76‑1 (MQ=255) tcCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATgg < 1:1535506/76‑1 (MQ=255) tATGGCAGGGGAGGGGTCGCTCATTTTTTTTTCTCTAATGGGGTTTTTGTAGGGTACGACAAGGGTTATTGCgggt > 2:1221341/1‑74 (MQ=255) aTGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATg > 2:1650413/1‑76 (MQ=255) ggCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATgcg < 1:1174882/76‑1 (MQ=255) ggCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATgcg < 2:760343/76‑1 (MQ=255) cAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATgcggc > 2:1410165/1‑76 (MQ=255) aggtgaggtgTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATgcggcg > 1:1304187/1‑76 (MQ=255) aggtgaggtgTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATgcggcg > 1:1158344/1‑76 (MQ=255) tgaggtgTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATGACGCGTaa > 2:199070/1‑76 (MQ=255) | CTTAATAACCGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATGCGGCGTAA > REL606/4574177‑4574327 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |