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_0009132,173,3630C.87.5% 49.1 / ‑1.1 16intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (11/3);  total (12/4)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3640C.93.8% 54.2 / ‑1.9 16intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (12/3);  total (12/4)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

CCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACC  >  NC_000913/2173130‑2173561
                                                                                                                                                                                                                                          ||                                                                                                                                                                                                     
ccgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAACCGATGGAGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAAAAGCAATTAAAATGGTGAGAGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTCTGTCCCagcagc                                                                                                                                                                                     >  2:347720/1‑251 (MQ=255)
                     cAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                                                                                                >  1:123233/1‑251 (MQ=255)
                       aTAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATccg                                                                                                                                                              >  1:90310/1‑249 (MQ=255)
                                        ccaaaGGGTGAGGCTCAAAATGATGCAACCCGAGATTAAGGGGCGGAACAGATTCCCACGTTGCAGGGCGCCCGCCACCGCAACAAAGAAGCCGAGGGGGGAAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGAAATTAAAAGGGTGAGTGGGATAAAAACCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGa                                                                                                                                             <  1:42265/247‑1 (MQ=255)
                                         cacaGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGCCCGCCGa                                                                                                                                             >  1:232691/1‑250 (MQ=255)
                                          acaGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTACGGAACAGCTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGAGGCAAGATCGCCCAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCAGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCAGGCTCAAGGCCAATCAGGAACTcccaacacccgcac                                                                                                                                           >  2:140039/1‑239 (MQ=255)
                                                             gaggaccccggagtttagggggcgaaaCAGATTCCCCCGATGCACGGCGCCCGCCATCGCCACGAAAAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTAc                                                                                                                        <  1:58314/226‑1 (MQ=255)
                                                                                      aaCAGATTGCCACGATGCAAGGCGACCGCCATCGCCACGAAGAAGCCGTTGTTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGACCACGAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCAAGCAGCAACGCAGGATCAAGGCCAATCGGGAACTCATGACAGCCGAACTTCGCCTGTAAACGCCTACAAGACTGCTTACCGACGGGCGtac                                                                                                >  2:265641/1‑248 (MQ=255)
                                                                                                         cGGCGACCGCCATCGACACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCACTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGCTGGGCGTTAAACCATCCATGATGGGttt                                                                             >  1:218420/1‑250 (MQ=255)
                                                                                                                                 ggCCGATGGTGGCAAGACCGCCAAACGGCTGCATAGTAGTCCCCGCCACACAAACAGAAATAAAAAGGTGGGGTGGTAAAAAAATCGGTCTTGCCGTACCCCCCGCCGTATGTCCCAGCAGCACCGCCGGATCAATGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATc                                                    <  2:325339/250‑1 (MQ=255)
                                                                                                                                                      cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCa                                                                        >  2:325526/1‑210 (MQ=255)
                                                                                                                                                                   cTGATTACCCGGCAAACAAACAGCAATTAAAATGGAGAGTGGGATAAAAATCAGGCTTGCCGATACACCCAACGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAACGCCGAACTTCGCCTGCAAACGACTACGAGACTTCTTAGCGATGGACATTAAACCATCCATGATGGCTTTAATAACACGTGGCATTAGCATCCTCACTGCCGCCGTGTTTACCGCCACCTGCAATacac                  >  2:321451/1‑251 (MQ=255)
                                                                                                                                                                      aTTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATAAGGCTTGCCGATACCACA‑‑GCCGTATGTCCCAACAGCAACGCCGGAGCAAGGCCAATCAGGAACTCCTGACCGCCTAACTTCGCCTGTAAACTACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTAGCATTACCAGCATCGCTCCCGCCGTTTTTACCACCAGCTACAATACACCtt               >  1:185666/1‑251 (MQ=255)
                                                                                                                                                                                acacaAACAGCGATTAAAATGGTGAGTGGGATAACAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCACCGACGGATCAAAGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACTACTACGAACCTGCTTAGCGATGAGCGTTGAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCACCTGCACTACACCttatccatcgt      >  2:296054/1‑241 (MQ=255)
                                                                                                                                                                                acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGt      >  1:229249/1‑250 (MQ=255)
                                                                                                                                                                                   aaacaagaaataaaaaTGGAGGGTGGGAAAAAAACCAGGTTCGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTACTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAcc  <  1:203762/239‑1 (MQ=255)
                                                                                                                                                                                                                                          ||                                                                                                                                                                                                     
CCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACC  >  NC_000913/2173130‑2173561

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: