| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 2,885,984 | C→T | 100% | intergenic (‑109/+250) | ygcB ← / ← cysH | conserved hypothetical protein, member of DEAD box family/3'‑phosphoadenosine 5'‑phosphosulfate reductase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 2,885,984 | 0 | C | T | 79.4% | 60.1 / 9.6 | 29 | intergenic (‑109/+250) | ygcB/cysH | conserved hypothetical protein, member of DEAD box family/3'‑phosphoadenosine 5'‑phosphosulfate reductase |
| Reads supporting (aligned to +/- strand): ref base C (1/5); new base T (17/6); total (18/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.85e-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. | |||||||||||
GCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGG > W3110S.gb/2885920‑2886054 | gCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAaagaa < 1:1476391/71‑1 (MQ=255) tttttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAagaaga < 1:1907785/69‑1 (MQ=255) tttttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAagaaga < 1:1380039/69‑1 (MQ=255) tttttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAagaaga < 1:1041806/69‑1 (MQ=255) tttttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAagaaga < 1:557147/69‑1 (MQ=255) tttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGAc > 1:803459/1‑70 (MQ=255) cATAAATATCGAATTATAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAAtt > 1:567344/1‑70 (MQ=255) cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTg > 1:735244/1‑71 (MQ=255) cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTg > 1:50087/1‑71 (MQ=255) cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTg > 1:502630/1‑71 (MQ=255) cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTg > 1:1878665/1‑71 (MQ=255) cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTg > 1:1563392/1‑71 (MQ=255) cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTg > 1:1562741/1‑71 (MQ=255) cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTg > 1:1540390/1‑71 (MQ=255) cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTg > 1:1402895/1‑71 (MQ=255) cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTg > 1:1312525/1‑71 (MQ=255) cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTg > 1:1123587/1‑71 (MQ=255) aaTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTg > 1:773796/1‑71 (MQ=255) aaTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTg > 1:927973/1‑71 (MQ=255) aaTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAGTTGTTTg > 1:1872525/1‑71 (MQ=255) ttACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAAc < 1:187115/71‑1 (MQ=255) ttACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAAc < 1:1559926/71‑1 (MQ=255) ttACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAAc < 1:1091615/71‑1 (MQ=255) agaGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGc > 1:1633181/3‑71 (MQ=255) accgTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACg > 1:841586/4‑70 (MQ=255) aCAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGc > 1:1063408/1‑71 (MQ=255) tccTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGaaa < 1:22401/69‑1 (MQ=255) gCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGaaa < 1:1444001/70‑1 (MQ=255) gCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGaaa < 1:39299/70‑1 (MQ=255) cAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgg > 1:65098/1‑71 (MQ=255) | GCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGG > W3110S.gb/2885920‑2886054 |
| 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: ‑ |