Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | CP009273 | 4,173,770 | T→G | D866E (GAT→GAG) | rpoB → | RNA polymerase, beta subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 4,173,770 | 0 | T | G | 92.9% | 129.9 / 1.1 | 42 | D866E (GAT→GAG) | rpoB | RNA polymerase, beta subunit |
Reads supporting (aligned to +/- strand): ref base T (2/1); new base G (14/25); total (16/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.47e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01 |
GTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGT > CP009273/4173709‑4173832 | gTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATcc < 1:3059392/68‑1 (MQ=255) gAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTAt < 1:1473965/68‑1 (MQ=255) gAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTAt < 1:964228/68‑1 (MQ=255) gAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTAt > 1:2640675/1‑68 (MQ=255) aaGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATc > 1:3512319/1‑68 (MQ=255) agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg < 1:3031939/68‑1 (MQ=255) agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg < 1:1261374/68‑1 (MQ=255) agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg < 1:4147609/68‑1 (MQ=255) agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg < 1:4106016/68‑1 (MQ=255) agagATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCg > 1:2536701/1‑68 (MQ=255) gagaTCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGt < 1:4243035/68‑1 (MQ=255) cGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACCTTg > 1:3867035/1‑68 (MQ=255) gCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACAttct > 1:696049/1‑66 (MQ=255) gCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTgg > 1:1289935/1‑68 (MQ=255) gCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTgg > 1:1494148/1‑68 (MQ=255) cTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGt < 1:582275/68‑1 (MQ=255) cTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGt < 1:4196903/68‑1 (MQ=255) tGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTg > 1:1848962/1‑68 (MQ=255) tGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTg > 1:1348744/1‑68 (MQ=255) cATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCgg < 1:3738442/68‑1 (MQ=255) aTCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGa < 1:3093803/68‑1 (MQ=255) cGTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCg > 1:1331528/1‑68 (MQ=255) gTGGGTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCgg < 1:634501/68‑1 (MQ=255) gTGAAGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGc < 1:3299794/68‑1 (MQ=255) aaGCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGAc < 1:3362168/68‑1 (MQ=255) gCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGTCCGGTGGCGACAt > 1:1769664/1‑68 (MQ=255) gCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAt < 1:1249760/68‑1 (MQ=255) gCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAt > 1:3478653/1‑68 (MQ=255) gCTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAt > 1:2298173/1‑68 (MQ=255) cTGCGCTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACAtt > 1:3880043/1‑68 (MQ=255) gcTCTCCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtgg < 1:3960427/68‑1 (MQ=255) ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg < 1:4115686/68‑1 (MQ=255) ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg < 1:204712/68‑1 (MQ=255) ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg < 1:82150/68‑1 (MQ=255) ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg < 1:2673227/68‑1 (MQ=255) ctcCAAACTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg < 1:1068958/68‑1 (MQ=255) tcCAATCTGGAGGAATCTGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttgg < 1:215965/68‑1 (MQ=255) aaCTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaag < 1:3580310/68‑1 (MQ=255) aaCTGGAGGAATCCGGTATCGTATACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaag < 1:2509106/68‑1 (MQ=255) aCTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaagg < 1:1237484/68‑1 (MQ=255) cTGGAGGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaaggt > 1:2561427/1‑68 (MQ=255) cTGGAGGAATCCGGTATCGTTTACAATGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtacggt > 1:3292607/1‑68 (MQ=255) | GTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGT > CP009273/4173709‑4173832 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |