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 | 96.5% | 94.9 / ‑5.8 | 35 | I204S (ATC→AGC) | wcaA | putative colanic acid biosynthesis glycosyl transferase |
Reads supporting (aligned to +/- strand): ref base A (2/0); new base C (17/15); total (19/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.92e-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: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACC > NC_000913/2132647‑2132921 | ggATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGAATTTAACCGTGATTGATATGCAGCATc > 2:124405/1‑143 (MQ=255) tAGAGGGTAAACAGTTGAGGTTTTTTGACGGAGTGCTCGAATTTGTGTTGGTGCTTGCGGTCAAAGTGGTAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc < 2:8029/143‑1 (MQ=255) aggcaaaaCGGCCGAGATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctc < 2:649954/138‑1 (MQ=255) agGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctc > 1:716293/1‑143 (MQ=255) ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAACCGTGGAAGTACCCGGAGCATTTTTTCGGTGACGAGGTGATCTGCCTTTCACCGTGATTGATCTGCAGGCCATGCCTCGCctctg > 2:529570/1‑142 (MQ=255) ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCcctt > 2:433912/1‑143 (MQ=255) aagcAGCTGATATTTTTTGCTGGCGCTGTCGAATGTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCTGTTATGAGTTGATCTGCATTTCACTGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa < 1:136369/140‑1 (MQ=255) gATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCAc > 1:429200/1‑143 (MQ=255) ttgtggtGGCGGGGTCGAATTTGTCTGTGTGCTGGTGGTAAAGGTGGAAGTAGCCGGAGAATTTTTTCGCTGACGAGTTGATCTGCATTTCACCGTGATTGATATGCGGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCg < 2:356249/137‑1 (MQ=255) tttGGTGGGGCGGTCGGATTTGTGTTTGTGCTTGCGGTAAGAGTGGAAGTACCCGGAGAATTTTTTAGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGc < 2:199194/143‑1 (MQ=255) cTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa < 1:564638/143‑1 (MQ=255) cTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGGGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa < 2:235726/143‑1 (MQ=255) tGGCGTGGTCGACTTTGTCTTGTTGCTGGCGGTAAAAGTGGACGTACCCGGGGAATTTTTTCGGTGACGAGGTGACCTGCATGTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAc < 2:840082/143‑1 (MQ=255) tcgGTCTAATTTGTCTTTGTGCTTGCGGTAACAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTcca < 1:519593/142‑1 (MQ=255) gcgGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTcca < 1:516761/143‑1 (MQ=255) gTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCCCCGTGCTccaccc > 2:801277/1‑142 (MQ=255) tCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTctaccaa > 1:451984/1‑142 (MQ=255) gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTACGCAGTACTccaccacca > 2:34298/1‑143 (MQ=255) gcgctccaaGCGGCAGTACCCGCAGACTTTTTTCGCCGACGACGTTATCTGCATTTCCCCGTGATTGCTATGCAGGCCCTGCGTCGCCTCTTCTACTTTCCACGGCTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGt < 2:397187/136‑1 (MQ=255) ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTaa > 1:208544/1‑143 (MQ=255) gaaaaaGTGGAAGTACCCGGAGAATTTTTTCGGGGACGAGGTGATCTGGATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACATTCCACGGTTCGGCGTACTCAACCACCATCCCCGGGAAGATGTCGTAAt > 2:78423/3‑143 (MQ=255) gTGGAAGTACCCGGAGAATTTTTTCGGTGGCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg < 2:291124/143‑1 (MQ=255) gTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTGCTCCGCCACCATCCGCAGGAAGATGTCGTAGTCCtgtg > 1:360057/1‑143 (MQ=255) tACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAATTTACCTTCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt < 1:804166/143‑1 (MQ=255) gagaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGAAGATGTCGTGATCTTGTGACGCTTTCAGTTc > 1:482000/1‑143 (MQ=255) aggatgttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACGTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg < 1:273779/137‑1 (MQ=255) ggTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGGGCCGCTTTCCGTTCCGTATCgcgcaa > 1:703435/1‑138 (MQ=255) gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACGACCATCCGCAGGAAGATGTCGTAATCCTGTGGCGCTTTCAGTTCGGTATCGaacacg > 1:556065/1‑141 (MQ=255) cctcgtctcCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCtt < 1:23736/138‑1 (MQ=255) tCTGCATTTCACCGAGATTGATAGGCCGGTTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGaaa < 2:591039/143‑1 (MQ=255) tGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGACGATGTCGTAATCCCGTGCACCCTTCAGTTCGGTCTCGAACCAACGCCGATTGAAAca > 2:672870/1‑142 (MQ=255) tGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTCATCCTGTGCCGCTTTCAGTTCGGTATCTACCCAACACTCTTTGAAACg > 1:668410/1‑143 (MQ=255) gTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGAAATCCACTACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGCACAAACACCCTTTGAAACGCCATGCCCAgg > 1:801328/1‑143 (MQ=255) gatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCACGTCAAGa > 1:546316/1‑143 (MQ=255) tatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAcc > 1:769381/1‑143 (MQ=255) | GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACC > NC_000913/2132647‑2132921 |
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: ‑ |