| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 1,472,567 | C→G | pseudogene (2049/3495 nt) | ydbA → | putative outer membrane protein, N‑terminal fragment |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,472,567 | 0 | C | G | 97.0% | 116.3 / ‑3.1 | 33 | pseudogene (2049/3495 nt) | ydbA | putative outer membrane protein, N‑terminal fragment |
| Reads supporting (aligned to +/- strand): ref base C (1/0); new base G (17/15); total (18/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01 | |||||||||||
AGGTCATGTGGTTAACAACGGCACGGTGGTGATTGCCGATGGGGTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGCATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTATGGCGACTATACGTTGCCGGATGT > NC_000913/1472489‑1472644 | aGGTCATGTGGTTAACAACGGCACGGTGGTGATTGCCGATGGGGTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGGAt > 2:258206/1‑81 (MQ=255) cATGTGGTTAACAACGGCACGGTGGTGATTGCCGATGGGGTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGGATCAAt > 2:258199/1‑81 (MQ=255) cATGTGGTTAACAACGGCACGGTGGTGATTGCCGATGGGGTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGGATCAAt > 1:258209/1‑81 (MQ=255) aTGTGGTTAACAACGGCACGGTGGTGATTGCCGATGGGGTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGGATCAAt < 1:258172/80‑1 (MQ=255) tgGTTAACAACGGCACGGTGGTGATTGCCGATGGGGTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGGATCAATGTTg < 2:258173/81‑1 (MQ=255) gTTAACAACGGCACGGTGGTGATTGCCGATGGGGTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGGATCAATGTTGa < 1:258174/80‑1 (MQ=255) ttAACAACGGCACGGTGGTGATTGCCGATGGGGTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGGATCAATGTTGAAg < 1:258175/81‑1 (MQ=255) aacGGCACGGTGGTGATTGCCGATGGGGTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGGATCAATGTTGAAGGTAt > 2:258214/1‑80 (MQ=255) aacGGCACGGTGGTGATTGCCGATGGGGTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGGATCAATGTTGAAGGTAt > 2:258204/1‑80 (MQ=255) acGGCACGGTGGTGATTGCCGATGGGGTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGGATCAATGTTGAAGGTATg > 2:258202/1‑80 (MQ=255) acGGCACGGTGGTGATTGCCGATGGGGTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGGATCAATGTTGAAGGTATg > 1:258225/1‑80 (MQ=255) tgATTGCCGATGGGGTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGGATCAATGTTGAAGGTATGAACGGTAACAACg < 2:258176/81‑1 (MQ=255) aTTGCCGATGGGGTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGGATCAATGTTGAAGGTATGAACGGTAACAACGGt < 1:258177/81‑1 (MQ=255) ccGATGGGGTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGGATCAATGTTGAAGGTATGAACGGTAACAACGGTAAt < 2:258178/80‑1 (MQ=255) cGATGGGGTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGGATCAATGTTGAAGGTATGAACGGTAACAACGGTAATag > 1:258211/1‑81 (MQ=255) aTGGGGTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGGATCAATGTTGAAGGTATGAACGGTAACAACGGTAATag > 2:258210/1‑79 (MQ=255) aTGGGGTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGGATCAATGTTGAAGGTATGAACGGTAACAACGGTAATag > 1:258215/1‑79 (MQ=255) gggTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGGATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCg < 1:258179/81‑1 (MQ=255) gTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGGATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGaa < 1:258180/81‑1 (MQ=255) gTTCTGGACTGATCAAGCAGGGCGACAGGATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCAtt < 1:258181/81‑1 (MQ=255) gTTCTGGACTGATCAAGCAGGGCGACAGGATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCAt > 1:258223/1‑80 (MQ=255) ttCTGGACTGATCAAGCAGGGCGACAGGATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTa < 1:258182/81‑1 (MQ=255) cTGGACTGATCAAGCAGGGCGACAGGATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTAt < 2:258183/80‑1 (MQ=255) ggACTGATCAAGCAGGGCGACAGGATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTATGGc > 2:258216/1‑81 (MQ=255) gACTGATCAAGCAGGGCGACAGGATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTATGGc > 2:258220/1‑80 (MQ=255) aaGCAGGGCGACAGCATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTATGGCGACTATACg > 2:258231/1‑81 (MQ=255) cAGGGCGACAGGATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTATGGCGACTATACGTTg < 2:258184/81‑1 (MQ=255) aGGGCGACAGGATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTATGGCGACTATACGTTGc > 1:258217/1‑81 (MQ=255) aGGGCGACAGGATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTATGGCGACTATACGTTGc > 2:258222/1‑81 (MQ=255) ggCGACAGGATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTATGGCGACTATACGTTGCCg > 1:258235/1‑81 (MQ=255) gCGACAGGATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTATGGCGACTATACGTTGCCgg < 2:258185/81‑1 (MQ=255) gACAGGATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTATGGCGACTATACGTTGCCGGAt < 1:258186/81‑1 (MQ=255) cAGGATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTATGGCGACTATACGTTGCCGGAtgt > 2:258221/1‑81 (MQ=255) | AGGTCATGTGGTTAACAACGGCACGGTGGTGATTGCCGATGGGGTTACGGGTTCTGGACTGATCAAGCAGGGCGACAGCATCAATGTTGAAGGTATGAACGGTAACAACGGTAATAGCAGCGAAGTGCATTATGGCGACTATACGTTGCCGGATGT > NC_000913/1472489‑1472644 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |