| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,509,579 | G→A | 30.9% | A42T (GCA→ACA) | ackA → | acetate kinase A and propionate kinase 2 |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,509,579 | 0 | G | A | 30.9% | 79.8 / 42.6 | 55 | A42T (GCA→ACA) | ackA | acetate kinase A and propionate kinase 2 |
| Reads supporting (aligned to +/- strand): ref base G (16/22); new base A (5/12); total (21/34) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.49e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ATCATCGATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCTTTAGGTGCAGGCGCCGCTCA > minE/1509513‑1509643 | atcatcGATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACg < 1:74760/71‑1 (MQ=255) tcatcGATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGt < 1:1997282/71‑1 (MQ=255) cGATGCAGTAAATGGTGAAGATTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATc > 1:185271/1‑70 (MQ=255) cGATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATc > 1:503030/1‑70 (MQ=255) gATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaa < 1:1972831/71‑1 (MQ=255) gATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaa < 1:1282013/71‑1 (MQ=255) gATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaa < 1:1284660/71‑1 (MQ=255) gATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCa > 1:737569/1‑70 (MQ=255) gATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCa < 1:2061452/70‑1 (MQ=255) gATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCa < 1:851385/70‑1 (MQ=255) gTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCaaatgga > 1:1342000/1‑70 (MQ=255) gTAAATGGTGAAGAGTACCTTTCTGGTTTAGACGAATGTTTCCACCTGCCCGAAGCACGTATCa < 1:1343242/64‑1 (MQ=255) aaTGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaaatggaaaa < 1:1982101/70‑1 (MQ=255) aaTGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCaaatggaaaat > 1:1354607/1‑71 (MQ=255) aaTGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCaaatggaaaat > 1:156113/1‑71 (MQ=255) aTGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaaatggaaaatg < 1:1095843/71‑1 (MQ=255) aTGGGGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaaatggaaaatg < 1:467874/71‑1 (MQ=255) ggTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaaatggaaaatgga > 1:42599/1‑71 (MQ=255) ggTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaaatggaaaatgga > 1:1893927/1‑71 (MQ=255) gTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaaatggaaaatgg > 1:303430/1‑69 (MQ=255) aaGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaaatgg < 1:393136/59‑1 (MQ=255) gagTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaaatggaa > 1:474362/1‑59 (MQ=255) gagTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCaa < 1:2200792/71‑1 (MQ=255) agTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCaa > 1:1305914/1‑70 (MQ=255) gTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCAAATGGAAAATGGACGGCAATa < 1:388835/71‑1 (MQ=255) tACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCa < 1:869092/49‑1 (MQ=255) tACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCa < 1:10670/49‑1 (MQ=255) tACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCa < 1:1111094/49‑1 (MQ=255) tACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATa > 1:1060067/1‑70 (MQ=255) tACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATa > 1:1431783/1‑70 (MQ=255) tACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCa < 1:1693450/49‑1 (MQ=255) tACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCa < 1:232851/49‑1 (MQ=255) ccTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCAAATGGAAAATGGACGGCAATaaa < 1:50033/70‑1 (MQ=255) ccTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCAAATGGAAAATGGACGGCAATaaa < 1:306416/70‑1 (MQ=255) tCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCaaatggaaaa > 1:41597/1‑52 (MQ=255) cTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGa < 1:686138/70‑1 (MQ=255) cTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGa < 1:1073730/70‑1 (MQ=255) cTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGa < 1:234917/70‑1 (MQ=255) cTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCAAATGGAAAATGGACGGCAATAAACAGGa < 1:1227913/70‑1 (MQ=255) gTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATaaa > 1:232628/1‑62 (MQ=255) gCCGAATGTTTCCACCTGCCCGAAACACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGAGGCttt < 1:370385/71‑1 (MQ=255) gAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCTTTAgg > 1:28689/1‑71 (MQ=255) gAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCTTTAgg > 1:26102/1‑71 (MQ=255) aTGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCTTTAgg > 1:1359644/1‑69 (MQ=255) aTGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCTTTAg > 1:1637691/1‑68 (MQ=255) tttCCACCTGCCCGAAACACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCttt < 1:1972934/63‑1 (MQ=255) aCCTGCCCGAAACACGTATCAAATGGAAAATGGACGGCAATa < 1:1933316/42‑1 (MQ=255) ccTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCttt < 1:2148543/57‑1 (MQ=255) ccTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCttt < 1:498338/57‑1 (MQ=255) cTGCCCGAAACACGTATCAAATGGAAAATGGACGGCAATaaa > 1:1586534/1‑42 (MQ=255) cTGCCCGAAACACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCTTTAGGt < 1:1178269/60‑1 (MQ=255) cccGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCTTTAGGTGCAGGCGCCGCt < 1:1993433/69‑1 (MQ=255) cccGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCTTTAGGTGCAGGCGCCGCTca < 1:63066/71‑1 (MQ=255) ccGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCttt < 1:1806343/52‑1 (MQ=255) aaGCACGTATCAAATGGAAAATGGACGGCAATaaa > 1:291289/1‑35 (MQ=255) gCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCTTTa < 1:1231779/48‑1 (MQ=255) | ATCATCGATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCTTTAGGTGCAGGCGCCGCTCA > minE/1509513‑1509643 |
| 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: ‑ |