| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 4,235,478 | T→C | 31.8% | intergenic (+71/‑428) | pgi → / → yjbE | glucosephosphate isomerase/extracellular polysaccharide production threonine‑rich protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,235,478 | 0 | T | C | 31.8% | 32.2 / 16.1 | 22 | intergenic (+71/‑428) | pgi/yjbE | glucosephosphate isomerase/extracellular polysaccharide production threonine‑rich protein |
| Reads supporting (aligned to +/- strand): ref base T (7/8); new base C (3/4); total (10/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01 | |||||||||||
GCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTA > NC_000913/4235349‑4235564 | gcagcCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTa < 1:204687/148‑1 (MQ=255) aTTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAt > 4:1115850/1‑142 (MQ=255) aTTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGACGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACaa < 3:1115850/142‑2 (MQ=255) tAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc < 4:37133/114‑1 (MQ=255) tAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc > 3:37133/1‑114 (MQ=255) gcgcGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat > 4:1021192/1‑114 (MQ=255) gcgcGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat < 3:1021192/114‑1 (MQ=255) cgcgGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTAc < 1:1000266/119‑1 (MQ=255) cgcgGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTAc > 2:1000266/1‑119 (MQ=255) ggCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt > 1:126864/1‑103 (MQ=255) ggCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAAt < 2:126864/103‑1 (MQ=255) cGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTCGCATCAg < 3:185247/56‑1 (MQ=11) cGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTCGCATCAg > 4:185247/1‑56 (MQ=2) ttCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAt > 2:830067/1‑51 (MQ=255) ttCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTAt < 1:830067/51‑1 (MQ=255) gtaggCCTGATAAGACGCGACCGCGTCGCATCAGGCATACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTta < 3:1040467/103‑1 (MQ=21) gtaggCCTGATAAGACGCGACCGCGTCGCATCAGGCATACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTta > 4:1040467/5‑107 (MQ=21) ataagACGCGACCGCGTCGCATCAGGCATACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAAcc < 3:232393/91‑1 (MQ=37) aagACGCGACCGCGTCGCATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgattga < 2:548857/81‑1 (MQ=38) aagACGCGACCGCGTCGCATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACgattga > 1:548857/3‑83 (MQ=37) cGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAAcc > 3:550911/1‑92 (MQ=255) cGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAAcc < 4:550911/92‑1 (MQ=255) | GCAGCCACGATAGCTCGACCAATGGTCTGATTAACCGCTATAAAGCGTGGCGCGGTTAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTA > NC_000913/4235349‑4235564 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |