| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 4,528,968 | 0 | A | G | 73.0% | 61.4 / 22.0 | 37 | L242P (CTA→CCA) | yjhH | predicted lyase/synthase |
| Reads supporting (aligned to +/- strand): ref base A (0/10); new base G (12/15); total (12/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.51e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACG > W3110S.gb/4528913‑4529033 | tttGTTCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCg > 1:2612832/1‑71 (MQ=255)tttGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCg > 1:3647783/1‑71 (MQ=255) agaCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTc > 1:957107/1‑71 (MQ=255) agaCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTc > 1:3548612/1‑71 (MQ=255) agaCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTc > 1:2256650/1‑71 (MQ=255) agaCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGATTTTTATTCAGCGTCGCAGCGGTc > 1:3893962/1‑71 (MQ=255) ggTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTa < 1:759580/45‑1 (MQ=255) ggTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTa < 1:594686/45‑1 (MQ=255) ggTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTa < 1:1933312/45‑1 (MQ=255) gtgtTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCg < 1:76156/71‑1 (MQ=255) gtgtTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCg < 1:360307/71‑1 (MQ=255) gtgtTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTa < 1:884426/44‑1 (MQ=255) gtgtTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTa < 1:874215/44‑1 (MQ=255) aaTAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGgcgc < 1:3473516/71‑1 (MQ=255) aaTAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGgcgc < 1:65511/71‑1 (MQ=255) aTAGCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGgcgc < 1:191291/70‑1 (MQ=255) gCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGAt > 1:1263112/1‑49 (MQ=255) gCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGAt > 1:1300200/1‑49 (MQ=255) gCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGAt > 1:4103195/1‑49 (MQ=255) gCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTa < 1:1140141/71‑1 (MQ=255) gCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTa < 1:2071018/71‑1 (MQ=255) gCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTa < 1:616315/71‑1 (MQ=255) gCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTa < 1:1317209/71‑1 (MQ=255) gCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTa < 1:4056643/71‑1 (MQ=255) gCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTa < 1:3913881/71‑1 (MQ=255) gCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTa < 1:1440440/71‑1 (MQ=255) gCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTa < 1:1066243/71‑1 (MQ=255) gCGGGCAGTTGTAGTGGTTTTTTATTCAGCGTAGCAGCGGTCGCCAGAt > 1:3289627/1‑49 (MQ=39) cGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTa < 1:2589108/70‑1 (MQ=255) cGGGCAGTTGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTa < 1:2977809/70‑1 (MQ=255) ttGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAg < 1:1535057/39‑1 (MQ=255) ttGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAg < 1:1292225/39‑1 (MQ=255) ttGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAg < 1:4137391/39‑1 (MQ=255) ttGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAg < 1:1484906/39‑1 (MQ=255) ttGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAg < 1:3561134/39‑1 (MQ=255) tGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGAc > 1:2949118/1‑71 (MQ=255) tGTAGTGGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGAc > 1:642250/1‑71 (MQ=255) gtagtgGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCgg > 1:269840/1‑59 (MQ=255) gtagtgGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCgg > 1:2409273/1‑59 (MQ=255) gtagtgGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCgg > 1:942654/1‑59 (MQ=255) gtagtgGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCgg > 1:1204309/1‑59 (MQ=255) gtagtgGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCgg > 1:3518688/1‑59 (MQ=255) gtagtgGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCgg > 1:1293815/1‑59 (MQ=255) gtagtgGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCgg > 1:1595494/1‑59 (MQ=255) gtagtgGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCgg > 1:44773/1‑59 (MQ=255) gtagtgGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCgg > 1:3763709/1‑59 (MQ=255) gtagtgGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACg > 1:1409367/1‑71 (MQ=255) gtagtgGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACg > 1:3061021/1‑71 (MQ=255) gtagtgGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACg > 1:391661/1‑71 (MQ=255) gtagtgGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACg > 1:3763751/1‑71 (MQ=255) | TTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACG > W3110S.gb/4528913‑4529033 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |