Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,419,494 | G→A | 40.3% | Y173Y (TAC→TAT) | wzyE ← | predicted Wzy protein involved in ECA polysaccharide chain elongation |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,419,494 | 0 | G | A | 40.3% | 23.8 / 37.5 | 57 | Y173Y (TAC→TAT) | wzyE | predicted Wzy protein involved in ECA polysaccharide chain elongation |
Reads supporting (aligned to +/- strand): ref base G (19/15); new base A (12/11); total (31/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.94e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01 |
TCAGCAAGCCAAAGGCGACCGTGCTGACGAGGAAAAACAGCCACGCTTTGCTGTCCTGGCGCAGAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGTTTTAACGCCACGCCGGAGACTT > minE/2419427‑2419558 | tCAGCAAGCCAAAGGCGACCGTGCTGACGAGGAAAAACAGCCACGCTTTGCTGTCCTGGCGCAGAAAg < 1:2345353/68‑1 (MQ=255) cAAGCCAAAGGCGACCGTGCTGACGAGGAAAAACAGCCACGCTTTGCTGTCCTGGCGCAGAAAGTAGac < 1:460383/69‑1 (MQ=255) cGTGCTGACGAGGAAAAACAGCCACGCTTTGCTGTCCTGGCGCAGAAAGTAGACCACCAGCATCGCCgg > 1:1507282/1‑69 (MQ=255) cGTGCTGACGAGGAAAAACAGCCACGCTTTGCTGTCCTGGCGCAGAAAGTAGACCACCAGCATCGCCgg > 1:1865662/1‑69 (MQ=255) cGTGCTGACGAGGAAAAACAGCCACGCTTTGCTGTCCTGGCGCAGAAAGTAGACCACCAGCATCGCCgg > 1:1963460/1‑69 (MQ=255) cGTGCTGACGAGGAAAAACAGCCACGCTTTGCTGTCCTGGCGCAGAAAGTAGACCACCAGCATCGCCgg > 1:3177878/1‑69 (MQ=255) gAGGAAAAACAGCCACGCTTTGCTGTCCTGGCGCAGAAAGTAGACCACCAGCATCGCCGGGATGaaaa > 1:1036230/1‑68 (MQ=255) gAGGAAAAACAGCCACGCTTTGCTGTCCTGGCGCAGAAAGTAGACCACCAGCATCGCCGGGATGaaaa > 1:797264/1‑68 (MQ=255) gAGGAAAAACAGCCACGCTTTGCTGTCCTGGCGCAGAAAGTAGACCACCAGCATCGCCGGGATGaaaa > 1:453863/1‑68 (MQ=255) gAGGAAAAACAGCCACGCTTTGCTGTCCTGGCGCAGAAAGTAGACCACCAGCATCGCCGGGATGaaaa > 1:3213481/1‑68 (MQ=255) gAGGAAAAACAGCCACGCTTTGCTGTCCTGGCGCAGAAAGTAGACCACCAGCATCGCCGGGATGaaaa > 1:2487969/1‑68 (MQ=255) gAGGAAAAACAGCCACGCTTTGCTGTCCTGGCGCAGAAAGTAGACCACCAGCATCGCAGGGATGaaaa > 1:89522/1‑68 (MQ=255) gAGGAAAAACAGCCACGCTTTGCTGTCCTGGCGCAGAAAGTAGACCACCAGCATAGCCGGGATGaaaa > 1:972055/1‑68 (MQ=255) ggAAAAACAGCCACGCTTTGCTGTCCTGGCGCAGAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGt < 1:6813/69‑1 (MQ=255) aGCCACGCTTTGCTGTCCTGGCGCAGAAAGTATACCACCAGCATCGCCGGGATGAAAAAGTAAaagaag < 1:1976381/69‑1 (MQ=255) aGCCACGCTTTGCTGTCCTGGCGCAGAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAaagaag < 1:1146231/69‑1 (MQ=255) aGCCACGCTTTGCTGTCCTGGCGCAGAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAaagaag < 1:2567495/69‑1 (MQ=255) gCTTTTCTATCCTGGCGCAGAAAATAAACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGtttt < 1:936344/69‑1 (MQ=255) gCTTTGCTGTCCTGGCGCAGAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGtttt < 1:1497213/69‑1 (MQ=255) gCTTTGCTATCCTGGCGCAGAAAATAAAGCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGtttt < 1:2588551/69‑1 (MQ=255) gCTTTGCTATCCTGGCGCAGAAAATAAACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGtttt < 1:3130670/69‑1 (MQ=255) gCTTTGCTATCCTGGCGCAGAAAATAAACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGtttt < 1:2574332/69‑1 (MQ=255) gCTTTGCTATCCTGGCGCAGAAAATAAACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGtttt < 1:329675/69‑1 (MQ=255) gCTTTGCTATCCTGGCGCAGAAAATAAACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGtttt < 1:451728/69‑1 (MQ=255) gCTTTGCTATCCTGGCGCAGAAAATAAACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGtttt < 1:2391559/69‑1 (MQ=255) gCTTTGCTATCCTGGCGCAGAAAATAAACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGtttt < 1:599083/69‑1 (MQ=255) gCTTTGCTATCCTGGCGCAGAAAATAAACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGtttt < 1:2247613/69‑1 (MQ=255) gCTTTGCTATCCTGGCGCAGAAAATAAACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGtttt < 1:102138/69‑1 (MQ=255) cTTTGCTATCCTGGCGCAGAAAATAAACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGtttt < 1:2758709/68‑1 (MQ=255) cTGTCCTGGCGCAGAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGTTTTAacgc < 1:2128684/68‑1 (MQ=255) cTGTCCTGGCGCAGAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGTTTTAacgc < 1:2317662/68‑1 (MQ=255) cTGGCGCAGAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGTTTTAacgccacgc > 1:3059591/1‑68 (MQ=255) cTGGCGCAGAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGTTTTAacgccacgc > 1:2210070/1‑68 (MQ=255) cTGGCGCAGAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGTTTTAacgccacgc > 1:2563120/1‑68 (MQ=255) cTGGCGCAGAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGTTTTAacgccacgc > 1:3213398/1‑68 (MQ=255) cTGGCGCAGAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGTTTTAacgccacgc > 1:870611/1‑68 (MQ=255) cTGGCGCAGAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGTTTTAacgccacgc > 1:3301518/1‑68 (MQ=255) cTGGCGCAGAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGTTTTAacgccacgc > 1:2460550/1‑68 (MQ=255) cTGGAGCAGAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGTTTTAacgccacgc > 1:3306806/1‑68 (MQ=255) tGGCGCAGAAAATAAGCCACCAGCATCGCCGGGATGAAAAAGTaa > 1:2846839/1‑45 (MQ=37) tGGCGCAGAAAATAAACCACCAGCATCGCCGGGATGAAAAAGTaa > 1:302878/1‑45 (MQ=38) tGGCGCAGAAAATAAACCACCAGCATCGCCGGGATGAAAAAGTaa > 1:904959/1‑45 (MQ=38) tGGCGCAGAAAATAAACCACCAGCATCGCCGGGATGAAAAAGTaa > 1:2229694/1‑45 (MQ=38) tGGCGCAGAAAATAAACCACCAGCATCGCCGGGATGAAAAAGTaa > 1:1278913/1‑45 (MQ=38) tGGCGCAGAAAATAAACCACCAGCATCGCCGGGATGAAAAAGTaa > 1:713985/1‑45 (MQ=38) tGGCGCAGAAAATAAACCACCAGCATCGCCGGGATGAAAAAGTaa > 1:545357/1‑45 (MQ=38) tGGCGCAGAAAATAAACCACCAGCATCGCCGGGATGAAAAAGTaa > 1:137440/1‑45 (MQ=38) tGGCGCAGAAAATAAACCACCAGCATCGCCGGGATGAAAAAGTaa > 1:412530/1‑45 (MQ=38) tGGCGCAGAAAATAAACCACCAGCATCGCCGGGATGAAAAAGTaa > 1:2011250/1‑45 (MQ=38) tGGCGCAGAAAATAAACCACCAGCATCGCCGGGATGAAAAAGTaa > 1:3242923/1‑45 (MQ=38) tGGCGCAGAAAATAAACCACCAGCACCGCCGGGATGAAAAAGTaa > 1:2031193/1‑45 (MQ=37) gAAAGTAGACCACCAGCATCTCCGGTATGAAAAAGTAAAAGAAGCGTTTTAACGCCACGCCGGAGACtt < 1:2439260/69‑1 (MQ=255) gAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGTTTTAACGCCACGCCGGAGACtt < 1:552119/69‑1 (MQ=255) gAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGTTTTAACGCCACGCCGGAGACtt < 1:3425415/69‑1 (MQ=255) gAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGTTTTAACGCCACGCCGGAGACtt < 1:1777500/69‑1 (MQ=255) gAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGTTTTAACGCCACGCCGGAGACtt < 1:2808268/69‑1 (MQ=255) gAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGTTTTAACGCCACGCCGGAGACtt < 1:2065792/69‑1 (MQ=255) aaaGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGTTTTAACGCCACGCCGGAGACtt < 1:436933/68‑1 (MQ=255) | TCAGCAAGCCAAAGGCGACCGTGCTGACGAGGAAAAACAGCCACGCTTTGCTGTCCTGGCGCAGAAAGTAGACCACCAGCATCGCCGGGATGAAAAAGTAAAAGAAGCGTTTTAACGCCACGCCGGAGACTT > minE/2419427‑2419558 |
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: ‑ |