| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 579,680 | G→A | 100% | intergenic (+371/‑377) | ybcW → / → nohB | hypothetical protein/DNA packaging protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 579,680 | 0 | G | A | 80.7% | 53.9 / 10.9 | 31 | intergenic (+371/‑377) | ybcW/nohB | hypothetical protein/DNA packaging protein |
| Reads supporting (aligned to +/- strand): ref base G (1/5); new base A (15/10); total (16/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.29e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
AAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACG > W3110S.gb/579618‑579723 | aaaGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCa > 1:927234/1‑69 (MQ=25)aaaGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCa > 1:3122373/1‑69 (MQ=25)aaaGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCa < 1:1669934/69‑1 (MQ=25)aaaGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCa > 1:2119847/1‑69 (MQ=25) aaGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAg > 1:1194607/1‑69 (MQ=25) gCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCTCATATTCAgtg > 1:2518159/1‑69 (MQ=255) ggTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTg < 1:2629773/68‑1 (MQ=255) tgttgttCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTAtt < 1:3485715/69‑1 (MQ=255) gttgttCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTg > 1:2513837/1‑69 (MQ=35) cGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTAtt < 1:2092519/59‑1 (MQ=255) tctgttctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGtttt < 1:1439732/69‑1 (MQ=25) tctggtctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGtttt > 1:2103562/1‑69 (MQ=255) tctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACg > 1:143039/1‑68 (MQ=25) ctTAGTTGACATGAGGTTGGCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAgtt > 1:289154/1‑58 (MQ=11) ctTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTa < 1:2147377/65‑1 (MQ=35) ctTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTa < 1:861329/65‑1 (MQ=35) ctTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTa < 1:53148/65‑1 (MQ=35) ctTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTa < 1:49667/65‑1 (MQ=35) ctTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTa < 1:3082181/65‑1 (MQ=35) ctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAgtt > 1:2544834/1‑58 (MQ=255) ctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAg < 1:324843/56‑1 (MQ=255) ctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAg < 1:2242115/56‑1 (MQ=255) ctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAg > 1:1504895/1‑56 (MQ=255) ctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGtttt > 1:3592922/1‑63 (MQ=25) ctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGtttt > 1:176160/1‑63 (MQ=25) ctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTAc < 1:2368068/66‑1 (MQ=25) ctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTAc < 1:391052/66‑1 (MQ=25) ctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTAc > 1:2867563/1‑66 (MQ=25) ctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTAc < 1:2834372/66‑1 (MQ=25) ctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACg > 1:378792/1‑67 (MQ=25) ctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACg > 1:539132/1‑67 (MQ=25) | AAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACG > W3110S.gb/579618‑579723 |
| 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: ‑ |