Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,711,531 | T→G | 7.1% | Q4P (CAG→CCG) | gpsA ← | glycerol‑3‑phosphate dehydrogenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,711,531 | 0 | T | G | 7.1% | 168.9 / 4.5 | 56 | Q4P (CAG→CCG) | gpsA | glycerol‑3‑phosphate dehydrogenase |
Reads supporting (aligned to +/- strand): ref base T (22/30); new base G (2/2); total (24/32) | |||||||||||
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.71e-01 |
AGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTCGAGATCGAGCACCGCCACGC > NC_002947/4711384‑4711676 | agaTAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACGGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTg < 8:96864/149‑1 (MQ=255) gaTAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGt < 2:301369/149‑1 (MQ=255) aTAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGtt < 1:250292/149‑1 (MQ=255) cgcgGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACAGCAACCTGTTGCGGTTctgg < 3:241696/149‑4 (MQ=255) cgGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCa > 5:338891/1‑149 (MQ=255) ggATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCa > 7:165103/1‑148 (MQ=255) ggATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACAGCAACCTGTTGCGGTTctgggcg < 6:393854/149‑7 (MQ=255) tcGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATatccaa > 5:214055/1‑149 (MQ=255) gcgATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATatccaatc < 1:302137/149‑1 (MQ=255) tCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCAtt > 5:349742/1‑148 (MQ=255) cacaCGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTa < 7:24405/148‑1 (MQ=255) acacGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTaaa < 1:326790/149‑1 (MQ=255) cacGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTaaa < 1:260991/148‑1 (MQ=255) cacGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTa > 2:260991/1‑146 (MQ=255) acGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATatccaat > 5:202361/1‑139 (MQ=255) acGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATatccaat < 6:202361/139‑1 (MQ=255) acGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGc > 1:238677/1‑149 (MQ=255) gCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGcc > 7:360339/1‑148 (MQ=255) aTGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCaa > 7:21867/1‑148 (MQ=255) tGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCaaaa < 5:381120/149‑1 (MQ=255) gCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAg < 6:134623/148‑1 (MQ=255) ccTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGtt < 6:164589/149‑1 (MQ=255) ccTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGtt < 5:368596/149‑1 (MQ=255) cTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGtt > 4:23352/1‑148 (MQ=255) ggTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGc < 6:349742/149‑1 (MQ=255) tCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCa < 8:114467/149‑1 (MQ=255) gcgcATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCAt < 1:38690/148‑1 (MQ=255) gcgcATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCAt > 3:285745/1‑148 (MQ=255) aTTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCt > 7:362992/1‑149 (MQ=255) gcgcACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGcc > 2:139526/1‑149 (MQ=255) cGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGCCCCGGCCATTATACGGCTCGCCTGCTa > 7:186501/1‑148 (MQ=255) cGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTAc < 3:23352/149‑1 (MQ=255) cacaCCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGa < 8:165103/149‑1 (MQ=255) acCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAAcc < 6:43190/149‑1 (MQ=255) ttCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAgtgt < 5:240472/148‑1 (MQ=255) cGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTTTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTg < 8:360339/148‑1 (MQ=255) cGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTg < 3:84743/148‑1 (MQ=255) ccAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTaaa > 4:321979/1‑85 (MQ=255) agcagGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAgt > 2:252029/1‑139 (MQ=255) agcagGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAgt < 1:252029/139‑1 (MQ=255) agcagGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTgg < 3:251362/145‑1 (MQ=255) agcagGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTgg > 4:251362/1‑145 (MQ=255) ggTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGc < 6:338891/149‑1 (MQ=255) tttGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCgg < 8:21867/149‑1 (MQ=255) ggcggTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGcc < 4:285745/149‑1 (MQ=255) gCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGCCCCGGCCTTTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCa > 3:279159/1‑149 (MQ=255) cgccCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTc < 7:9368/149‑1 (MQ=255) cgccCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTc < 1:139526/149‑1 (MQ=255) gactgAACAGCAACCTGTTGCGGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTc > 5:219096/5‑147 (MQ=255) cccAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTc > 3:111350/1‑147 (MQ=255) cAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAg > 1:69565/1‑140 (MQ=255) cAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAg < 2:69565/140‑1 (MQ=255) tgAACAGCAACCTGTTGCGGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTc > 1:50663/2‑144 (MQ=255) ccGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTCGAGATCGAg > 4:369956/1‑149 (MQ=255) cGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGACATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTCGAGATagag > 3:374837/1‑148 (MQ=255) gCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTCGAGATCGAGCACCGCCAcgc < 7:124288/148‑1 (MQ=255) | AGATAGCGCGGATTCTCGCGATTCACACGCATGGCCTCGGCCTGCGCCGGGTCGCGCATCCATTGGCGCACCGGCACACCGTTTTCCGCCAGCAGGTTTGCCACGGCGGTGCCGAAGCTGCCGCCGCCCAGAACCGCAACAGGTTGCTGTTCAGTCATATCCAATCCATTAAAGCCAAAAGTTAAGTGGCGACCCGGCCATTATACGGCTCGCCTGCTACAGAACCAGTGTCTGGACGTTTGCGGTACTTGCCATGGCATTTTCGGAGCCGTCGAGATCGAGCACCGCCACGC > NC_002947/4711384‑4711676 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |