| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 2,132,787 | A→C | I204S (ATC→AGC) | 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,787 | 0 | A | C | 96.3% | 62.5 / NA | 28 | I204S (ATC→AGC) | wcaA | putative colanic acid biosynthesis glycosyl transferase |
| Reads supporting (aligned to +/- strand): ref base A (1/0); new base C (13/13); total (14/14) | |||||||||||
| 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.75e-01 | |||||||||||
AACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCC > NC_000913/2132664‑2132904 | aaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc < 1:69491/129‑1 (MQ=255) gCTGATATTTTTTGCTGGCGCTGTCGAATTAGTCTTAAAAATTACAAAAAAAGTGAAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCg > 1:162533/1‑129 (MQ=255) atTTTTTGCTGGCGAGGACGGATTTGTCTTTGTGCTTGCGGTAAAAGGGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGACGGCGCTGcgccgcagctt > 2:563506/1‑124 (MQ=255) cgtgcccgGTCTAATTTGTCTTTGTGCTTGCGGTAAATGTGGAAGTATCCGGAGAAGTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTc < 2:488487/123‑1 (MQ=255) tGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCAGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTcc > 1:407264/1‑129 (MQ=255) ggCGCGGTCGCATTTGTCTTTGTGCTTGCGGTAAAAGTAAAAGTACCCGGAGAATATTTTCGATGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGGCGCGTCTGCTACTTTCCa > 1:487246/1‑129 (MQ=255) gcgGTCGATTTTGTCTTTGTGCTTGCAATAAAAGACGAAATAAAAGAAAAATTTTTTCGGTAAAGAGGTGATCTGACTTTCATCGTTATATATATGCAGGCTCTGCGTCGCCTGTTCTACTTTCGACgg > 1:575177/1‑129 (MQ=255) ggTCGCATTTGTCTTTGTGCTTGCGGTAAAAAAAGAAGTACCAGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGtt > 1:1211/1‑129 (MQ=255) ttGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGGGATTGATATGAAGACTCTGCGTCGCGCCCTCCCCTTTCCAAGGTTGGCCGtgc > 2:154710/1‑127 (MQ=255) gTCTTTGTGCTTGCAGTCAACGTGGTAGTACTCGGAGTATTTTTCCGGTGACGAGGTGATCTGCTTTTCACTTTTATTTTTTTGTTTTTTTTGTGTCGCTTCTTCTACTTGCCACGGTTCGCCGTACTc < 1:598676/129‑1 (MQ=255) gTCTTTGCTCTGTCGGTAAACGTGGACGTACCCGGAGCATTTTTTCGGTGACGAGGGTATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc < 2:153447/129‑1 (MQ=255) ttGCGGTATAGGTGGAAGTGCCCGGTGAATTTGTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCGACTTTCCACGGTTCGCCGTACTCCACCACCATc < 2:373062/129‑1 (MQ=255) ttGCGGTAACCGTGGAAGTACCCGGAGAATTTTTTAGATGAAGAAATAAACTACATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATc > 1:278489/1‑129 (MQ=255) tAAAAGTTGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACTGTGATTGATATGCAGGCTCTGCGTTGCTTTTTCTACTTTCTACGGTTCGCCGTACTCCACCACCATCCGCAgg < 1:617270/129‑1 (MQ=255) tcgggtaGTGCCCGGAGAGTTGTTTCGGTGGCGGGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGGACTCCACCACCATCCGCAGGAAGa < 2:449237/123‑1 (MQ=255) aGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGTCTCTGCGTCGTCTTTTTTTTTTTCTTTGTTTCGCCGTACTCCACCACCATCCGCAGGAAGa < 1:21857/129‑1 (MQ=255) ggAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTAGTCAACTATCCACGGTACATCGGCCGCCACCACCATCCTCAGGaacaagt > 2:221502/1‑124 (MQ=255) agaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCTTCTTTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg < 1:86353/129‑1 (MQ=255) aTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGGCGCCACTTCTACTTTCCACGGTTCGACGTACTCCACCACCATCCGAAAGAAGAAGTCGAAATACCGTGGCg > 2:422785/1‑129 (MQ=255) ttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCtt < 2:1211/129‑1 (MQ=255) tCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAAACCTGTGCCGCTTCCa > 2:509553/1‑129 (MQ=255) cGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATACAAAATATAAAAAACATATTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAg > 1:633032/1‑129 (MQ=255) aGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGACGTCGTACTCCTGCGCCGCTTCCAGGTCGGTAt > 2:567830/1‑129 (MQ=255) gCATTTCACCGTGATTGATATACAAGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCTCCTTCTGTAGGAAGATGTTGTTATTCTGTGCCGCTTGCAGTGCGGTATCGaacaaac < 1:376142/129‑1 (MQ=255) gTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGc > 2:259886/1‑50 (MQ=255) gTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGc < 1:259886/50‑1 (MQ=255) tggcgctGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACGCCACCGCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACTAACACTCTTTGAAACGc < 2:204901/123‑1 (MQ=255) tgatATGCAGGCTCTGCGTCGCCTCTTCTCCTTCCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTTGGTATCGAACAAACACTCTGTGAAACGcc < 1:488976/129‑1 (MQ=255) | AACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCC > NC_000913/2132664‑2132904 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |