| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,126,319 | A→C | 75.7% | G105G (GGT→GGG) | pheS ← | phenylalanine tRNA synthetase, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,126,319 | 0 | A | C | 75.7% | 64.4 / 16.1 | 41 | G105G (GGT→GGG) | pheS | phenylalanine tRNA synthetase, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base A (3/7); new base C (20/11); total (23/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.53e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01 | |||||||||||
GTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGA > minE/1126257‑1126376 | gTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGAccgccg < 1:321473/69‑1 (MQ=255) tAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGAccgccg < 1:1397440/68‑1 (MQ=255) gCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTc < 1:1860737/69‑1 (MQ=255) ccAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCa > 1:2132630/1‑69 (MQ=255) ccAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCa > 1:3205970/1‑69 (MQ=255) cTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGc > 1:2066082/1‑68 (MQ=255) cGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGcgacga < 1:2864777/68‑1 (MQ=255) gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa < 1:2194106/68‑1 (MQ=255) gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa < 1:561889/68‑1 (MQ=255) gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa < 1:171666/68‑1 (MQ=255) gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa < 1:75833/68‑1 (MQ=255) gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa < 1:1544558/68‑1 (MQ=255) gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa < 1:1379135/68‑1 (MQ=255) gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa < 1:879628/68‑1 (MQ=255) gTCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa < 1:1332349/68‑1 (MQ=255) tCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa < 1:1073963/67‑1 (MQ=255) tCGATGGTACGGGTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa < 1:882632/67‑1 (MQ=255) gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGCCGACCCGGCAGGGAGACATCAATCGTTTccgcc > 1:2966534/1‑68 (MQ=255) gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATc > 1:731398/1‑55 (MQ=255) gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGttt > 1:2339030/1‑63 (MQ=255) gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc > 1:526925/1‑68 (MQ=255) gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc > 1:3403063/1‑68 (MQ=255) gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc > 1:3329/1‑68 (MQ=255) gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc > 1:740394/1‑68 (MQ=255) gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc > 1:812173/1‑68 (MQ=255) gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc > 1:2599318/1‑68 (MQ=255) gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc > 1:2726898/1‑68 (MQ=255) gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc > 1:2604776/1‑68 (MQ=255) gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc > 1:2498043/1‑68 (MQ=255) gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc > 1:2267239/1‑68 (MQ=255) gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc > 1:1998917/1‑68 (MQ=255) gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc > 1:1926931/1‑68 (MQ=255) gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc > 1:1840568/1‑68 (MQ=255) gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc > 1:1643913/1‑68 (MQ=255) gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc > 1:1489509/1‑68 (MQ=255) gggTCACCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc > 1:1262032/1‑68 (MQ=255) gggTCACCGGATGCAGCCAGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCAATCGTTTccgcc > 1:3399510/1‑68 (MQ=255) tcaCCGGATGCAGCCCGCCGTTTTCAATGCGACGACCCGGCAGGGAGACATCaa < 1:2658158/52‑1 (MQ=255) aCCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGa < 1:3052475/69‑1 (MQ=255) aCCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGa < 1:3216484/69‑1 (MQ=255) aCCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGa < 1:1344406/69‑1 (MQ=255) | GTAAAGCCAAGCTCACCGAAGAAACTTTCGATACGGTCGATGGTACGGGTAACCGGATGCAGACCGCCGTTTTCAATGCGACGACCTGGCAGAGAGACATCAATCGTTTCCGCCGCCAGA > minE/1126257‑1126376 |
| 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: ‑ |