| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 1,569,965 | T→C | 100% | Q191Q (CAA→CAG) | yddW ← | predicted liprotein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 1,569,965 | 0 | T | C | 87.2% | 106.1 / 6.9 | 39 | Q191Q (CAA→CAG) | yddW | predicted liprotein |
| Reads supporting (aligned to +/- strand): ref base T (4/1); new base C (13/21); total (17/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.47e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
CTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGTACCGG > W3110S.gb/1569899‑1570026 | cTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGca < 1:964244/71‑1 (MQ=255)cTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGca < 1:2343839/71‑1 (MQ=255)cTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGca < 1:1782262/71‑1 (MQ=255)cTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGca > 1:1612080/1‑71 (MQ=255)cTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGca < 1:1595142/71‑1 (MQ=255) cAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACg > 1:1123858/1‑71 (MQ=255) gAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGgttgtt > 1:3531241/1‑68 (MQ=255) gAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGgttgtt > 1:2400153/1‑68 (MQ=255) gAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGgttgtt > 1:1858762/1‑68 (MQ=255) gAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTga < 1:490994/70‑1 (MQ=255) gAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTga < 1:1387658/70‑1 (MQ=255) aaaGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGc > 1:3219728/1‑40 (MQ=255) aaaGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACagag > 1:1017013/1‑71 (MQ=255) gCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTc < 1:1512100/70‑1 (MQ=255) gCCAGACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTc < 1:1112728/70‑1 (MQ=255) gACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTc < 1:834796/71‑1 (MQ=255) gACGTTCTGATCCAGTCGCGGTGCTGCACATAGCCGCTCGCCGGTTGTTgaga > 1:1500983/1‑53 (MQ=39) gACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAg > 1:1645000/1‑68 (MQ=255) gACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAg > 1:979562/1‑68 (MQ=255) gACGTTCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAg > 1:3029770/1‑68 (MQ=255) ttCTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTg < 1:3677671/71‑1 (MQ=255) cTGATCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCt > 1:2958128/1‑68 (MQ=255) tCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGTa < 1:747051/71‑1 (MQ=255) tCCAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGTa < 1:667659/71‑1 (MQ=255) ccAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACaga > 1:639032/1‑46 (MQ=255) ccAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACaga > 1:1144869/1‑46 (MQ=255) ccAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACaga > 1:2636788/1‑46 (MQ=255) ccAGTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATa > 1:1388870/1‑67 (MQ=255) gtagcgGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCt < 1:2253702/57‑1 (MQ=255) gTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCt < 1:3105277/60‑1 (MQ=255) gTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCt < 1:2817970/60‑1 (MQ=255) gTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCt < 1:135976/60‑1 (MQ=255) gTCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCt < 1:1117494/60‑1 (MQ=255) tCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGc < 1:3147540/46‑1 (MQ=255) tCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGc < 1:2221648/46‑1 (MQ=255) tCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGc < 1:1489841/46‑1 (MQ=255) tCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGa < 1:2297974/61‑1 (MQ=255) tCGCGGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGa < 1:218729/61‑1 (MQ=255) gcgGTGCTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTa > 1:1022513/1‑46 (MQ=255) gtgcTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGTACCgg > 1:155600/1‑65 (MQ=255) | CTGAACCTCAGGGATGCCCGGGTCGAGGACAAAGCGATCGCCAGACGTTCTGATCCAGTCGCGGTGTTGCACATAGACGCTCGCCGGTTGTTGAGACAGAGTGCTATTCAGTTCCCTGATAGTACCGG > W3110S.gb/1569899‑1570026 |
| 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: ‑ |