Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,480,611 | T→G | 11.4% | C146G (TGC→GGC) | PP_3972 → | short‑chain dehydrogenase/reductase family oxidoreductase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,480,611 | 0 | T | G | 11.4% | 83.1 / 5.5 | 35 | C146G (TGC→GGC) | PP_3972 | short‑chain dehydrogenase/reductase family oxidoreductase |
Reads supporting (aligned to +/- strand): ref base T (16/15); new base G (2/2); total (18/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.78e-01 |
TCATCGCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGGCCGC > NC_002947/4480462‑4480753 | tcatcGCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCt < 1:603701/150‑1 (MQ=255) cccTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGTCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGgc < 1:105601/150‑1 (MQ=255) cccTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCAGGCAAAGgc < 1:585146/150‑1 (MQ=255) acacCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCtt < 2:673973/150‑1 (MQ=255) cacCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATACAGGGGCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGTCTGCAAAGGCGCACGCTTg < 1:520604/150‑1 (MQ=255) gcGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAg > 1:774598/1‑150 (MQ=255) gCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGccc < 2:484189/150‑1 (MQ=255) agagCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGa < 1:193503/150‑1 (MQ=255) gagCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGAt > 1:802376/1‑150 (MQ=255) cAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATgaa < 1:742486/150‑1 (MQ=255) agaagCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATgaaga > 2:841859/1‑150 (MQ=255) gaagCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATgaagaa > 2:511053/1‑150 (MQ=255) gaagCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGC‑ATAGACCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATgaagaa < 2:422025/150‑1 (MQ=255) cccGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGAt < 1:841859/150‑1 (MQ=255) gggtCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCAc < 1:511053/150‑1 (MQ=255) gggCCGCCGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTgg > 1:693083/1‑97 (MQ=255) gggCCGCCGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTgg < 2:693083/97‑1 (MQ=255) cGCTGCGTGGCAAAGTTGTGGCAATCACAGGGGCATCCACCGGT‑ATAGGCCGCGCCGCCGCGCACGCTTTTGCCTGCCAAGGCGCCCGCTTGGTTCTCGCAGCCCGAGATGAAGAAGCATTATTCCACGTTCTGGCTGAATGCACCCAtt > 2:274463/1‑150 (MQ=255) tGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTAc > 1:511553/1‑150 (MQ=255) aaaGTTGTGGTAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGAt > 1:286123/1‑150 (MQ=255) aGTTGTGGTAATCACAGGGGCATCCACCCGCAACCTGACGC‑CGGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGc < 2:496279/150‑1 (MQ=255) gTAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCg > 2:459703/1‑150 (MQ=255) tAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGTAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCgg > 2:381009/1‑150 (MQ=255) tAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCgg > 2:549504/1‑150 (MQ=255) tAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCAGGCAAAGgc > 1:285904/1‑62 (MQ=38) tAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCAGGCAAAGgc < 2:285904/62‑1 (MQ=38) aTCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGAACGGATGGTGGTACGGATGCCATAGCGGTc > 1:217600/1‑150 (MQ=255) aTCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCgg < 1:145756/137‑1 (MQ=255) aTCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCgg > 2:145756/1‑137 (MQ=255) cAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCAc > 1:78151/1‑50 (MQ=255) cAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCAc < 2:78151/50‑1 (MQ=255) cAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCa < 2:218325/82‑1 (MQ=255) cAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCa > 1:218325/1‑82 (MQ=255) cgcgccgcCGCCCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAg < 1:365509/150‑1 (MQ=255) gtgcGCCTTTGCAGGCAAAGGCGTGCGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCt > 2:532462/4‑150 (MQ=255) cTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCGGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGgccgc > 1:689090/1‑150 (MQ=255) tGCAAAGGCGCACGCTTGGTGCTCACAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTgg > 1:804007/1‑80 (MQ=255) tGCAAAGGCGCACGCTTGGTGCTCACAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTgg < 2:804007/80‑1 (MQ=255) | TCATCGCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGC‑ATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGGCCGC > NC_002947/4480462‑4480753 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |