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.80.7% 75.4 / 21.0 31intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/5);  new base . (17/8);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 5.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01
*NC_0009132,173,3640C.80.7% 77.1 / 19.3 31intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/5);  new base . (17/8);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 5.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

CATTAGCCGCCAGTTGGGTGTGC‑AAACCAATAGTTTGCGTCGCGATCC‑ACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATC  >  NC_000913/2173124‑2173589
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                 
cATTAGCCGCCAGTTGGGTGTGC‑AAACCAATAGTTTGCGTCGCGATCC‑ACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGCCCGCCATCGCCACGCAGAAGCCGATGGTGGAAAGATCGCCAAAC‑GGCAGCACCTGCTTACCCGGAACACAAACAGCAATTAACATGGTGCGTGGGATAAAAATCAAG‑CTTGCCTATACCACC‑‑GCCGTATGTccc                                                                                                                                                                                                                       >  2:185117/1‑251 (MQ=255)
    aGCCGCCAGTTGGGTGTGC‑AAACAAATAGTTTGCGTCGCGATCC‑ACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTCCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCcagc                                                                                                                                                                                                                    <  1:7256/250‑1 (MQ=255)
       cgccGGTTGGGT‑TGCCAAACCAATAGTTTGCGTGGC‑ATCCCACAGGGTAATGCTCATAATGATGCCCCCCGAGATTAAGGGGCGGACAAGATTCCCCCGAGGCCCGGCGCCCGCCATCGCCACGAAAAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAgcagca                                                                                                                                                                                                                <  2:308417/251‑1 (MQ=255)
              tGGGTGTGC‑AAACCAATAGTTTGCGTCGCGATCC‑ACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCgg                                                                                                                                                                                                         >  1:99614/1‑251 (MQ=255)
                      c‑AAACCAATAGATTGCGTCGCGATCC‑ACAGGGTGATGCTCATAATGCTGACACCCGAGGTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACACACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAGCGCCGGATCAAGGc                                                                                                                                                                                                 >  2:97294/1‑251 (MQ=255)
                        aaaCCAATAGTTTGCGTCGCGATCC‑ACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                                                                                                                >  2:71425/1‑251 (MQ=255)
                           aaaaTAGTTGGCGCCGCGACCC‑ACGGGGTAATGCCAAAAATGATGACCCCCGAGATAAGGGTGGAGAACAAATTCCCCCGAGGCCGGGCCACCGCCCTCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTCCCCGGCACACAAACAGCAATAAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCGGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                                                                                                                                                                                             <  2:99614/249‑1 (MQ=255)
                            cAATAGTTTGCGTCGCGATCC‑ACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                                                                                                                            >  1:211489/1‑251 (MQ=255)
                             aaTAGTTTGCGCCGCGACCC‑ACAGGGGGATGCTAAAAATGATGCCACCCGAGATTAAGGGGCGGAACAAATTACCACGATGCACGGCGACCGCCAGCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCGGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCAGT‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCa                                                                                                                                                                                           <  2:53022/251‑1 (MQ=255)
                                    tGCGTCGCGATCC‑ACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCATGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTc                                                                                                                                                                                    >  1:374829/1‑251 (MQ=255)
                                                                   atgatgACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                                                                      >  1:155312/1‑251 (MQ=255)
                                                                     aatgaCGCCCGAGATTAAGGGGCGGAACAGATTTCCCCGAGGCGCGGCTACCGCCATCGCCACGAAGAAGCGGATGGTGGCAAGACCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGAAAAAAAATCAGG‑TTGGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCTGATCAAGGCCAATCAGGAACGCCTGACCGCCGAACGTCGCCTGTAAACGACTAc                                                                                                                                                    <  1:170259/250‑1 (MQ=255)
                                                                      aggaCACCCGAGCTTAAGGTGCGGAACAGATTTACACGATGCACGGCGACCGCCATCGGCACGAAGAAGCCGATGGTGGCAAGATCGCCACAC‑GGCATCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGCGTGAGATAAACATCAGA‑CTTACCGCTA‑‑ACACCGCCGGATGTCCCAGCAGCAACGCCGAATCAAGGCCAATCAGGCAATCCTGCCCGCCGAACTTCGCCTGTaacaggctcc                                                                                                                                                    >  2:153515/3‑242 (MQ=255)
                                                                       tgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                   >  1:285395/1‑250 (MQ=255)
                                                                                agaTTAAGGTGCGGAACAGATTTCCAAGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGAATGTCCCAGCAGAAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAACCGACTACGAgcttgctt                                                                                                                                          >  1:166008/1‑250 (MQ=255)
                                                                                                        cactatccccgcagaccgcatccGCCACGAAGAACCCAATGGTGGACAGTCCCCCAAAC‑GGCCGCACCTGATTCCCCGGCACACAAACAGCATTTAAAATGGTGAGGGGGATAAAAACCAGG‑TTTGCCGA‑ACC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCatg                                                                                                                 <  2:64349/229‑1 (MQ=255)
                                                                                                                                 cACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAtt                                                                                         >  1:58896/1‑250 (MQ=255)
                                                                                                                                       aaccccgggggtgcaaaGATC‑CCAAACCGGCAGCACCGTATTCCCCGCCACACAAACAGAAATAAAAAGGGGGGGTGG‑AAAAAAAACAGGGTTTCCCGAT‑CC‑CCCCGCCGTATGTCCCAGAAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAgcggca                                                                                  <  2:155312/236‑1 (MQ=255)
                                                                                                                                                       cctcttccgatcgggCAGCCCCTGATTCCCCGGCACACAAACAGAAATAAAAAGGGTGAGGGGGATAAAAATCAGG‑CTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGttt                                                                   <  1:345837/238‑1 (MQ=255)
                                                                                                                                                        gagcccaaaaC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATAAAAATGGGGTGTGGGAAAAAAACCAGG‑CTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGttttt                                                                 <  1:255302/244‑1 (MQ=255)
                                                                                                                                                          ccGCCAAAC‑GGAAGCCCCGGATTCCCCGCCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑TTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAc                                                               <  1:199207/250‑1 (MQ=255)
                                                                                                                                                             ccAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCCTGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGc                                                            >  2:312117/1‑251 (MQ=255)
                                                                                                                                                                    ggCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGttt                                                                   >  2:345837/1‑238 (MQ=255)
                                                                                                                                                                      cagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGc                                                     >  1:86537/1‑250 (MQ=255)
                                                                                                                                                                        ccaCCTAATTCCCCGGCACAAAAACAGCAATTAAAATGGTGAGTGGTAAAAAAATCAGT‑CTTGCCGTCAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGCCCAATCAGGAACTCCTGACCGCCGAACTTCTCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAAt                                                  <  1:19029/250‑1 (MQ=255)
                                                                                                                                                                                          cacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTCCCGCCAGCTGCAATACACCTTTGACATCGTa                                 >  2:368161/1‑250 (MQ=255)
                                                                                                                                                                                                          aaTGGTGAGTGGGCTAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCg                 >  1:184365/1‑250 (MQ=255)
                                                                                                                                                                                                                gAGTGGGATAAAAATCAGC‑CTGGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCa          <  1:93612/251‑1 (MQ=255)
                                                                                                                                                                                                                 aGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACCCCTTTGACATCGGAACCCGCGAGGATGCCGCTAATc           >  2:312201/1‑249 (MQ=255)
                                                                                                                                                                                                                    gggATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGccc      >  1:151835/1‑251 (MQ=255)
                                                                                                                                                                                                                        aaaaaaCCAGG‑CTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATc  <  1:71425/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                 
CATTAGCCGCCAGTTGGGTGTGC‑AAACCAATAGTTTGCGTCGCGATCC‑ACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAAC‑GGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATC  >  NC_000913/2173124‑2173589

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