Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,533,003 | G→A | 57.1% | G230G (GGC→GGT) | prmB ← | N5‑glutamine methyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,533,003 | 0 | G | A | 57.1% | 8.0 / 35.6 | 49 | G230G (GGC→GGT) | prmB | N5‑glutamine methyltransferase |
Reads supporting (aligned to +/- strand): ref base G (3/18); new base A (22/6); total (25/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.55e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01 |
CTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGA > minE/1532945‑1533061 | cTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCg > 1:1377236/1‑69 (MQ=255) tGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCg > 1:81175/1‑68 (MQ=255) gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa < 1:1780471/69‑1 (MQ=255) gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa < 1:2600999/69‑1 (MQ=255) gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa < 1:3078580/69‑1 (MQ=255) gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa < 1:359206/69‑1 (MQ=255) gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa < 1:428468/69‑1 (MQ=255) gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa < 1:471664/69‑1 (MQ=255) gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa < 1:1863726/69‑1 (MQ=255) gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa < 1:1698573/69‑1 (MQ=255) gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa < 1:830593/69‑1 (MQ=255) gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa < 1:1186644/69‑1 (MQ=255) gagaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa < 1:1347696/69‑1 (MQ=255) agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt > 1:1713197/1‑69 (MQ=255) agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt > 1:833071/1‑69 (MQ=255) agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt > 1:69819/1‑69 (MQ=255) agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt > 1:1824874/1‑69 (MQ=255) agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt > 1:865044/1‑69 (MQ=255) agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt > 1:629469/1‑69 (MQ=255) agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt > 1:1434617/1‑69 (MQ=255) agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt > 1:2158315/1‑69 (MQ=255) agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt > 1:1327621/1‑69 (MQ=255) agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt > 1:2622736/1‑69 (MQ=255) agaATGCGACGCGTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTAt > 1:1216770/1‑69 (MQ=255) agaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa < 1:918297/68‑1 (MQ=255) agaATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa < 1:3392177/68‑1 (MQ=255) cgTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGTGTTCCGGCTCGTGGCGGTATTCGTTTGGCAg < 1:2454019/69‑1 (MQ=255) cgTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAg < 1:1363975/69‑1 (MQ=255) cgTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAg < 1:146431/69‑1 (MQ=255) cgTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAg < 1:168046/69‑1 (MQ=255) cgTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAg < 1:1926338/69‑1 (MQ=255) cgTCAGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAg < 1:2079523/69‑1 (MQ=255) gTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTa < 1:345114/56‑1 (MQ=255) aGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGtg > 1:259405/1‑67 (MQ=255) aGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTc > 1:2980474/1‑68 (MQ=255) aGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCg > 1:2585475/1‑69 (MQ=255) aGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCg > 1:1673822/1‑69 (MQ=255) aGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCg > 1:756519/1‑69 (MQ=255) aGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCg > 1:1801012/1‑69 (MQ=255) aGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCg > 1:2300387/1‑69 (MQ=255) aGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCg > 1:3392044/1‑69 (MQ=255) aGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCg > 1:3191484/1‑69 (MQ=255) aGTTTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCg > 1:2729403/1‑69 (MQ=255) aGTCTCAGGCCGTCGGTGCCAGATGCCAGACCGAGTTCCGGCTCGTGTCGGTATTCGTTTGGCAGGTCg > 1:244740/1‑69 (MQ=255) tgccagatgccagaACGAGTTCCGGCTCGTGGCGGTATTCGTTTGTCAg < 1:694455/49‑1 (MQ=37) ccAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATc > 1:385840/1‑68 (MQ=255) aGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGa < 1:2318447/68‑1 (MQ=255) aGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGa < 1:2775895/68‑1 (MQ=255) aGATGCCAGGCCCAGTTCCGGCTCGTGACGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGa < 1:1731235/68‑1 (MQ=255) gATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTCGGCAGGTCGGACATATCTTCCGCATCGa < 1:2356513/67‑1 (MQ=255) | CTGCCGCGTTACCGAGAATGCGACGCGTCAGTTTCAGGCCGTCAGTGCCAGATGCCAGGCCCAGTTCCGGCTCGTGGCGGTATTCGTTTGGCAGGTCGGACATATCTTCCGCATCGA > minE/1532945‑1533061 |
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: ‑ |