| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 2,132,787 | A→C | 100% | 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 | 92.3% | 65.8 / ‑5.5 | 26 | I204S (ATC→AGC) | wcaA | putative colanic acid biosynthesis glycosyl transferase |
| Reads supporting (aligned to +/- strand): ref base A (2/1); new base C (8/15); total (10/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.38e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
GGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTC‑GCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGG > NC_000913/2132653‑2132924 | ggTAGAGGGTAGACTGCTGATATTTTTCGCGGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTc < 2:569031/143‑1 (MQ=255) gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTATCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGATGAGGTGATATGCATTGCACCGTGATTGATATGAAGTATCCGCGTC‑g > 2:131522/1‑143 (MQ=255) tggagGGTAAACAGCTGTGATTTTTCGCAGGCGCGGTCGGAGTGGGGTTTGTGCGTGCGGTAAAAGTGGAAGTCCCCGGAGAATTTTTTCGGTGACGAGGTGACCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTC‑Gc < 2:353263/141‑1 (MQ=255) tAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGTATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTC‑Gc < 1:315083/143‑1 (MQ=255) gggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGGCGAGGTGTTCTGCATTTCACCGCGATTGATATGCAGGATATGCGtaggctat > 1:275157/1‑136 (MQ=255) gggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACACGGAGAATTTTTTCGGTGAGGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTATGCGTC‑GCctct > 1:444861/1‑143 (MQ=255) aaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTC‑GCCTCTTCTAc < 1:155365/143‑1 (MQ=255) accgcagATATAGTTTGCGGGCGCCCTCGAATTTGCCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTC‑GCCTCTTCTACt < 1:221695/137‑1 (MQ=255) tatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTATAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCACGCCCTGCGTC‑GCCTCTTCTTCTTTCCACgg > 2:529567/1‑143 (MQ=255) gTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTC‑GCCTCTTCTACTTTCCACGGTTCGCCGTACTccacca < 1:620963/143‑1 (MQ=255) aTTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTC‑GCCTCTTCTACTTTCCACGGTTCGCCGTACTCCCCCACCATc > 2:185381/1‑143 (MQ=255) ttCTTTGTGCTTGCGGTAACAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGTCTCTGCGTC‑GCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCa < 1:485478/142‑1 (MQ=255) ctCTTTGTTCTTGCTGTATAAGTGGAAGTACCCCGAGACTTTTTTCGGTGACGACGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTC‑GCCTCTTCTACTTTCCACTGTTCGCCGTACTCCACCACCATCCGCa < 2:270067/142‑1 (MQ=255) tgtgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTC‑GCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAg < 1:435587/143‑1 (MQ=255) gcgtGCGGTAAAAGTGTAAGTACCCGGACAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTC‑GCCTCTTCTACGTCCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATg < 2:539911/140‑1 (MQ=255) cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTC‑GCCTCTTCTACTTTCCACGCTTCGCCGTACTCCACCACCATCCGCAGGAAGATGt > 1:141440/1‑143 (MQ=255) tGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTC‑GCCTCTTCTACTTTCTACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGAGGTCg > 1:118510/1‑143 (MQ=255) gaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGGTGTGC‑CGACGATTGTCTGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCgg < 1:150170/143‑1 (MQ=255) ttttttCTGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTC‑GCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTAt < 2:303755/143‑1 (MQ=255) ggTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTC‑GCCTCTTCCACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTGTCAGTTCGGTATCGAACAAACActcgc > 2:64828/1‑141 (MQ=255) tttCACCGTGATTGATATGCAGGCTCTGCGTC‑GCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGCTATCCAACAACCACTCTTTGACCCGccct > 1:43860/1‑141 (MQ=255) ttgatATGCAGGCTCTGCGTC‑GCCTCTACTACATTCCACGGTTCGCCGTACTCCACCACCATCCGACGGAAGAAATTTTATTCCTGTGCCGCCTTCGGGTTGGTATCGAACAAACACTCTTTGACACGCCAGGCCCAGGTaaa > 2:198914/1‑143 (MQ=255) gatATGCAGGCTCTGCGTC‑GCCTCTTCTACTTTCCACGGTTCGTCGTTTTCTATCTCCGTTCGCAGGAAGATGGCGTAAGCCTGGGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGa < 1:566997/143‑1 (MQ=255) tGCTGGCTCTGCGTC‑GCCTCTTCTACGTTCCACGGTTCGCCGTACTCCACCACCGTCCGCAGGAAGATGTCGTAATTTTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTg < 2:320187/143‑1 (MQ=255) tGCAGTCTCTGCGTC‑GCCTCTTCTACTTTCCACGGTTCGCCGTACTCCTGCACCATCCTCATGAAGATGTCTTAATTCTGTGTCTCTTGCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTg < 2:264797/143‑1 (MQ=255) gCAGGCTCTGCGTC‑GCCTCTTCTACTTTGCATGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTGCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTgg < 2:68314/143‑1 (MQ=255) | GGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTC‑GCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGG > NC_000913/2132653‑2132924 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |