Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,253,967 | G→A | 14.3% | K81K (AAG→AAA) | leuB → | 3‑isopropylmalate dehydrogenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,253,967 | 0 | G | A | 14.3% | 62.2 / 3.6 | 28 | K81K (AAG→AAA) | leuB | 3‑isopropylmalate dehydrogenase |
Reads supporting (aligned to +/- strand): ref base G (12/12); new base A (2/2); total (14/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01 |
CTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAAGATCCGCTCGCAACTG‑GGCCTGTTCGCCAAC‑CTGCGCCCGGCCATCCTTTACCCGCAACTGGCCGATGCCTCGTCGCTGAAGCCGGAAATCGTGTC > NC_002947/2253824‑2254101 | atCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATcgag > 2:289880/2‑150 (MQ=255) gCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATcgagcgag < 2:465897/150‑1 (MQ=255) gcACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGa > 2:335245/1‑150 (MQ=255) gacgTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAgcgc > 1:406370/1‑150 (MQ=255) gCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATcga > 1:193065/1‑118 (MQ=255) gCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATcga < 2:193065/118‑1 (MQ=255) ggTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAAg > 1:214299/1‑150 (MQ=255) ccgccATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAAGATcc < 1:335245/150‑1 (MQ=255) aCGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAAGATCCGCTCGCAACTG‑GGcc < 2:395715/150‑1 (MQ=255) tGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGaa < 1:112125/118‑1 (MQ=255) tGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGaa > 2:112125/1‑118 (MQ=255) tGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAAGATCCGCTCGCAACTG‑GGCCTGTTCGCCAAc‑c < 2:222930/150‑1 (MQ=255) ccTGGAGCGTGCGCGCAACGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGAACGAGCGAGACATTCGCCCGGAGCGCCGCCTGCTGAAGATCTGCTCGCAACTG‑GGCCTGTTCTCCAAC‑CTGCGCCCGGCCAt > 2:354292/1‑150 (MQ=255) cTGGAGCGTGCGCGCAAGGTCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAAGATCCGCTCGCAACTG‑GGCCTGTTCGCCAAC‑CTGCGCCCGGCCATc < 2:406370/150‑1 (MQ=255) gcgcAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAAGATc > 1:413316/1‑97 (MQ=255) gcgcAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAAGATc < 2:413316/97‑1 (MQ=255) aaGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAAGATCCGCTCGCAACTG‑GGCCTGTTCGCCAAC‑CTGCGCCCGGCCATCCTTTACCCGCAACta > 1:93883/1‑149 (MQ=255) cGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAAATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAAGATCCGCTCGCAACTG‑GGCCTGTTCGCCAAC‑CTGCGCCCCGccctt > 1:285601/1‑127 (MQ=255) cGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAAATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAAGATCCGCTCGCAACTG‑GGCCTGTTCGCCAAC‑CTGCGCCCCGccctt < 2:285464/130‑4 (MQ=255) cGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAAATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAAGATCCGCTCGCAACTG‑GGCCTGTTCGCCAAC‑CTGCGCCCCGccctt < 2:285601/130‑4 (MQ=255) cGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAAATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAAGATCCGCTCGCAACTG‑GGCCTGTTCGCCAAC‑CTGCGCCCCGccctt > 1:285464/1‑127 (MQ=255) tgGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAAGATCCGCTCGCAACTG‑GGCCTGTTCGCCAAC‑CTGCGCCCGGCCATCCTTTACCCGCAACTGGCCGATGCCTCGTCGCTGa > 1:519036/1‑150 (MQ=255) gggCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAAGATCCGCTCGCAACTG‑GGCCTGTTCGCCAAC‑CTGCGCCCGGCCATCCTTTACCCGCAACTGGCCGATGCCTCGTCGCTGAAGc < 2:214299/150‑1 (MQ=255) cTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAAGATCCGCTCGCAACTG‑GGCCTGTTCGCCAAC‑CTGCGCCCGGCCATCCTTTACCCGCAACTGGCCGATGCCTCGTCGCTGAAGCCgg < 2:551531/150‑1 (MQ=255) gggcgggcCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAAGATCCGCTCGCAACTG‑CGCCTGTTCGCCAACGCCGC‑CCCGGCCAGCCGTTACCCGCAACTGGCCGATGCCTCGTCGCTGAAGCCGGAAAt < 1:537723/150‑1 (MQ=255) gcCGAAGTGGGACCCGATCGAGCTAAACATTCTCCCGGAGCGCGGCCTGCTGAAGATCCGCTCG‑AATTGCGTCCTGTTCGCCAAC‑CTGCGCCCGGCCATCCTTTACCCGCAACTGGCCGATGCCTCGTCGCTGAAGCCGGAAATCGTGTc < 1:289880/150‑1 (MQ=255) gcCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAAGATCCGCTCGCAACTG‑GGCCTGTTCGCCAAC‑CTGCGCCCGGCCATCCTTTACCCGCAACTGGCCGATGCCTCGTCGCTGAAGCCGGAAATCGTGTc > 1:535097/1‑150 (MQ=255) gcCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAAGATCCGCTCGCAACTG‑GGCCTGTTCGCCAAC‑CTGCGCCCGGCCATCCTTTACCCGCAACTGGCCGATGCCTCGTCGCTGAAGCCGGAAATCGTGTc > 2:607588/1‑150 (MQ=255) | CTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAAGATCCGCTCGCAACTG‑GGCCTGTTCGCCAAC‑CTGCGCCCGGCCATCCTTTACCCGCAACTGGCCGATGCCTCGTCGCTGAAGCCGGAAATCGTGTC > NC_002947/2253824‑2254101 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |