Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,718,528 | T→C | T359A (ACG→GCG) | anmK ← | anhydro‑N‑acetylmuramic acid kinase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,718,528 | 0 | T | C | 97.2% | 118.7 / ‑4.1 | 36 | T359A (ACG→GCG) | anmK | anhydro‑N‑acetylmuramic acid kinase |
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (23/12); total (23/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.61e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01 |
AAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTC > NC_000913/1718402‑1718660 | aaGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTgcg > 2:247742/1‑140 (MQ=255) aaGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTgcg > 2:828187/1‑140 (MQ=255) ttCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTg > 1:786674/1‑140 (MQ=255) gAGGGCAATACGCCCTCCCTAACGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAg > 1:177535/1‑140 (MQ=255) tACGCCCTCCCTAACGTTCCAAGTGTAACGGCATACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTg < 1:238310/140‑1 (MQ=255) ccctAACGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATccc > 2:50368/1‑140 (MQ=255) aaCGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCa > 1:63496/1‑140 (MQ=255) aaCGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCa > 2:47012/1‑140 (MQ=255) cGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGc > 1:92275/1‑140 (MQ=255) tCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTc > 1:166682/1‑140 (MQ=255) ccAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTcc > 1:589857/1‑140 (MQ=255) cAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCg < 1:277723/140‑1 (MQ=255) cAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCg < 1:724022/140‑1 (MQ=255) cAAGTGTAACGGCAGACGCAGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCg > 1:1058518/1‑140 (MQ=255) gtAACGGCAGGTGGGGAAAAAAAAATTCAGTTAACTCTGATATCCCGGGTAAGCGGGGAAAAAAGCCCCCGGTACCCCCTCCGGGCTTGCGCCAGTGACGGAAGGGAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAgg < 2:177535/140‑1 (MQ=255) aCGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCaa > 2:123156/1‑140 (MQ=255) gACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAgg > 2:594940/1‑140 (MQ=255) gACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAgg > 2:848363/1‑140 (MQ=255) gACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAgg > 2:86260/1‑140 (MQ=255) aaaaaaTCCAGTTAACTCTGATACCACGGGTTGGCGGGGAAAATCGCCCCCGGTCCCGCCTCCTGGGTTGCGCCAGTGACGGAGGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCaa < 2:589857/139‑1 (MQ=255) gTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATg < 1:437815/140‑1 (MQ=255) tAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTc < 1:594940/140‑1 (MQ=255) ctctGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATc < 1:687809/140‑1 (MQ=255) ctctGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATc < 2:249762/140‑1 (MQ=255) tctGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCg < 1:123156/140‑1 (MQ=255) tctGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCg < 1:50368/140‑1 (MQ=255) tctGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCg < 2:63496/140‑1 (MQ=255) cGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGc > 2:154126/1‑140 (MQ=255) cGGGTTAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGc > 2:825929/1‑140 (MQ=255) ttAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTc > 2:451662/1‑122 (MQ=255) ttAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTc < 1:451662/122‑1 (MQ=255) ttAGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAAc > 1:197056/1‑140 (MQ=255) aGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACgg > 1:810399/1‑140 (MQ=255) aGCGGGGAAAATAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACgg > 2:380876/1‑140 (MQ=255) aaaaTAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCATGCGAAAGCCAATGCTTCCATGTCATCGCCACTAAt > 2:836980/1‑125 (MQ=255) aaaaTAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCATGCGAAAGCCAATGCTTCCATGTCATCGCCACTAAt < 1:836980/125‑1 (MQ=255) aaaTAGCCCCCAGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCggtg < 2:92275/140‑1 (MQ=255) aGTACCGCCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTc > 1:1093617/1‑140 (MQ=255) | AAGACTTTCCTGGTCTGGAGGGCAATACGCCCTCCCTAACGTTCCAAGTGTAACGGCAGACGCGGTAAGAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTC > NC_000913/1718402‑1718660 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |