Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,655,626:1 | +A | 29.0% | intergenic (‑81/+71) | epd ← / ← tktA | D‑erythrose 4‑phosphate dehydrogenase/transketolase A |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,655,626 | 1 | . | A | 29.0% | 65.1 / 38.0 | 31 | intergenic (‑81/+71) | epd/tktA | D‑erythrose 4‑phosphate dehydrogenase/transketolase A |
Reads supporting (aligned to +/- strand): ref base . (8/14); new base A (5/4); total (13/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.33e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01 |
CGCGCAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCG > NC_002947/5655479‑5655770 | cgcgcAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CGTTCGCGGGCAAGCCCGCTCCCACAg‑g < 2:272998/149‑1 (MQ=255) cgcgcAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAg‑g < 5:382783/149‑1 (MQ=255) cgcgcAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGGGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAg‑g < 7:198481/149‑1 (MQ=255) cAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCGCCCCCCg‑ggggt > 4:225679/1‑146 (MQ=255) aGGACAGAGCTGCCGATAGGGACATAACCGTTGAGTGCTAATTTGTAGGGACGGTGGGGCATTCGGGGCTCAGGAGTCAGTTCAATTGTCGTCAGCCTGCGGAATTTTCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTc < 8:384529/149‑1 (MQ=255) gACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCAt < 8:3446/149‑1 (MQ=255) cGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTCGCCTTCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATg > 6:237384/1‑148 (MQ=37) gaaaCGGCCATAAGCGTTGAGGCCAACTTTGTATGGACGGTGGGGCATTCGTGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCTACCTCtgtg < 8:175445/146‑1 (MQ=255) aTACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGTCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTgg < 6:201673/149‑1 (MQ=255) gCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTCGCCTTCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACC‑CTGTggg < 5:237384/149‑1 (MQ=21) ccATACCCGTGGCCTGCAACTTGGTCGGGACGGTGGGGCATTCGGGGCTCACGCGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCgg < 4:360793/149‑1 (MQ=38) ccATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGt > 2:78072/1‑125 (MQ=255) ccATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGt < 1:78072/125‑1 (MQ=255) ccATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCgg < 8:380716/149‑1 (MQ=255) cATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCggg < 3:166342/149‑1 (MQ=255) attaCCGTTGAGTGCAACTTTGTAGGGACTGTGGGTCATTTTTGTCTCACGAGACAGTTCAATTGTAGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTAGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCtgt < 4:301999/136‑1 (MQ=255) aTAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCtgt > 3:301999/1‑139 (MQ=255) aTAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGc < 7:56443/149‑1 (MQ=255) aCCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTATCATGCAACCTCTGTGGGAGCGGGCtt < 1:12280/148‑1 (MQ=255) cAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGccc < 1:342270/141‑1 (MQ=255) cAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGccc > 2:342270/1‑141 (MQ=255) aaCTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACctct > 5:46603/1‑122 (MQ=255) aaCTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACctct < 6:46603/122‑1 (MQ=255) gTAGGGACGGTGGGGCATTCGGGCCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGGGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGa < 3:225679/148‑1 (MQ=255) aGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGc < 4:361512/122‑1 (MQ=25) aGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGc > 3:361512/1‑122 (MQ=25) ggACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAg > 6:47710/1‑95 (MQ=255) gCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCg > 5:54545/1‑148 (MQ=255) cccGCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAgg > 5:216240/1‑129 (MQ=25) cccGCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAgg < 6:216240/129‑1 (MQ=25) cccGCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGaa < 2:37281/134‑1 (MQ=37) cccGCTCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGaa > 1:37281/1‑134 (MQ=37) tCCCACAGAGGTTGCATGAAACC‑CTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATg > 7:177619/1‑149 (MQ=38) tCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATg > 3:181991/1‑149 (MQ=255) cacaG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGcc > 6:370376/1‑148 (MQ=255) acaG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCAGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg > 4:314560/1‑148 (MQ=255) | CGCGCAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGG‑‑‑‑CATTCGCGGGCAAGCCCGCTCCCACAG‑GGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCG > NC_002947/5655479‑5655770 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |