| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 1,068,568 | T→A | 100% | G22G (GGT→GGA) | ymdF → | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 1,068,568 | 0 | T | A | 76.9% | 21.3 / 3.4 | 13 | G22G (GGT→GGA) | ymdF | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base T (0/3); new base A (6/4); total (6/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.92e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
GCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATC > W3110S.gb/1068506‑1068624 | gCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAgg > 1:1098740/1‑71 (MQ=255)gCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAgg > 1:1158506/1‑71 (MQ=255)gCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAgg > 1:1182455/1‑71 (MQ=255)gCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAgg > 1:1673121/1‑71 (MQ=255)gCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAgg > 1:1952709/1‑71 (MQ=255)gCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAgg > 1:730400/1‑71 (MQ=255) cGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTa < 1:1685225/71‑1 (MQ=255) cGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTa < 1:414055/71‑1 (MQ=255) cGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTa < 1:761039/71‑1 (MQ=255) tGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAg < 1:1419091/70‑1 (MQ=255) gagaATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCa < 1:641708/71‑1 (MQ=255) gagCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCa < 1:1078912/71‑1 (MQ=255) gCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATc < 1:261071/71‑1 (MQ=255) | GCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATC > W3110S.gb/1068506‑1068624 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |