| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,236,903 | G→A | 31.4% | 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 | 31.4% | 46.8 / 25.6 | 35 | A150T (GCT→ACT) | purD | phosphoribosylglycinamide synthetase phosphoribosylamine‑glycine ligase |
| Reads supporting (aligned to +/- strand): ref base G (10/14); new base A (8/3); total (18/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.46e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AGAACCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCT > minE/2236839‑2236970 | aGAACCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCg > 1:271203/1‑71 (MQ=255) gAACCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCgg > 1:502610/1‑71 (MQ=255) gAACCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCgg > 1:2007309/1‑71 (MQ=255) cgcTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGc < 1:1669060/71‑1 (MQ=255) gcTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCg > 1:1917589/1‑71 (MQ=255) gCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTAt < 1:1467024/71‑1 (MQ=255) tATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGa < 1:540010/58‑1 (MQ=255) aTCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTg > 1:820191/1‑71 (MQ=255) tGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGc > 1:787701/1‑70 (MQ=255) cGTGAGAAATGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGAt < 1:1303832/71‑1 (MQ=255) cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGAt < 1:744128/71‑1 (MQ=255) cGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGAt < 1:1157528/71‑1 (MQ=255) gTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGAt > 1:473324/1‑70 (MQ=255) aaaGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAgg < 1:1272340/50‑1 (MQ=255) gcgcCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaag > 1:1851649/1‑70 (MQ=255) ccAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATc < 1:1946935/48‑1 (MQ=255) ccAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGa < 1:1952773/37‑1 (MQ=255) cAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTAt > 1:1865072/1‑46 (MQ=255) aTCGTCATTATAGCGGACGGGCTGACTGCCGGGAAAGGCGTTATCg < 1:1494594/46‑1 (MQ=38) aTCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGtt < 1:464847/42‑1 (MQ=255) aTCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcgg > 1:431164/1‑70 (MQ=255) tCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGc < 1:1840497/49‑1 (MQ=255) tCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcgg > 1:2165399/1‑69 (MQ=255) gTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagc > 1:872723/1‑71 (MQ=255) ggACGGTCTGGCTGCCGGGAAAGGCGTTATAGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCAc > 1:996815/1‑70 (MQ=255) gACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGa < 1:1194140/71‑1 (MQ=255) gACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaa > 1:1329693/1‑57 (MQ=255) gACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGa > 1:694342/1‑71 (MQ=255) gACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGa > 1:291573/1‑71 (MQ=255) gACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGa > 1:766653/1‑71 (MQ=255) cGGTCTGGCTGCCGGGAAAGGCGTTATCGTGTCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGata < 1:1507771/71‑1 (MQ=255) cGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGata < 1:2112813/71‑1 (MQ=255) cGGGCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaag < 1:1378105/47‑1 (MQ=255) gTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATg < 1:553251/71‑1 (MQ=255) cTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCt > 1:465639/1‑71 (MQ=255) | AGAACCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCGGCTGTTCACGATATGCT > minE/2236839‑2236970 |
| 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: ‑ |