Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,974,485 | A→G | 30.2% | A262A (GCA→GCG) | metK → | methionine adenosyltransferase 1 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,974,485 | 0 | A | G | 30.2% | 67.3 / 35.2 | 76 | A262A (GCA→GCG) | metK | methionine adenosyltransferase 1 |
Reads supporting (aligned to +/- strand): ref base A (36/17); new base G (14/9); total (50/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.04e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01 |
TGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGAAA > minE/1974417‑1974551 | tGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCa < 1:935139/69‑1 (MQ=255) gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTtctccg > 1:2060898/1‑61 (MQ=255) gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg > 1:3215994/1‑64 (MQ=255) gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg > 1:2960879/1‑64 (MQ=255) gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg > 1:3418355/1‑64 (MQ=255) gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg > 1:2548129/1‑64 (MQ=255) gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg > 1:628568/1‑64 (MQ=255) gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg > 1:715876/1‑64 (MQ=255) gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg > 1:2060899/1‑64 (MQ=255) gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg > 1:848540/1‑64 (MQ=255) gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg > 1:1804815/1‑64 (MQ=255) gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg > 1:1427821/1‑64 (MQ=255) gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg > 1:923293/1‑64 (MQ=255) gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg > 1:896607/1‑64 (MQ=255) gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg > 1:1354889/1‑64 (MQ=255) cGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAg < 1:320879/67‑1 (MQ=255) cGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAg < 1:1653393/67‑1 (MQ=255) cGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAg < 1:2528568/67‑1 (MQ=255) gTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATc > 1:1986681/1‑69 (MQ=255) tATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa < 1:3083886/69‑1 (MQ=255) tATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa > 1:1428852/1‑69 (MQ=255) tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa < 1:854794/63‑1 (MQ=255) tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa < 1:3046012/63‑1 (MQ=255) tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa < 1:3283788/63‑1 (MQ=255) tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa < 1:1418921/63‑1 (MQ=255) tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa < 1:1431015/63‑1 (MQ=255) tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa < 1:2132095/63‑1 (MQ=255) tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa < 1:718575/63‑1 (MQ=255) ccTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAGATCCATCAAAAGTGGACCGtt < 1:2567405/69‑1 (MQ=255) ccTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAGATCCATCAAAAGTGGACCGtt < 1:1824285/69‑1 (MQ=255) tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc > 1:539605/1‑68 (MQ=255) tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc > 1:3139853/1‑68 (MQ=255) tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc > 1:313293/1‑68 (MQ=255) tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc > 1:1276951/1‑68 (MQ=255) tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc > 1:42752/1‑68 (MQ=255) tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc > 1:1664341/1‑68 (MQ=255) tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc > 1:2630749/1‑68 (MQ=255) tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc > 1:3190460/1‑68 (MQ=255) tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc > 1:3229188/1‑68 (MQ=255) tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc > 1:2817366/1‑68 (MQ=255) tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc > 1:2945154/1‑68 (MQ=255) tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc > 1:2145291/1‑68 (MQ=255) tACGGCGGCATGGCGCGTCACGGAGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGAc > 1:2172804/1‑63 (MQ=255) gcggcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAg > 1:1214761/1‑69 (MQ=255) gcggcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAg > 1:734721/1‑69 (MQ=255) cggcATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGc < 1:53537/69‑1 (MQ=255) cggcATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGc < 1:654477/69‑1 (MQ=255) cggcATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGc < 1:222714/69‑1 (MQ=255) ggcATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAg < 1:2657417/40‑1 (MQ=25) gcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt > 1:8868/1‑69 (MQ=255) gcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt > 1:3334707/1‑69 (MQ=255) gcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt > 1:40918/1‑69 (MQ=255) gcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt > 1:110733/1‑69 (MQ=255) gcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt > 1:1378215/1‑69 (MQ=255) ggCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACgc < 1:2430509/69‑1 (MQ=255) ggCGCGTCACGGTGGCGGTGCACTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACgc < 1:3292727/69‑1 (MQ=255) cACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGt < 1:2756557/69‑1 (MQ=255) aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg > 1:3069050/1‑42 (MQ=255) aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg > 1:2103211/1‑42 (MQ=255) aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg > 1:2233488/1‑42 (MQ=255) aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg > 1:242123/1‑42 (MQ=255) aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg > 1:2582660/1‑42 (MQ=255) aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg > 1:121259/1‑42 (MQ=255) aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg > 1:2644193/1‑42 (MQ=255) aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGtt > 1:2103713/1‑69 (MQ=255) aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGt < 1:2078837/68‑1 (MQ=255) aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGt < 1:885863/68‑1 (MQ=255) gCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTc > 1:1296277/1‑68 (MQ=255) gCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTc > 1:1088949/1‑68 (MQ=255) gCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTc > 1:152024/1‑68 (MQ=255) gCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTc > 1:81371/1‑68 (MQ=255) gCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTc > 1:2502417/1‑68 (MQ=255) gCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTc > 1:2611672/1‑68 (MQ=255) gCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGaaa < 1:1890090/69‑1 (MQ=255) gCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGaaa < 1:791583/69‑1 (MQ=255) gCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGaaa < 1:662400/69‑1 (MQ=255) gCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGaaa < 1:2982599/69‑1 (MQ=255) | TGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGAAA > minE/1974417‑1974551 |
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: ‑ |