| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,650,420 | A→T | 100% | L215L (CTA→CTT) | glgX → | glycogen debranching enzyme |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,650,420 | 0 | A | T | 82.3% | 146.9 / 25.4 | 62 | L215L (CTA→CTT) | glgX | glycogen debranching enzyme |
| Reads supporting (aligned to +/- strand): ref base A (1/10); new base T (51/0); total (52/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.02e-10 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
ACCGCGCTGGAACTGCTGCCAGTGGCGCAGTTTGCCAGTGAACCACGTCTGCAACGCATGGGGCTAAGTAACTACTGGGGTTACAACCCGG > minE/2650355‑2650445 | aCCGCGCTGGAACTGCTGCCAGTGGCGCAGTTTGCCAGTGAACCACGTCTGCAACGCATGGGGCTaa < 1:574864/67‑1 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:422813/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1987574/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:2047580/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:2058342/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:2164752/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:2231461/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:240288/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:26790/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:269413/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:30057/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:308207/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:338901/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:377325/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:386688/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:395791/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1052314/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:431602/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:48724/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:540632/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:558206/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:572105/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:597804/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:653458/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:705429/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:79926/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:940952/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:946603/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:95398/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1710519/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1033968/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1144185/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1215194/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1250934/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1261832/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1272197/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1300737/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1327655/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1373084/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1387892/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1437534/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1497388/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1561273/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1782895/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1708024/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1699376/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1646748/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1653406/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1683411/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1692533/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1698600/3‑67 (MQ=255) gttgctgcCGGTGGCACAATTTGCCAGCGAACCACGTCTGCAACGCATGGGGCTTAGTAACTACTgg > 1:1695832/3‑67 (MQ=255) gcCAGTGGCGCAGTTTGCCAGTGAACCACGTCTGCAACGCATGGGGCTAAGTAACTACTGGGGtt < 1:40112/65‑1 (MQ=255) gcCAGTGGCGCAGTTTGCCAGTGAACCACGTCTGCAACGCATGGGGCTAAGTAACTACTGGGGtt < 1:95033/65‑1 (MQ=255) gcCAGTGGCGCAGTTTGCCAGTGAACCACGTCTGCAACGCATGGGGCTAAGTAACTACTGGGGtt < 1:1104069/65‑1 (MQ=255) gcCAGTGGCGCAGTTTGCCAGTGAACCACGTCTGCAACGCATGGGGCTAAGTAACTACTGGGGtt < 1:248442/65‑1 (MQ=255) gcCAGTGGCGCAGTTTGCCAGTGAACCACGTCTGCAACGCATGGGGCTAAGTAACTACTGGGGtt < 1:500883/65‑1 (MQ=255) gcCAGTGGCGCAGTTTGCCAGTGAACCACGTCTGCAACGCATGGGGCTAAGTAACTACTGGGGtt < 1:1397979/65‑1 (MQ=255) gcCAGTGGCGCAGTTTGCCAGTGAACCACGTCTGCAACGCATGGGGCTAAGTAACTACTGGGGtt < 1:1645301/65‑1 (MQ=255) gcCAGTGGCGCAGTTTGCCAGTGAACCACGTCTGCAACGCATGGGGCTAAGTAACTACTGGGGtt < 1:1611831/65‑1 (MQ=255) gcCAGTGGCGCAGTTTGCCAGTGAACCACGTCTGCAACGCATGGGGCTAAGTAACTACTGGGGtt < 1:1562038/65‑1 (MQ=255) ccAGTGGCGCAGTTTGCCAGTGAACCACGTCTGCAACGCATGGGGCTAAGTAACTACTGGGGtt < 1:907324/64‑1 (MQ=255) cgcAGTTTGCCAGTGAACCACGTCTGCAACGCATGGGGCTAAGTAACTACTGGGGTTACAACCCgg > 1:997778/1‑66 (MQ=255) | ACCGCGCTGGAACTGCTGCCAGTGGCGCAGTTTGCCAGTGAACCACGTCTGCAACGCATGGGGCTAAGTAACTACTGGGGTTACAACCCGG > minE/2650355‑2650445 |
| 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: ‑ |