Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 370,911 | C→T | 6.9% | P131L (CCT→CTT) | dgcA → | dimethylglycine dehydrogenase subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 370,911 | 0 | C | T | 6.9% | 168.4 / 2.2 | 58 | P131L (CCT→CTT) | dgcA | dimethylglycine dehydrogenase subunit |
Reads supporting (aligned to +/- strand): ref base C (26/28); new base T (2/2); total (28/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01 |
CATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGAC > NC_002947/370768‑371053 | catcCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCCACGCTGATGTCGCCATCGGGCATTCGCGAACCtgtg > 8:61130/1‑148 (MQ=255) cGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCCGGCATTCGCGAACCTGTGCACcgc > 4:50700/1‑149 (MQ=255) gCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACcgcg > 8:135411/1‑149 (MQ=255) ccAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCt > 8:90130/1‑149 (MQ=255) cAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTg > 6:317782/1‑149 (MQ=255) aCCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCaaa < 3:209308/149‑1 (MQ=255) ctggCCGACGCCCTGCACAAGCTTGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAAc < 5:275933/146‑1 (MQ=255) ccGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGATCCGGTGCACCGCGCCACCTGCAAAACCATc > 6:235139/1‑149 (MQ=255) ccGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATc < 8:121499/149‑1 (MQ=255) cgccGCCATGCCCAAGCATGGCGCCAAGATCATGATCCAGCTTACCCACACCGCCCGTCCCGCGCGCTGGGACGGCTTCAACTGGCCGACGCCGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCg < 8:12956/146‑1 (MQ=255) gACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCa < 3:225356/133‑1 (MQ=255) gACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCa > 4:225356/1‑133 (MQ=255) gACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGa < 3:7898/149‑1 (MQ=255) aCGCCATGCACAAGCATGGCGCCAAGATCATGGTCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAg < 7:260119/149‑1 (MQ=255) gCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACTTGTGCACCGCGCCACCTGCAAAACCATCGAgg > 4:255088/1‑148 (MQ=255) gCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACTTGTGCACCGCGCCACCTGCAAAACCATCGAgg > 4:250772/1‑148 (MQ=255) gCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACTTGTGCACCGCGCCACCTGCAAAACCATCGAgg < 3:255088/148‑1 (MQ=255) gCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACTTGTGCACCGCGCCACCTGCAAAACCATCGAgg < 3:250772/148‑1 (MQ=255) gCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGt > 5:159495/1‑149 (MQ=255) cacaAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAaga < 7:61130/149‑1 (MQ=255) gCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCt > 2:83063/1‑149 (MQ=255) cATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTg > 3:260985/1‑149 (MQ=255) gcgcCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCgg < 6:57069/149‑1 (MQ=255) cgcCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCggg < 5:317782/149‑1 (MQ=255) aTCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCgg < 6:238058/149‑1 (MQ=255) tCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTAc < 4:136269/149‑1 (MQ=255) ccACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACgcgc < 2:42094/143‑1 (MQ=255) ccACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACgcgc > 1:42094/1‑143 (MQ=255) ccACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAAACATAGAGGTGGAAGAGATCTGGCGGCTAATCGGCAACTAAGCGCAGCCTg > 4:42281/1‑149 (MQ=255) acaTGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCa > 2:278185/1‑99 (MQ=255) ggCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTc < 7:90130/148‑1 (MQ=255) cgctcgcGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGACGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCaa < 3:9425/149‑1 (MQ=255) gctcgcGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCa > 5:198631/1‑147 (MQ=255) tcgcGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAACACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCGaaa > 4:277532/1‑147 (MQ=255) gcgcTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAgg > 6:265698/1‑149 (MQ=255) cTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGggcg < 7:135411/149‑1 (MQ=255) aCTGGCCGACGCTTATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTg < 3:42281/149‑1 (MQ=255) aCTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTg < 4:260985/149‑1 (MQ=255) aCTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTg > 3:78102/1‑149 (MQ=255) cTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTg > 5:272737/1‑148 (MQ=255) gCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAAt > 8:35170/1‑149 (MQ=255) ccGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTATGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAAtt < 1:83063/149‑1 (MQ=255) cGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAg > 3:272830/1‑104 (MQ=255) cGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAg < 4:272830/104‑1 (MQ=255) cGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTg > 3:228517/1‑149 (MQ=255) gACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCGGGACGGCGTGGAATTGt < 7:303745/149‑1 (MQ=255) gCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCgg < 7:35170/149‑1 (MQ=255) tGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTg < 3:44895/149‑1 (MQ=255) ccATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCa < 6:198631/149‑1 (MQ=255) aTCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGc < 2:299102/149‑1 (MQ=255) tCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGca < 4:228517/149‑1 (MQ=255) gggCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGcaca > 1:59611/1‑149 (MQ=255) tCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCg > 4:132279/1‑149 (MQ=255) gcgAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTg < 1:63449/116‑1 (MQ=255) gcgAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTg > 2:63449/1‑116 (MQ=255) gcgAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCa > 2:189273/1‑128 (MQ=255) gcgAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCa < 1:189273/128‑1 (MQ=255) gcgAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGAc > 1:308603/1‑149 (MQ=255) | CATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGAC > NC_002947/370768‑371053 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |