| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 4,591,445 | G→A | 42.1% | T18I (ACA→ATA) | hsdR ← | endonuclease R |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 4,591,445 | 0 | G | A | 42.1% | 6.7 / 19.2 | 19 | T18I (ACA→ATA) | hsdR | endonuclease R |
| Reads supporting (aligned to +/- strand): ref base G (7/4); new base A (2/6); total (9/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.70e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01 | |||||||||||
ACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTT > W3110S.gb/4591378‑4591510 | acaGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGt > 1:1073393/1‑69 (MQ=255) cGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAAt > 1:1103142/1‑71 (MQ=255) tGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTg > 1:1038623/1‑60 (MQ=255) gCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGtta > 1:280840/1‑61 (MQ=255) gCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCAtt < 1:1645816/71‑1 (MQ=255) cATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTg < 1:1623151/71‑1 (MQ=255) aGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCtgt > 1:91328/1‑71 (MQ=255) gTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGttat > 1:962548/1‑56 (MQ=255) gTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCtgtg < 1:355519/71‑1 (MQ=255) tCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATc < 1:1639864/55‑1 (MQ=255) cGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATc < 1:1255460/54‑1 (MQ=255) gACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTGTc < 1:1515387/70‑1 (MQ=255) aCGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCgg < 1:187440/62‑1 (MQ=255) cccTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGc > 1:1321161/1‑57 (MQ=255) aTTCAAAATTGGATTTATTCATCATTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGcc > 1:1548090/1‑71 (MQ=255) cAAAATTGGATTTATTCATAATTATTATTAATCCATTg > 1:384659/1‑38 (MQ=25) ttATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAcatctcatc > 1:285619/1‑71 (MQ=255) aTTCATCATTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAcatctcatct < 1:1348997/70‑1 (MQ=255) aTTCATCATTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAAcatctcatct < 1:1346919/70‑1 (MQ=255) tcatTATTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCtt < 1:797301/71‑1 (MQ=255) | ACAGGCGATGGCATAAGTGAAGTCGTTGACGCCCTTCAGGAATTCAAAATTGGATTTATTCATCATTGTTATTAATCCATTGCTGTGCGGGCCTGTCCAAATATTTAAGGCCCATAACATCTCATCTTAGCTT > W3110S.gb/4591378‑4591510 |
| 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: ‑ |