Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,087,768 | G→C | 7.4% | R259R (CGC→CGG) ‡ | metR‑II ← | homocysteine‑binding DNA‑binding transcriptional regulator |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,087,768 | 0 | G | C | 7.4% | 168.6 / 3.6 | 54 | R259R (CGC→CGG) ‡ | metR‑II | homocysteine‑binding DNA‑binding transcriptional regulator |
Reads supporting (aligned to +/- strand): ref base G (24/26); new base C (2/2); total (26/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01 |
AACCGGGCCTGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTCGGCGGG > NC_002947/3087621‑3087914 | aaCCGGGCCTGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGgc > 7:283173/1‑149 (MQ=255) cGGGCCTGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCt < 4:43618/149‑1 (MQ=255) gCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCtt < 6:29623/140‑1 (MQ=255) gCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCtt > 5:29623/1‑140 (MQ=255) cACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACa > 8:33857/1‑148 (MQ=255) cACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACa > 7:254922/1‑148 (MQ=255) aCGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATa < 5:157381/149‑1 (MQ=255) tAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCcg > 3:136522/1‑149 (MQ=255) acCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCcgcgc > 2:351634/1‑148 (MQ=255) acCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGCGCCTTGGCgagc > 1:22620/1‑133 (MQ=255) acCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGCGCCTTGGCgagc < 2:22620/136‑4 (MQ=255) cAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCACCAGCATGTCACTGCGCAGTGCAGCAAAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAAcgg > 2:320620/1‑146 (MQ=255) cGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTAtt < 5:358193/149‑1 (MQ=255) gtgAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGt < 6:300052/148‑1 (MQ=255) cGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACgg > 3:148795/1‑149 (MQ=255) gTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGt > 2:270981/1‑149 (MQ=255) cTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTgg < 8:254922/149‑1 (MQ=255) gCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGgcagca < 4:136522/148‑1 (MQ=255) gCAGAAAATCACGCATGTACGGCACCTCCAGCATGCCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGgcaggcag > 2:106626/1‑146 (MQ=255) cAGAAAATCACGCATGTGCGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGa > 5:281279/1‑149 (MQ=255) cAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGa > 6:199943/1‑149 (MQ=255) aaaTCACGCATGTACGGCACCTCCAGCATGCCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGAcgcc < 1:106626/149‑1 (MQ=255) aaTCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCCGCAGAcgccg > 7:330932/1‑149 (MQ=255) aTGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACt > 8:101794/1‑101 (MQ=255) aTGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCgg > 6:364300/1‑148 (MQ=255) aTGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTGGGCGGTTACATAGCCGCGCGAACt < 7:101794/101‑1 (MQ=255) tGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGaa < 5:199943/149‑1 (MQ=255) tGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGaa < 1:351634/149‑1 (MQ=255) tGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGaa < 5:364300/149‑1 (MQ=255) gTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGaa < 8:13055/148‑1 (MQ=255) ccTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGt < 8:283173/149‑1 (MQ=255) tCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGcatcatcat < 4:148795/149‑1 (MQ=255) gcAGTGCAGCATAGAGGGTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTc < 3:299087/149‑1 (MQ=255) gTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGt > 7:33861/1‑149 (MQ=255) tGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCTGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGtt < 6:281279/149‑1 (MQ=255) tAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAgt < 1:270981/149‑1 (MQ=255) ngagGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAgt > 4:315161/2‑148 (MQ=255) gTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAATCAACCa > 3:51/1‑113 (MQ=255) gTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCa < 4:51/113‑1 (MQ=255) tttGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAgtacgtac < 7:33857/149‑1 (MQ=255) gCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCACGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCt < 6:180165/148‑1 (MQ=255) tGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTg > 6:33961/1‑133 (MQ=255) tGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTg < 5:33961/133‑1 (MQ=255) tGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGc > 8:284760/1‑149 (MQ=255) gCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCg > 4:7713/1‑149 (MQ=255) cAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGggcagg > 6:41496/1‑76 (MQ=255) gCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGcatc < 7:276841/109‑1 (MQ=255) gCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGcatc > 8:276841/1‑109 (MQ=255) ccGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTc < 3:7713/149‑1 (MQ=255) ccGCCAAGCGCCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGAc < 7:50103/145‑1 (MQ=255) ccGCCAAGCGCCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGAc > 8:50103/2‑146 (MQ=255) gCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTcgg > 2:166346/1‑149 (MQ=255) cAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCt < 3:220763/116‑1 (MQ=255) cAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCt > 4:220763/1‑116 (MQ=255) aGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTCGGCggg < 3:315161/149‑1 (MQ=255) | AACCGGGCCTGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTCGGCGGG > NC_002947/3087621‑3087914 |
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: ‑ |