Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 1,430,136 | C→T | 23.8% | pseudogene (965/1527 nt) R102H (CGC→CAC) |
lomR → insH ← |
ECK1366:JW5884+JW5904:b4570; Rac prophage region; hypothetical protein IS5 transposase and trans‑activator |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 1,430,136 | 0 | C | T | 23.8% | 78.8 / 22.5 | 42 | pseudogene (965/1527 nt) R102H (CGC→CAC) | lomR insH | ECK1366:JW5884+JW5904:b4570; Rac prophage region; hypothetical protein IS5 transposase and trans‑activator |
Reads supporting (aligned to +/- strand): ref base C (10/22); new base T (8/2); total (18/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.04e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01 |
GGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTA > W3110S.gb/1430069‑1430197 | ggTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGAt > 1:3702095/1‑71 (MQ=255) ttGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCgg > 1:1699079/1‑71 (MQ=255) gAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGc < 1:1456991/70‑1 (MQ=255) gAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGc < 1:3732947/70‑1 (MQ=255) aTTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGgcg < 1:710188/71‑1 (MQ=255) ttGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGa < 1:1852234/58‑1 (MQ=255) tGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGa < 1:1869542/57‑1 (MQ=255) gggCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGTGATCCGGCAGGGCGCTATcc < 1:3847468/71‑1 (MQ=255) gggCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATcc < 1:826829/71‑1 (MQ=255) ttGATGCTGCTCGAGGAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATcc < 1:417210/64‑1 (MQ=255) ttGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGTGATCCGGCAGGGCGCTATCCa > 1:4112795/1‑65 (MQ=255) ttGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGTGATCCGGCAGGGCGCTATCCa > 1:1569449/1‑65 (MQ=255) tGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCg > 1:5106398/1‑49 (MQ=255) tGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCTGGGCGCTATCCAGGGa > 1:296196/1‑68 (MQ=255) tGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATaa > 1:363701/1‑71 (MQ=255) gATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCt > 1:1395635/1‑58 (MQ=255) gATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAAt > 1:1197852/1‑71 (MQ=255) aTGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATc < 1:126602/71‑1 (MQ=255) gctCGAGCAGGTGGCGGAAATTCATGATGGTGGTGTGATCCGGCAGGGCGCTATCCAGGGa > 1:1956090/1‑61 (MQ=255) gctCGAGCAGGTGGCGGAAATTCATGATGGTGGTGTGATCCGGCAGGGCGCTATCCAGGGAt > 1:3750387/1‑62 (MQ=255) ctCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCa < 1:3809067/70‑1 (MQ=255) gAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCaaa < 1:1187399/69‑1 (MQ=255) gAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCaaa < 1:749821/69‑1 (MQ=255) gAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCaaa < 1:2499032/69‑1 (MQ=255) gAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCaaa < 1:4130611/69‑1 (MQ=255) tGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCa < 1:147786/70‑1 (MQ=255) ggCGGAAATTCATGATGGTGGTGTGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATg > 1:2922513/1‑71 (MQ=255) gCGGAAATTCATGATGGTGGTGTGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATg > 1:3984114/1‑70 (MQ=255) gCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCa < 1:2287258/68‑1 (MQ=255) cGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATgg < 1:3571987/70‑1 (MQ=255) aaaTTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATg < 1:1532003/66‑1 (MQ=255) aTTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGa > 1:4114107/1‑71 (MQ=255) cATGATGGTGGTGTGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAgg > 1:2240305/1‑55 (MQ=255) cATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGAttt > 1:2445668/1‑71 (MQ=255) tgatgGTGGTGTGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCg < 1:4818995/71‑1 (MQ=255) tgatgGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCg < 1:2483508/71‑1 (MQ=255) tggtggTGTGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGAATGGAGGCGATTTCGTa > 1:677199/1‑70 (MQ=255) gtgtGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCg < 1:2650496/63‑1 (MQ=255) gtgCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGc < 1:4069969/56‑1 (MQ=255) gtgCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGc < 1:4383714/56‑1 (MQ=255) tgCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGgcgc < 1:1265850/47‑1 (MQ=255) tgCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGgcgc < 1:2787648/47‑1 (MQ=255) gCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGa > 1:3004385/1‑56 (MQ=255) | GGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTA > W3110S.gb/1430069‑1430197 |
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: ‑ |