| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 4,574,252 | 0 | T | G | 45.7% | 38.1 / 17.1 | 35 | S5S (TCA→TCC) | hpaX | hypothetical 4‑hydroxyphenylacetate permease |
| Reads supporting (aligned to +/- strand): ref base T (2/17); new base G (16/0); total (18/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.19e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTAATAACCGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATGCGGCGTA > REL606/4574178‑4574326 | ttAATAGCCCCCTGTTGCCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCaggtg < 2:85452/76‑1 (MQ=255)ttAATAACCGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGACGGATCGATACTCTCCGGTATGGCaggtg < 1:745995/76‑1 (MQ=255) taataaCCGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCaggtga < 2:1637565/76‑1 (MQ=255) aataaCCGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCaggggag > 1:110001/1‑76 (MQ=255) aCCGCCTGGTGGCCGGCAGTCAGCGCTTTATGCTGATTCGCGGGATCGATACTCTCCGGTATGGCAGGTGAGgtgt < 2:1189750/76‑1 (MQ=255) gCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGTGgcgc > 2:1298759/1‑76 (MQ=255) ccTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCTGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCt < 1:1484203/76‑1 (MQ=255) cTGTTGTCCGGCAGTCAGCGCTTTGTGCNTATTCGCCGGCTCGATACTCTNCGATATGGCAGGGGGGGTGTCGCTc > 1:258947/1‑76 (MQ=255) cTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTc < 1:1339955/76‑1 (MQ=255) tgttgtCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCa < 2:1655424/76‑1 (MQ=255) tgttggCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCa < 2:1397542/76‑1 (MQ=255) tgtCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCAttt < 1:1396739/76‑1 (MQ=255) tgtCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCAttt < 1:959337/76‑1 (MQ=255) gtCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGGGTCGCTTAtttt > 1:191359/1‑76 (MQ=255) cGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGAGGGGGCGCTCAttttttt > 1:327646/1‑74 (MQ=255) aGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGGGTCGCTCCTTTTTTTTTct > 1:205572/1‑76 (MQ=255) gTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGGGGGGGTCATTTTTTTTTctt > 2:433053/1‑75 (MQ=255) cAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGAGGCGTGGCTCCTTTTTCTTCctctc > 1:285758/1‑75 (MQ=255) aGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGTGTCGCTCCTTTTTTTTTctctta > 2:1169091/1‑74 (MQ=255) cTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATggg < 1:537965/76‑1 (MQ=255) ttATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGGGGCGCTCCTTTTTCTTCCTCTCATGGGag > 2:887800/1‑75 (MQ=255) tATGCTGATTCGCCTGATCGATACTCTCCGGTATGGCAGGGGGGGGGTCGCTCCTTTTTTTTCCTCTTATGGGGtt > 1:920583/1‑76 (MQ=255) tATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGGGTCGCTCATTTTTTTTCCTCTCATGGGGtt > 2:1504726/1‑76 (MQ=255) tGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGAtttt < 1:192254/76‑1 (MQ=255) cTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGAGGGGTCGCTCATTTTTTTTCCTCTCATGGGAtttttt > 2:1074752/1‑75 (MQ=255) gATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGTGTCGCTCATTTTTTTTCCTCTAATGGGGTTTTTGTa > 2:1475089/1‑76 (MQ=255) gATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGAGGGGTCGCTCCTTTTTTTTTCTCTTATGGGGTTTTTGTa > 1:1337907/1‑76 (MQ=255) tCGCCGGATCTATACTCTCCGGTATGGCAGGGGGGGGGGCGCTCCTTTTTTTTTCTCTTTTTGGGTTTTTTTTggg > 1:1104744/1‑76 (MQ=255) tCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAggg < 1:664150/76‑1 (MQ=255) gCCGGATCGATACTNTCCGGTATGGCAGGTGAGGTGTCGCTCATTTATCTTCCTCTAATGGGATTTGTGGNGTGTa < 2:1766505/76‑1 (MQ=255) ccGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTAc < 1:911940/76‑1 (MQ=255) cGGATCGATACTCTCCGGTATGGCAGGGGAGGGGTCCCTCATTTTTTTTTCTCTTGTGGGGTTTTTTTAGGGGACg > 2:1410073/1‑76 (MQ=255) cGGATCGATACTCTCCGGTATGGCAGGGGAGGGGGCGCTCATTTTTTTTTCTCTTATGGGGTTTTTTTAGGGGACg > 1:957598/1‑76 (MQ=255) cGGATCGATACTCTCCAGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGTGTACg < 1:112180/76‑1 (MQ=255) ggATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGa < 2:1075041/76‑1 (MQ=255) gATCGATACTCTCCGGTATGGCAGGGGGGGGGGCGCTCATTTTTTTTCCTCTTATGGGGTTTTTTTAGGGgacgac > 1:1334039/1‑76 (MQ=255) tCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACaa < 1:1053486/76‑1 (MQ=255) cGATACTCTCCGGTATGGCAGGGGAGGGGTCGCTCCTTTTTTTTCCTCTAATGGGGTTTTTTTAGGGGACGACAAg > 2:933871/1‑76 (MQ=255) tACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGtt < 1:461533/76‑1 (MQ=255) tctcCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTAt < 2:614063/76‑1 (MQ=255) cctcCGGTATGGCAGGTGAGGTGTCGCTCATTGTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTAt < 2:861125/75‑1 (MQ=255) ccGGTATGGCAGGGGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGGACGACAAGGGTTATGGc > 2:1385250/1‑76 (MQ=255) cAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATgcggc > 1:1227517/1‑76 (MQ=255) cAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATgcggc > 2:627511/1‑76 (MQ=255) aggtgaggtgTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATgcggcg < 2:1671740/76‑1 (MQ=255) aggtgaggtgTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATgcggcg > 2:1011571/1‑76 (MQ=255) aggtgaggtgTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATgcggcg < 1:316243/76‑1 (MQ=255) aggtgaggtgACGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTGTGGCGGATgcggcg > 1:142603/1‑76 (MQ=255) gtgaggtgTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATGCGGCGTa < 1:1241665/76‑1 (MQ=255) gtgaggcgACGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATGCGGCGTa > 1:1636910/1‑76 (MQ=255) | TTAATAACCGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGCTCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATGCGGCGTA > REL606/4574178‑4574326 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |