Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 2,872,022 | G→A | 29.2% | R157C (CGC→TGC) | ispD ← | 2‑C‑methyl‑D‑erythritol 4‑phosphate cytidylyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,872,022 | 0 | G | A | 29.2% | 38.3 / 18.8 | 34 | R157C (CGC→TGC) | ispD | 2‑C‑methyl‑D‑erythritol 4‑phosphate cytidylyltransferase |
Reads supporting (aligned to +/- strand): ref base G (11/13); new base A (6/4); total (17/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.08e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01 |
CTGAGGATGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACGCCAACAATCGCGCGAGGTCATCCTGATGCAAACA > NC_000913/2871882‑2872160 | cTGAGGATGGAATCCGCAATAGTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCGAt < 1:420874/143‑1 (MQ=255) gATGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTCCGATCAACg > 1:91983/1‑143 (MQ=255) tctCCGCAATATGCCAGCGCCGAGGCTTCGCCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCGATCAACGGTATg < 1:7761/141‑1 (MQ=255) cGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCGATCAACGGTATGAGCa < 1:617410/143‑1 (MQ=255) gCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGACAGACAGTCATGTAGCAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGACGTTGCGTTCAACGGTGTGGGcga > 1:337453/1‑141 (MQ=255) aaTATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGGGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCGATCAACGGTATGAGCAATg < 2:91088/143‑1 (MQ=255) cgcgGAGGTTTGCTCGGTAATAGTCGCGCCTGCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGc < 1:644561/143‑1 (MQ=255) gcggAGGCTTCGTCGGTAGTAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCAATCAACGGTATGAGCAATGGCATTTTTGcc < 2:564289/140‑1 (MQ=255) cTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAACAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCGATCAACGGTATGAGCCATGGCATTTTTGCCCGGTTCg > 2:68804/1‑143 (MQ=255) tcgtcgGTAGTAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAGATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGc < 2:270687/143‑1 (MQ=255) tcgGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTGTTGCCCGGTTCGGAACg > 2:51571/1‑143 (MQ=255) tcgGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGAGCCAGCGCGTGACACCCAGCC‑CGGCAATCGACGGTATGAGCACTGGCATTTTTGCCCGGTTCGGCACg > 1:123014/1‑143 (MQ=255) ggTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCAATCAACGGTATGAGCAATGGCATATTTGCCCGGTAAGGAACGtt > 2:48010/1‑143 (MQ=255) aatctgcgcgACTTGATTGAGAGCGCGCGTCAGCGAGTCATGTACCAGCTCACGAGGGAAACATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACCTTTCa < 2:546210/137‑1 (MQ=255) aaTAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTCCAATCAACGGTATGAGCAATGGCATTTTTGCCGGGTTCGGCACGTTTCa > 1:621966/1‑143 (MQ=255) tCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCGATCAACGGTATGACCCATGGCATCTTTGACCGCTTCCCCACGGCTCCTAGTa > 2:256812/1‑143 (MQ=255) gcgcCTTCATTGTGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAATATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCAATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATc < 2:448504/143‑1 (MQ=255) cTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCAATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCgcgc > 1:337145/1‑143 (MQ=255) ggcAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGAt < 2:169361/141‑1 (MQ=255) cAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATcc < 2:647114/143‑1 (MQ=255) cAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATcc > 2:513347/1‑143 (MQ=255) aCAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCAATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGGGCGGCGTGGATCCCCCCCGTGCGGCTg > 1:508887/1‑143 (MQ=255) aGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGACACGTTTCCTAGTATCGCGCACTCGTGCGGCGAGGATCCCCCCCGTGCGGCTggg > 1:471998/1‑142 (MQ=255) aGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGACACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTggg > 1:40890/1‑142 (MQ=255) cccACGAGGGAAAAATTGGGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCAAGCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGttt < 1:220484/141‑1 (MQ=255) aCGAGGGAAAGAGGGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTGTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGc < 1:200259/143‑1 (MQ=255) gggAAAGATTGCGGCGTCAGCGCGTGCCA‑TCAGCCGGTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGCTCaa < 2:417642/143‑1 (MQ=255) aaaaaTTGCGGCGTCAGCGGGTGCCA‑TACGCCGTTGCAATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACGc < 2:445433/143‑1 (MQ=255) aaTTGCGGCGTCAGCGCGTGCCA‑TACGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACGCcaa < 2:183114/143‑1 (MQ=255) aTTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTCTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACGCcaac > 2:150668/1‑143 (MQ=255) gcgTCAGCGCGTGCCA‑TAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCCCTCAACGCCGACAAgcgc > 1:75598/1‑143 (MQ=255) aGCGCGTTCCA‑TAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGGAGCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACGCCAACAATCGCGCGAg < 2:356673/143‑1 (MQ=255) gTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACGCCAACAATCGCGCGAGGTCATCCTGATGCAAAc > 1:556983/1‑143 (MQ=255) ttGCAATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACGCCAACAATCGCGCGAGGTCATCCTGATGCAAACa > 1:35989/1‑143 (MQ=255) | CTGAGGATGGAATCCGCAATATTCCAGCGCCGAGGCTTCGTCGGTAATAGTCGCGCCTTCATTTAGAGCGCGCGTCAGACAGTCATGTAACAGCTCACGAGGGAAAAATTGCGGCGTCAGCGCGTGCCA‑TAAGCCGTTGCGATCAACGGTATGAGCAATGGCATTTTTGCCCGGTTCGGCACGTTTCATAGTATCGCGCACTGGTGCGGCGAGGATCCCCCCCGTGCGGCTGGTTTCGCTCAACGCCAACAATCGCGCGAGGTCATCCTGATGCAAACA > NC_000913/2871882‑2872160 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |