| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 2,173,364 | Δ1 bp | intergenic (‑491/+918) | gatD ← / ← gatB | galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,173,364 | 0 | C | . | 80.0% | 46.1 / 11.7 | 20 | intergenic (‑491/+918) | gatD/gatB | galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component |
| Reads supporting (aligned to +/- strand): ref base C (0/4); new base . (10/6); total (10/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.67e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01 | |||||||||||
GGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGG > NC_000913/2173139‑2173609 | gggTGTGCAAACCAAAATTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACCCCCGAGATTAAGGTGCGAAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGa < 1:31367/251‑1 (MQ=255) cgATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg > 2:195128/1‑251 (MQ=255) aTCCACAGGGAGATGCTCATAATGATGCCACCCGAGATTAAGGTGCGGAACAGATTTCAACGAAGCACGGCCACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCACACGGCACCACCTGATTACCCGGAACACAAAAAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCATGCCGATACCACCC‑‑CCGTATGTCCCAGCAGCAACGCCGGACAAAGACCAATCAGACACTCCTGAccgcc > 2:172718/1‑251 (MQ=255) atatttgtgaCACCCGAGATTAAGGGGGGGAAAGGATTTCCGTGATCCAGGGCGACCGCCATCGCCACTAATAAGCTGAGGGTGGAAAGACCGCCAAAGGGCAGCACCTGATCGCCCGGCGCAGTAACACCATTTAAAAGGGGGGGTGGGATAGAAACCAGTTTTGCGGTTCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCATGGCCTGTCAGGAGCTCCTGACCGCCGAACTTCGCCTGTAAACg < 1:250913/244‑1 (MQ=255) acccgaccttacggTGCGGAACAGACTTCCACGCTGCACGGCGACCGCCATCGCCACGACGAAGCCGATGGTGGCAAGCTCGCCAAACGCCAGCACCTGATTACCCGGCCCACCCACAGCAATTAAAATGGTGCGTGGGATACAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAACGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACCAGCCt > 2:41521/13‑251 (MQ=255) cGGAACGGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAACATGGTGAGTGGGATACAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGt > 1:166002/1‑251 (MQ=255) ccGATGCACGGCAACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGACCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGGTAAACCATCCatga < 1:62536/250‑1 (MQ=255) aCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATCAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTAAGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAtt > 1:125669/1‑249 (MQ=255) agCCGATGGGGGCAAGATCGCAAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATc < 1:162722/251‑1 (MQ=255) cGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACt < 1:195128/251‑1 (MQ=255) gATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACATCAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTg < 1:37572/251‑1 (MQ=255) gATCGACAAACGGCAGCACCTGATTACCCGGAACACACACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGACGACACCACC‑‑GCCGTATGTCCCAGCAGCAAAGCCGGGTCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTCGCGATGGGCGTTAAACCAGCCATGATGGATGTACTCACCCGTGGCATTAGCAGCATCACTGCAGCCGttttt > 2:246555/1‑251 (MQ=255) cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc < 1:186784/212‑1 (MQ=255) cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc > 2:186784/1‑212 (MQ=255) cGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCt > 1:275503/1‑251 (MQ=255) acaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTATCGATGGGCGTTAAACCATCCATGAGGGGTTTAATCACCCGTGGCATTAGCAGCATCACTCCCGCCGTTTTTACCGCCAGCGGCAATACACCTTTGACATCGTAAc > 2:143973/1‑251 (MQ=255) aaaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAc > 1:40620/1‑170 (MQ=255) aaaaGGGTGAGTGGGATAAAAACCAGGCTTGCCGATAC‑‑CCCCGCCGTAGGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGcc < 1:174811/251‑1 (MQ=255) gggTGAGTGGGATAAAATTCGGGCTGCCCGAT‑CC‑CTCCGCCGTATGTCCCATCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTGCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGGATCGCTGCCGCCGTTTTTGGCGCCAGCTGCAGTACACCTTTGACATCGTAACCCGCGAGGATGCCGATa < 1:41521/250‑1 (MQ=255) gCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACgg < 1:156145/250‑1 (MQ=255) cgatacccccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGTCGATAATCAGCCCCATCACAAAACCCACGGTGa < 1:16672/245‑1 (MQ=255) gataccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAc > 2:232406/8‑251 (MQ=255) cccaccgccgTTTCTCCCAGCAGCAACTCCGAATAAAGGCCAATCAGGAACTCTTGCCCGCCGAACTTCGCCTGTAAACGACTACGGGCTTGCTTAGCGTTGGGCGTTAAACCATCCATGATGGGTTTAACCACCCGTTGCATTAGCAGCATCACTGTCGCAGTTTTTACCGCCAGCTGCAATACACCTTTGACAGCGTAACCCTCGAGGATGCCGATAATCAGCCCCATCACAAAAGCCACGGTTACAgg < 2:180675/247‑1 (MQ=255) | GGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGG > NC_000913/2173139‑2173609 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |