| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 2,960,928 | T→C | 51.3% | K53E (AAG→GAG) | recC ← | exonuclease V (RecBCD complex), gamma chain |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 2,960,928 | 0 | T | C | 51.3% | ‑2.3 / 51.3 | 37 | K53E (AAG→GAG) | recC | exonuclease V (RecBCD complex), gamma chain |
| Reads supporting (aligned to +/- strand): ref base T (12/6); new base C (6/13); total (18/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.02e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01 | |||||||||||
CCCGGACGAACATATCCCAGATAAAGCTCGCTGGCAGCGGAAAATCAATGTTTGCCGCAATACCAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGT > W3110S.gb/2960859‑2960992 | cccGGACGAACATATCCAAGATAAAGCTCGCTGGCAGCGGAAAATCAATGTTTGCCGCAATACCAAACttt < 1:1608894/71‑1 (MQ=255) ggACGAACATATCCCAGATAAAGCTCGCTGGCAGCGGAAAATCAATGTTTGCCGCAATACCAAACTCTTGc < 1:27947/71‑1 (MQ=255) aCGAACATATCCCAGATAAAGCTCGCTGGCAGCGGAAAATCAATGTTTGCCTCAATACCAAACTTTTGCGa > 1:1800762/1‑71 (MQ=255) cGAACATATCCCAGATAAAGCTCGCTGGCAGCGGAAAATCAATGTTTGCCGCAATACCAAACTTTTGCGAc < 1:2633715/71‑1 (MQ=255) cGAACATATCCCAGATAAAGCTCGCTGGCAGCGGAAAATCAATGTTTGCCGCAATACCAAACTCTTGCGa > 1:1102678/1‑70 (MQ=255) atatCCCAGATAAAGCTCGCTGGCAGCGGAAAATCAATGTTTGCCGCAATACCAAACTCTTGCGACAGGGt < 1:844326/71‑1 (MQ=255) ccAGATAAAGCTCGCTGGCAGCGGAAAATCAATGTTTGCCGCAATACCAAACTTTTGCGACAGGGTCAtt < 1:996077/70‑1 (MQ=255) ccAGATAAAGCTCGCTGGCAGCGGAAAATCAATGTTTGCCGCAATACCAAACTTTTGCGACAGGGTCAtt < 1:1801721/70‑1 (MQ=255) ccAGATAAAGCTCGCTGGCAGCGGAAAATCAATGTTTGCCGCAATACCAAACTCTTGCGACAGGGTCa > 1:1401368/1‑68 (MQ=255) gATAAAGCTCGCTGGCAGCGGAAAATCAATGTTTGCCGCAATACCAAACTCTTGCGACAGGGTCATTTGCa < 1:1933959/71‑1 (MQ=255) tAAAGCTCGCTGGCAGCGGAAAATCAATGTTTGCCGCAATACCAAACTTTTGCGACAGGGTCATTTGCAGc < 1:481224/71‑1 (MQ=255) tAAAGCTCGCTGGCAGCGGAAAATCAATGTTTGCCGCAATACCAAACTCTTGCGACAGGGTCATTTGCAGc < 1:934617/71‑1 (MQ=255) aGCTCGCTGGCAGCGGAAAATCAATGTTTGCCGCAATACCAAACTCTTGCGACAGGGTCATTTGCAGCCAc < 1:266245/71‑1 (MQ=255) tCGCTGGCAGCGGAAAATCAATGTTTGCCGCAATACCAAACTTTTGCGACAGGGTCATTTGCAGCCACtgt > 1:1954931/1‑71 (MQ=255) cTGGCAGCGGAAAATCAATGTTTGCCGCAATACCAAACTCTTGCGACAGGGTCATTTGCAGCCACTGTGc < 1:600570/70‑1 (MQ=255) tGGCAGCGGAAAATCAATGTTTGCCGCAATACCAAACTCTTGCGACAGGGTc < 1:798878/52‑1 (MQ=255) gCAGCGGAAAATCAATGTTTGCCGCAATACCAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCAt < 1:1900020/70‑1 (MQ=255) aGCGGAAAATCAATGTTTGCCGCAATACCAAACTCTTGCGACAGGGTCATTTGCAGCCACTGTGCCATAcc < 1:1583750/71‑1 (MQ=255) gAAAATCAATGTTTGCCGCAATACCAAACTCTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTa > 1:600028/1‑71 (MQ=255) gAAAATCAATGTTTGCCGCAATACCAAACTCTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTa < 1:1400787/71‑1 (MQ=255) aaaaTCAATGTTTGCCGCAATACCAAACTTTTGCGAc < 1:2419168/37‑1 (MQ=255) aaTCAATGTTTGCCGCAATACCAAACTCTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTAc < 1:1534077/69‑1 (MQ=255) aTCAATGTTTGCCGCAATACCAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACttt > 1:1612038/1‑71 (MQ=255) tttGCCGCAATACCAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCAc > 1:1865433/1‑68 (MQ=255) ccGCAATACCAAACTCTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAat < 1:467259/70‑1 (MQ=255) tACCAAACTCTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCAcc < 1:1122311/59‑1 (MQ=255) tACCAAACTCTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGa < 1:1713602/62‑1 (MQ=255) aCCAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCa > 1:2078420/1‑41 (MQ=255) aCCAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCa > 1:628749/1‑41 (MQ=255) aCCAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTg > 1:1377772/1‑71 (MQ=255) aCCAAACTCTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTg > 1:15276/1‑71 (MQ=255) ccAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCa > 1:1557167/1‑55 (MQ=255) ccAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCAc > 1:2349155/1‑56 (MQ=255) ccAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTgg > 1:276621/1‑71 (MQ=255) ccAAACTCTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTgg > 1:547544/1‑71 (MQ=255) cAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCa > 1:1707569/1‑39 (MQ=255) cAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCa > 1:798222/1‑39 (MQ=255) aaaCTCTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGt > 1:341550/1‑70 (MQ=255) tttGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAatcat > 1:2144688/1‑58 (MQ=255) tttGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAatcat > 1:1388876/1‑58 (MQ=255) | CCCGGACGAACATATCCCAGATAAAGCTCGCTGGCAGCGGAAAATCAATGTTTGCCGCAATACCAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGT > W3110S.gb/2960859‑2960992 |
| 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: ‑ |