| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 2,132,568 | (C)5→6 | coding (830/840 nt) | wcaA ← | putative colanic acid biosynthesis glycosyl transferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,132,563 | 1 | . | C | 95.5% | 64.8 / ‑3.2 | 22 | coding (835/840 nt) | wcaA | putative colanic acid biosynthesis glycosyl transferase |
| Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (9/12); total (9/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01 | |||||||||||
GCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCG‑C‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGG > NC_000913/2132434‑2132690 | gccCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCG‑c‑cc > 1:121187/1‑133 (MQ=255) cccACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCG‑c‑ccc > 2:419440/1‑133 (MQ=255) cccACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCG‑c‑ccc < 2:404601/133‑1 (MQ=255) ggACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCC‑CCCCGAATACCATCAGc > 1:100577/1‑133 (MQ=255) tttGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCC‑CCCCGAATACCATCAGCCAGACGTTTg < 1:508739/133‑1 (MQ=255) ccGAACAAAATTGAGAAACACGATCGGCAAGATCCATGCAGTACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCC‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCAc < 2:210269/133‑1 (MQ=255) aaCAAAAATGAGCAACACGATAGGCAAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCC‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGa > 2:5841/1‑133 (MQ=255) aaaaTGAGCACCACGCTAGGCAAGAACCCTGCCGCACGGGCCTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCC‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGAcagc < 2:130837/133‑1 (MQ=255) aTGAGCAACACGATAGGCAAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCC‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACagcagc < 2:16829/133‑1 (MQ=255) tGAGCAACACGATAGGCAAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCC‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCg > 1:393690/1‑133 (MQ=255) acGATAGGCAAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCC‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCaaa > 2:322101/1‑133 (MQ=255) aGGCAAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCC‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTgcg < 1:328008/133‑1 (MQ=255) gCAAGAACCATGCAGAACGGGCGTTAACTCCAGCTGTTGGCTCGCAGATCTTCCAGCATCTTTAGCGCC‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGcc < 2:215506/133‑1 (MQ=255) gAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCC‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGt < 1:419440/133‑1 (MQ=255) aTGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCC‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACTGCAGCGTCAGCAAAGTGCGCCCGGTCATAc > 1:175283/1‑133 (MQ=255) tGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCC‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACATCAGCGTCTGCACAGGGCGCCATGTCATACg > 1:179333/1‑133 (MQ=255) gcagcaCGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCC‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGt < 2:121187/133‑1 (MQ=255) gggCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCC‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGc < 2:100577/133‑1 (MQ=255) gggCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCC‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGc < 2:393690/133‑1 (MQ=255) gTAAACTCCAGCTGTAGGCGCGCACATCTTCCCGCATCTTTCGCGCC‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGAt < 2:179333/133‑1 (MQ=255) caccAGCTGCTGGCGTGCAGATCTTCCAACATTTTTAGCG‑CTCCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGt < 2:253761/131‑1 (MQ=255) tCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCC‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTa < 2:310249/133‑1 (MQ=255) cAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCGCC‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTaga > 2:304237/1‑133 (MQ=255) aTCTTTAGCGCC‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTgg > 1:331227/1‑133 (MQ=255) tCTTTAGCGCC‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGgc > 1:430850/1‑133 (MQ=255) tAGCGCC‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCgg < 1:283202/133‑1 (MQ=255) | GCCCACAGATTGTTGAGGACGTTCTTTTTGCGCCACACCGAACAAAAATGAGCAACACGATAGGCAAGAACCATGCAGCACGGGCGTAAACTCCAGCTGTTGGCGCGCAGATCTTCCAGCATCTTTAGCG‑C‑CCCCGAATACCATCAGCCAGACGTTTGCCGTTACGCACCGACAGCAGCGTCAGCAAAGTGCGCCAGGTCATACGTTTATTGCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGG > NC_000913/2132434‑2132690 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |