| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 880,072 | C→T | 59.0% | A187A (GCC→GCT) | yciM → | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 880,072 | 0 | C | T | 59.0% | 9.7 / 25.1 | 39 | A187A (GCC→GCT) | yciM | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base C (3/13); new base T (11/12); total (14/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.30e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01 | |||||||||||
GGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGC > minE/880012‑880134 | ggTGAAGCTTGGTAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCCCTgcagca < 1:452480/69‑1 (MQ=255)ggTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCCCTgcagca < 1:585026/69‑1 (MQ=255)ggTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCCCTgcagca < 1:1770882/69‑1 (MQ=255)ggTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCCCTgcagca < 1:2514956/69‑1 (MQ=255)ggTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCCCTgcagca < 1:1040941/69‑1 (MQ=255)ggTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCCCTgcagca < 1:1888192/69‑1 (MQ=255)ggTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCCCTgcagca < 1:2358594/69‑1 (MQ=255)ggTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCCCTgcagca < 1:23169/69‑1 (MQ=255)ggTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCCCTgcagca < 1:2160961/69‑1 (MQ=255) tAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTGCAGCATATGGCCAGc < 1:2419426/67‑1 (MQ=255) tAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTGCAGCATATGGCCAGc < 1:995707/67‑1 (MQ=255) tAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTGCAGCATATGGCCAGc < 1:936385/67‑1 (MQ=255) tAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTGCAGCATATGGCCAGc < 1:627916/67‑1 (MQ=255) tAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTGCAGCATATGGCCAGc < 1:3155922/67‑1 (MQ=255) tAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTGCAGCATATGGCCAGc < 1:3095434/67‑1 (MQ=255) tAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTGCAGCATATGGCCAGc < 1:2718285/67‑1 (MQ=255) tAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTGCAGCATATGGCCAGc < 1:2113283/67‑1 (MQ=255) tAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTGCAGCATATGGCCAGc < 1:1122283/67‑1 (MQ=255) tAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTGCAGCATATGGCCAGc < 1:1786300/67‑1 (MQ=255) tAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTGCAGCATATGGCCAGc < 1:1422715/67‑1 (MQ=255) aaaGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTGCAGCATATGGCCAGc < 1:595686/66‑1 (MQ=255) cgTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTgca > 1:2040484/1‑39 (MQ=255) cgTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTgca > 1:681834/1‑39 (MQ=255) cgTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTgca > 1:1231627/1‑39 (MQ=255) cgTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTgca > 1:453673/1‑39 (MQ=255) cgTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTgca > 1:1296539/1‑39 (MQ=255) cgTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTgca > 1:3106415/1‑39 (MQ=255) cgTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTgca > 1:286747/1‑39 (MQ=255) cgTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTgca > 1:15557/1‑39 (MQ=255) cgTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTgca > 1:2659753/1‑39 (MQ=255) cgTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTgca > 1:1809185/1‑39 (MQ=255) cgTCGAAATTGCCCATTTCTACTGTGAGTTAGCTCTgca > 1:2386435/1‑39 (MQ=255) gAAATTGCCCATTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATg < 1:3302772/69‑1 (MQ=255) gtgAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAgggg > 1:3327714/1‑69 (MQ=255) gtgAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAggg < 1:1411910/68‑1 (MQ=255) gtgAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAggg < 1:662739/68‑1 (MQ=255) gtgAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAggg < 1:859424/68‑1 (MQ=255) gtgAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAggg > 1:1131728/1‑68 (MQ=255) tAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGggcggc > 1:1160526/1‑67 (MQ=255) | GGTGAAGCTGGGTAAAGATAAACAGCGCGTCGAAATTGCCCATTTCTACTGTGAGTTAGCCCTGCAGCATATGGCCAGCGACGATCTCGATCGTGCCATGACCTTGCTAAAAAAAGGGGCGGC > minE/880012‑880134 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |