Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 3,085,073 | G→C | G194A (GGC→GCC) | rsmE → | 16S rRNA (uracil(1498)‑N(3))‑methyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,085,073 | 0 | G | C | 80.0% | 19.6 / 1.1 | 10 | G194A (GGC→GCC) | rsmE | 16S rRNA (uracil(1498)‑N(3))‑methyltransferase |
Reads supporting (aligned to +/- strand): ref base G (1/1); new base C (3/5); total (4/6) | |||||||||||
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.51e-01 |
GAATCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCC > NZ_CP009273/3084999‑3085138 | gaaTCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGCCCCGGAAGGCGGttt > 1:376687/1‑90 (MQ=255) gcgcCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGCCCCGGAAGGCGGTTTATCGGCAGATGaaa < 1:133704/90‑1 (MQ=255) aGTAACAGCATCAATCCGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGCCCCGGAAGGCGGTTTATCGGCAGATGAAATTGcc < 1:404000/90‑1 (MQ=255) ttGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCaa > 1:87571/1‑90 (MQ=255) gTTACCGGTTGAACGCGTCCGCCTGCTGATTGCCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTAc < 1:218729/90‑1 (MQ=255) cGGTTGAACGCGTCCGCCTGCTGATTGCCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGata < 2:354325/90‑1 (MQ=255) gTTGAACGCGTCCGCCTGCTGATTGCCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATc > 1:37922/1‑90 (MQ=255) ttGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATcc < 2:402544/90‑1 (MQ=255) gATTGCCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATc < 1:365911/70‑1 (MQ=255) gATTGCCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATc > 2:365911/1‑70 (MQ=255) | GAATCTTCACCCGCGCGCCAGTAACAGCATCAATACGTTGCCGTTACCGGTTGAACGCGTCCGCCTGCTGATTGGCCCGGAAGGCGGTTTATCGGCAGATGAAATTGCCATGACTGCCCGCTATCAATTTACTGATATCC > NZ_CP009273/3084999‑3085138 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |