| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 155,463 | 0 | A | T | 73.2% | 123.0 / 53.0 | 71 | *247K (TAA→AAA) | ecpD | predicted periplasmic pilin chaperone |
| Reads supporting (aligned to +/- strand): ref base A (0/19); new base T (48/4); total (48/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.02e-13 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGC > W3110S.gb/155422‑155518 | gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa > 1:2525476/1‑62 (MQ=255)gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa > 1:2752550/1‑62 (MQ=255)gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa > 1:2952428/1‑62 (MQ=255)gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa > 1:1584004/1‑62 (MQ=255)gTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCaa > 1:2009037/1‑62 (MQ=255) gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg > 1:629215/1‑40 (MQ=255) gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg > 1:163214/1‑40 (MQ=255) gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg > 1:296513/1‑40 (MQ=255) gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg > 1:2693008/1‑40 (MQ=255) gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg > 1:1245213/1‑40 (MQ=255) gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg > 1:814563/1‑40 (MQ=255) gCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg > 1:3105383/1‑40 (MQ=255) tCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaa < 1:37877/61‑1 (MQ=255) tGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTc > 1:2118847/1‑62 (MQ=255) tGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTc > 1:2619712/1‑62 (MQ=255) tGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTc > 1:898911/1‑62 (MQ=255) tGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTc > 1:598157/1‑62 (MQ=255) tGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTc > 1:1842160/1‑62 (MQ=255) tGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTc > 1:1005332/1‑62 (MQ=255) tGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTc > 1:3025451/1‑62 (MQ=255) tGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTc > 1:1269347/1‑62 (MQ=255) tGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTc > 1:307320/1‑62 (MQ=255) tGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTc > 1:3099378/1‑62 (MQ=255) tGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTc > 1:1128530/1‑62 (MQ=255) cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTg > 1:740561/1‑47 (MQ=255) cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa > 1:2359400/1‑61 (MQ=255) cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa > 1:273898/1‑61 (MQ=255) cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa > 1:2537864/1‑61 (MQ=255) cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa > 1:2534426/1‑61 (MQ=255) cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa > 1:518645/1‑61 (MQ=255) cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa > 1:213511/1‑61 (MQ=255) cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa > 1:212864/1‑61 (MQ=255) cGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCa > 1:2437568/1‑61 (MQ=255) ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg < 1:1264575/62‑1 (MQ=255) ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg < 1:3066711/62‑1 (MQ=255) atatCCTTGCCAGATTAAAGGTTGGCGTTGCCTTCAATAGATCCACCAAAATCATTAATGGc < 1:2853834/62‑1 (MQ=255) atatCCTTGCCAGATTAAAGGTTGGCGTTGCCTTCAATAGATCCACCAAAATCATTAATGGc < 1:1051166/62‑1 (MQ=255) atatCCTTGCCAGATTAAAGGTTGGCGTTGCCTTCAATAGATCCACCAAAATCATTAATGGc < 1:420463/62‑1 (MQ=255) atatCCTTGCCAGATTAAAGGTTGGCGTTGCCTTCAATAGATCCACCAAAATCATTAATGGc < 1:1080205/62‑1 (MQ=255) atatCCTTGCCAGATTAAAGGTTGGCGTTGCCTTCAATAGATCCACCAAAATCATTAATGGc < 1:110329/62‑1 (MQ=255) atatCCTTGCCAGATTAAAGGTTGGCGTTGCCTTCAATAGATCCACCAAAATCATTAATGGc < 1:1226801/62‑1 (MQ=255) atatCCTTGCCAGATTAAAGGTTGGCGTTGCCTTCAATAGATCCACCAAAATCATTAATGGc < 1:1424373/62‑1 (MQ=255) atatCCTTGCCAGATTAAAGGTTGGCGTTGCCTTCAATAGATCCACCAAAATCATTAATGGc < 1:1459250/62‑1 (MQ=255) atatCCTTGCCAGATTAAAGGTTGGCGTTGCCTTCAATAGATCCACCAAAATCATTAATGGc < 1:1695561/62‑1 (MQ=255) atatCCTTGCCAGATTAAAGGTTGGCGTTGCCTTCAATAGATCCACCAAAATCATTAATGGc < 1:901979/62‑1 (MQ=255) atatCCTTGCCAGATTAAAGGTTGGCGTTGCCTTCAATAGATCCACCAAAATCATTAATGGc < 1:2849673/62‑1 (MQ=255) atatCCTTGCCAGATTAAAGGTTGGCGTTGCCTTCAATAGATCCACCAAAATCATTAATGGc < 1:711417/62‑1 (MQ=255) atatCCTTGCCAGATTAAAGGTTGGCGTTGCCTTCAATAGATCCACCAAAATCATTAATGGc < 1:1989009/62‑1 (MQ=255) atatCCTTGCCAGATTAAAGGTTGGCGTTGCCTTCAATAGATCCACCAAAATCATTAATGGc < 1:2009356/62‑1 (MQ=255) atatCCTTGCCAGATTAAAGGTTGGCGTTGCCTTCAATAGATCCACCAAAATCATTAATGGc < 1:2124371/62‑1 (MQ=255) atatCCTTGCCAGATTAAAGGTTGGCGTTGCCTTCAATAGATCCACCAAAATCATTAATGGc < 1:816318/62‑1 (MQ=255) atatCCTTGCCAGATTAAAGGTTGGCGTTGCCTTCAATAGATCCACCAAAATCATTAATGGc < 1:2485423/62‑1 (MQ=255) tatCCTTGCCAGATTAAAGGTTGGCGTTGCCTTCAATAGATCCACCAAAATCATTAATGGc < 1:2888294/61‑1 (MQ=255) tatCCTTGCCAGATTAAAGGTTGGCGTTGCCTTCAATAGATCCACCAAAATCATTAATGGc < 1:590451/61‑1 (MQ=255) tGCCTGATTTCCGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa > 1:1929423/1‑62 (MQ=255) tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCTCCAAAGTCATTAATGGCGTaaaaa > 1:1825447/1‑62 (MQ=255) tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa > 1:786928/1‑62 (MQ=255) tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa > 1:559134/1‑62 (MQ=255) tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa > 1:980006/1‑62 (MQ=255) tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa > 1:362365/1‑62 (MQ=255) tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa > 1:328906/1‑62 (MQ=255) tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa > 1:2848822/1‑62 (MQ=255) tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa > 1:2785328/1‑62 (MQ=255) tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa > 1:2736282/1‑62 (MQ=255) tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa > 1:2564451/1‑62 (MQ=255) tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa > 1:2271616/1‑62 (MQ=255) tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa > 1:2244118/1‑62 (MQ=255) tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa > 1:2143625/1‑62 (MQ=255) tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaaa > 1:1527990/1‑62 (MQ=255) tGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaaa > 1:1708097/1‑61 (MQ=255) tGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGc < 1:2921289/61‑1 (MQ=255) | GTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGC > W3110S.gb/155422‑155518 |
| 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: ‑ |