| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 1,533,715 | T→C | 100% | G62G (GGT→GGC) | ydcC → | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 1,533,715 | 0 | T | C | 79.2% | 41.0 / 8.0 | 24 | G62G (GGT→GGC) | ydcC | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base T (4/1); new base C (11/8); total (15/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.15e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
TTGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGACACCATTGCCAGAGTTGTATCCTGTATCAGTCCTGCAA > W3110S.gb/1533678‑1533782 | ttGGGGAAACACATCCCGATTTTTTGAAGCAATATGGCGATTTTGAAAATGGTATTCCTGTTCACGacac > 1:159947/1‑70 (MQ=25)ttGGGGAAACACATCCCGATTTTTTGAAGCAATATGGCGATTTTGAAAATGGTATTCCTGTTCACGacac < 1:798305/70‑1 (MQ=25) tGGGGAAACACATCCCGATTTTTTGAAGCAATATGGCGATTTTGAAAATGGTATTCCTGTTCACGACAcc < 1:79533/70‑1 (MQ=25) tGGGGAAACACATCCCGATTTTTTGAAGCAATATGGCGATTTTGAAAATGGTATTCCTGTTCACGACAcc < 1:92163/70‑1 (MQ=25) aaCACATCCCGATTTTTTGAAGCAATATGGCGATTTTGAAAATGGTATTCCTGTTCACGaca > 1:1903484/1‑62 (MQ=18) acaTCCCGATTTTTTGAAGCAATATGGCGATTTTGAAAATGGTATTCCTGTTCACGACACCATTGCCagag > 1:1879461/1‑71 (MQ=255) tCCCGATTTTTTGAAGCAATATGGCGATTTTGAAAATGGTATTCCTGTTCACGACACCATTGCCAGAGTTg > 1:1752215/1‑71 (MQ=255) tCCCGATTTTTTGAAGCAATATGGCGATTTTGAAAATGGTATTCCTGTTCACGACACCATTGCCAGAGTTg > 1:1905834/1‑71 (MQ=255) ccGATTTTTTGAAGCAATATGGCGATTTTGAAAATGGTATTCCTGTTCACGACACCATTGCCAGAGTTGTa > 1:2810312/1‑71 (MQ=255) ttttttGAATCAATATGGCTATTTTGAAAATGGTATTCCTGTTCACGacac < 1:1953008/51‑1 (MQ=11) tttttGAAGCAATATGGCGATTTTGAAAATGGTATTCCTGTTCACGACACCATTGCCAGAGTTGTATCCTg > 1:2715736/1‑71 (MQ=255) ttttGAAGCAATATGGCGATTTTGAAAATGGTATTCCTGTTCACGACACCATTGCCAGAGTTGTATCCTg < 1:958296/70‑1 (MQ=255) ttttGAAGCAATATGGCGATTTTGAAAATGGTATTCCTGTTCACGACACCATTGCCAGAGTTGTATCCTg < 1:1320711/70‑1 (MQ=255) gCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGACACCATTGCCAGAGTTGTATCCTGTATCAGTc > 1:171729/1‑71 (MQ=34) gCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGACACCATTGCCAGAGTTGTATCCTGTATCAGTc > 1:377641/1‑71 (MQ=34) tatGGTGATTTTGAAAATGGTATTCCTGTTCACGACACCATTGCCAGAGTTGTATCCTGTa < 1:1998779/61‑1 (MQ=35) tatGGCGATTTTGAAAATGGTATTCCTGTTCACGACACCATTGCCAGAGTTGTATCCTGTATc < 1:104880/63‑1 (MQ=255) atGGCGATTTTGAAAATGGTATTCCTGTTCACGACACCa > 1:296688/1‑39 (MQ=255) atGGCGATTTTGAAAATGGTATTCCTGTTCACGACACCa > 1:747095/1‑39 (MQ=255) tGGTGATTTTGAAAATGGTATTCCTGTTCACGACACCATTGCCAGAGTTGTATCCTGTATCAGTCCTGCa > 1:2755764/1‑70 (MQ=35) tGGTGATTTTGAAAATGGTATTCCTGTTCACGACACCATTGCCAGAGTTGTATCCTGTATCAGTCCTGCa > 1:693193/1‑70 (MQ=35) tGGCGATTTTGAAAATGGTATTCCTGTTCACGACACCa > 1:719832/1‑38 (MQ=255) tGGCGATTTTGAAAATGGTATTCCTGTTCACGACACCATTGCCAGAGTTGTATCCTGTa > 1:682165/1‑59 (MQ=255) tGGCGATTTTGAAAATGGTATTCCTGTTCACGACACCATTGCCAGAGTTGTATCCTGTATCAGTCCTGCaa < 1:1318931/71‑1 (MQ=21) | TTGGGGAAACACATCCCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGACACCATTGCCAGAGTTGTATCCTGTATCAGTCCTGCAA > W3110S.gb/1533678‑1533782 |
| 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: ‑ |