| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 4,574,252 | 0 | T | G | 38.5% | 28.1 / 15.7 | 26 | S5S (TCA→TCC) | hpaX | hypothetical 4‑hydroxyphenylacetate permease |
| Reads supporting (aligned to +/- strand): ref base T (3/13); new base G (10/0); total (13/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.08e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTAATAACCGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGC‑TCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATGCGGCG > REL606/4574178‑4574324 | ttAATAACCGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCaggtg < 1:293974/76‑1 (MQ=255) taataaCCGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCaggtga < 2:117692/76‑1 (MQ=255) ataaCCGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCaggtgtgg > 1:1669154/1‑76 (MQ=255) gttgtCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGAGGGGTCGC‑TCAt > 2:1538685/1‑76 (MQ=255) gtCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATTCTCTCCGGTATGGCAGGGGCGGTGTCGCTTCAttt > 2:654209/1‑76 (MQ=255) cGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGGGGGGC‑TCATTTttct > 1:1063731/1‑76 (MQ=255) cGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGAGGGGCGGC‑TGCTTTttct > 1:637428/1‑76 (MQ=255) aGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGGGTCGC‑TCCTTTTTTTTCct > 2:921178/1‑76 (MQ=255) gTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGCGGGGGGGGCC‑CCCTTTTTCTTCctc > 2:636449/1‑76 (MQ=255) gcTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGC‑TCATTTTTCTTCCTCTAATgg < 1:1049244/76‑1 (MQ=255) tttATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGGGGCGG‑TCCTTTTTTTTCCTCTAATgggg > 2:1838907/1‑75 (MQ=255) tATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGC‑TCATTTTTCTTCCTCTAATGGGAtt < 2:686019/76‑1 (MQ=255) tGATTCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGGGGGGC‑TCATTTTTTTTTCTCTTATGGGGttttttt > 1:1124563/1‑74 (MQ=255) tCGCCGGATCGATACTCTCCGGTATGGCAGGGGGGGGGTCGC‑TCCTTTTTTTTCCTCTAATGGGGTTTTTTTAggg > 2:486822/1‑76 (MQ=255) cGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGC‑TCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGt < 2:816893/76‑1 (MQ=255) cGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGC‑TCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGt < 2:847896/76‑1 (MQ=255) gCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGC‑TCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTa < 1:1318318/76‑1 (MQ=255) gCCGGATCGATACTCTCCGGTATGGCAGGGGGGGGGTCGC‑TCCTTTTTTTTCCTCTAATGGGGTTTTTTTAgggga > 1:341852/1‑74 (MQ=255) gCCGGATCGATACTCTCCGGTATGGCAGGGGAGGGGGCCC‑TTCATTTTCTTCCTCTAATTGGATTTTTTTTGGGtg > 1:275500/1‑75 (MQ=255) ccGGATCGATACTCTCCGGTATGGCAGGGGAGGGGTCGC‑TCCTTTTTCTTTCTCTACTGGGGTTTTTGTTGGGGAc > 2:831875/1‑76 (MQ=255) cGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGC‑TCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACg < 2:959983/76‑1 (MQ=255) cGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGC‑TCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACg < 1:448517/76‑1 (MQ=255) cGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGC‑TCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACg < 1:362069/76‑1 (MQ=255) cGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGC‑TCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACg < 1:362051/76‑1 (MQ=255) gATCGATACTCTCCGGTATGGCAGGGGGGGGGTCGC‑TCATTTTTTTTTCTCTTATTGGGTTTTTTTAGGGgacgac > 1:335534/1‑76 (MQ=255) aTCGATACTCTCCGGTATGGCAGGTGATGTGTCGC‑TCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACa < 1:1292128/76‑1 (MQ=255) cGATACTCTCCGGTATGGCAGGGGGGGGGTCGC‑TCATTTTTTTTCCTCTTATGGGGTTTTTGTAGGGGACGACAAg > 2:535690/1‑76 (MQ=255) gATACTCTCCGGTATGGCAGGGGAGGGGGCGC‑TCCTTTTTTTTCCTCTTATGGGATTTTTGTAGGGTACGACAAgg > 2:74965/1‑76 (MQ=255) aTACTCTCCGGTATGGCAGGTGGGGGGGCGC‑TCCTTTTTCTTTCTCTTATGGGATTTTTTTTGGGTACGACAAGGt > 2:1802469/1‑76 (MQ=255) tACTCTCCGGTATGGCAGGTGAGGTGTCGC‑TCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGtt < 2:950017/76‑1 (MQ=255) tctcCGGTATGGCAGGTGAGGTGTCGC‑TCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTAt < 2:590392/76‑1 (MQ=255) ccGGTATGGCAGGGGAGGGGTCGC‑TCATTTTTTTTCCTCTAATGGGGTTTTTTTTGGGGACGACAAGGGTTTTggg > 2:642078/1‑75 (MQ=255) cGGTATGGCAGGTGAGGTGTCGC‑TCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCg < 2:375191/76‑1 (MQ=255) gTATGGCAGGTGAGGTGTCGC‑TCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGa > 2:1465542/1‑76 (MQ=255) tGGCAGGTGAGGTGTCGC‑TCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATgc > 1:200337/1‑76 (MQ=255) aggtgaggtgTCGC‑TCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATgcggcg < 2:1377/76‑1 (MQ=255) aggtgaggtgTCGC‑TCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATgcggcg > 1:399641/1‑76 (MQ=255) | TTAATAACCGCCTGTTGTCCGGCAGTCAGCGCTTTATGCTGATTCGCCGGATCGATACTCTCCGGTATGGCAGGTGAGGTGTCGC‑TCATTTTTCTTCCTCTAATGGGATTTTTGTAGGGTACGACAAGGTTTATGGCGGATGCGGCG > REL606/4574178‑4574324 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |