Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,952,864 | G→A | 22.5% | intergenic (‑173/+95) | leuQ ← / ← rsmC | tRNA‑Leu/16S RNA m2G1207 methylase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,952,864 | 0 | G | A | 22.5% | 64.6 / 9.6 | 31 | intergenic (‑173/+95) | leuQ/rsmC | tRNA‑Leu/16S RNA m2G1207 methylase |
Reads supporting (aligned to +/- strand): ref base G (14/10); new base A (3/4); total (17/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.71e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCA > minE/2952794‑2952933 | aaaaTGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTgg < 1:3679536/71‑1 (MQ=255) aaaaTGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTgg < 1:2176706/71‑1 (MQ=255) aaaTGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGgc < 1:3566406/71‑1 (MQ=255) aTGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGgcgc < 1:2897069/71‑1 (MQ=255) cATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCt > 1:1130310/1‑70 (MQ=255) tCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTAt > 1:2422417/1‑71 (MQ=255) aGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATc < 1:424187/70‑1 (MQ=255) aaCTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTgg > 1:1989121/1‑55 (MQ=255) aaCTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGcc > 1:3613752/1‑71 (MQ=255) cTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTAt > 1:2795973/1‑63 (MQ=255) gCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTAc < 1:513333/69‑1 (MQ=255) gCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACg > 1:2330371/1‑70 (MQ=255) aGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCt < 1:2009551/71‑1 (MQ=255) tttGGTATGCATGCGGGGGCAGATGCCAGAAGCGACGCTGgcgc < 1:2307398/44‑1 (MQ=255) tGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGgcgc < 1:3035790/42‑1 (MQ=255) tGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAActt > 1:2109310/1‑69 (MQ=255) gTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGc < 1:2008351/71‑1 (MQ=255) atgcatgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAg > 1:2079663/1‑71 (MQ=255) atgcatgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAg > 1:1465149/1‑71 (MQ=255) tgcattcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAgg < 1:2309817/58‑1 (MQ=255) catgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGtt > 1:464906/1‑71 (MQ=255) ggggCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTa < 1:2750459/70‑1 (MQ=255) cagatgccagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAg < 1:950931/67‑1 (MQ=255) tgccagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGAt > 1:1144213/1‑70 (MQ=255) gccagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGAt > 1:1406971/1‑69 (MQ=255) gccagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGAt > 1:1388348/1‑69 (MQ=255) cagatgcGACGCTGACGCGTCTTATCTGGCCTACAAAGGGCTAACGTGCAGGttttt > 1:2483765/1‑56 (MQ=18) cagatgcGACGCTGACGCGTCTTATCTGGCCTACAAAGGGCTAACGTGCAGGttttt > 1:2479270/1‑56 (MQ=18) agatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGc < 1:3136668/71‑1 (MQ=255) atgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGt > 1:3801982/1‑71 (MQ=255) gACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCa > 1:1922309/1‑70 (MQ=255) gACGCTGACGCGTCTTATCTGGCCTACAAAGGGCTa < 1:3618509/36‑1 (MQ=14) gACGCTGACGCGTCTTATCTGGCCTACAAAGGGCTAACGt < 1:2633690/40‑1 (MQ=17) gACGCTGACGCGTCTTATCTGGCCTACAAAGGGCTAACGt > 1:3460426/1‑40 (MQ=17) cGCTGACGCGTCTTATCTGGCCTACAAAGGGCTAACGt < 1:380235/38‑1 (MQ=17) gCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCAcgcc < 1:1581746/71‑1 (MQ=255) gCTGACGCGTCTTATCTGGCCTACAAAGGGCTAACGTGCAGGttt < 1:2356655/45‑1 (MQ=25) cTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGt < 1:1876783/47‑1 (MQ=255) ggCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAg > 1:3599125/1‑47 (MQ=255) gcgcGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCa < 1:428120/70‑1 (MQ=255) | AAAATGTCCATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCA > minE/2952794‑2952933 |
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: ‑ |