Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,668,066 | T→C | 13.9% | intergenic (+84/‑83) | prpF → / → prpD | aconitate isomerase/2‑methylcitrate dehydratase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,668,066 | 0 | T | C | 13.9% | 85.4 / 6.3 | 36 | intergenic (+84/‑83) | prpF/prpD | aconitate isomerase/2‑methylcitrate dehydratase |
Reads supporting (aligned to +/- strand): ref base T (17/14); new base C (3/2); total (20/16) | |||||||||||
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.45e-01 |
TCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTA > NC_002947/2667918‑2668199 | tCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGttt > 3:69783/1‑149 (MQ=255) aTGAGCCGTAGTGCTTGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGa < 5:383413/149‑1 (MQ=255) gCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGttt > 1:232604/1‑143 (MQ=255) gCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGttt < 2:232604/143‑1 (MQ=255) cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCtg < 6:381921/149‑1 (MQ=255) tAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCtgtg < 4:69783/149‑1 (MQ=255) tCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATcc > 5:341566/1‑149 (MQ=255) cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGatat > 7:181381/1‑140 (MQ=255) cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGatat < 8:181381/140‑1 (MQ=255) tGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGAT‑TACCTGTgg > 7:355327/1‑147 (MQ=18) gATGGCGGGCTGGGTTCGGGTGCCTGGCTATAGCTTCTAAACCACCAC‑TGGCGCTCGGCTTCAAGTGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCtgt < 2:190065/133‑1 (MQ=38) gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCtgt > 7:173407/1‑133 (MQ=255) gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCtgt < 8:173407/133‑1 (MQ=255) gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCtgt > 1:190065/1‑133 (MQ=255) gAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGAT‑TACCTGTGGGAGCgg < 8:355327/148‑1 (MQ=18) gcgTGCCTGGCGATAGCTTCTAAAGCACCACGTGGCGCGCGGCTTCAAGAGGGGCCGCTTTGCGCCCCTTCCCGGGCTTGCCCCCTCCCACAGGTA‑ATAACAGGGTTTGAATGGTGTTATATCCCGGTGGGCGCGGGCAAGCCcgcg > 1:39018/1‑147 (MQ=255) gCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACCCGAA‑ATAACAGGGGTTGTATGGGGTGGTATGCCGTGGGGAGCGGGCAAGGCCGCGCAGGGCt > 8:291950/1‑149 (MQ=37) tGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGAT‑TACCTGTgg > 3:64917/1‑122 (MQ=14) tGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGAT‑TACCTGTgg < 4:64917/122‑1 (MQ=14) gCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGca < 6:183962/147‑1 (MQ=255) ttCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATGACAGGGTTTTTATGCTGTGATATGCCTGTGGGTGCGGGCAAGCCCGCGCAGGGCTGGACAGCCGcccc > 3:117795/1‑149 (MQ=37) aaGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGAGGCCCTTCGCGGGCATGCCCCCTCCCACAGGTA‑ATCACACGGTTTGAATGGTGTGATAACCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGGACAGCAGCCCCGGCaa > 6:164039/1‑149 (MQ=38) ggCGGTCGGCTTCAAGTGGGGCCGCTTTGCGCCCCTTGGCGGGGTTCCCCGCTCCCACAGTTA‑ATCCCCGGTTTTTAAGGCTGTTATATCCTTGTGGGAGGGGGAAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAggg < 6:341566/149‑1 (MQ=255) gcgcgcGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGCAA‑ATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAgg < 7:291950/147‑1 (MQ=255) tCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTTCCCGCTCCCACAGGTA‑ATAACAGGGGTTGAATGCTGTGATATACCTGTGGGAGCGGGCAACCCCGCGAAGGGCTGCACAGCAGCCCCCGCAATTGTGAAGGGAGcag > 8:8035/1‑147 (MQ=37) cTTCAAGAGGGGACGCTTTGCGGCCCTTCGCGGGCTTTCCCGCTCCCACACGTA‑ATCAAACGGTTTTAAAGCTGTGATATACCTGAGGGAGCGGGCAACCCCCCCAAGGGCTGGCCAGCAGCCCCCGCAATTAAGAAGTGAGCAGGat > 6:88947/1‑148 (MQ=255) cAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGGTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGCAGGGCTGGACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGatat > 3:61682/1‑147 (MQ=255) agGGGCCGCTTGGGGGCCGTGCGCGGGGTGGCCGGCTCCCACTGGTA‑ATCACTGGGTTTGAATGCTGTTATATCCGTGGGGGAGCGGGAAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAgcgc < 7:111457/149‑1 (MQ=255) tggCGGCCCTTCGCGGGCTTCCCCGCTCCCCCCGTTA‑ATCACCGGTTTTTAATGCTGTCATATCCCGGGGGGTGGGGGAAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGa < 4:61682/146‑1 (MQ=255) ggCCCTTCGCGGGCTTGCCCCCTACCACAGGTA‑ATAACAAGGTTTTAATTGTGTGATATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCAGCCCAGCAGCCCCGGCCAATCTGAAGGGAGCTGGATATTAGCGCCAACGTAGACCtg > 2:315859/1‑147 (MQ=25) cTTCGCGGGCTTGCCCCCTCCCAAAGGTA‑ATCACAGGGTTTGAATGCTGTGATATACCAGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTAAGAAAGGAGCTGGAAATGAGCGCCAACGTAGACCTAAACg > 2:111337/1‑148 (MQ=37) ccccGCTCCCACCGTTA‑ATCCCCGGGTTTCAATGCTGTTATAACCCTGGGGGTGCGGGAAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAtt < 2:39018/147‑1 (MQ=255) cccGCTCCCACAGGTA‑ATCACAGGGTTTGATTGCTGTCATATCCCTGTGGGAGCGGGAAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc < 7:8035/149‑1 (MQ=255) ccGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc > 4:5565/1‑148 (MQ=255) ccGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc > 8:339574/1‑148 (MQ=255) tCCCACAGGGATATCACAGCATTCAAACCCTGTGAT‑TACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAAAGCGGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCg > 5:150592/1‑149 (MQ=17) gTA‑ACCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTa > 1:212624/1‑148 (MQ=255) acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGa > 7:286749/1‑90 (MQ=255) acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGa < 8:286749/90‑1 (MQ=255) | TCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTA > NC_002947/2667918‑2668199 |
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: ‑ |