| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 1,489,608 | 1 | . | G | 59.5% | 75.2 / 53.8 | 37 | coding (660/801 nt) | ydcF | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base . (15/0); new base G (0/22); total (15/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.07e-10 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.68e-02 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTTATCTCTCACTACTCACTGGCGAGCTGCCGCGTTTACGC‑‑‑GATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGC > W3110S.gb/1489568‑1489670 | gTTATCTCTCACTACTCACTGGCGAGCTGCCGCGTTTACGC‑‑‑GATGATAGCGATGGCTACGGTccc > 1:2898880/1‑65 (MQ=255)gTTATCTCTCACTACTCACTGGCGAGCTGCCGCGTTTACGC‑‑‑GATGATAGCGATGGCTACGGTccc > 1:2259042/1‑65 (MQ=255) ctcACTACTCACTGGCGAGCTGCCGCGTTTACGC‑‑‑GATGATAGCGATGGCTACGGTCCCCGCgg > 1:1067335/1‑63 (MQ=255) ctcACTACTCACTGGCGAGCTGCCGCGTTTACGC‑‑‑GATGATAGCGATGGCTACGGTCCCCGCgg > 1:492447/1‑63 (MQ=255) ctcACTACTCACTGGCGAGCTGCCGCGTTTACGC‑‑‑GATGATAGCGATGGCTACGGTCCCCGCgg > 1:2704500/1‑63 (MQ=255) ctcACTACTCACTGGCGAGCTGCCGCGTTTACGC‑‑‑GATGATAGCGATGGCTACGGTCCCCGCgg > 1:260677/1‑63 (MQ=255) ctcACTACTCACTGGCGAGCTGCCGCGTTTACGC‑‑‑GATGATAGCGATGGCTACGGTCCCCGCgg > 1:1086854/1‑63 (MQ=255) aCTGGCGAGCTGCCGCGTTTACGC‑‑‑GATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTAt > 1:1438317/1‑65 (MQ=255) aCTGGCGAGCTGCCGCGTTTACGC‑‑‑GATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTAt > 1:1617970/1‑65 (MQ=255) aCTGGCGAGCTGCCGCGTTTACGC‑‑‑GATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTAt > 1:1045122/1‑65 (MQ=255) aCTGGCGAGCTGCCGCGTTTACGC‑‑‑GATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTAt > 1:2068236/1‑65 (MQ=255) aCTGGCGAGCTGCCGCGTTTACGC‑‑‑GATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTAt > 1:1811736/1‑65 (MQ=255) aCTGGCGAGCTGCCGCGTTTACGC‑‑‑GATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTAt > 1:1827547/1‑65 (MQ=255) aCTGGAGAGCTGCCGCGTTTACGC‑‑‑GATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTAt > 1:1949561/1‑65 (MQ=255) aGCTGCCGCGTTTACGC‑‑‑GATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCAc > 1:2677154/1‑65 (MQ=255) gCCGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:2620436/65‑1 (MQ=255) gCCGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:706323/65‑1 (MQ=255) gCCGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:2830546/65‑1 (MQ=255) gCCGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:2514124/65‑1 (MQ=255) gCCGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:2924298/65‑1 (MQ=255) gCCGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:3010691/65‑1 (MQ=255) gCCGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:366295/65‑1 (MQ=255) gCCGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:426465/65‑1 (MQ=255) gCCGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:503627/65‑1 (MQ=255) gCCGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:533738/65‑1 (MQ=255) gCCGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:603030/65‑1 (MQ=255) gCCGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:629395/65‑1 (MQ=255) gCCGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:635320/65‑1 (MQ=255) gCCGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:123820/65‑1 (MQ=255) gCCGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:2277507/65‑1 (MQ=255) gCCGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:1489653/65‑1 (MQ=255) gCCGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:2045098/65‑1 (MQ=255) gCCGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:1498855/65‑1 (MQ=255) gCCGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:1651961/65‑1 (MQ=255) ccGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:175609/64‑1 (MQ=255) ccGCGTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:373597/64‑1 (MQ=255) cgcgTTTACGCGATGATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACg < 1:1799823/63‑1 (MQ=255) c‑‑‑gATTATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGc > 1:2450763/1‑63 (MQ=255) c‑‑‑gATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGc > 1:313162/1‑63 (MQ=255) c‑‑‑gATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGc > 1:1281473/1‑63 (MQ=255) c‑‑‑gATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGc > 1:1556898/1‑63 (MQ=255) c‑‑‑gATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGc > 1:1702597/1‑63 (MQ=255) c‑‑‑gATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGc > 1:1707620/1‑63 (MQ=255) c‑‑‑gATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGc > 1:2464327/1‑63 (MQ=255) c‑‑‑gATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGc > 1:1903039/1‑63 (MQ=255) c‑‑‑gATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGc > 1:192299/1‑63 (MQ=255) c‑‑‑gATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGc > 1:2276293/1‑63 (MQ=255) c‑‑‑gATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGc > 1:2456378/1‑63 (MQ=255) c‑‑‑gATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGc > 1:2579888/1‑63 (MQ=255) | GTTATCTCTCACTACTCACTGGCGAGCTGCCGCGTTTACGC‑‑‑GATGATAGCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGC > W3110S.gb/1489568‑1489670 |
| 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: ‑ |