Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 1,718,597 | A→T | 10.4% | W336R (TGG→AGG) | 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,597 | 0 | A | T | 10.4% | 161.6 / 7.1 | 48 | W336R (TGG→AGG) | anmK | anhydro‑N‑acetylmuramic acid kinase |
Reads supporting (aligned to +/- strand): ref base A (23/20); new base T (2/3); total (25/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.60e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAATCGTTCGCAGCC > NC_000913/1718470‑1718744 | gAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCa > 1:334387/1‑148 (MQ=255) aaTTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTc < 2:334387/148‑1 (MQ=255) tGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGaa > 4:180427/1‑116 (MQ=255) tGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGaa < 3:180427/116‑1 (MQ=255) cGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACg > 4:503054/1‑145 (MQ=255) cGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACg < 3:503054/145‑1 (MQ=255) gAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCTGGCTTGCCTGGCgga > 3:286039/1‑94 (MQ=255) gAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCTGGCTTGCCTGGCgga < 4:286039/96‑3 (MQ=255) cGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCggtggt < 2:385539/126‑1 (MQ=255) cGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCggtggt > 1:385539/1‑126 (MQ=255) ggCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATg > 2:957786/1‑86 (MQ=255) ggCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATg < 1:957786/86‑1 (MQ=255) ggCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGgtggtggtg > 3:1052436/1‑122 (MQ=255) ggCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGgtggtggtg < 4:1052436/122‑1 (MQ=255) tGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAg > 4:501212/1‑135 (MQ=255) tGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAg < 3:501212/135‑1 (MQ=255) gTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGgtggtggtg > 1:7523/1‑111 (MQ=255) gTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGgtggtggtg < 2:7523/111‑1 (MQ=255) gAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTa > 1:60615/1‑148 (MQ=255) gCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCgg < 2:60615/148‑1 (MQ=255) aTTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTcc > 2:992183/1‑58 (MQ=255) aTTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTcc < 1:992183/58‑1 (MQ=255) aTTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATgccg > 3:245054/1‑122 (MQ=255) aTTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATgccg < 4:245054/122‑1 (MQ=255) ttCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATgc > 2:792311/1‑118 (MQ=255) ttCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATgc < 1:792311/118‑1 (MQ=255) ggTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCggtggtggt > 2:764626/1‑89 (MQ=255) ggTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCggtggtggt < 1:764626/89‑1 (MQ=255) ggTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCtgtg < 2:424863/99‑1 (MQ=255) ggTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCtgtg > 1:424863/1‑99 (MQ=255) gTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGa > 1:324364/1‑137 (MQ=255) gTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGa < 2:324364/137‑1 (MQ=255) tAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCa > 3:171584/1‑125 (MQ=255) tAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCa < 4:171584/125‑1 (MQ=255) aaTCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATgccg > 4:170515/1‑112 (MQ=255) aaTCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATgccg < 3:170515/112‑1 (MQ=255) ccAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCggt < 1:387728/73‑1 (MQ=255) ccAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCggt > 2:387728/1‑73 (MQ=255) ttcGCCTGGCTTGCCTGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTa < 4:299919/129‑1 (MQ=255) ttcGCCTGGCTTGCCTGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTa > 3:299919/3‑131 (MQ=255) ccAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTcca > 4:125751/1‑133 (MQ=255) ccAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTcca < 3:125751/133‑1 (MQ=255) ccAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAAt > 2:695920/1‑148 (MQ=255) tggcttgCCTGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCa < 4:1142363/135‑1 (MQ=255) gCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTc > 3:176611/1‑71 (MQ=255) gCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTc < 4:176611/71‑1 (MQ=255) cAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAATCGTTc > 3:9042/1‑148 (MQ=255) ccAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAATCGTTCGCAg > 2:722502/1‑148 (MQ=255) aGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAATCGTTCGCAGcc > 3:753531/1‑148 (MQ=255) | GAAAAATTCAGTTAACTCTGATATCACGGGTTAGCGGGGAAAATAGCCCCCAGTACCGTCTCCTGGCTTGCGCCAGTGACGGAAGGCAGATTTCCTGGTAATCCCGCCAGCGTCCGCCAGGCAAGCCAGGCGAAAGCCAATGCTTCCATGTCATCGCCACTAATGCCAACGGCATCGGTGGTGGTGACTTCTGTGCCTGGCAGTAATGCCGCCAGACGCGCCATGAGTAGCGGATTACGACTACCTCCACCACATACCATCAATCGTTCGCAGCC > NC_000913/1718470‑1718744 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |