| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 4,643,405 | T→C | 29.0% | H182H (CAT→CAC) | creD → | inner membrane protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 4,643,405 | 0 | T | C | 29.0% | 41.5 / 22.0 | 31 | H182H (CAT→CAC) | creD | inner membrane protein |
| Reads supporting (aligned to +/- strand): ref base T (12/10); new base C (4/5); total (16/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.04e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGC > W3110S.gb/4643337‑4643471 | gAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCatat > 1:4679617/1‑71 (MQ=255) tAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCtt > 1:1896474/1‑71 (MQ=255) tAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCtt > 1:4540754/1‑71 (MQ=255) tAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCtt > 1:4494255/1‑71 (MQ=255) tAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCtt > 1:377454/1‑71 (MQ=255) acggaacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTa > 1:867800/1‑71 (MQ=255) gaacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCt > 1:2514621/1‑71 (MQ=255) aacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTACCTg < 1:1606168/71‑1 (MQ=255) cggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCttt < 1:1066390/64‑1 (MQ=255) cggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTACCTGaa < 1:3920987/71‑1 (MQ=255) cTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGa < 1:4900555/71‑1 (MQ=255) cATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGc > 1:4159193/1‑71 (MQ=255) gAACCCGGCACCGGGTTAGAGCAAGGCTGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGa < 1:1424880/70‑1 (MQ=255) gAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTgg > 1:310758/1‑66 (MQ=255) gAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCttt > 1:3297114/1‑50 (MQ=255) gAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCACATCCCTTTa > 1:3291404/1‑51 (MQ=255) aaCCCGGCACCGGGTTAGAGCAAGGCTGGCAGGGCGTGCACATCCCTTTACCTGAAGGGGACTGGCGGAAg < 1:4596182/71‑1 (MQ=255) cGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTa > 1:4267992/1‑46 (MQ=255) cACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAAc > 1:640468/1‑70 (MQ=255) ccGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCTGAAGCAGAACCTg < 1:304517/71‑1 (MQ=255) gggTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGaa < 1:4546154/71‑1 (MQ=255) agCAAGGCGGGCAGGGCGTGCACATCCCTTTACCTGAAGGGGACTGGCGGAAg > 1:5064711/1‑53 (MQ=255) cGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGccc > 1:4415272/1‑71 (MQ=255) ggCAGGGCGTGCACATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTg < 1:3620553/71‑1 (MQ=255) gCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTg < 1:4441208/58‑1 (MQ=255) gCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGa < 1:2551789/71‑1 (MQ=255) cAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGaa > 1:2710152/1‑71 (MQ=255) ggCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGa < 1:3697698/43‑1 (MQ=255) ggCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGTCCCTGAAttt < 1:156429/71‑1 (MQ=255) cGTGCACATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAAttt > 1:2430359/1‑69 (MQ=255) tGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGc < 1:751985/71‑1 (MQ=255) tGCACATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCt < 1:484299/43‑1 (MQ=255) tatCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGaa < 1:3690248/48‑1 (MQ=255) tatCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGaa < 1:629282/48‑1 (MQ=255) | GAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGC > W3110S.gb/4643337‑4643471 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |