Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 2,811,134 | A→G | H122R (CAC→CGC) | mprA → | DNA‑binding transcriptional repressor MprA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,811,134 | 0 | A | G | 96.8% | 96.1 / ‑4.6 | 31 | H122R (CAC→CGC) | mprA | DNA‑binding transcriptional repressor MprA |
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (14/16); total (15/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.84e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01 |
TGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGA > NC_000913/2811006‑2811262 | tgCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAg > 2:120389/1‑133 (MQ=255) tgCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAg > 2:487761/1‑133 (MQ=255) atcatcCCGTGCCAACGCGACGCGTATTGCCGGTGAACTGGTAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTgcgc < 2:260718/133‑1 (MQ=255) atcatcCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTgcgc < 1:467243/133‑1 (MQ=255) tcatcCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTgcgc < 1:315014/132‑1 (MQ=255) tcatcCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTgcgc < 2:458163/132‑1 (MQ=255) tcatcCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTgcgc > 2:315014/1‑132 (MQ=255) tcatcCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGcgcg > 2:598072/1‑133 (MQ=255) tcatcCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTgcgc > 1:458163/1‑132 (MQ=255) catcCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGa < 1:74996/133‑1 (MQ=255) ccGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGtt > 2:18115/1‑133 (MQ=255) ccGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGtt > 2:310135/1‑133 (MQ=255) tACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTa < 1:310135/133‑1 (MQ=255) tACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTa < 1:18115/133‑1 (MQ=255) acgcgacgcgTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTaccacc > 1:406274/1‑133 (MQ=255) cgcgacgcgTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCg < 1:598072/133‑1 (MQ=255) cgcgacgcgTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCg > 1:132555/1‑133 (MQ=255) cgcgTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCgcagca > 2:84416/1‑133 (MQ=255) gCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAActgcc < 1:469155/133‑1 (MQ=255) aaaaCGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTAGCACCGCATCATAACTGGCTGCATCAACTCTgg > 2:190311/1‑133 (MQ=255) aaaCGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGt < 1:84416/133‑1 (MQ=255) cgGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCcgcg < 2:406274/133‑1 (MQ=255) aaCGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCacaac < 2:323715/133‑1 (MQ=255) gccgcTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAAt < 1:342561/133‑1 (MQ=255) cATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCACCACAACAGCAAAGGCTTAGCTCGAGAAAATGACCCGCaact > 2:279348/1‑131 (MQ=255) tGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGct < 1:487761/133‑1 (MQ=255) cAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGctct > 1:617824/1‑133 (MQ=255) cAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGctct > 1:34776/1‑133 (MQ=255) cAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTATGGTCAGCGCTCAGCACAAAAGAAAAAGATCAGCTCTAGCAAATCTCCCGGAAATTGctct > 2:206897/1‑133 (MQ=255) aaaaaGGTCGCGAGTTTTCGCGCGAGGTTTTACCACCGCAGCATAACTGCTTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGa < 1:190311/133‑1 (MQ=255) ggTCGCGAGTTTTTGTGCGAGGTTTTACTACCTCAGCATAATTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTGGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGa < 2:216382/133‑1 (MQ=255) | TGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGA > NC_000913/2811006‑2811262 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |