| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 4,574,252 | 0 | T | G | 36.7% | 48.4 / 11.9 | 30 | S5S (TCA→TCC) | hpaX | hypothetical 4‑hydroxyphenylacetate permease |
| Reads supporting (aligned to +/- strand): ref base T (2/17); new base G (11/0); total (13/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.43e-06 | |||||||||||
| 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 | ctgaataaCCGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCaggt < 1:49418/73‑1 (MQ=255) ccGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTc < 1:1448069/76‑1 (MQ=255) cGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCg < 1:1120335/76‑1 (MQ=255) ccTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCt < 1:654537/76‑1 (MQ=255) gtCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGAGGTGTCGCTCAtttt > 1:28420/1‑76 (MQ=255) ggCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTttctt < 2:182260/76‑1 (MQ=255) cAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGAGGGGTTGCTCCTTTTTTTTcc > 2:783918/1‑76 (MQ=255) aGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTaa < 1:634539/76‑1 (MQ=255) aGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGCGGTGGCGCTTATTTTTTTTCctctca > 2:1560800/1‑74 (MQ=255) cgcTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATg < 1:438635/76‑1 (MQ=255) gcTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGTGTCGCTCATTTTTTTTCCTCTAATgg > 1:1861154/1‑76 (MQ=255) cTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATggg < 2:510083/76‑1 (MQ=255) tATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGGGGGGCTCATTTTTTTTCTTCTAATTGGAtt > 2:371559/1‑76 (MQ=255) aTGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGAttt < 2:1625219/76‑1 (MQ=255) gCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGGGTGGCTCATTTTTTTTTCTCTTATGGGGttttt > 1:1220501/1‑76 (MQ=255) gCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGAGGGGTCGCTCATTTTTTTTCCTCTTATGGGGttttt > 1:135850/1‑76 (MQ=255) cTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGGGTCGCTCATTTTTTTTTCTCTAATGGGGtttttt > 1:1403190/1‑75 (MQ=255) gATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGTGTCGCCCCTTTTTCTTTTTTTTATGGGAtttttttt > 1:732681/1‑73 (MQ=255) aTTCGCCGGATCGATACTCTCCGGTATGGCAGGGGAGGGGTCGCTCCTTTTTCTTTCTCTTATTGGATTTTTTTAg > 1:276435/1‑76 (MQ=255) tCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGGGTCGCTCATTTTTTTTTCTCTTTTGGGGTTTTTTTTggg > 1:1212274/1‑76 (MQ=255) tCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGGGTCGCTCATTTTTCTTCCTCTTATGGGGTTTTTGTTggg > 2:346671/1‑76 (MQ=255) gCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTa < 2:1584786/76‑1 (MQ=255) ccGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTAc < 1:1813405/76‑1 (MQ=255) ggATCGATACTCTCCGGTATGGCAGGGGGGGTGTCGCTCCTTTTTTTTTCTCTTATGGGGTTTTTTTAGGGGACGa > 1:245996/1‑76 (MQ=255) gATCGATACTCTCCGGTATGGCAGGGGAGGTGTCGCTCCTTTTTTTTCCTCTAATGGGATTTTTTTTGGGTACGAc > 2:1736331/1‑76 (MQ=255) aTCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACa < 1:624805/76‑1 (MQ=255) cGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAg < 2:408642/76‑1 (MQ=255) gATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAgg < 2:1841492/76‑1 (MQ=255) aTACTCTCCGGTATGGCAGGGGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGGTTTTTGTAGGGGACGACAAGGt > 2:1011326/1‑76 (MQ=255) tACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGtt < 1:793652/76‑1 (MQ=255) ccGGTATGGCAGGGGAGGTGTCGCTCATTTTTCTTTCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGc > 2:465279/1‑76 (MQ=255) gTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGa > 2:1392521/1‑76 (MQ=255) gTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGa > 1:206363/1‑76 (MQ=255) tATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGAt < 2:488938/76‑1 (MQ=255) aTGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATGTTTGTAGGGTACGACAAGGTTTATGGCGGATg < 2:1786266/76‑1 (MQ=255) gCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATgcgg < 1:303150/76‑1 (MQ=255) aggtgaggtgTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATgcggcg > 2:968084/1‑76 (MQ=255) gtgaggtgTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATGCGGCGTa < 1:213734/76‑1 (MQ=255) gtgaggtgTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATGCGGCGTa < 1:219065/76‑1 (MQ=255) gtgaggtgTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATGCGGCGTa < 1:1793819/76‑1 (MQ=255) tgaggtgTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATGCGGCGTaa > 2:1423726/1‑76 (MQ=255) tgaggtgTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATGCGGCGTaa > 1:1727219/1‑76 (MQ=255) | CTTAATAACCGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATGCGGCGTAA > REL606/4574177‑4574327 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |