Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 85,077 | T→G | 7.6% | intergenic (+55/+75) | aroE → / ← PP_0075 | shikimate 5‑dehydrogenase/choline sulfate transporter |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 85,077 | 0 | T | G | 7.6% | 154.1 / 3.7 | 53 | intergenic (+55/+75) | aroE/PP_0075 | shikimate 5‑dehydrogenase/choline sulfate transporter |
Reads supporting (aligned to +/- strand): ref base T (25/24); new base G (2/2); total (27/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01 |
AAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTG > NC_002947/84940‑85224 | aaGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGc < 2:285237/149‑1 (MQ=255) gCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCg < 6:343744/149‑1 (MQ=255) tCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCtgt < 8:84725/148‑1 (MQ=255) tCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCtgt < 8:79743/148‑1 (MQ=255) gcgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCtgtg < 1:360221/145‑1 (MQ=255) gcgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAg > 6:285450/1‑149 (MQ=255) gtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCC‑ATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACG‑TACTTTGTGGGAACgg < 1:238685/149‑1 (MQ=18) gCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTc > 7:272881/1‑149 (MQ=255) ggATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTg < 6:142578/142‑1 (MQ=255) ggATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTg > 5:142578/1‑142 (MQ=255) ggATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTcgc > 8:275297/1‑148 (MQ=255) ggTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCaaa < 5:285450/149‑1 (MQ=255) ggTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCa > 6:331222/1‑147 (MQ=255) cTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTAGGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAc‑‑cc < 6:176582/149‑1 (MQ=255) aCTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCTTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGaa < 7:1164/149‑1 (MQ=255) tGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATc > 4:40118/1‑149 (MQ=255) gcgcCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGGGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCa > 6:162708/1‑137 (MQ=255) gcgcCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGAGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCa < 5:162708/137‑1 (MQ=255) cgcCGGCAACTGTCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑GGTCGTACGCCGATCCTGTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATc < 4:222366/147‑1 (MQ=255) cgcCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATcc > 5:66270/1‑148 (MQ=255) aCTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCAt < 5:269369/149‑1 (MQ=255) cgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCt < 8:256806/143‑1 (MQ=255) cgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCt > 7:256806/1‑143 (MQ=255) ggggCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCaaa < 6:66270/149‑1 (MQ=255) gCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCggg < 3:40118/147‑1 (MQ=255) ccGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGca > 7:146094/1‑148 (MQ=255) tagggCCGCCAAGCGGCCCTATCGCGACGCAACGCCGTTACCACAAAGTA‑CGTCGTACGCCGCTCCACTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAACCAGAC‑‑CCAGAATCCCCACCTCAATCCTCAACCCCGATCGGGcaca > 8:139354/3‑149 (MQ=37) cGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGc > 8:212209/1‑85 (MQ=255) gggCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAAc < 8:272881/149‑1 (MQ=255) ccGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTcc < 8:146094/149‑1 (MQ=255) cTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTcc > 7:243762/1‑122 (MQ=255) cTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTcc < 8:243762/122‑1 (MQ=255) gcccaTCGCGACGCAAGGCCGCGCCCACAGAGGATAGGCGTAAGCCG‑TACGTTGTGGGAACGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGCATCCCAACCGAATTCCTCAAACCGGATCGGGCACCACTCCACCCCGTcccg > 4:88614/5‑147 (MQ=9) ccTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAAc < 7:214385/119‑1 (MQ=255) ccTATCGCGACGCAAGGCCGTTACCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAAc > 8:214385/1‑119 (MQ=255) cTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCt > 3:260588/1‑149 (MQ=255) tCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCt > 7:204603/1‑149 (MQ=255) gcgACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTg > 4:124887/1‑148 (MQ=255) cgACGCAAGGCCGTACCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGTGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCACAATCCACAAACCCGATCGGGCACAACTCCACCCCTTCCAGCGTCTGc > 8:176599/1‑148 (MQ=39) cgctCCCACAGAGGATCGGCGTACGACG‑TACTTTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCAc < 6:143270/118‑1 (MQ=18) cgctCCCACAGAGGATCGGCGTACGACG‑TACTTTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCAc > 5:143270/4‑121 (MQ=18) gTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGc < 3:185698/136‑1 (MQ=255) gTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGc > 4:185698/1‑136 (MQ=255) ttCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCAc > 3:293917/1‑149 (MQ=255) cccACAAAGTA‑CGTCGTACGCCGATCCTCTGGGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCAcc < 7:139354/148‑1 (MQ=255) ccACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGa < 3:106905/142‑1 (MQ=255) ccACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGa > 4:106905/1‑142 (MQ=255) ccACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTg > 5:326430/1‑149 (MQ=255) cacaAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGAGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGc > 1:44496/1‑149 (MQ=255) acaAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCg < 8:204603/149‑1 (MQ=255) acaAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCg > 6:29826/1‑149 (MQ=255) gTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACTACCAGAATCCCAACCTCATTCCTCAAACCGGAACGGGCACACCTCCACCCATGCCACGTGGTGCACCACCACGATCACATGCGGca > 1:103305/1‑148 (MQ=25) gTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTg < 4:260588/149‑1 (MQ=255) | AAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTA‑CGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAC‑‑CCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTG > NC_002947/84940‑85224 |
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: ‑ |