Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 3,969,337 | G→T | 35.0% | E278* (GAG→TAG) | wecE → | TDP‑4‑oxo‑6‑deoxy‑D‑glucose transaminase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 3,969,337 | 0 | G | T | 35.0% | 15.5 / 7.6 | 19 | E278* (GAG→TAG) | wecE | TDP‑4‑oxo‑6‑deoxy‑D‑glucose transaminase |
Reads supporting (aligned to +/- strand): ref base G (7/5); new base T (2/5); total (9/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.50e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGAAAGAAGCGGAAATCATGGCGGTGTTTCATTACATTCCGCTGCACG > CP009273/3969198‑3969481 | cAGCTATTTGATGTCCGATCTGCAAGCGGCGTACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCCGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGt < 2:14791/149‑1 (MQ=255) ttGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTcc < 2:173771/149‑1 (MQ=255) tCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCTAGCTGCCGTCGATTCCCGATgg > 1:329554/1‑149 (MQ=255) tCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCATTTGCCGAAGTCGTGACGTATCTAGCTGCCGTCGATTCCCGATgg > 1:137815/1‑149 (MQ=255) cGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGCTCATATCAACCCGCACCGTCTGGGGCTGAGCCAAAACCACTAAGACGCGGGGGCGCGTCTGGGGAAAGCCAGGCGTGGCGAGCGGCCGGCGATTCCCGATGGCt > 2:75175/1‑149 (MQ=255) cAAGCTGCATACCCGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCa < 2:183680/149‑1 (MQ=255) ggggGCAAGTGGACGCAGCGGATCGTATCAACCAGGAAGGTCTGGCGCGGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGGATATGTTc < 2:57381/149‑1 (MQ=255) gggCGCAGTGGGAAGCAGCGCATCGTATCAACCCGCAACCTCGGGCGCTGTGGCAAGACTACTACGATGTGTTAGCGCCGCTGGCGAAAGGCGGGCGGATCTAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTc < 2:243595/149‑1 (MQ=255) gcgcAACTGGAACCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCTAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTa < 2:309973/149‑1 (MQ=255) tGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTaaa > 1:176280/1‑149 (MQ=255) agcGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCCAATGTTCTACATTAAACTCCggg > 2:158240/1‑149 (MQ=255) agcGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGAGCATATGTTCTACATTAAACTGCggg > 2:117690/1‑149 (MQ=255) aGCAACGTCGGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCTGTCGTTTTTCGCCATCTTCGTGGTGCACACACATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAgc < 1:75175/149‑1 (MQ=255) gCAACGTCTGGCGCTGTGGCAAAACTACTCCGATGCGTTAGCGCCTCTGGCGAAAGCCCGGCGTATCTAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGGACGCGCATAGGTTCTACATTAAACTGCGGGATATTGATGACCGGAgcg < 2:95308/149‑1 (MQ=255) gCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGa > 2:139929/1‑149 (MQ=255) gCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGa > 1:220222/1‑149 (MQ=255) aaaaCTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCTAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGaaa < 1:324160/149‑1 (MQ=255) tctGGCGAAAGCCGGGCGTATCTAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTCCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGAAAGAAGCGGAAATCATGGCGGTGTTTc < 2:137815/149‑1 (MQ=255) tATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGAAAGAAGCGGAAATCATGGCGGTGTTTCATTACATTCCGCTGCACg > 1:320200/1‑149 (MQ=255) | CAGCTATTTGATGTCCGATCTGCAAGCTGCATACCTGTGGGCGCAACTGGAAGCAGCGGATCGTATCAACCAGCAACGTCTGGCGCTGTGGCAAAACTACTACGATGCGTTAGCGCCTCTGGCGAAAGCCGGGCGTATCGAGCTGCCGTCGATTCCCGATGGCTGCGTGCAGAACGCGCATATGTTCTACATTAAACTGCGGGATATTGATGACCGGAGCGCGTTGATTAACTTTCTGAAAGAAGCGGAAATCATGGCGGTGTTTCATTACATTCCGCTGCACG > CP009273/3969198‑3969481 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |