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.89.7% 101.3 / 6.6 29intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (16/10);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 7.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01
*NC_0009132,173,3640C.89.7% 103.4 / 4.0 29intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (16/10);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 7.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGC‑GACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAA‑TCAGGCTTGCCG‑ATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGC  >  NC_000913/2173128‑2173614
                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                                          
aGCCGCCAGTTGGGTGAGCAAACCAATAGTTTGCGTCGCGCTCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGC‑GACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATACAAA‑TCAGGCTGCCCG‑ATACCACA‑‑GCCGTATGTCCcagc                                                                                                                                                                                                                                             >  2:165388/1‑250 (MQ=255)
   cgccAGCTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGC‑GACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAA‑TCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCagcagc                                                                                                                                                                                                                                          >  1:122053/1‑250 (MQ=255)
                              tggCGTCGCGATCCACAGGGTGATGCTCATAAGGATGCCCCCCGAGATTAAGGTGCGAAACAGATTTCCACGATGCACGGC‑GACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAA‑TCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAAc                                                                                                                                                                                                               <  1:1126/248‑1 (MQ=255)
                                gCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGC‑GCCCGCCAGCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATGACCCGGCACACAAACAGCAATGAAAATGGTG‑AGTGGGAAAAAAA‑TCAGGCTTGCCG‑AAACCACC‑‑GCCTTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTc                                                                                                                                                                                                             <  1:244468/250‑1 (MQ=255)
                                               gggTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGC‑GACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAA‑TCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                                                                                                                                             >  1:158160/1‑251 (MQ=255)
                                               gggTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGC‑GACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAA‑TCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                                                                                                                                             >  1:370187/1‑251 (MQ=255)
                                                ggtgcggCTCATAAGGATGACTCCCGAGATTAGGGTGGGGAACAGATGTCCCCGGTCCCCGGG‑GACCGCCACCGCCACGAAGAAGCCAATGGGGGCAAGAGCGCAAACCGGCAGCACCTGATTACCCGGCACCCAAACAGCAATAAAAAGGGTG‑AGTGGGATAAAAA‑TCGGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGAACTTCg                                                                                                                                                                                            <  2:349984/245‑1 (MQ=255)
                                                ggTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGC‑GACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAA‑TCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                                                                                                                                            >  1:126864/1‑251 (MQ=255)
                                                          aTAATGATGCCACCCGAGGTTACGGAGCGGAACCGATTTCCACGATGCACGGC‑GACCGCCATCGCCACGACGAAGCCCATGGAGCCCAGCTCGCCACACGGCAGCCCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGAGGGATAAAAA‑TCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGAAACGCCGGAACAAGGCCAATCAGGAACACCTGACCGCCGAACTTCGCCTGTAAAcc                                                                                                                                                                                  >  2:429552/1‑250 (MQ=255)
                                                            aaTGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGC‑GACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAA‑TCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCCGGAACTCCTGACCGCCGAACATCGCCTGTaacagac                                                                                                                                                                                >  2:85834/1‑246 (MQ=255)
                                                             atgcggccccccGAGATTAGGGGGCGGAACAGATTCCCCCGATG‑ACGGCCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAGGGTG‑AGTGGGATAAAAA‑TCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                                                                                               <  2:126864/242‑1 (MQ=255)
                                                                                   ggCGGAACAGATTTCCACGATGCCCGGC‑GACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAA‑TCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACGCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                                                                                                                                         <  1:254493/250‑1 (MQ=255)
                                                                                         aCAGATTCCCACGGTCCACGGC‑GACCGCCACCGCCCCGAAGAAGCCGATGGTGGAAAGAACGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAAGGGTG‑AGTGGGATAAAAA‑TCAGGCTTGCCG‑ATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                                                                                                                                   <  1:301688/251‑1 (MQ=255)
                                                                                            gATTTCCACGATGCACGGC‑GACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAA‑TCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcc                                                                                                                                                 >  1:61788/1‑250 (MQ=255)
                                                                                                       aggcccggagccggcaaTCGCCACGAAAAAGCGGATGGTGGAAAGAGCGCCAAACGGCAGCACCTGATTCCCCGGCACACAAACAGAAATTAAAATGGTG‑AGGGGGATAAAAA‑TCAGGCTTGCCG‑ATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg                                                                                                                                      <  2:122053/235‑1 (MQ=255)
                                                                                                           cGGC‑GACAGCCATCGCCCCGAAGAAGCCGCTGGTGGCAAGACCGCAAAACGGCAGCACCTGATAACCCGGAACCCAAACAGCAATTAAACAGGTG‑AGTGGGATCAAAA‑TCAGGCTAGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGAAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt                                                                                                                                  >  2:184350/1‑250 (MQ=255)
                                                                                                                      cATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAA‑TCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTg                                                                                                                       >  2:450542/1‑251 (MQ=255)
                                                                                                                          gCCACGAAGAAACCGATGGTGGCAAGATCGCCAAACGGCCGCACCTGATTACCCGACACAAAAACAGAAATTACAATGGTG‑CGTGGGAAAAAAC‑TCAGGCTTGCCG‑ATACCACC‑‑GCAGTATGTCCCAGCAGCAACGACGGATCAAGGCAAATCAGGAACTCCTGACCGCCGAACTTAGCCTGTAAACGACTACGAGCCTGCTTCGCGATGGGCGTTAAACCATCCATGATGGGATTAATAACCCGAGGCAt                                                                                                                   >  2:255031/1‑251 (MQ=255)
                                                                                                                                                 aaGATCGCCAAACGGCAGCACCTGATTAACCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAA‑TCCGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAACAGCATCACTGCCGCCGttt                                                                                            >  2:10260/1‑251 (MQ=255)
                                                                                                                                                  gaacctcaaatgggaagccccggtttACCCGGCACACAA‑CAACAATTAAAAAGGGG‑GGGGGGGTAAAAAATAAGGTTTGCGG‑TTACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGtttt                                                                                           <  2:317540/227‑1 (MQ=255)
                                                                                                                                                  agaccgcaaaccGGCAGCTCCTGATTACCCGCCACAAAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAA‑TCAGGCGTGCCG‑ATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGtttt                                                                                           <  1:85834/240‑1 (MQ=255)
                                                                                                                                                                 agcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAA‑TCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCa                                                                             >  1:113637/1‑250 (MQ=255)
                                                                                                                                                                            ccccgaccaacaaaagaaaataaaaaGGGTG‑AGTGGGAAAAAAA‑CCGGGCT‑GCCGGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAAt                                                                           <  2:167211/219‑1 (MQ=255)
                                                                                                                                                                            ccccGGCACACAAACAG‑AAAAAAAAAGGTGGCGTGGGAAAAAAA‑TCGGGCTT‑‑CG‑CTATCCCCCCGCCGTATTCCCCAGCAGCACCGCCGGATCAGGGCCAATCAGGAACCCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTg                                                                  <  2:61788/249‑1 (MQ=255)
                                                                                                                                                                            aCCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAA‑TCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAAt                                                                           >  1:167211/1‑241 (MQ=255)
                                                                                                                                                                              ccGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAA‑TCAGGCTTGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaa                                                                                                                                                    >  2:257094/1‑166 (MQ=255)
                                                                                                                                                                              ccGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAA‑TCAGGCTGGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaa                                                                                                                                                    <  1:257094/166‑1 (MQ=255)
                                                                                                                                                                                                    aaggggg‑AGTGGGATAAAAA‑TCAGGTTTGCCG‑ATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGa                                         <  1:95258/245‑1 (MQ=255)
                                                                                                                                                                                                                       aaaaaaacggggtgccgtatccccccGCCTTAGGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTTTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAGCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATc                           <  2:271888/229‑1 (MQ=255)
                                                                                                                                                                                                                                ttGCCG‑ATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCAGCTTAGCGATGGGCGTTAAACCATCCATGATGGGCTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACCTCGTAACCCGCGAGGATGCCGCTAATCAGCCCCATCACAACACCCACg               >  2:337922/1‑251 (MQ=255)
                                                                                                                                                                                                                                 cctggccgaccccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTGAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTGAATCACCCGTGGCATTGGCAGCATTACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACAGCGTACCCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCa                 <  1:429552/237‑1 (MQ=255)
                                                                                                                                                                                                                                                cgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGc  >  2:408146/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                                          
AGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGC‑GACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTG‑AGTGGGATAAAAA‑TCAGGCTTGCCG‑ATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGC  >  NC_000913/2173128‑2173614

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: ‑