Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,528,342 | G→C | 22.4% | G214A (GGC→GCC) | ftsA → | cell division protein FtsA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,528,342 | 0 | G | C | 22.4% | 24.6 / 7.2 | 18 | G214A (GGC→GCC) | ftsA | cell division protein FtsA |
Reads supporting (aligned to +/- strand): ref base G (7/7); new base C (2/2); total (9/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-01 |
GCCGCCCAGAACATCGAGAAGTGTGTGCGCCGCTGCGGCCTGGAAATCGACGACATCATCCTCGAGCAGCTGGCTTCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGACCAGGTGACCAACGACATCGCCATGGCCCTGCGTACCCCGACCCAGTACGCCGAAGAAATCAAGATC > NC_002947/1528200‑1528478 | gccgccCAGAACATCGAGAAGTGTGTGCGCCGCTGCGGCCTGGAAATCGACGACATCATCCTCGAGCAGCTGGCTTCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCaccacc > 1:790142/1‑150 (MQ=255) cgccCAGAACATCGAGAAGTGTGTGCGCCGCTGCGGCCTGGAAATCGACGACATCATCCTCGAGCAGCTGGCTTCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGa < 1:567039/150‑1 (MQ=255) ggccAGAACATCGAGAAGTGTGTGCGCCGCTGCGGCCTGGAAATCGACGACATCATCCTCGAGCAGCTGGCTTCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGTGCCGCCACCGat < 2:219664/148‑2 (MQ=255) cAGAACATCGAGAAGTGTGTGCGCCGCTGCGGCCTGGAAATCGACGACATCATCCTCGAGCAGCTGGCTTCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGACATc > 2:621596/1‑150 (MQ=255) gatgacATCATCCTCGAGCAGCTGGCTTCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATc < 2:759481/147‑1 (MQ=255) atcatcCTCGAGCAGCTGGCTTCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGACACCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGAtt < 2:790142/150‑1 (MQ=255) cGAGCAGCTGGCTTCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGa < 1:46095/150‑1 (MQ=255) ttCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGACCAGGTGACCaa > 1:368765/1‑150 (MQ=255) ggCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGTGCCGCCACCGat < 1:632295/76‑2 (MQ=255) ggCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGTGCCGCCACCGat > 2:632295/1‑75 (MQ=255) aaaaaGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCACGGAAATTGGCGACCAGGTGACCAACGACATCGCCATGGCCCTGCGTACCCCGa > 1:250726/1‑150 (MQ=255) aaaaaGAGCTGGGCGTGTGCCTGGTGGACATCCGTAGCCGGCA‑GTTCGAGCTGACCGGCTTCACCGAGGGGTGCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGACCAGGTTACGCACCACCTCGCCATGGCCCTGCGTACCCCGa > 1:337503/1‑150 (MQ=255) aaaaGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCTACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGACCAGGAATGCAACGACATCGCCATGGCCCTGCGTACCCCGAc < 1:585651/150‑1 (MQ=255) aaaGAGCTGGGCGTTTGCCTGTTGGACATCGGT‑GGCGGCACCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGACCAGGTGACCAACGACATCGCCATGGCCCTGCGTACCCCGAcc < 2:337503/150‑1 (MQ=255) cTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGACCAGGTGACCAACGACATCGCCATGGCCCTGCGTACCCCGACCCAGTAc > 1:730802/1‑150 (MQ=255) tGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGACCAGGTGACCAACGACATCGCCATGGCCCTGCGTACCCCGACCCAGTACg > 1:272352/1‑150 (MQ=255) gCCTGGTGGACATCGGT‑GGCGGCACCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGACCAGGTGACCAACGACATCGCCATGGCCCTGCGTACCCCGACCCAGTACGCCGAAGaaa < 2:730802/150‑1 (MQ=255) gatgtCGGT‑GGTGCCGCCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGACCAGGTGACCAACGACATCGCCATGGCCCTGCGTACCCCGACCCAGTACGCCGAAGAAATCAAGATc > 1:298754/5‑150 (MQ=255) | GCCGCCCAGAACATCGAGAAGTGTGTGCGCCGCTGCGGCCTGGAAATCGACGACATCATCCTCGAGCAGCTGGCTTCGGCCTACTCGGTACTGACCGATGATGAAAAAGAGCTGGGCGTGTGCCTGGTGGACATCGGT‑GGCGGCACCACCGACATCGCCATCTTCACCGAGGGTGCCATCCGCCACACGGCGGTCATCCCGATTGCTGGCGACCAGGTGACCAACGACATCGCCATGGCCCTGCGTACCCCGACCCAGTACGCCGAAGAAATCAAGATC > NC_002947/1528200‑1528478 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |