Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 190,065 | G→A | 11.7% | 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 | 11.7% | 147.0 / 8.6 | 51 | R521H (CGT→CAT) | ldcC | lysine decarboxylase 2, constitutive |
Reads supporting (aligned to +/- strand): ref base G (28/17); new base A (4/2); total (32/19) | |||||||||||
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.99e-01 |
TGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAAC > minE/189996‑190134 | tGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGt < 1:2251199/71‑1 (MQ=255) gACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGttg > 1:1545284/1‑69 (MQ=255) gACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTg < 1:3456136/70‑1 (MQ=255) gAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGa > 1:550116/1‑70 (MQ=255) gCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCg < 1:893249/71‑1 (MQ=255) gggCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgtag > 1:1280975/1‑71 (MQ=255) gCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgtagta < 1:2844800/71‑1 (MQ=255) gCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCg < 1:1913072/66‑1 (MQ=255) tatGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTagag > 1:3342698/1‑70 (MQ=255) tGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgta > 1:995718/1‑62 (MQ=255) tGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgta > 1:2987016/1‑62 (MQ=255) tGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaa > 1:927767/1‑70 (MQ=255) aggGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaaaa > 1:2771031/1‑63 (MQ=255) gggATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaaa < 1:2156612/60‑1 (MQ=255) aTCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATaa > 1:3132968/1‑71 (MQ=255) cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCgta > 1:3784601/1‑46 (MQ=255) cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAAcc > 1:3721888/1‑58 (MQ=255) cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAcc > 1:1751251/1‑71 (MQ=255) cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAcc > 1:397263/1‑71 (MQ=255) cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAc > 1:3751403/1‑70 (MQ=255) ccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAAc > 1:1664714/1‑56 (MQ=255) ggcggcGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgc < 1:3767392/71‑1 (MQ=255) gcggagCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgct > 1:2416242/1‑71 (MQ=255) gcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtc < 1:2875556/71‑1 (MQ=255) ggTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAc > 1:3534410/1‑58 (MQ=255) gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTaga < 1:628541/38‑1 (MQ=255) gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc < 1:4049635/69‑1 (MQ=255) gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt > 1:2480545/1‑71 (MQ=255) gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt < 1:1680610/71‑1 (MQ=255) gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt < 1:2393807/71‑1 (MQ=255) gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt < 1:1466335/71‑1 (MQ=255) aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGaa > 1:1306369/1‑39 (MQ=255) aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt > 1:1335423/1‑69 (MQ=255) gCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAg < 1:738377/71‑1 (MQ=255) aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCAtc > 1:989448/1‑70 (MQ=255) ataccTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCAt < 1:1156892/66‑1 (MQ=255) ttCCTCGACGAACGTGGGATCGTAGTAGAGAAAAcc > 1:3688156/1‑36 (MQ=255) ttCCTCGACGAACGTGGGATCGTAGTAGAGAAAAcc > 1:1837572/1‑36 (MQ=255) accTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTcctg < 1:2034573/58‑2 (MQ=19) ccTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCa > 1:2442605/1‑71 (MQ=255) cTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCAt > 1:3512026/1‑71 (MQ=255) cTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCAt > 1:3564507/1‑71 (MQ=255) tCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATc > 1:3450749/1‑71 (MQ=255) gacgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTta > 1:2721051/1‑38 (MQ=255) gacgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATc < 1:1918597/69‑1 (MQ=255) gacgaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGa > 1:1369840/1‑71 (MQ=255) gaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAg < 1:178726/56‑1 (MQ=255) gaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTcct > 1:2411829/1‑50 (MQ=16) gaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTcct > 1:1916936/1‑50 (MQ=16) cGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAc > 1:1728748/1‑36 (MQ=255) cGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaaaaa < 1:2821963/71‑2 (MQ=255) gTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGt < 1:2450675/53‑1 (MQ=255) | TGGACGAGCAGGGCAATATGAGCGAGGAGGGGATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAAC > minE/189996‑190134 |
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: ‑ |