Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,091,305 | A→G | S58G (AGC→GGC) | rsmE → | 16S rRNA m(3)U1498 methyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,091,305 | 0 | A | G | 87.6% | 362.8 / 38.1 | 121 | S58G (AGC→GGC) | rsmE | 16S rRNA m(3)U1498 methyltransferase |
Reads supporting (aligned to +/- strand): ref base A (7/8); new base G (58/48); total (65/56) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.91e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01 |
Only 100 of 122 total aligned reads displayed.
ATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGA > NC_000913/3091256‑3091351 | aTTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCa < 1:84106/50‑1 (MQ=255) ttGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAg < 2:3703190/50‑1 (MQ=255) gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc < 2:5374836/49‑1 (MQ=255) gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc > 1:3888300/1‑49 (MQ=255) gTTTGACGGAAGCAACCAGGTCTTTGACGCCGACATTACCAGCGCCGGc > 2:3703326/1‑49 (MQ=38) tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa > 1:557445/1‑49 (MQ=255) tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa < 1:2350718/49‑1 (MQ=255) tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCaaa < 2:4398021/49‑1 (MQ=255) aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa < 1:3753338/49‑1 (MQ=255) aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa > 1:5382530/1‑49 (MQ=255) aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa < 2:1041653/49‑1 (MQ=255) cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaaa < 1:4368400/50‑1 (MQ=255) cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaaa < 1:3152660/50‑1 (MQ=255) cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa > 2:4083509/1‑48 (MQ=255) cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa > 2:2399948/1‑48 (MQ=255) cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa > 1:3460044/1‑48 (MQ=255) cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa > 1:4582075/1‑48 (MQ=255) cGGTAGCAACCAGGTCTTTGACGCAGAAATTACCAGCGCCGGCaaaaa > 2:1461778/1‑48 (MQ=38) ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg < 2:4578382/50‑1 (MQ=255) ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg < 2:3035670/50‑1 (MQ=255) ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg > 2:4929008/1‑50 (MQ=255) gTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc < 2:3811783/50‑1 (MQ=255) gTAGCAACCAGGTATTTGCCGCCGAAATTACCAGCGCCGGCAAAAAAAGc < 1:3919910/50‑1 (MQ=255) tAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCg > 1:5011585/1‑50 (MQ=255) aGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt > 2:4945200/1‑50 (MQ=255) gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt < 1:2578313/49‑1 (MQ=255) gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg > 1:1398859/1‑50 (MQ=255) gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg > 1:1595759/1‑50 (MQ=255) cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg > 1:1548486/1‑49 (MQ=255) aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg > 2:2688049/1‑49 (MQ=255) aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGa < 2:4788840/50‑1 (MQ=255) aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGa < 2:559854/50‑1 (MQ=255) cccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg > 2:3614357/2‑48 (MQ=255) aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGaa < 2:2783688/50‑1 (MQ=255) aCCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTgg > 1:2328986/1‑48 (MQ=255) aCCAGGTCTTTGACGCCGAAAGTACCAGCGCCGGCAAAAAAAGCGTGGaa < 1:2060755/50‑1 (MQ=38) ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg < 1:1674216/50‑1 (MQ=255) ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg < 2:2036910/50‑1 (MQ=255) ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg > 1:2988985/1‑50 (MQ=255) ccAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAg < 1:3615174/50‑1 (MQ=255) cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt < 2:1675711/50‑1 (MQ=255) cAGGCCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGCGGAAGt > 1:2446451/1‑50 (MQ=39) aGGTCTTTGACGCCGAAATTACCAGCGCCGGC‑AAAAAAGCGTGGAAGTg > 2:191516/1‑49 (MQ=38) ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg > 2:5249183/1‑49 (MQ=255) ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa < 1:2202967/50‑1 (MQ=255) ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa < 2:4917500/50‑1 (MQ=255) tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAGAGCGTGGAAGTGAAg > 2:2284692/1‑50 (MQ=39) tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg > 1:4367631/1‑50 (MQ=255) cTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg > 2:3642779/1‑50 (MQ=255) tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt > 1:2346862/1‑50 (MQ=255) tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt > 1:5362593/1‑50 (MQ=255) tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt > 2:3877358/1‑50 (MQ=255) tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt > 2:1215642/1‑50 (MQ=255) tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt > 1:1163114/1‑50 (MQ=255) ttGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTg > 2:1893706/1‑50 (MQ=255) tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg < 1:1568639/47‑1 (MQ=255) tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt < 1:594942/48‑1 (MQ=255) tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt < 2:5280449/48‑1 (MQ=255) tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg < 2:3989031/49‑1 (MQ=255) tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg < 1:2502245/49‑1 (MQ=255) tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc > 2:1274190/1‑50 (MQ=255) tGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAgg < 1:1252942/47‑1 (MQ=255) tGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTg < 2:4528036/49‑1 (MQ=255) gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc < 1:1408848/49‑1 (MQ=255) gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc < 1:2921515/49‑1 (MQ=255) aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg > 1:5147252/1‑50 (MQ=255) cGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg < 1:3962045/50‑1 (MQ=255) gCCGAATTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg > 2:3178475/1‑49 (MQ=38) gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg > 2:1197880/1‑49 (MQ=255) gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg > 1:1988179/1‑49 (MQ=255) gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGa < 1:5182451/50‑1 (MQ=255) gCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCGGGa < 2:442028/50‑1 (MQ=255) ccGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGaa < 2:2273862/50‑1 (MQ=255) cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg > 2:3611718/1‑50 (MQ=255) cGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg > 1:3196964/1‑50 (MQ=255) gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg < 2:4471970/50‑1 (MQ=255) gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg < 1:1907788/50‑1 (MQ=255) gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg > 1:5225806/1‑50 (MQ=255) aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc > 2:17598/1‑50 (MQ=255) aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc < 1:2532334/50‑1 (MQ=255) aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc > 1:4750294/1‑50 (MQ=255) aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc < 1:1831861/50‑1 (MQ=255) aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc > 1:3031480/1‑49 (MQ=255) aaTTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc > 1:4085270/1‑49 (MQ=255) aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc > 2:3290938/1‑48 (MQ=255) aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc > 2:1785931/1‑48 (MQ=255) aTTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa < 1:1106948/50‑1 (MQ=255) ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc > 2:2619659/1‑47 (MQ=255) ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAg < 1:2094877/50‑1 (MQ=255) tACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc > 2:1896194/1‑46 (MQ=255) aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGa < 2:5106026/49‑1 (MQ=255) aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGa < 1:4209492/49‑1 (MQ=255) ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc > 1:5094205/1‑44 (MQ=255) ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc < 1:3460902/50‑1 (MQ=255) ccAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc > 2:2070820/1‑44 (MQ=255) cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg < 2:4702724/50‑1 (MQ=255) cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg < 2:4939352/50‑1 (MQ=255) aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga > 2:1086679/1‑50 (MQ=255) gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac < 2:2476790/50‑1 (MQ=255) gcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacga > 1:2626634/1‑50 (MQ=255) | ATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGA > NC_000913/3091256‑3091351 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |