Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,668,079 | Δ1 bp | 17.2% | intergenic (+97/‑70) | 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,079 | 0 | A | . | 17.2% | 60.9 / 12.7 | 29 | intergenic (+97/‑70) | prpF/prpD | aconitate isomerase/2‑methylcitrate dehydratase |
Reads supporting (aligned to +/- strand): ref base A (16/8); new base . (3/2); total (19/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-01 |
TCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCGGC > NC_002947/2667934‑2668226 | tCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGATATcc > 5:341566/1‑149 (MQ=255) cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGatat < 8:181381/140‑1 (MQ=255) cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCAC‑TGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTA‑ATCACAGGGTTTGAATGCTGTGatat > 7:181381/1‑140 (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: ‑ |