Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,236,903 | G→A | 45.9% | A150T (GCT→ACT) | purD → | phosphoribosylglycinamide synthetase phosphoribosylamine‑glycine ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,236,903 | 0 | G | A | 45.9% | 19.2 / 81.4 | 61 | A150T (GCT→ACT) | purD | phosphoribosylglycinamide synthetase phosphoribosylamine‑glycine ligase |
Reads supporting (aligned to +/- strand): ref base G (20/13); new base A (18/10); total (38/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.97e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTAGAACCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCT > minE/2236837‑2236970 | gTAGAACCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGc > 1:2264255/1‑71 (MQ=255) aGAACCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCg > 1:3185880/1‑71 (MQ=255) aCCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGa < 1:777273/71‑1 (MQ=255) ccTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGaa < 1:3185941/71‑1 (MQ=255) ccTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGa > 1:153489/1‑70 (MQ=255) cgcTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAg < 1:4009243/69‑1 (MQ=255) gCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGtt < 1:1725880/69‑1 (MQ=255) cGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATc < 1:1248767/71‑1 (MQ=255) cGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATc < 1:2476028/71‑1 (MQ=255) aTCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGAGTGCCGGGa < 1:366115/57‑1 (MQ=255) tCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAgg < 1:2672240/60‑1 (MQ=255) tCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTgg > 1:3117701/1‑71 (MQ=255) tCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTgg > 1:3839979/1‑71 (MQ=255) tCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTgg > 1:1775064/1‑71 (MQ=255) tCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGa < 1:670178/56‑1 (MQ=255) tGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCg > 1:1497393/1‑66 (MQ=255) gCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGa < 1:747367/53‑1 (MQ=255) cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGc < 1:1016544/68‑1 (MQ=255) cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGc < 1:2485555/68‑1 (MQ=255) agaAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGa > 1:3731541/1‑66 (MQ=255) aaaGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCt > 1:790524/1‑71 (MQ=255) aaaGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTTCCGGGAAAGGCg > 1:1987108/1‑52 (MQ=255) aaaGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCg > 1:1086079/1‑52 (MQ=255) gcgcgcCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGaaa > 1:812703/1‑44 (MQ=255) aaTCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcg > 1:2270006/1‑70 (MQ=255) aTCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTAt > 1:1954187/1‑44 (MQ=255) aTCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTAt > 1:1426525/1‑44 (MQ=255) aTCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTAt > 1:3061269/1‑44 (MQ=255) tCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGAt < 1:1535155/52‑1 (MQ=255) tCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGAt < 1:1530242/52‑1 (MQ=255) tCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTg < 1:4051871/59‑1 (MQ=255) tCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTg < 1:3760827/59‑1 (MQ=255) cGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGa > 1:536630/1‑53 (MQ=255) cGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGa > 1:2039352/1‑53 (MQ=255) tCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACg < 1:516141/53‑1 (MQ=255) tCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACg < 1:967498/53‑1 (MQ=255) tCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACg < 1:2863338/53‑1 (MQ=255) tCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg < 1:2685062/71‑1 (MQ=255) cATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGaagcggaagcgg > 1:2017247/1‑71 (MQ=255) cATTAAAGCGGAAGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGaagcggaagcgg > 1:1598905/1‑71 (MQ=255) aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGc > 1:741351/1‑71 (MQ=255) tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaagaa < 1:858603/58‑1 (MQ=255) tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTg > 1:1039613/1‑71 (MQ=255) tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTg > 1:1710888/1‑71 (MQ=255) tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTg > 1:617218/1‑71 (MQ=255) tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTg > 1:3245023/1‑71 (MQ=255) tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTg > 1:914242/1‑71 (MQ=255) aaaGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGa > 1:4250096/1‑46 (MQ=255) aaaGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGa > 1:269238/1‑46 (MQ=255) aaaGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGa > 1:1806436/1‑46 (MQ=255) aaaGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGa > 1:3829756/1‑46 (MQ=255) aaaGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGt > 1:2567528/1‑71 (MQ=255) aaaGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGt > 1:3874952/1‑71 (MQ=255) aaGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGtt < 1:690621/71‑1 (MQ=255) gACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGa > 1:524953/1‑71 (MQ=255) gACGGTCTGGCTGCCGGGAAAGGCGTTATCGTAGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGa > 1:649557/1‑71 (MQ=255) tCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGc > 1:2859478/1‑71 (MQ=255) tCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGc > 1:1514422/1‑71 (MQ=255) gcTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGc < 1:2267709/70‑1 (MQ=255) cTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCt > 1:4005696/1‑71 (MQ=255) cTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCt > 1:3107326/1‑71 (MQ=255) | GTAGAACCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCT > minE/2236837‑2236970 |
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: ‑ |