| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,908,113 | C→A | 46.4% | L451L (CTG→CTT) | 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,113 | 0 | C | A | 46.4% | 16.5 / 54.8 | 56 | L451L (CTG→CTT) | recC | exonuclease V (RecBCD complex), gamma chain |
| Reads supporting (aligned to +/- strand): ref base C (7/23); new base A (16/10); total (23/33) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.16e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01 | |||||||||||
GCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGG > minE/1908046‑1908178 | gCAACGCCAGCACATCCTCTGACACAACACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTcc < 1:448177/68‑1 (MQ=255) 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:128012/1‑57 (MQ=255) cacaAAGCGGCTGTCTGGCAGCGATAACAGGCTGATAAACGCTTCAAGTACAGGAtgc > 1:1310669/1‑57 (MQ=255) cacaAAGCGGCTGTCTGGCAGCGATAACAGGCTGATAAACGCTTCAAGTACAGGAtgc > 1:1341921/1‑57 (MQ=255) cacaAAGCGGCTGTCTGGCAGCGATAACAGGCTGATAAACGCTTCAAGTACAGGAtgc > 1:1470919/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:2666681/63‑1 (MQ=255) agcgATAACAGGCTGATAAACGCTTCAAGTACAGGATGCGACTGCCGCGCCCGACGGTCGGAAATg < 1:1005413/63‑1 (MQ=255) agcgATAACAGGCTGATAAACGCTTCAAGTACAGGATGCGACTGCCGCGCCCGACGGTCGGAAATg < 1:1051941/63‑1 (MQ=255) agcgATAACAGGCTGATAAACGCTTCAAGTACAGGATGCGACTGCCGCGCCCGACGGTCGGAAATg < 1:455133/63‑1 (MQ=255) agcgATAACAGGCTGATAAACGCTTCAAGTACAGGATGCGACTGCCGCGCCCGACGGTCGGAAATg < 1:1226905/63‑1 (MQ=255) agcgATAACAGGCTGATAAACGCTTCAAGTACAGGATGCGACTGCCGCGCCCGACGGTCGGAAATg < 1:2936981/63‑1 (MQ=255) agcgATAACAGGCTGATAAACGCTTCAAGTACAGGATGCGACTGCCGCGCCCGACGGTCGGAAATg < 1:2826029/63‑1 (MQ=255) agcgATAACAGGCTGATAAACGCTTCAAGTACAGGATGCGACTGCCGCGCCCGACGGTCGGAAATg < 1:1354482/63‑1 (MQ=255) agcgATAACAGGCTGATAAACGCTTCAAGTACAGGATGCGACTGCCGCGCCCGACGGTCGGAAATg < 1:1709483/63‑1 (MQ=255) agcgATAACAGGCTGATAAACGCTTCAAGTACAGGATGCGACTGCCCCGCCCGACGGTCGGAAATg < 1:2615080/63‑1 (MQ=255) gATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACg < 1:2150344/69‑1 (MQ=255) gATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACg < 1:1035944/69‑1 (MQ=255) gATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACg < 1:76064/69‑1 (MQ=255) gATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACg < 1:492828/69‑1 (MQ=255) gATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACg < 1:1635877/69‑1 (MQ=255) gATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACg < 1:1089417/69‑1 (MQ=255) gATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACg < 1:1255964/69‑1 (MQ=255) gATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACg < 1:1710319/69‑1 (MQ=255) tAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGaa > 1:3270436/1‑47 (MQ=255) tAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGAt > 1:1399658/1‑69 (MQ=255) tAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGAt > 1:2556589/1‑69 (MQ=255) tAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACCACGGTCGGAAATGGCGTAAGGTAGGTAACGAt > 1:2983654/1‑69 (MQ=255) aaaCGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATc < 1:2888054/69‑1 (MQ=255) aaaCGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATc < 1:162639/69‑1 (MQ=255) aaaCGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATc < 1:2319893/69‑1 (MQ=255) aaCGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATc < 1:137411/68‑1 (MQ=255) aaCGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATc < 1:2801869/68‑1 (MQ=255) aaCGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATc < 1:809567/68‑1 (MQ=255) ttCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAgg < 1:2781366/69‑1 (MQ=255) ttCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAgg < 1:1678655/69‑1 (MQ=255) ttCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAgg < 1:2057057/69‑1 (MQ=255) ttCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAgg < 1:2565698/69‑1 (MQ=255) ttCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAgg < 1:2444881/69‑1 (MQ=255) ttCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAgg < 1:1384795/69‑1 (MQ=255) tCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGGAAGGTAGGTAACGATCCGCAgg < 1:2076431/68‑1 (MQ=255) | GCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTCCAGTACCGGATGTGACTGCCGCGCACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGG > minE/1908046‑1908178 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |