| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,821,953 | T→A | 15.2% | L173Q (CTG→CAG) | metW → | methionine biosynthesis protein MetW |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 5,821,953 | 0 | T | A | 15.2% | 150.3 / 13.6 | 58 | L173Q (CTG→CAG) | metW | methionine biosynthesis protein MetW |
| Reads supporting (aligned to +/- strand): ref base T (25/24); new base A (4/5); total (29/29) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TGGTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGCACCATGC > NC_002947/5821811‑5822100 | tGGTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACAGCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTc < 1:623341/150‑1 (MQ=255) gTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTCGa < 2:1017630/150‑1 (MQ=255) ccTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTCGagta < 1:1037148/150‑4 (MQ=255) cTGGCGACCAAAGGCCGCATGCCGGTATCGGACAACATTCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCAc > 1:1112406/1‑150 (MQ=255) gACNAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTc > 2:700202/1‑139 (MQ=255) gACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTc < 1:700202/139‑1 (MQ=255) gACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCa < 2:804752/150‑1 (MQ=255) gACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCa > 1:372311/1‑150 (MQ=255) cAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCg > 1:718240/1‑150 (MQ=255) gCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACgggt < 2:766696/150‑1 (MQ=255) cATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCTCTTGCACCGCAACgggtg < 2:99088/150‑1 (MQ=255) tGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACCTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACgggtggg < 2:638549/150‑1 (MQ=255) ccGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACgggtggggt < 1:470878/150‑1 (MQ=255) aTCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGCCGACCACTTGCACCGCAACGGGTGGGGTGGCCg > 2:222018/1‑150 (MQ=255) tCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCgg > 1:895773/1‑150 (MQ=255) gACTTCATGCCGTATACCGGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTa > 2:742442/1‑150 (MQ=255) ctgcATGCCGTATACCTGGTACAACACGCCGAACATCCGCTGCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATg < 2:474017/147‑1 (MQ=255) ttCATGCCGTATACCTGGTACAACACGCCCAACATCCACTGTTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATgg < 2:1112406/150‑1 (MQ=255) tCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGc < 1:441245/150‑1 (MQ=255) cATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGcc > 2:687009/1‑150 (MQ=255) aTGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACgggtg < 2:924289/131‑1 (MQ=255) aTGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACgggtg > 1:924289/1‑131 (MQ=255) tataCCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAAtctt < 2:718240/150‑1 (MQ=255) tataCCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAAtctt < 2:895773/150‑1 (MQ=255) cTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTc > 2:719657/1‑97 (MQ=255) cTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCAGGCGGTc < 1:719657/97‑1 (MQ=255) tGGTACAACACGCCGAACCTCCACTTCGGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGt < 2:255851/150‑1 (MQ=255) tGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGTCCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTa > 2:606475/1‑147 (MQ=255) tGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCg < 2:550810/93‑1 (MQ=255) tGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCg > 1:550810/1‑93 (MQ=255) tGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTa < 1:606475/147‑1 (MQ=255) ggTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTg > 2:920181/1‑150 (MQ=255) gCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACgggtgg > 2:459192/1‑109 (MQ=255) gCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACgggtgg < 1:459192/109‑1 (MQ=255) gCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTa > 1:417839/1‑136 (MQ=255) gCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTa < 2:417839/136‑1 (MQ=255) ccACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGt < 2:74286/150‑1 (MQ=255) cTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCGAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTcag < 2:171064/150‑1 (MQ=255) tCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTcagca < 1:687009/150‑1 (MQ=255) gCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGAt > 1:96314/1‑127 (MQ=255) gCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGAt < 2:96314/127‑1 (MQ=255) cACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTGTAGGTGAGATCGGCATCTACCGCGTCAGCAGTcc > 2:580018/1‑150 (MQ=255) ttCGAAGAACTTGTGCACGAACGCCGAGCTAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGGGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAagag > 1:29767/1‑150 (MQ=255) gaagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCaa < 1:947464/68‑1 (MQ=255) gaagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCaa > 2:947464/1‑68 (MQ=255) aagaaCTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATc > 1:1112055/1‑150 (MQ=255) gcTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGCCCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCt > 1:435451/2‑150 (MQ=255) cACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATc < 2:81280/139‑1 (MQ=255) cACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATc > 1:81280/1‑139 (MQ=255) cAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCgga > 1:903175/1‑150 (MQ=255) caagtggtCGACCGCCAGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCNGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCggag < 1:494682/143‑1 (MQ=255) aaatggtCGACCGAGAGGCGGTCGACCACTTACACCGTAAAGGGAGGGGTCGCCGGCTAGGGCAGAATCGTCTAGGGGAGATAAGGGTCAACCGCGTCAGCAGGGAGGACCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCggagg > 2:494682/7‑150 (MQ=255) cTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCg > 1:176473/1‑150 (MQ=255) tCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGc > 1:406677/1‑150 (MQ=255) tCGACCGCCAGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGc > 1:647121/1‑150 (MQ=255) gACCGCCAGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGATCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATc > 2:888622/1‑146 (MQ=255) gACCGCCAGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGATCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATc < 1:888622/146‑1 (MQ=255) ccTGGCGGACGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGCACCATGc > 1:98445/1‑150 (MQ=255) | TGGTACCTGGCGACCAAAGGCCGCATGCCGGTATCGGACTTCATGCCGTATACCTGGTACAACACGCCGAACATCCACTTCTGCACCTTCGCCGACTTCGAAGAACTTGTGCACGAACGCCGGGCCAAGGTACTCGACCGCCTGGCGGTCGACCACTTGCACCGCAACGGGTGGGGTGGCCGGCTATGGCCTAATCTTCTAGGTGAGATCGGCATCTACCGCGTCAGCAGTCCGGGCCTGCAAGAGCATCAGCTCGCGGTCTGACGCCCACCGGAGGAATCGCACCATGC > NC_002947/5821811‑5822100 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |