Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,236,903 | G→A | 27.7% | 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 | 27.7% | 108.4 / 45.9 | 65 | A150T (GCT→ACT) | purD | phosphoribosylglycinamide synthetase phosphoribosylamine‑glycine ligase |
Reads supporting (aligned to +/- strand): ref base G (27/20); new base A (7/11); total (34/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.67e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01 |
ACCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCTGGC > minE/2236842‑2236973 | aCCTGCGCTGGCGTATCTGCGTGAGAAAGTCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGa < 1:3042397/71‑1 (MQ=255) aCCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGa < 1:1098481/71‑1 (MQ=255) ccTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGa > 1:2167761/1‑70 (MQ=255) ccTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGa > 1:787939/1‑70 (MQ=255) gCGTATCTGCGTGAGAAAGGCGGGCCAATCGTCATTAAGGCGGACGGTCTGGCTGCCGGGaaa > 1:669681/1‑63 (MQ=255) gCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGaaa > 1:3328702/1‑63 (MQ=255) gCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTAt > 1:3436195/1‑71 (MQ=255) gCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTAt > 1:2042696/1‑71 (MQ=255) gCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTAt > 1:3482935/1‑71 (MQ=255) tGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGtt < 1:331035/62‑1 (MQ=255) tGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGtt < 1:412297/62‑1 (MQ=255) tGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGc > 1:2671843/1‑70 (MQ=255) tGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGc > 1:435384/1‑70 (MQ=255) cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGc < 1:841690/68‑1 (MQ=255) cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGc < 1:1143134/68‑1 (MQ=255) cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGc < 1:2311075/68‑1 (MQ=255) cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGc < 1:2202056/68‑1 (MQ=255) cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGc < 1:2165099/68‑1 (MQ=255) cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGAt < 1:2239992/71‑1 (MQ=255) cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGAt < 1:2155433/71‑1 (MQ=255) cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGAt < 1:3714704/71‑1 (MQ=255) cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGAt < 1:3935053/71‑1 (MQ=255) tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGa > 1:1958875/1‑71 (MQ=255) tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGa > 1:2657215/1‑71 (MQ=255) tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGa > 1:1436898/1‑71 (MQ=255) gagaAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGa < 1:2479587/70‑1 (MQ=255) aGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTgg < 1:3385566/71‑1 (MQ=255) aGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTgg < 1:1800926/71‑1 (MQ=255) cgcgcCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaag < 1:2965709/71‑1 (MQ=255) cgcCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaaga < 1:1701084/70‑1 (MQ=255) aTCGTCATTAAAGCGGACGGTCTGGCTGCCGGGaaa < 1:2567473/36‑1 (MQ=255) tCGTCATTAAAGCTGACGGTCTGGCTGCCGGGaaa < 1:3482463/35‑1 (MQ=255) tCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTAt < 1:3051856/43‑1 (MQ=255) cGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaag > 1:3292206/1‑71 (MQ=255) cGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaa > 1:926289/1‑70 (MQ=255) gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTAt < 1:1046668/41‑1 (MQ=255) gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCg > 1:3634237/1‑43 (MQ=255) gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCg > 1:1544179/1‑43 (MQ=255) gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATg > 1:3217465/1‑51 (MQ=255) gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATg > 1:1468003/1‑51 (MQ=255) gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATg > 1:235601/1‑51 (MQ=255) gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATg > 1:2405440/1‑51 (MQ=255) gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagc > 1:279236/1‑71 (MQ=255) gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagc > 1:2046206/1‑71 (MQ=255) gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagc > 1:2338066/1‑71 (MQ=255) gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGAGTTATCg > 1:2035120/1‑43 (MQ=255) gTCATAAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATg > 1:2108892/1‑51 (MQ=255) tCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGtt < 1:1215150/38‑1 (MQ=255) tCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGtt < 1:868098/38‑1 (MQ=255) tCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGtt < 1:2899263/38‑1 (MQ=255) tCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGtt < 1:4021/38‑1 (MQ=255) tCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGtt < 1:37305/38‑1 (MQ=255) tCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGtt < 1:3814224/38‑1 (MQ=255) cATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg > 1:4118916/1‑70 (MQ=255) cATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGaagcggaagcgg > 1:3002319/1‑71 (MQ=255) aaaGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGt > 1:2496986/1‑71 (MQ=255) gACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGa < 1:2042139/71‑1 (MQ=255) gACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGa < 1:2487637/71‑1 (MQ=255) gACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGa < 1:3126265/71‑1 (MQ=255) gACGGTATGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGa > 1:531738/1‑71 (MQ=255) cGGTCTGGCTGCCGGGAAAGGCGTTATCGTGTCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGata > 1:3164264/1‑71 (MQ=255) tGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCTg > 1:1865823/1‑71 (MQ=255) tGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCTg > 1:448438/1‑71 (MQ=255) tGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCTg > 1:3237202/1‑71 (MQ=255) tGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCTg > 1:2600266/1‑71 (MQ=255) gCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCt < 1:808172/68‑1 (MQ=255) gCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCTGGc > 1:1515742/1‑71 (MQ=255) | ACCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCTGGC > minE/2236842‑2236973 |
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: ‑ |