| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,908,119 | C→A | 41.9% | P449P (CCG→CCT) | recC ← | exonuclease V (RecBCD complex), gamma chain |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,908,119 | 0 | C | A | 41.9% | 33.0 / 53.4 | 62 | P449P (CCG→CCT) | recC | exonuclease V (RecBCD complex), gamma chain |
| Reads supporting (aligned to +/- strand): ref base C (13/23); new base A (16/10); total (29/33) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.11e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.05e-01 | |||||||||||
ACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTA > minE/1908057‑1908185 | acaTCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTCCAGTACCGGAtgt < 1:1797800/69‑1 (MQ=255) ccTCTGACACAAAACGGCTGTCAGGCAGTGATAACAGGCTGATAAACGCTTCCAGTACCGGATGTGACt > 1:1806252/1‑69 (MQ=255) cacaAAGCGGCTGTCTGGCAGCGATAACAGGCTGATAAACGCTTGAAGTACAGGAtgc > 1:150155/1‑57 (MQ=255) cacaAAGCGGCTGTCTGGCAGCGATAACAGGCTGATAAACGCTTCAAGTACAGGAtgc > 1:1005876/1‑57 (MQ=255) cacaAAGCGGCTGTCTGGCAGCGATAACAGGCTGATAAACGCTTCAAGTACAGGAtgc > 1:937102/1‑57 (MQ=255) cacaAAGCGGCTGTCTGGCAGCGATAACAGGCTGATAAACGCTTCAAGTACAGGAtgc > 1:924196/1‑57 (MQ=255) cacaAAGCGGCTGTCTGGCAGCGATAACAGGCTGATAAACGCTTCAAGTACAGGAtgc > 1:856731/1‑57 (MQ=255) cacaAAGCGGCTGTCTGGCAGCGATAACAGGCTGATAAACGCTTCAAGTACAGGAtgc > 1:562903/1‑57 (MQ=255) cacaAAGCGGCTGTCTGGCAGCGATAACAGGCTGATAAACGCTTCAAGTACAGGAtgc > 1:3394662/1‑57 (MQ=255) cacaAAGCGGCTGTCTGGCAGCGATAACAGGCTGATAAACGCTTCAAGTACAGGAtgc > 1:3382878/1‑57 (MQ=255) cacaAAGCGGCTGTCTGGCAGCGATAACAGGCTGATAAACGCTTCAAGTACAGGAtgc > 1:3312802/1‑57 (MQ=255) cacaAAGCGGCTGTCTGGCAGCGATAACAGGCTGATAAACGCTTCAAGTACAGGAtgc > 1:2796162/1‑57 (MQ=255) cacaAAGCGGCTGTCTGGCAGCGATAACAGGCTGATAAACGCTTCAAGTACAGGAtgc > 1:1884097/1‑57 (MQ=255) cacaAAGCGGCTGTCTGGCAGCGATAACAGGCTGATAAACGCTTCAAGTACAGGAtgc > 1:1485565/1‑57 (MQ=255) cacaAAGCGGCTGTCTGGCAGCGATAACAGGCTGATAAACGCTTCAAGTACAGGAtgc > 1:1470919/1‑57 (MQ=255) cacaAAGCGGCTGTCTGGCAGCGATAACAGGCTGATAAACGCTTCAAGTACAGGAtgc > 1:128012/1‑57 (MQ=255) cacaAAGCGGCTGTCTGGCAGCGATAACAGGCTGATAAACGCTTCAAGTACAGGAtgc > 1:1310669/1‑57 (MQ=255) cacaAAGCGGCTGTCTGGCAGCGATAACAGGCTGATAAACGCTTCAAGTACAGGAtgc > 1:1341921/1‑57 (MQ=255) caAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCac < 1:1068663/69‑1 (MQ=255) tGTCAGGCAGTGATAACAGGCTGATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCg > 1:1496157/1‑68 (MQ=255) tGTCAGGCAGTGATAACAGGCTGATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCg > 1:1528405/1‑68 (MQ=255) agcgATAACAGGCTGATAAACGCTTCAAGTACAGGATGCGACTGCCGCGCCCGACGGTCGGAAATg < 1:1005413/63‑1 (MQ=255) agcgATAACAGGCTGATAAACGCTTCAAGTACAGGATGCGACTGCCGCGCCCGACGGTCGGAAATg < 1:2666681/63‑1 (MQ=255) agcgATAACAGGCTGATAAACGCTTCAAGTACAGGATGCGACTGCCGCGCCCGACGGTCGGAAATg < 1:2826029/63‑1 (MQ=255) agcgATAACAGGCTGATAAACGCTTCAAGTACAGGATGCGACTGCCGCGCCCGACGGTCGGAAATg < 1:2936981/63‑1 (MQ=255) agcgATAACAGGCTGATAAACGCTTCAAGTACAGGATGCGACTGCCGCGCCCGACGGTCGGAAATg < 1:1226905/63‑1 (MQ=255) agcgATAACAGGCTGATAAACGCTTCAAGTACAGGATGCGACTGCCGCGCCCGACGGTCGGAAATg < 1:455133/63‑1 (MQ=255) agcgATAACAGGCTGATAAACGCTTCAAGTACAGGATGCGACTGCCGCGCCCGACGGTCGGAAATg < 1:1354482/63‑1 (MQ=255) agcgATAACAGGCTGATAAACGCTTCAAGTACAGGATGCGACTGCCGCGCCCGACGGTCGGAAATg < 1:1709483/63‑1 (MQ=255) agcgATAACAGGCTGATAAACGCTTCAAGTACAGGATGCGACTGCCGCGCCCGACGGTCGGAAATg < 1:1051941/63‑1 (MQ=255) agcgATAACAGGCTGATAAACGCTTCAAGTACAGGATGCGACTGCCCCGCCCGACGGTCGGAAATg < 1:2615080/63‑1 (MQ=255) gATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACg < 1:1710319/69‑1 (MQ=255) gATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACg < 1:492828/69‑1 (MQ=255) gATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACg < 1:2150344/69‑1 (MQ=255) gATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACg < 1:1089417/69‑1 (MQ=255) gATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACg < 1:1635877/69‑1 (MQ=255) gATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACg < 1:1255964/69‑1 (MQ=255) gATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACg < 1:76064/69‑1 (MQ=255) gATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACg < 1:1035944/69‑1 (MQ=255) tAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGaa > 1:3270436/1‑47 (MQ=255) tAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGAt > 1:2556589/1‑69 (MQ=255) tAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGAt > 1:1399658/1‑69 (MQ=255) tAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACCACGGTCGGAAATGGCGTAAGGTAGGTAACGAt > 1:2983654/1‑69 (MQ=255) aaaCGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATc < 1:2888054/69‑1 (MQ=255) aaaCGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATc < 1:2319893/69‑1 (MQ=255) aaaCGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATc < 1:162639/69‑1 (MQ=255) aaCGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATc < 1:2801869/68‑1 (MQ=255) aaCGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATc < 1:137411/68‑1 (MQ=255) aaCGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATc < 1:809567/68‑1 (MQ=255) ttCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAgg < 1:2565698/69‑1 (MQ=255) ttCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAgg < 1:2781366/69‑1 (MQ=255) ttCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAgg < 1:2444881/69‑1 (MQ=255) ttCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAgg < 1:1384795/69‑1 (MQ=255) ttCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAgg < 1:1678655/69‑1 (MQ=255) ttCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAgg < 1:2057057/69‑1 (MQ=255) tCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGGAAGGTAGGTAACGATCCGCAgg < 1:2076431/68‑1 (MQ=255) tACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTATGGTAGGTAACGATCCGCAGGTGCGCt > 1:1405575/1‑69 (MQ=255) tACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCt > 1:2233407/1‑69 (MQ=255) tACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCt > 1:2245814/1‑69 (MQ=255) tACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCt > 1:2292333/1‑69 (MQ=255) tACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCt > 1:2361799/1‑69 (MQ=255) tACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCt > 1:2418969/1‑69 (MQ=255) ccGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTa > 1:325243/1‑68 (MQ=255) | ACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTA > minE/1908057‑1908185 |
| 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: ‑ |