Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,668,081 | C→A | 26.8% | intergenic (+99/‑68) | 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,081 | 0 | C | A | 26.8% | 55.8 / 11.4 | 33 | intergenic (+99/‑68) | prpF/prpD | aconitate isomerase/2‑methylcitrate dehydratase |
Reads supporting (aligned to +/- strand): ref base C (11/11); new base A (7/2); total (20/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.37e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01 |
TCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCGGC > NC_002947/2667934‑2668226 | tCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATcc > 5:341566/1‑149 (MQ=255) tGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGAT‑TACCTGTgg > 7:355327/1‑147 (MQ=18) 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 > 8:339574/1‑148 (MQ=255) ccGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc > 4:5565/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) tgtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGccc > 7:329924/1‑149 (MQ=255) tgtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGcc > 8:38468/1‑148 (MQ=255) taatatCCCTGTGGGAGGGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCt < 5:88947/146‑1 (MQ=255) tatatCCCTGTGGGTGGGGGAAAGCCCGCGAAGGGCTGCTTCTCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTc < 1:315859/147‑1 (MQ=255) atCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCg > 6:26098/1‑146 (MQ=255) atCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCg < 5:26098/146‑1 (MQ=255) atCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCGGc > 2:308496/1‑148 (MQ=255) | TCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCGGC > NC_002947/2667934‑2668226 |
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: ‑ |