| 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 | 97.1% | 93.8 / NA | 35 | I204S (ATC→AGC) | wcaA | putative colanic acid biosynthesis glycosyl transferase |
| Reads supporting (aligned to +/- strand): ref base A (1/0); new base C (13/21); total (14/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.00e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT > NC_000913/2132647‑2132925 | ggATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGGTCTGCATTTCACCGTGATTGATATGCAGGCTc > 2:39395/1‑143 (MQ=255) agcTGGTAGAGGGGGAACAGCGGATATGTTGTGCGGGCGCGGGCGAGTTTGTCTTTGTGCTTGCGGTAAAAGTGGAGGTAGCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTg < 1:397785/141‑1 (MQ=255) cTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTcc > 1:44232/1‑143 (MQ=255) tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa > 1:380051/1‑143 (MQ=255) tatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCGCCGTGATTGATCTGCTGGCTCTGCGTCGCCCCTTCTCCTTTcctccg > 1:536938/1‑139 (MQ=255) ttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACTCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATCTGCAGGCTTTGCGTGTGCTCTTCTACATTCCAGGGttgg > 2:294679/1‑141 (MQ=255) tttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTATGCGTCGCCTCTTCTACTTTCCACGGTTcac > 1:520806/1‑141 (MQ=255) gtGGGGCGGGCGAATTTGTGTTTGTGGTTGGGGTAGAGGTGTAGGTCCCCGGAGAATTTTTTCGGTGACGTGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTTCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa < 2:195862/142‑1 (MQ=255) cTGGCGCAGTCAAATTGGTCTGTGTGCTGGCTGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa < 2:810268/143‑1 (MQ=255) tGGCGCGGTCGAATTTGTCTTTGTGATTGAGGTAAAAGTGGAAGTACCCGGAGAATGTTTTCGGTGAGGAGGTGATCTACATTGCACCGAGATTGCTGTGCAGGACCTGCGTGGCATCTTCTACTTTACGCGGTGGGCCgcag > 2:228888/1‑140 (MQ=255) gTCGGATTTCTGTGTGTGCTGGCGGTTAAGGTGGAAGTTCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacca < 2:699593/143‑1 (MQ=255) tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTGCGTTGCACGGGTAGCCGTAGTCCACCACCCCCCGCAGGa > 2:252176/1‑143 (MQ=255) tgtgCTTGCGGTAAACGTGGAAGTCCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAg < 2:368607/143‑1 (MQ=255) tgtgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCCCGGTTCGCCGTACTCCACCACCATCCCCAGGAAg > 2:163137/1‑143 (MQ=255) tgtgCTTGCGGTAAAAGTGCAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAg < 1:24011/143‑1 (MQ=255) tgtgCGTGCGGGAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAg < 1:257012/143‑1 (MQ=255) cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGt > 1:514970/1‑143 (MQ=255) cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGt > 1:511683/1‑143 (MQ=255) aaaGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTACGTCGCCTCTTCTACTTTCCACGGTTCTCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCt > 2:837554/1‑143 (MQ=255) tACCCGGAGAATTTTCTCGCTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt < 1:67989/143‑1 (MQ=255) ctctgataaTTCTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa < 1:233521/136‑1 (MQ=255) cccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa < 1:327474/143‑1 (MQ=255) cccGGAGAATTTGTTCGGTGACGAGGTGATCTGCATTTCACTGTGATTGATATTCAGGCTCTGCGTCGCCTCTTCTACTGTCCACGGTTCGCCGGACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa < 1:343515/143‑1 (MQ=255) tttttCTGTGACGAGGTTATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATc < 1:83415/143‑1 (MQ=255) gatctCGGAGACAGGGAGAACTGCAACTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATc < 1:209042/139‑1 (MQ=255) gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa < 1:244568/143‑1 (MQ=255) aCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAacac < 1:34750/143‑1 (MQ=255) aCGAGGCGCTCGGCATGTCACCGTGATGGCTCTGCCGGCTCTGCGTCGGCTCCTCTCCTTTCCACGGTTCGCCGTACTCCACCGCCATCTGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAacac < 1:358469/143‑1 (MQ=255) gTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCCGTTCTGTATCGAACAAACACTCttt > 1:65049/1‑143 (MQ=255) aCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGccc < 2:511683/143‑1 (MQ=255) aCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGccc < 2:514970/143‑1 (MQ=255) atatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAc < 2:456913/143‑1 (MQ=255) gCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTgg < 2:660293/143‑1 (MQ=255) gCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGGCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTgg < 1:360604/143‑1 (MQ=255) cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTCATCCTGTGACACTCTCACTTCGGTATCGAACACACACTCTTTGACACGCCATGCCCAGGTAAAGACCTGGt > 2:772324/1‑143 (MQ=255) | GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT > NC_000913/2132647‑2132925 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |