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.90.6% 118.9 / 4.8 32intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base . (18/11);  total (19/13)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01
*NC_0009132,173,3640C.90.6% 119.5 / 4.2 32intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base . (18/11);  total (19/13)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACA‑AACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTG  >  NC_000913/2173124‑2173604
                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                
cATTAGCCGCACGTTGGGTGTGCCAACCAATAGTTTGCGTCGCGATCCACAGGGCGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGCTGGTGGCAAGATCGCCAAACGG‑CAGCAACTGATTACCCGGCACACACAAC‑GCAATTAACATGGTGAGTGGGCTAAAAATCAGGCTTGACGATA‑‑ACACCGCCGTATGTccc                                                                                                                                                                                                                                      >  2:368944/1‑251 (MQ=255)
            gTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACA‑AACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGcc                                                                                                                                                                                                                          >  2:91780/1‑251 (MQ=255)
                 atatgCACACCAATAGTTTGCGTCGAGATCCACCGGGTGCTGCTCATAATGATGCCACCCGCGATTAAGGTACGGAACAGATTTCCACGATGCACGGCGACCGACAACGACACGAAGAAGACGATGGTGGAAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACA‑AACAGCAATTAAAATGGGGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGGCCCAGCAGCAACGCCGGATc                                                                                                                                                                                                                     >  2:233473/4‑251 (MQ=255)
                            aaTAGTTTGCGTCGCGATCCACAGGGTGATGCTCAAAATGATGACACCCGAGATTAAGGTGCGGAACAGATTCCCACGATGCACGGCGCCCGCCAGCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACA‑AACAGCAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCa                                                                                                                                                                                                          <  2:159968/251‑1 (MQ=255)
                               aGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACA‑AACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                                                                                                                                        >  1:401090/1‑250 (MQ=255)
                                                            ccATAATGATGACCCCCGAAATTAAGGTGGGGAACAGATTTCCCCGATGCCCGGCGCCCGCCAGCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACA‑AACAGCAATGAAAATGGTGAGGGGGATAAAAATCAGTCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaaa                                                                                                                                                                          <  1:271493/250‑1 (MQ=255)
                                                                             gtgctctcccaatctaaCAGATTCCCACGATGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACA‑AACAGCAATTAAAATGGTGAGTGGGATAAAAACCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAgcttgctt                                                                                                                                                         <  1:366234/236‑1 (MQ=255)
                                                                             gagaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACA‑AACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAgcttgctt                                                                                                                                                         >  2:114906/1‑251 (MQ=255)
                                                                               gaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACA‑AACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAgcttgctt                                                                                                                                                         >  1:228165/1‑249 (MQ=255)
                                                                                        gcgaaaCAGATTCCCCCGATGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTCCCCGGCACACA‑AACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCg                                                                                                                                              <  1:11451/247‑1 (MQ=255)
                                                                                          ggAACAGTTTCCCACGAGGCACGGCGACCGCCACCGCCACGAAGAAGCCGATGGGGGCAAGATCGC‑AAACGGGCAGCCCCTGATTACCCGGCACACA‑AACAGCAATTAAAATGGTGAGGGGGATAAAAACCAGGCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGTCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGtt                                                                                                                                            <  2:350357/251‑1 (MQ=255)
                                                                                            aaCAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAACCGG‑CAGCACCTGATTACCCGGCACACA‑CACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAgcttgctt                                                                                                                                                         >  2:366234/1‑236 (MQ=255)
                                                                                                     ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACA‑AACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTGAACGACTACGAGCTTGCTTAGCGATGGGCGTTCAAccatccat                                                                                                                                 >  2:266553/1‑251 (MQ=255)
                                                                                                                        ccATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACA‑AACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGt                                                                                                              <  1:91780/251‑1 (MQ=255)
                                                                                                                         caccGCCCCGAAGAGGCCGATGGGGGCAAGATCGCAAAAGGC‑CAGCCCCTGATTACCCGGCACACA‑AACAGAAATAAAAAGGGGGAGTGGGATAAAAATCAGGTTTGCCAATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTg                                                                                                             <  2:334811/248‑1 (MQ=255)
                                                                                                                          aTCGCCACGCAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACA‑AACAGCAATTAAAATGGTGCGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTgg                                                                                                            >  2:240965/1‑251 (MQ=255)
                                                                                                                                                      gATCGCCAAACGG‑CAGAACCTGATTACCCGGCACACA‑AACAGCAATTAAAATGGTGAGAGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTCAACGACTACGAGCTTGCTTAGCGATGGTCGTTCAACCATCAATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGtttt                                                                                 >  2:215698/1‑250 (MQ=255)
                                                                                                                                                        tCGCCAAACGG‑CAGCACCTGATTACCCGGCACACA‑AACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAc                                                                              >  2:195057/1‑251 (MQ=255)
                                                                                                                                                           ccAAACGG‑CAGCACCTGATTACCCGGCACACA‑AACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGc                                                                           >  1:326752/1‑251 (MQ=255)
                                                                                                                                                                 gg‑CAGCACCTGATTACCCGGCACACA‑AACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTg                                                                     <  1:240965/251‑1 (MQ=255)
                                                                                                                                                                 gg‑CAGCACCTGATTACCCGGCACACA‑AACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTg                                                                     >  1:89753/1‑251 (MQ=255)
                                                                                                                                                                          cTGATTACCCGGCACACA‑AACAGCAATTCAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCACTCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGCCTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATacac                                                             >  2:45203/1‑251 (MQ=255)
                                                                                                                                                                                          ca‑AACAGCAATTAAAATGGTGAGTGGTATAAAAATCAGTCTTGCCTTTCCCACC‑‑GCCGTTTGTCCCAGCAGCAACGCCGGATCAATGCCAATTAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCGTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCTTCACTGCCGCCGTTGTTACCGCCAGCTGCAATACACCTTTGACATCGTAAcc                                             <  1:364133/251‑1 (MQ=255)
                                                                                                                                                                                             aaCAGCAATTAAACTGGCGAGAGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGAATGTCACAGCAGCCACGACGGATCAAGGCCAACCAGGAACTCCAGACCGCCGAACTACGCCTGAAAACGACTACGAACTTGCTTAGCAATGGGCATTAAACCATCCATGATGAGTTTAATCCCCCGTTACATTAGCAGCATCACTGCCGCCTTTTTTACCGCCAGCTGCAATACACCTTTGACATCGAAAcacc                                           >  2:414545/1‑248 (MQ=255)
                                                                                                                                                                                             aaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGttttt                                                                                <  1:37087/214‑1 (MQ=255)
                                                                                                                                                                                             aaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGttttt                                                                                >  2:37087/1‑214 (MQ=255)
                                                                                                                                                                                                 gCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcgcg                                         >  1:245830/1‑249 (MQ=255)
                                                                                                                                                                                                                    ggATAAAAACCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGAAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCTTTAAACCAGCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACTCCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGcccc                    <  1:190738/251‑1 (MQ=255)
                                                                                                                                                                                                                         aaaaTCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcac               >  1:78912/1‑251 (MQ=255)
                                                                                                                                                                                                                          aaaccAGGCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATTACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcaca              <  1:114906/247‑1 (MQ=255)
                                                                                                                                                                                                                                     gCCGATACCACC‑‑GACGTCTGTACCAGCAGCACCGCCGGATCCAGGCCAATCAGGAACTCATGACCGCCGAACTTCGCCTGTAAACGCCTACGAGCTTGCTTCGCGATGGGCGTTAAACCATACATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACggg   >  2:392419/1‑250 (MQ=255)
                                                                                                                                                                                                                                      ggcccaacccgccgccgTAGGTCCAAGCAGCAACGCCGGAACAAGGCCAAACAGGACCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGTTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAc      <  2:84828/238‑1 (MQ=255)
                                                                                                                                                                                                                                         gatccaccCC‑CCTTATCTCCCAGCAGCAACGCCGGTTCAAGGCCAATCATGACCACCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAACCCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTg  <  2:83565/246‑1 (MQ=255)
                                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                
CATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACA‑AACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTG  >  NC_000913/2173124‑2173604

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