Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,070,763 | G→T | R65S (CGT→AGT) | rutE ← | putative malonic semialdehyde reductase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,070,763 | 0 | G | T | 96.7% | 99.9 / ‑2.8 | 30 | R65S (CGT→AGT) | rutE | putative malonic semialdehyde reductase |
Reads supporting (aligned to +/- strand): ref base G (1/0); new base T (17/12); total (18/12) | |||||||||||
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.64e-01 |
GTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCAT > NC_000913/1070686‑1070838 | gTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTc > 1:217451/1‑81 (MQ=255) gTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTc > 1:217444/1‑81 (MQ=255) aGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTc > 1:217445/1‑81 (MQ=255) gacgaTAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTc > 2:217439/1‑81 (MQ=255) gaTAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGc > 1:217461/1‑81 (MQ=255) aTAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGcc > 1:217446/1‑81 (MQ=255) tAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCg < 1:217408/81‑1 (MQ=255) cGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGc < 1:217409/81‑1 (MQ=255) cGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGc > 2:217440/1‑81 (MQ=255) tGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCggg > 2:217458/1‑81 (MQ=255) aCGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTa < 2:217410/81‑1 (MQ=255) cGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTaa < 2:217411/81‑1 (MQ=255) ggCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTaaa > 1:217449/1‑80 (MQ=255) gcgcGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAaca > 1:217462/1‑81 (MQ=255) ggTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGAt < 1:217412/81‑1 (MQ=255) tCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATcc > 1:217454/1‑81 (MQ=255) ggTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGc < 1:217413/81‑1 (MQ=255) tttGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCg > 1:217456/1‑78 (MQ=255) cAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGaa > 1:217448/1‑79 (MQ=255) cGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCt > 2:217460/1‑80 (MQ=255) cGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTg < 2:217414/81‑1 (MQ=255) ggAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGt > 2:217468/1‑81 (MQ=255) ggAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGt > 2:217472/1‑81 (MQ=255) gAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTc > 1:217463/1‑81 (MQ=255) gTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGcc < 2:217415/81‑1 (MQ=255) cGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCa < 1:217416/81‑1 (MQ=255) gggCGCAGACTTTCTGTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCAt < 2:217417/81‑1 (MQ=255) gcgcAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCAttt > 1:217457/1‑81 (MQ=255) cAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTtcat < 1:217418/81‑1 (MQ=255) cAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATACGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTtcat < 1:217419/81‑1 (MQ=255) | GTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCAT > NC_000913/1070686‑1070838 |
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: ‑ |