Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 2,173,363 | Δ2 bp | intergenic (‑490/+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,363 | 0 | C | . | 95.2% | 76.7 / ‑3.9 | 21 | intergenic (‑490/+919) | gatD/gatB | galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component |
Reads supporting (aligned to +/- strand): ref base C (0/1); new base . (17/3); total (17/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.90e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
* | NC_000913 | 2,173,364 | 0 | C | . | 95.2% | 77.3 / ‑1.9 | 21 | 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/1); new base . (17/3); total (17/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.90e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01 |
GCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGG > NC_000913/2173123‑2173609 || gCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTcc > 1:253111/1‑251 (MQ=255) gTTGGGTGGGCAACCCAAAGGTTTGCTGCGGGATCCACAGGGTGATGCTCAAAATGAGGACCCCCGAGCTTAAGGTGCGAAAAAGATTTCCACGATGCACGGCGACCGCCTCCGCCCCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCCCACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTGGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGcc < 2:90089/251‑1 (MQ=255) aTGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATATCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACACACAGAAATTAACATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTACTGACCGCCGAACTTCGcact > 2:347579/1‑248 (MQ=255) agaTTCAGGTGCGGAACAGAATTCAACGATGCACGGCGACCGCCATAGCCAAGAAGACGCCGATGGTGGCCAGATCGCCAAAAGGCAGCACCTGATTACCCGGCACACAAACACCCATTCAAATGGGGAGTGGGAAAAAAATACGGCTTGCCGATACAACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCCGTCAACGCCCACGAACCTGCtt > 2:268826/1‑250 (MQ=255) attaaggcgaggcacccaTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGACGGTGGCAAGCTCGCCAACCGGCAGCACCTGCTTAACCGGCACCCAACCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGCTCACGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGc > 2:237127/18‑251 (MQ=255) cgcGATGCACGGCGCCCGCCTTCGCTACGAAGATGCCGTTGGTGCAAAGAACGCAAAACGGAGGCTCTTGTTCACCCGGCACACACCCAGCAATTAAAATGGTGAGTGGGATAAAAATAAGCCTTGCAGATACCACC‑‑GCCGTATGTCCTAGCAGCAACTCCGTATCAAGGCCAATCAGGACCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATTGGCGTTAAACCATCCatg < 2:284876/249‑1 (MQ=255) aCGGCGACCGCCATCTCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACCAACCGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt > 2:136468/1‑251 (MQ=255) cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGACAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCAGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc > 1:304143/1‑251 (MQ=255) aTCGCCACGAAGAAGCCGATGGTGGAAAGACCGAAAAAGGGCAGCACCTGATTACAAGGCACACAAACAGAAATTAAAATGGCGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCAAGCAGAAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTgg > 1:405820/1‑251 (MQ=255) aCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCa > 1:56308/1‑216 (MQ=255) agcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCa > 1:340538/1‑250 (MQ=255) aCCCGGCACACAAACAGCCATTAAAATGGTGAGTGGAATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGGGCCTGATTAGCGAGGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCCTAACTGACGCCGTTGTAACAGACAGCTGCAATACACCTTTg > 2:276136/1‑250 (MQ=255) cccGGCACACAAAAAGCAGTTAAAATGGTGAGTGGGAGAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCACAGCAGCAACGCAGGATCAAGGCAAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATACATGATGGGTTTAATCACCCGTTGCATTAGCAGCCTCACTGCCGCCGTTTTTACCGCCAGCTGCAATACAccgatgac > 2:215770/1‑245 (MQ=255) ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGAc > 1:86446/1‑250 (MQ=255) aaaaTTAAATTGGTGAGTGGAATAAAACCCTGTCCTCCCGATTCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAGGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGCCTACGAGCCTGCTTAGCGATGTGCGTTAAACCAGCCATGTTGGGTTTAATCACCCGTGGCATTAGCAGCATGACTTCCGCCGTTTTTACCTCCAGCTGCAATACACCTTGGACTTCGTAACCCGCGAg < 1:135201/249‑1 (MQ=255) aaaaaGGGTGAGGGGGATAAAAATAAGGCTGGCCGAT‑CC‑CCCCGCCGTATGGCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATg < 1:136468/249‑1 (MQ=255) ggtaTAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTACCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCCCCCGTGGCATTAGCAGAATCACTGCCGCCGTTTTTACCGCCAGCTGCAATCCACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGccc > 2:81897/4‑251 (MQ=255) gggATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGccc > 1:15074/1‑251 (MQ=255) ggATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGcc > 1:100944/1‑238 (MQ=255) cAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAc > 1:380654/1‑251 (MQ=255) gCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCCGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCCGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGCTACTCAGCCCCATCACAAAAccc > 2:348813/1‑250 (MQ=255) accaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAgg > 1:8236/5‑251 (MQ=255) ccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGCCTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgc > 1:157361/4‑157 (MQ=255) || GCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGG > NC_000913/2173123‑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: ‑ |