| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 190,065 | G→A | 30.5% | R521H (CGT→CAT) | ldcC → | lysine decarboxylase 2, constitutive |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 190,065 | 0 | G | A | 30.5% | 48.8 / 25.1 | 36 | R521H (CGT→CAT) | ldcC | lysine decarboxylase 2, constitutive |
| Reads supporting (aligned to +/- strand): ref base G (11/14); new base A (2/9); total (13/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.59e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01 | |||||||||||
TGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAA > minE/189996‑190133 | tGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGt < 1:1253571/71‑1 (MQ=255) gACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTgg < 1:1690917/71‑1 (MQ=255) gACGAGCAGGGAAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTgg > 1:3765164/1‑71 (MQ=255) aGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgt < 1:316631/70‑1 (MQ=255) gCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgtag > 1:33399/1‑69 (MQ=255) atatGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgt < 1:3482976/64‑1 (MQ=255) cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAcc > 1:1743702/1‑71 (MQ=255) ccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACct < 1:2744473/71‑1 (MQ=255) ggcggcGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGcc < 1:2081397/59‑1 (MQ=255) ggcggcGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgc > 1:728347/1‑71 (MQ=255) gcTGGTAGCAAAATTCCTCGACGAACGTGGGATCgta < 1:2506143/37‑1 (MQ=255) cTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCg > 1:1643744/1‑49 (MQ=255) gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt < 1:77913/71‑1 (MQ=255) gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt > 1:3645604/1‑71 (MQ=255) gaaaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACctgctg < 1:742703/57‑1 (MQ=19) gaaaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACctgctg < 1:3118408/57‑1 (MQ=19) aaaaTTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTa < 1:3716269/71‑1 (MQ=255) aaaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGtata < 1:722563/49‑1 (MQ=12) aaaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCAt < 1:2288050/71‑1 (MQ=255) aaTTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTAt > 1:1495476/1‑70 (MQ=255) aaTACCTCGACGATCATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCAtc > 1:3036022/1‑70 (MQ=255) aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCAt < 1:1999686/70‑1 (MQ=255) aTTCCTCGACGAACGTGGGATCGTAGTAGAGaaaa > 1:504756/1‑35 (MQ=255) aTTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtct < 1:460844/60‑1 (MQ=255) aTTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTg < 1:3616792/71‑1 (MQ=255) cTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTAt < 1:2701062/71‑1 (MQ=255) cTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTAt < 1:2073456/71‑1 (MQ=255) cgacgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATc > 1:2661880/1‑70 (MQ=255) gacgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCg < 1:3073030/70‑1 (MQ=255) gacgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGa > 1:1524719/1‑71 (MQ=255) gacgaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCg < 1:498848/70‑1 (MQ=255) acgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAc > 1:2580940/1‑41 (MQ=255) acgaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGAt < 1:742748/71‑1 (MQ=255) cgaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGAt > 1:2295930/1‑70 (MQ=255) aaCGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctg < 1:1476725/41‑1 (MQ=255) aCGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaaaa < 1:3567069/71‑1 (MQ=255) gTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctg > 1:2180636/1‑38 (MQ=255) | TGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAA > minE/189996‑190133 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |