| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 615,101 | T→C | 100% | G250G (GGT→GGC) | serC → | 3‑phosphoserine/phosphohydroxythreonine aminotransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 615,101 | 0 | T | C | 76.4% | 22.3 / 2.4 | 21 | G250G (GGT→GGC) | serC | 3‑phosphoserine/phosphohydroxythreonine aminotransferase |
| Reads supporting (aligned to +/- strand): ref base T (4/1); new base C (0/16); total (4/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.35e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
CCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTATCTGGTCTGGTCTTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAG > minE/615037‑615161 | ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg > 1:2079359/1‑69 (MQ=255)ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg > 1:2927475/1‑69 (MQ=255)ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg > 1:3232092/1‑69 (MQ=255)ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg > 1:3372302/1‑69 (MQ=255)ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg > 1:1859209/1‑69 (MQ=255)ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg > 1:1784204/1‑69 (MQ=255)ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg > 1:1621876/1‑69 (MQ=255)ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg > 1:1470328/1‑69 (MQ=255)ccATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTAtctggtctgg > 1:581050/1‑69 (MQ=255) tCCTCAACGATAACGATTCCATGTTCAACACGCCGCCGACATTTGCCTGGTACTTATctggcctggt < 1:601711/67‑1 (MQ=255) tCCTCAACGATAACGATTCCATGTTCAACACGCCGCCGACATTTGCCTGGTACTTATctggcctggt < 1:2080678/67‑1 (MQ=255) tCCTCAACGATAACGATTCCATGTTCAACACGCCGCCGACATTTGCCTGGTACTTATctggcctggt < 1:169106/67‑1 (MQ=255) tCCTCAACGATAACGATTCCATGTTCAACACGCCGCCGACATTTGCCTGGTACTTATctggcctggt < 1:2652213/67‑1 (MQ=255) tCCTCAACGATAACGATTCCATGTTCAACACGCCGCCGACATTTGCCTGGTACTTATctggcctggt < 1:268536/67‑1 (MQ=255) tCCTCAACGATAACGATTCCATGTTCAACACGCCGCCGACATTTGCCTGGTACTTATctggcctggt < 1:272242/67‑1 (MQ=255) ccTCAACGATAACGATTCCATGTTCAACACGCCGCCGACATTTGCCTGGTACTTATctggcctggt < 1:2123783/66‑1 (MQ=255) cTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTATCTGGTCTGGTCTTTa < 1:2952560/56‑1 (MQ=255) ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg < 1:511090/67‑1 (MQ=255) ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg < 1:383505/67‑1 (MQ=255) ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg < 1:716339/67‑1 (MQ=255) ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg < 1:1186500/67‑1 (MQ=255) ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg < 1:3231803/67‑1 (MQ=255) ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg < 1:2939064/67‑1 (MQ=255) ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg < 1:576087/67‑1 (MQ=255) ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg < 1:2922007/67‑1 (MQ=255) ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg < 1:2765107/67‑1 (MQ=255) ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg < 1:2321859/67‑1 (MQ=255) ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg < 1:2136972/67‑1 (MQ=255) ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg < 1:2012570/67‑1 (MQ=255) ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg < 1:164962/67‑1 (MQ=255) ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg < 1:1649539/67‑1 (MQ=255) ccgACATTTGCCTGGTACTTATCTGGCCTGGTCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg < 1:1592519/67‑1 (MQ=255) ccgACATTTGCCTGGTACTTATCTGGCCTGGGCTTCAAATGGCTGAAAGCGAACGGCGGTGTGGCTg < 1:451959/67‑1 (MQ=255) tatctGGTCTGGTCTTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATg > 1:201560/1‑53 (MQ=255) tatctGGTCTGGTCTTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAatcaatca > 1:2363037/1‑67 (MQ=255) tatctGGTCTGGTCTTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAg > 1:2506641/1‑68 (MQ=255) tatctGGTCTGGTCTTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAg > 1:1719596/1‑68 (MQ=255) | CCATCCTCAACGATAACGGCTCCATGTTTAACACGCCGCCGACATTTGCCTGGTATCTATCTGGTCTGGTCTTTAAATGGCTGAAAGCGAACGGCGGTGTAGCTGAAATGGATAAAATCAATCAG > minE/615037‑615161 |
| 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: ‑ |