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 | 93.5% | 81.0 / ‑4.8 | 31 | I204S (ATC→AGC) | wcaA | putative colanic acid biosynthesis glycosyl transferase |
Reads supporting (aligned to +/- strand): ref base A (2/0); new base C (11/18); total (13/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.68e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01 |
GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT > NC_000913/2132647‑2132925 | gtaTCTGGTAGAGAGTAAACAGCTTATATTTTTTTCTGGCGCGGTTTAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTc < 1:352738/141‑1 (MQ=255) tCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGGAGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTGTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc < 1:81573/143‑1 (MQ=255) tCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGAGGTTGAATTTGTATTTGTGCATGCGGGAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc < 1:117303/143‑1 (MQ=255) tAGAGGGTAAACAGCTGATATTTTTTGCTTTTGTTGTTGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc < 1:666370/143‑1 (MQ=255) agagGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGAGATCTGCAAATCACTGTGATTGATATGCAGGCTCTGCGTCGcc > 1:783921/1‑143 (MQ=255) gggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctct < 1:659773/143‑1 (MQ=255) gCTGATATTTTTTGTTGGCGCGGTCTATTTTGTCTTTGTTCTTGCTGTAAAAGTGGAAGTACCCGGAGAATGGTGACGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTc < 1:89324/143‑1 (MQ=255) tggttttttttGCTGGCGCTGTCTTATTTGTCTTTGTGCTTTCGGTAAAAGTGGAAGTACTCGGTGAATTTTTTCGGTGACGAGGTTATCTGCATTTCACCGTGATTGATATGCAGGCTCTGTGTCGCCTCTTCTACTTTCCa < 2:14770/138‑1 (MQ=255) tatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACgg < 1:393724/143‑1 (MQ=255) ttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCAACGTGATTGATATGCAGGCTCTGCGTCACCTCTTCTACTTTCCACGGTTCg > 2:106574/1‑143 (MQ=255) tGGCGCGGTCGAATTTGTGTTTGTACTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGAGCTGCGTCGCCTCTTGTACTTTCCACGGTGCGCAGtgc > 2:285189/1‑141 (MQ=255) gTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGATGTACCCGGATAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacca < 1:203012/143‑1 (MQ=255) tCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGTCGAGGTGATCTGCATTTCACAGTGATTGATATGAAGACTCTGCGTCGACTCTTCTACTTTCCACGGTTCGCCGGACTccaccac > 1:212918/1‑143 (MQ=255) aTTTGTCTTTGTGCTTCCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGCGATCTGCATATCACCGTGATTGACACGCACGCTCTGCGTCGCCTCTCCAACTTCCCCCGGTTCACCGTACTccaccaccacc > 2:329788/1‑141 (MQ=255) tGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCTGGCTCTGCGTCGCTTCTTCTACTTTCCACGGTTCGCCGTACTTCACCACCATCCGc > 1:546020/1‑143 (MQ=255) tgtgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAg < 2:596930/143‑1 (MQ=255) tGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACAGTTCGCCGCACTCCACCACCATCCGCAAGaaaaagaag > 2:470909/1‑136 (MQ=255) gCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCCCAGGAAGCTGTCGt > 2:224416/1‑143 (MQ=255) ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCAGCGTGATTAATATACAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTaa < 1:600077/143‑1 (MQ=255) tGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACACCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGc > 1:140462/1‑143 (MQ=255) gAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATACAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCTTCCGCAGGAAGATGTCGTAATCCTGTGCCg < 2:516341/143‑1 (MQ=255) aCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGGACCCCACCACCATCCGAAGGAAGATGTCGGAAGCCCGTGACGCGTTc > 2:784492/1‑143 (MQ=255) agaATTTTTTCGGTGACGAGGTGATCTGCATTTCATCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg < 1:429923/143‑1 (MQ=255) aaattgttCTGTGCCGAGGTGATCTACATGTCGGCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCTCTGTTCGCGGTGCTCCACCCCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTGTCAGTTCGGt < 2:825097/137‑1 (MQ=255) tgttCGGTGTGGGGGTGTTCTGCATGTCACCGTTATTGTTATGTGGGCTCTGCGTCGCCTCTTCTACTTGCCACGGTGCGCCGTACTCCTCGACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCg < 2:689845/141‑1 (MQ=255) ttCGGTGCCGAGGTGGTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCATGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaa < 2:783921/143‑1 (MQ=255) gAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTATAAGTACGATATCGAACAAACActc > 1:430267/1‑143 (MQ=255) gCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTTCACGGTTTGTTGTATTCCATCACCAGCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGc < 1:735107/143‑1 (MQ=255) cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCACCCCCCGGAAGATGTCGTAATCCTGTACCGCTATCAAATCGATAACGAACAAAAACACTTAGAAACGCAAtaacaat > 1:45988/1‑137 (MQ=255) ttgacATGCAGGCTCTGCTTTGCCTCTTCTACTTTCCACAGTTCGCCGTTCTCCACCACCATCCGTAGGAAGATCGCGAAAGCCTGGGCCGATTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTaaa < 1:108574/143‑1 (MQ=255) cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGCTGATGTCTTAATCCTCTACCGCATTCAAATCAGTATCGAACAAAAACTCTTTGAAACGCAATGCCCAGGTAAAAACCTGGt > 1:749151/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 < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |