| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,345,435 | 1 | . | C | 32.4% | 29.4 / 52.8 | 34 | coding (415/1140 nt) | wza | lipoprotein required for capsular polysaccharide translocation through the outer membrane |
| Reads supporting (aligned to +/- strand): ref base . (16/7); new base C (0/11); total (16/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.27e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTAAGCGCCCGGTAATATCACTGCGGATTTCAGCGAGCGTTTT‑CCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTGCCGGT > minE/1345393‑1345504 | cTAAGCGCCCGGTAATATCACTGCGGATTTCAGCGAGCGTTTT‑CCCGACTACGTGGACCTTGcc < 1:1937933/64‑1 (MQ=255)cTAAGCGCCCGGTAATATCACTGCGGATTTCAGCGAGCGTTTT‑CCCGACTACGTGGACCTTGcc < 1:3233289/64‑1 (MQ=255) gcCCGGTAATATCACTGCGGATTTCAGCGAGCGTTTT‑CCCGACTACGTGGACCTTGCCGATATACGGGTaa < 1:108929/71‑1 (MQ=255) gcCCGGTAATATCACTGCGGATTTCAGCGAGCGTTTT‑CCCGACTACGTGGACCTTGCCGATATACGGGTaa < 1:845418/71‑1 (MQ=255) gcCCGGTAATATCACTGCGGATTTCAGCGAGCGTTTT‑CCCGACTACGTGGACCTTGCCGATATACGGGTaa < 1:303990/71‑1 (MQ=255) gATTTCAGCGAGCGTTTTCCCCGACTACGTGGACCTTGCCGATATACTGGTAAAACATAGTGCCGTCAGGc < 1:4029677/71‑1 (MQ=255) gATTTCAGCGAGCGTTTTCCCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGc < 1:2503899/71‑1 (MQ=255) gATTTCAGCGAGCGTTTTCCCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGc < 1:2268043/71‑1 (MQ=255) gATTTCAGCGAGCGTTTTCCCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGc < 1:1974099/71‑1 (MQ=255) gATTTCAGCGAGCGTTTTCCCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGc < 1:3153512/71‑1 (MQ=255) gATTTCAGCGAGCGTTTTCCCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGc < 1:1864541/71‑1 (MQ=255) gATTTCAGCGAGCGTTTTCCCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGc < 1:3435134/71‑1 (MQ=255) gATTTCAGCGAGCGTTTTCCCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGc < 1:3547324/71‑1 (MQ=255) gATTTCAGCGAGCGTTTTCCCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGc < 1:484851/71‑1 (MQ=255) gATTTCAGCGAGCGTTTTCCCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGc < 1:1234366/71‑1 (MQ=255) gATTTCAGCGAGCGTTTTCCCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGc < 1:1234332/71‑1 (MQ=255) gATTTCAGCGAGCGTTTT‑CCCGACTACGTGTACCTTGCCGATATACGGGTaaaa > 1:1150838/1‑54 (MQ=255) gATTTCAGCGAGCGTTTT‑CCCGACTACGTGGACCTTGCCGATATACGGGTaaaa > 1:526552/1‑54 (MQ=255) gATTTCAGCGAGCGTTTT‑CCCGACTACGTGGACCTTGCCGATATACGGGTaaaa > 1:1399230/1‑54 (MQ=255) gATTTCAGCGAGCGTTTT‑CCCGACTACGTGGACCTTGCCGATATACGGGTaaaa > 1:1419444/1‑54 (MQ=255) gATTTCAGCGAGCGTTTT‑CCCGACTACGTGGACCTTGCCGATATACGGGTaaaa > 1:1545985/1‑54 (MQ=255) gATTTCAGCGAGCGTTTT‑CCCGACTACGTGGACCTTGCCGATATACGGGTaaaa > 1:1782426/1‑54 (MQ=255) gATTTCAGCGAGCGTTTT‑CCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGc > 1:3156527/1‑62 (MQ=255) gATTTCAGCGAGCGTTTT‑CCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGCt > 1:1220134/1‑71 (MQ=255) gATTTCAGCGAGCGTTTT‑CCCGACTACGTGGACCTTGCCGATATACCGGTaaaa > 1:2046287/1‑54 (MQ=255) gATTTCAGCGAGCGATTT‑CCCGACTACGTGGACCTTGCCGATATACGGGTaaa > 1:3715684/1‑53 (MQ=255) cgagcgTTTT‑CCCGACTACGTGGACCTTGCCGATATAc > 1:2361484/1‑38 (MQ=255) cgagcgTTTT‑CCCGACTACGTGGACCTTGCCGATATAc > 1:775311/1‑38 (MQ=255) cgagcgTTTT‑CCCGACTACGTGGACCTTGCCGATATAc > 1:3391276/1‑38 (MQ=255) cgagcgTTTT‑CCCGACTACGTGGACCTTGCCGATATAc > 1:1919619/1‑38 (MQ=255) cgagcgTTTT‑CCCGACTACGTGGACCTTGCCGATATAc > 1:3147885/1‑38 (MQ=255) cgagcgTTTT‑CCCGACTACGTGGACCTTGCCGATATAc > 1:2618359/1‑38 (MQ=255) gcgTTTT‑CCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTg < 1:1956147/71‑1 (MQ=255) gcgTTTT‑CCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTg < 1:2473473/71‑1 (MQ=255) tt‑CCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTGCCGgt > 1:2526070/1‑71 (MQ=255) tt‑CCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTGCCGgt > 1:915972/1‑71 (MQ=255) | CTAAGCGCCCGGTAATATCACTGCGGATTTCAGCGAGCGTTTT‑CCCGACTACGTGGACCTTGCCGATATACGGGTAAAACATAGTGCCGTCAGGCTGTACCCAGTTGCCGGT > minE/1345393‑1345504 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |