Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 941,442 | T→G | S162A (TCT→GCT) | dmsA → | dimethyl sulfoxide reductase subunit A |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 941,442 | 0 | T | G | 96.9% | 105.4 / ‑4.0 | 32 | S162A (TCT→GCT) | dmsA | dimethyl sulfoxide reductase subunit A |
Reads supporting (aligned to +/- strand): ref base T (0/1); new base G (12/19); total (12/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01 |
GAGCGCATTAGCTGGGAAGAAGCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGTCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAA > NC_000913/941364‑941509 | gAGCGCATTAGCTGGGAAGAAGCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGgct < 1:194998/81‑1 (MQ=255) gcgcATTAGCTGGGAAGAAGCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGgctat < 2:194999/81‑1 (MQ=255) aTTAGCTGGGAAGAAGCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGctatctat < 1:195000/81‑1 (MQ=255) aTTAGCTGGGAAGAAGCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGctatctat < 1:195001/81‑1 (MQ=255) aGCTGGGAAGAAGCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTatctatct < 1:195002/80‑1 (MQ=255) cTGGGAAGAAGCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGaa > 2:195044/1‑81 (MQ=255) ggAAGAAGCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACt < 2:195003/80‑1 (MQ=255) gaagCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGc < 1:195004/81‑1 (MQ=255) gaagCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGc < 2:195005/81‑1 (MQ=255) cGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGt < 2:195006/79‑1 (MQ=255) gACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACg > 1:195060/1‑81 (MQ=255) aCATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGc < 1:195007/81‑1 (MQ=255) atcGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGggc < 2:195008/81‑1 (MQ=255) tcGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGggcg < 2:195009/81‑1 (MQ=255) gcgACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGTCTATCTATCTGAACTATGGCACCGGTACGCTGggcggc < 2:195010/81‑1 (MQ=255) gACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCAc > 1:195050/1‑81 (MQ=255) ccAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCa < 2:195011/81‑1 (MQ=255) cAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCAt < 2:195013/81‑1 (MQ=255) cAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCAt < 1:195012/81‑1 (MQ=255) aaTATGCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATg < 2:195014/81‑1 (MQ=255) gCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCg > 2:195045/1‑81 (MQ=255) gCAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCg > 2:195051/1‑81 (MQ=255) cAGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCg < 1:195015/80‑1 (MQ=255) aGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCt > 1:195056/1‑81 (MQ=255) aGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCt > 1:195054/1‑81 (MQ=255) aGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCt < 2:195016/81‑1 (MQ=255) aGCGCCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCt > 2:195055/1‑81 (MQ=255) cgcCTGATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTcc < 2:195017/81‑1 (MQ=255) gATCAAAGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGgcc > 2:195049/1‑81 (MQ=255) aaGAGTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACTATGACCCGCTCCTGGCCGCCgg > 1:195047/1‑81 (MQ=255) agTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGaa < 1:195018/81‑1 (MQ=255) agTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGa > 1:195065/1‑80 (MQ=255) gTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGaaa < 1:195019/81‑1 (MQ=255) gTACGGCAACGAGGCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGaaa > 2:195064/1‑81 (MQ=255) | GAGCGCATTAGCTGGGAAGAAGCCTACGACATCATCGCGACCAATATGCAGCGCCTGATCAAAGAGTACGGCAACGAGTCTATCTATCTGAACTATGGCACCGGTACGCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAA > NC_000913/941364‑941509 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |